Brinkman's term for the steady-state Navier-Stokes equations with prescribed flux rate or pressure drop in perforated pipes

技術要約：混合境界条件を持つ穿孔管内における定常ナビエ・ストークス方程式のホモジナイゼーション

問題設定
本研究は、多数の小さな固体粒子を含む歪んだ管の中における、粘性非圧縮性流体の定常運動の漸近的挙動を調査するものである。領域 $\Omega_\varepsilon$ は、直径 $\varepsilon^3$ の粒子 $N(\varepsilon)$ 個を互いに距離 $\varepsilon$ で配置して取り除いた、有界で許容可能な領域 $\Omega \subset \mathbb{R}^3$ の穿孔された部分集合である。この構成は、ホモジナイゼーション理論における臨界領域（$\alpha = 3$）に対応しており、粒子のサイズが粒子間距離の3乗のスケールで変化している。

流体流動は、混合境界条件を持つ定常ナビエ・ストークス方程式によって支配される。具体的には：

1. 物理的な管壁（$\Gamma_W$）および粒子の表面（$\partial K_n^\varepsilon$）に対して、**滑りなし条件（no-slip conditions）**が課される。ここで、粒子は一定の速度 $\mu_n^\varepsilon$ で移動している場合がある。
2. 入口（$\Gamma_I$）および出口（$\Gamma_O$）には、ベルヌーイ圧力条件が適用される。これらは、速度の法線微分とベルヌーイ圧力（$\Phi = p + \frac{1}{2}|u|^2$）を含む。
3. 以下の2つの異なるシナリオが分析される：
   • 正味フラックス（Net Flux）が規定されている場合： 出口を通る全フラックス $F$ が固定されており、圧力降下は解によって決定される未知の定数となる。
   • 圧力降下が規定されている場合： 入口と出口の間の圧力降下が固定されており、その結果として生じるフラックスは未知の量となる。

主要な数学的課題は、外部力、フラックス率、または圧力降下に微小性の仮定を置くことなく、$\varepsilon \to 0$ の極限における有効（ホモジナイズされた）方程式を導出することである。これは、混合境界値問題においては、追加の制約なしではナビエ・ストークス方程式の標準的なエネルギー評価がしばしば成立しないという事実によって困難を伴う。

手法
著者らは、タルタールの一般的な戦略に従い、ホモジナイゼーション理論の枠組みの中でエネルギー法を用いている。証明は主に3つの段階で進行する：

1. 一様有界性（$\varepsilon$ に依存しない評価）：
   解 $(u_\varepsilon, \Phi_\varepsilon)$ に対する一様有界性を確立することが、最も重要かつ斬新なステップである。混合境界条件や未知の圧力降下の存在により、単純なエネルギー評価のみでは不十分であるため、標準的なケースとは異なり、著者らは単に解自体で運動方程式をテストすることに依存できない。

   • 代わりに、定常オイラー方程式の解に対するベルヌーイの法則に基づく背理法を利用する。
   • 解のノルムが $\varepsilon \to 0$ で発散すると仮定し、それらを正規化し、その極限がオイラー方程式を満たすことを示す。
   • 管壁における極限ベルヌーイ圧力の境界挙動を分析することで（局所的な正則性の結果と入口/出口の特定の幾何学的構造を用いて）、極限圧力が壁上で一定であることを示す。これにより、フラックスまたは圧力降下の条件との矛盾が生じ、解の列の一様有界性が証明される。
   • 規定フラックスの問題については、境界分析を容易にするために、粒子が入口および出口のシリンダー付近に蓄積しないという追加の幾何学的仮定がなされている。

2. 補正項（Correctors）および制限作用素（Restriction Operators）の構成：
   穿孔領域を扱うために、著者らは粒子の外部および粒子を取り囲む環状領域におけるストークス補正項 $(\omega_{n,1}^\varepsilon, \omega_{n,2}^\varepsilon)$ を構成する。これらの補正項は、制限作用素 $R_\varepsilon$ を定義するために使用される。

   • 粒子群の相対的なキャパシタンスが非ゼロとなる劣臨界領域における先行研究では、スカラーのキャパシタリティ補正項で十分であったが、臨界領域では、粒子群の非ゼロの相対キャパシタンスを考慮するために、これらのベクトル値のストークス補正項が必要となる。
   • 制限作用素 $R_\varepsilon(\phi)$ は、テスト関数を、粒子付近でゼロの境界条件を満たすように修正しつつ、発散フリー（divergence-free）の性質を保持しながら、ホモジナイズされた領域 $\Omega$ から穿孔された領域 $\Omega_\varepsilon$ へ写像する。

3. 極限への移行：
   一様有界性と制限作用素を用いて、ナビエ・ストークス方程式の弱形式において極限を通過させる。

   • 非線形対流項は、一様有界性とコンパクトな埋め込みによって保証される $L^4$ における強収束を介して収束する。
   • 流体と粒子の相互作用は、極限方程式における追加の項を生成する。この項は、ストークス抵抗行列 $G_n$ によって特徴付けられる、粒子によって及ぼされる抗力の和の極限的な挙動から生じる。

主な貢献と結果
本論文は、系のホモジナイズされた極限を記述する2つの主要な定理（定理2.4および2.5）を確立している：

• ブリンマン項（Brinkman Term）の出現： ホモジナイズされた領域 $\Omega$ における有効方程式は、ブリンマン項によって拡張されたナビエ・ストークス方程式となる：
  $$-\Delta u + (\nabla \times u) \times u + \nabla \Phi + Gu - J = f$$
  ここで：

  • $Gu$は、流速$u$に比例する**摩擦項**を表し、粒子の集団的な抗力（$G$ はスケーリングされたストークス抵抗行列の和の極限）から生じる。
  • $J$ は、移動する粒子自体の極限フラックスから生じるソース項である。
  • 入口および出口における境界条件は、ベルヌーイ圧力を含む混合型として維持される。

• 非ゼロのデータを持つ臨界領域への拡張： これらの結果は、劣臨界領域に関する先行研究（特に[59]）の結果を、臨界領域（$\alpha=3$）へと拡張するものである。決定的なことに、これらの結果は、外部力、フラックス率、または圧力降度に対する微小性の仮定なしに成立する。これは、一意性や有界性を確保するために微小データへの仮定に依存することが多い、多くのナビエ・ストークス・ホモジナイゼーションの結果と比較して、重要な進展である。

• 移動する粒子の取り扱い： 本フレームワークは、粒子の運動エネルギーが一様に有界であるという条件下（仮定2.11）において、非ゼロの一定速度 $\mu_n^\varepsilon$ を持つ粒子の扱いを許容する。これにより、固定された粒子モデルと比較して、より現実的な粒子・流体相互作用のモデリングが可能となる。

意義
本論文の最も重要な貢献は、無制限のデータの下での、定常ナビエ・ストークス方程式におけるブリンマンの法則の厳密な導出であると著者らは主張している。ベルヌーイの法則を用いた矛盾論法による一様有界性の導出は、特に、議論が $\varepsilon$ レベルで成立することを保証するために、穿孔の複雑な構造を考慮しなければならないという点で、非常に繊細な問題である。

粒子に対して消滅するテスト関数を修正することなく、真の制限作用素を構築し、非ゼロの粒子速度を持つ多分散の設定における圧力項の収束を扱うことにより、本研究は、粒子密度がブリンマン型の抗力を誘発するほど高く、かつダルシーの法則に至るほど高くはない、穿孔管内の流体流をモデリングするための堅牢な数学的基礎を提供する。これらの結果は、流体が固体の配列を通過する際の流体現象が基本的現象であるヘモダイナミクス、石油工学、航空力学などの分野に直接適用可能である。