Monolithic kinetic algorithm for heterogeneous porous media systems using a continuous one-domain approach

技術要約：不均一多孔質媒体システムのためのモノリシック・キネティック・アルゴリズム

問題提起
不均一な多孔質媒体における流体の流れ、特に自由流領域と多孔質構造を結合させたシステムのシミュレーションは、大きな課題を提示している。従来のアプローチでは、「2ドメイン」定式化に依存することが多く、これは自由流領域と多孔質領域に対して個別の支配方程式を解き、それらを結合するために明示的な界面境界条件（例：Beavers-Joseph条件）を必要とする。あるいは、「1ドメイン」アプローチを用いる手法もあるが、これらは歴史的に不均一媒体における一貫性の確保に苦慮してきた。

多孔質媒体のための既存の格子ボルツマン法（LBM）の定式化は、回収されるマクロな方程式における空隙率のスケーリングが不整合である場合が多い。具体的には、従来のモデルは、空隙率が空間的に変化する場合に、圧力勾配および粘性応力項を空隙率に対して正しくスケーリングできないことが多い。さらに、多くの定式化では、ブリンクマン粘性応力が（真の）固有速度ではなく、表観速度（ダルシー速度）の勾配を用いて表現されているため、不均一な領域において不正確さを招いている。加えて、標準的なLBMモデルは、有効体積粘性や音速を独立して調整する柔軟性に欠けることが多く、計算効率や安定性を制限している。

手法
著者らは、連続的な1ドメイン・アプローチに基づいた、モノリシックな単一ドメイン格子ボルツマンモデルを提案している。この手法は、以下の柱によって構築されている：

1. マクロモデル： 対象となる流体力学は、体積平均化法を用いて一般化ダルシー・フォルサイマー・ブリンクマンモデルへと導出される。このモデルは、全計算領域において有効であり、自由流と多孔質領域の間の遷移を物理特性（空隙率 $\epsilon$）の連続的な変化として扱う。支配方程式には以下が含まれる：

   • 体積平均化された密度（$\epsilon\rho$）に関する連続方程式。
   • ダルシー抵抗、フォルサイマー慣性抵抗、および有効粘性応力テンソル（$\sigma_{eff}$）を組み込んだ運動量方程式。
   • 極めて重要な点として、有効粘性応力は（固有速度 $u = v/\epsilon$ ではなく）表観速度 ($v$) の勾配を用いて表現されており、これにより可変空隙率におけるブリンクマン定式化との整合性が確保されている。

2. キネティック・ソルバー： モデルは、標準的な格子（例：D2Q9, D3Q19）上で単一の一般化されたキネティック方程式を利用する。

   • 平衡分布： 正しいモーメントとガリレイ不変性を保証するために、積形式の離散平衡分布が採用されている。これには、粘性テンソルの等方性を回復し、体積粘度を独立して調整することを可能にするための特定の補正項（$\Phi_{ext}$ および $\Phi_{\eta}$）が含まれる。
   • 強制項スキーム： ダルシー・フォルサイマー抵抗、外力、および必要な空隙率勾配の補正を組み込むために、一貫した一般化強制スキーム（Relax-to-Force法）が使用される。
   • 空隙率スケーリング： 局所的に保存される場は、体積平均化された密度と運動量として定義され、キネティック・モーメントに空隙率のスケーリングが自然に組み込まれている。
   • 速度評価： マクロな速度と強制項の間の暗黙的な関係を解くために反復アルゴリズムが提案されており、非線形なダルシー・フォルサイマー抵抗と粘性補正を処理する。

3. 境界条件： 著者らは、標準的なウェットノード境界条件（平衡スキーム、非平衡外挿、およびTamm-Mott-Smithスキーム）を多孔質媒体定式化に適応させている。質量連続性を考慮し、強制項の存在と空隙率の変化を反映させた上で、境界における固有流体密度を算出するための特定のプロセスが導出されている。

主な貢献

• 一貫した1ドメイン定式化： 本論文は、表観速度を用いて表現された、完全に一貫した空隙率スケーリングを持つブリンクマン粘性応力テンソルを回収する、初のLBM定式化を提示している。これにより、人工的な界面条件を設けることなく、空間的に変化する空隙率を扱うことが可能となる。
• 圧力および粘性項の補正： 著者らは、不均一媒体において生じる圧力勾配および粘性応力項の不整合を是正するために、特定の強制補正を導出し実装した。これにより、モデルが正しいマクロ極限を回収することを保証している。
• 独立したパラメータ調整： 本モデルは、音速と有効体積粘性を、剪断粘度とは独立して指定することができる。これにより、計算効率（高い音速を許容することによる）と数値的安定性（体積粘度の調整による）が向上する。
• モノリシック・アルゴリズム： 全領域に対して単一のキネティック方程式とソルバーを使用するため、自由流と多孔質領域が混在する複雑な形状の構築が簡素化される。

結果と検証
モデルは、一連のベンチマークと数値的研究を通じて検証された：

• 一貫性解析： 線形時間スペクトル解析により、モデルが正しい有効音速および剪断モードと音響モードの散逸率を回収することが確認された。解析の結果、提案された補正を伴う積形式の平衡分布は、従来の2次平衡分布とは異なり、ガリレイ不変性と等方性を保持することが示された。
• ベンチマーク流：
  • クーエット流： 重力下にある均質多孔質チャネルにおけるシミュレーションでは、様々なレイノルズ数およびダルシー数において、速度および密度プロファイルの両方で参照数値解と極めて良好な一致を示した。
  • ポアズイユ流： モデルは均質チャネルにおける流れを正確に捉えただけでなく、決定的なことに、不均一な空隙率プロファイルを持つチャネルにおいても正確に捉えた。流体が高空隙率領域を優先的に流れる「チャネリング効果」を正しく予測した。
• 数値デモンストレーション：
  • 透過性円柱の周囲の流れ： モデルは透過性円柱の周囲の拘束流をシミュレートし、レイノルズ数とコア空隙率に基づいた5つの明確な流れの領域を特定した。ウェイク（後流）の安定化、剥離した再循環、および遠方後流の不安定性などの現象を捉え、文献と良好な定性的一致を示した。
  • 多孔質障害物を含むリッド駆動キャビティ： 内部に多孔質障害物を持つキャビティをシミュレートし、透過性が渦形成、流れの浸透、および再循環パターンにどのように影響するかを示した。結果は、空隙率が増加するにつれて、固体的な挙動から透過的な挙動への遷移を浮き彫りにした。

意義と主張
著者らは、本研究が不均一な多孔質媒体における輸送をシミュレートするための、堅牢で統一されたフレームワークを提供すると主張している。明示的な界面境界条件を排除し、すべての流体力学項の一貫したスケーリングを保証することで、本モデルは以下のような用途に特に適している：

• 空隙率が空間的に変化する設計された多孔質構造の、初期段階のデバイススケールの設計研究。
• 自由流と透過性体（燃料電池、ろ過システム、多孔質燃焼など）の間の複雑な相互作用のシミュレーション。
• 自由流と多孔質媒体の間の遷移が連続的である、あるいは鋭い界面の解像が困難なシステムの輸送現象の捕捉。

本論文は、現在の研究が2次元の検証に焦点を当てているものの、定式化自体は本質的に多次元的であり、3次元幾何学への拡張が可能であることを結論付けている。今後の課題として、現実的なエンジニアリングシステムにおける結合輸送を完全に扱うための、エネルギー収支（熱伝達）、多成分混合物、および多相流の組み込みが必要であると述べている。