Navier-Stokes-Cahn-Hilliard system in a $3$D perforated domain with free slip and source term: Existence and homogenization

この論文「Navier-Stokes-Cahn-Hilliard system in a 3D perforated domain with free slip and source term: Existence and homogenization（自由滑り境界条件と源項を有する 3 次元穿孔領域における Navier-Stokes-Cahn-Hilliard 系：存在性と均質化）」の技術的概要を日本語でまとめます。

1. 問題設定 (Problem Setting)

本論文は、周期的に穿孔された多孔質媒質（3 次元）内の二成分非圧縮性混合流体の拡散界面モデルを解析対象としています。

• 支配方程式: 時間依存型の Navier-Stokes-Cahn-Hilliard (NSCH) 系。
  • 運動量方程式: 非圧縮性 Navier-Stokes 方程式に、Korteweg 型のカピラリ力項（$\lambda_\varepsilon \phi_\varepsilon \nabla \mu_\varepsilon$）と非保存的な源項 $G(\phi_\varepsilon)$ を含む Cahn-Hilliard 方程式が結合されています。
  • 相場方程式: 拡散界面の進化を記述する Cahn-Hilliard 型方程式。
• 領域: 固体インクルージョン（障壁）が周期的に配置された穿孔流体領域 $\Omega^\varepsilon_p$。
• 境界条件:
  • 外境界 $\partial \Omega$: ノー・スリップ条件（$u=0$）。
  • 固体インクルージョン表面 $\Gamma^\varepsilon_s$: **自由滑り（Free slip）**条件。法線方向速度はゼロ ($u \cdot n = 0$) だが、接線方向の剪断応力がゼロ ($[A^\varepsilon D(u)n] \cdot \tau = 0$) となります。これは Navier 滑りの一種です。
• パラメータ:
  • $\varepsilon$: 微視的な周期スケール。
  • $\lambda_\varepsilon$: カピラリ強度。$\varepsilon \to 0$ において極限値 $\lambda \in [0, +\infty)$ に収束します。
  • 外力 $g_\varepsilon$: カピラリ強度 $\sqrt{\lambda_\varepsilon}$ に比例してスケーリングされます。

2. 手法とアプローチ (Methodology)

解析は大きく 2 つの段階に分けて行われます。

A. 微視的解の存在性と一様な評価 (Existence and A Priori Estimates)

固定された $\varepsilon > 0$ に対して、弱解の存在と $\varepsilon$ 及び $\lambda_\varepsilon$ に対して一様な事前評価（a priori estimates）を導出します。

• 関数空間と拡張: 穿孔領域 $\Omega^\varepsilon_p$ 上で定義された関数を、固体領域を含めた全体領域 $\Omega$ へ拡張するための線形拡張作用素（Extension operators）を構成します。特に、質量保存則が成り立たない（源項 $G$ の存在）ため、相場 $\phi$ の空間平均が時間依存する点に注意し、双対レベルでの拡張を慎重に扱います。
• Bogovski˘ı 作用素: 穿孔領域における発散の逆作用素（Bogovski˘ı 作用素）を構成し、そのノルムが $\varepsilon$ に依存しないことを示します。これにより、圧力の存在と評価を確立します。
• エネルギー法: 自由エネルギー（運動エネルギー＋自由エネルギー密度）の時間変化を解析し、エネルギー散逸則を導出します。源項 $G$ による質量非保存を扱い、相場の空間平均の評価を別途行うことで、一様な評価を得ています。
• Galerkin 近似: Stokes 型の固有関数を基底として用いた Galerkin 近似により、弱解の存在を証明します。

B. 均質化 (Homogenization)

$\varepsilon \to 0$ の極限において、巨視的な有効モデルを導出します。

• 手法: 穿孔領域に適応された**周期的展開作用素（Periodic Unfolding Operator）**を用いた 2 尺度収束法（Two-scale convergence）を採用しています。
• ケース分け: カピラリ強度の極限 $\lambda$ の値に応じて、2 つの異なる巨視的モデルが得られます。
  1. $\lambda = 0$ (カピラリ消滅領域): 表面張力の効果が巨視的に無視できる場合。
  2. $\lambda \in (0, +\infty)$ (バランス領域): 表面張力と粘性、慣性力が拮抗する場合。

3. 主要な結果 (Key Results)

1. 弱解の存在と一様評価

穿孔領域における NSCH 系に対して、混合境界条件（外境界はノー・スリップ、固体表面は自由滑り）のもとで弱解の存在が証明されました。また、解のノルムが $\varepsilon$ と $\lambda_\varepsilon$ に対して一様に有界であることが示されました。これは、源項 $G$ が存在し質量保存則が破れる状況下での重要な結果です。

2. 巨視的有効モデルの導出

極限 $\varepsilon \to 0$ において、以下の 2 つの異なる巨視的方程式系が得られました。

• ケース (i): $\lambda = 0$ の場合

  • 得られる系は Stokes-Cahn-Hilliard (SCH) 型です。
  • 巨視的な対流項 $(u \cdot \nabla)u$ および移流項 $u \cdot \nabla \phi$ は消失します。
  • 運動量方程式は Stokes 方程式となり、カピラリ力のみが駆動力として残ります。
  • 有効テンソル $A^{hom}, B^{hom}, C^{hom}$ が細胞問題（Cell problems）を通じて定義されます。

• ケース (ii): $\lambda \in (0, +\infty)$ の場合

  • 得られる系は Navier-Stokes-Cahn-Hilliard (NSCH) 型です。
  • 巨視的な対流項 $(u \cdot \nabla)u$ と相場の移流項 $u \cdot \nabla \phi$ が現れます。
  • これらの項は $\sqrt{\lambda}$ に比例する係数を持ち、カピラリ力が流体の流れを駆動するメカニズムが巨視的に再現されます。
  • 有効テンソルは $\lambda=0$ の場合と同じですが、方程式の構造が異なります。

3. エネルギーの収束

微視的な全エネルギーが、巨視的な有効エネルギー汎関数に収束することが示されました。この巨視的エネルギーも、微視的な系と同様のエネルギー散逸則を満たします。

4. 論文の意義と新規性 (Significance and Novelty)

• 境界条件の一般化: 既存の多くの研究が「ノー・スリップ」条件を仮定しているのに対し、本論文は多孔質媒質の固体表面における「自由滑り（摩擦なしの滑り）」条件を NSCH 系に組み込み、その均質化を初めて厳密に扱いました。
• 質量非保存の扱い: 源項 $G(\phi)$ を含むことで質量保存則が成り立たないケースを扱い、相場の空間平均が時間変化する状況下でも、一様な評価と均質化が可能であることを示しました。
• カピラリ強度の役割の明確化: カピラリ強度 $\lambda$ の極限値によって、巨視的な流れの性質（Stokes 型か Navier-Stokes 型か）が根本的に変化することを数学的に証明しました。これは、界面張力が巨視的な輸送現象にどの程度寄与するかを定量的に評価する枠組みを提供します。
• 変数粘性と穿孔領域: 微視的な幾何構造に依存する変数粘性テンソルと、周期的穿孔領域という複雑な設定において、厳密な数学的解析（存在証明と均質化）を完結させた点で、多孔質媒質内の多相流の理論的基盤を強化しています。

結論

本論文は、多孔質媒質内の複雑な二相流を記述する NSCH 系について、自由滑り境界条件と非保存的な源項を考慮した厳密な数学的解析を行いました。その結果、微視的なパラメータ（特にカピラリ強度）に応じて、巨視的な流れが Stokes 型から Navier-Stokes 型へと遷移する現象を明らかにし、実用的な多孔質媒質モデルの基礎理論を確立しました。