Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media

本論文は、不均質多孔質媒体内の非圧縮性二相流について、ポアスケールのナビエ - ストークス・カahn・ヒルディル方程式から体積平均法を用いて厳密に導出したマクロモデルを提示し、閉鎖問題の解決による輸送係数の評価や濡れ性の形式的な取り込みを通じて、標準的な経験則に基づくモデルとの理論的差異を明確にし、数値シミュレーションによりその有効性を検証したものである。

要約

この論文は、**「多孔質媒体（スポンジのようなもの）の中を、混ざり合わない 2 つの液体が流れる様子」**を、ミクロな視点からマクロな視点へとつなぐ新しい数学モデルを開発したという研究です。

専門用語を並べると難しく聞こえますが、実は**「スポンジの中での油と水の戦い」**をより正確に予測するための「新しい地図（モデル）」を作ったお話です。

以下に、日常の例えを使ってわかりやすく解説します。

---

1. 研究の背景：なぜこんな難しいことをするの？

想像してください。地下の岩層（多孔質媒体）に、石油を回収しようとして水を注入するとします。

• ミクロな視点（ pore scale）： 岩の隙間（毛穴）のサイズで見ると、油と水は入り組んで複雑に動き、互いに押し合い、表面張力でくっついたり離れたりしています。これをすべてコンピュータでシミュレーションするには、「砂粒一つ一つ」を計算する必要があるため、時間がかかりすぎて現実的ではありません。
• マクロな視点（ Darcy scale）： 一方、油田全体のような大きなスケールでは、「平均してどれくらい流れているか」を知りたいだけです。

これまでの一般的なモデル（ダルシーの法則など）は、実験データに基づいた「経験則」で、「スポンジの性質」や「液体が壁にどうくっつくか（濡れ性）」を正確に反映しきれていないという問題がありました。

2. この研究のすごいところ：「アップスケーリング」という魔法

この研究チームは、「アップスケーリング（拡大縮小）」という手法を使いました。
これは、「ミクロな物理法則（ナヴィエ - ストークス方程式と Cahn-Hilliard 方程式）」を、数学的に丁寧に「平均化」して、マクロな世界でも使える新しい方程式を導き出したというものです。

重要なポイント：「濡れ性」の正体を解明

これまでのモデルでは、「水がスポンジにどうくっつくか（濡れ性）」は、実験で決めた「魔法の係数」で適当に調整していました。
しかし、この研究では、「壁との相互作用（濡れ性）」を、化学ポテンシャル（液体の「流れたい気持ち」を表す値）の中に、理論的に組み込みました。

• 例え話：
  • これまでのモデル： 「スポンジは水が好きだから、水はよく流れる」という**「経験則」**で流れる量を推測していた。
  • 今回のモデル： 「スポンジの表面が水分子をどう引き寄せるか」という**「分子レベルの仕組み」**を計算式の中に直接組み込んで、流れを予測する。

これにより、実験データに頼らずとも、スポンジの性質や液体の性質が変われば、自動的に正しい流れ方を予測できる「理にかなったモデル」が完成しました。

3. 具体的な成果：どんなテストをした？

開発した新しいモデルが本当に使えるか、いくつかのテストを行いました。

1. スポンジの中の単純な流れ（ポアズイユ流れ）：
   • 既知の答えがある単純な流れをシミュレーションし、モデルの計算結果が理論値と一致することを確認しました。「新しい地図は、既存の道も正しく描けている」ことを証明しました。
2. 油と水の境界（バクレー・レヴェレット問題）：
   • 油と水が入れ替わる際の「境界線」がどう動くかを計算しました。実験や従来の理論とよく一致し、**「境界がどう鋭く移動するか」**を捉えることができました。
3. 粘性指状現象（Viscous Fingering）：
   • 粘度の違う液体が混ざり合う時、指のような細い流れが生まれる現象です。
   • ここが面白い： 「スポンジが水を好む（親水性）」場合と「油を好む（疎水性）」場合で、その「指」の形がどう変わるかをシミュレーションしました。
   • 結果： 親水性だと、指が分かれて広がりやすくなり、疎水性だとまとまりやすくなるなど、**「濡れ性の違いが、流れの形を劇的に変える」**ことをモデルが再現できました。

4. まとめ：この研究は何を意味するのか？

この論文は、**「ミクロな世界の複雑なルールを、マクロな世界の予測にどう活かすか」**という、長年の課題に挑んだものです。

• 従来の方法： 「実験して、あてはめる」→ 条件が変わると使えなくなる。
• 今回の方法： 「物理法則から導き出す」→ 条件が変わっても、原理的に正しい答えを出せる。

「スポンジの中の油と水の動き」を、より正確に、より深く理解するための新しい「共通言語」を作ったと言えます。

このモデルは、石油の回収効率を上げる、地下に二酸化炭素を貯留する、あるいは地下水の浄化など、様々な実社会の問題を解決する際の、より信頼性の高い設計図として役立つでしょう。

技術概要

以下は、提示された論文「不均質多孔質媒体内の非圧縮性多相流のための Navier-Stokes-Cahn-Hilliard モデルのアップスケーリング（Upscaling the Navier-Stokes-Cahn-Hilliard model for incompressible multiphase flow in inhomogeneous porous media）」に関する詳細な技術的サマリーです。

1. 研究の背景と課題 (Problem)

多孔質媒体における二相流（例：地下水修復、石油回収、CO2 隔離など）は、工学的に極めて重要ですが、そのモデリングは困難を伴います。

• 課題: 多孔質媒体の幾何学的構造は非常に複雑であり、ポアスケール（微小スケール）での直接数値シミュレーション（DNS）は物理メカニズムの理解には有用ですが、時間・空間スケールが広範囲にわたるため、工学的なプロセス最適化には適用が困難です。
• 既存手法の限界: 従来のダルシー法則やその修正版（Brinkman-Forchheimer 方程式など）は、経験則に基づいており、相飽和度のみをパラメータとする相対透水性や毛管圧力の関数を用いています。しかし、これらは濡れ性の履歴効果（浸潤・排水の違い）や、ポアスケールの界面張力・濡れ境界条件の微視的詳細をマクロな式に体系的に反映できておらず、物理的な裏付けが不足しています。
• 目的: ポアスケールの物理法則（Navier-Stokes 方程式と Cahn-Hilliard 方程式）から、体積平均法を用いて厳密にマクロな（ダルシースケールの）モデルを導出すること。特に、濡れ性を化学ポテンシャルに形式的に組み込むことを目指します。

2. 手法 (Methodology)

本研究では、**体積平均法（Volume Averaging Method）**を採用し、ポアスケールの支配方程式をマクロスケールにアップスケーリングしました。

• ポアスケールモデル:
  • 拡散界面理論（Diffuse-interface theory）に基づき、秩序変数 $\phi$（0 から 1 の間）を用いて二相を記述。
  • 流れの支配方程式：非圧縮性 Navier-Stokes 方程式。
  • 界面の進化方程式：Cahn-Hilliard 方程式。
  • 自由エネルギー汎関数に、バルク自由エネルギー、界面エネルギー、および壁面自由エネルギー（濡れ性を記述）を含める。
• アップスケーリング過程:
  • 代表要素体積（REV）内で空間積分を行い、平均量と空間偏差（deviation）に分解（Gray's decomposition）する。
  • 空間平均定理と時間平均定理を適用し、ポアスケール方程式から平均化された方程式を導出。
  • クロージャ問題（Closure Problem）: 平均化によって生じる未決定項（空間偏差に起因する項）を、REV 内の局所的な問題として定義し、閉じるための係数（輸送係数など）を導出。
• 重要な導出:
  • 濡れ境界条件を直接平均化学ポテンシャルに組み込むことで、毛管力をマクロな運動量方程式の項として形式的に導出。
  • 単一流体モデル（Single-fluid model）として定式化し、従来の二流体モデル（各相に独立した速度と圧力を持つ）とは区別されるアプローチを確立。

3. 主要な貢献 (Key Contributions)

1. 厳密な単一流体マクロモデルの導出:
   経験則に頼らず、ポアスケールの NS-CH 方程式から体積平均法を用いて、熱力学的に整合性のあるマクロな二相流モデルを厳密に導出した。
2. 濡れ性の形式的な取り込み:
   従来の経験的な毛管圧力関数の代わりに、平均化学ポテンシャルの中に濡れ境界条件（接触角 $\theta$ や壁面自由エネルギー）を明示的に含めることで、濡れ性がマクロな流れに与える影響を理論的に記述可能にした。
3. 未閉項の体系的な処理:
   空間偏差に起因する未閉項（慣性項、粘性項、ドラッグ力など）を、REV 内の局所問題として定義し、有効粘度、相対透水性、分散テンソルなどの係数としてモデル化。
4. 格子ボルツマン法（LBM）による実装と検証:
   導出された平均化方程式を解くための新しい格子ボルツマン法を開発し、数値シミュレーションによりモデルの有効性を検証。

4. 数値結果と検証 (Results)

開発されたモデルとアルゴリズムの性能を、以下の 4 つのケースで検証しました。

1. 均質多孔質媒体内のポアズイユ流:
   • 位相場方程式を無視し、単相流として解析解と比較。
   • 異なるダルシー数（Da）において、数値解と解析解が良好に一致し、非線形抵抗項の挙動を正しく再現することを確認。
2. 流体 - 多孔質複合チャネル内のポアズイユ流:
   • 自由流体領域と多孔質領域の界面での流れをシミュレーション。
   • 界面での速度連続性と応力ジャンプ条件を正しく再現し、不均質多孔質媒体における流れの予測能力を確認。
3. Buckley-Leverett 問題（一相流置換）:
   • 二相置換におけるショックフロント（急激な飽和度変化）の伝播を解析。
   • 粘度比を変化させた場合のフロント飽和度を、理論解（Welge 接線法）と比較。最大相対誤差は 5.6% 以内であり、モデルが二相流の基本的なダイナミクスを捉えうることを示した。
4. 多孔質媒体内の粘性指状現象（Viscous Fingering）:
   • 濡れ性（接触角）が指状不安定性に与える影響を調査。
   • 親水性（接触角が小さい）の場合、指状構造が分散しやすく、指の形成が抑制されることを確認。これは実験結果と定性的に一致し、濡れ性の影響をモデルが適切に反映していることを示した。

5. 意義と結論 (Significance)

• 理論的統合: ポアスケールの微視的詳細（特に濡れ性と界面張力）とダルシースケールのマクロな記述を架橋する、統一された理論的枠組みを提供しました。
• 経験則からの脱却: 相対透水性や毛管圧力を単なる経験関数として扱うのではなく、微視的パラメータ（接触角、界面張力、孔隙構造）から導出可能な物理量として扱う道筋を開きました。
• 将来展望: 本研究で導出されたクロージャ係数（$a_\phi, b_\phi, c_\phi$ など）は、より複雑な幾何構造に対して厳密に評価する必要があります。また、従来の Leverett J 関数などの経験的モデルとの関係性を明確化し、実用的な多孔質媒体流のシミュレーションへの応用が期待されます。

総じて、この研究は多孔質媒体内の多相流モデリングにおいて、経験則に依存しない、物理的に厳密で拡張性の高い新しいアプローチを確立した点で画期的です。