[論文レビュー] A homogenized limit for the 2D Euler equations in a perforated domain

Matthieu Hillairet, Christophe Lacave+1 more|arXiv (Cornell University)|Jul 9, 2019

Advanced Mathematical Modeling in Engineering参考文献 28被引用数 2

ひとこと要約

本稿は、穴の大きさ a と穴間隔 ˜d が両方 0 に近づく臨界的状態 ˜d/a → ¯k > 0 のもとで、穴あき領域における2次元Euler方程式の最初の均質化極限を確立した。この状態では、多孔質媒体の体積分率が正で小さい値をとる。反射法を用いて、発散・回転問題に対する精密な推定を導出し、回転速度の結合を記述する有効行列項を含む均質化Euler型系を導いた。これは、従来の希薄（˜d/a → ∞）または稠密（˜d/a → 0）極限ではカバーされなかった、非ゼロで小さな多孔質体積分率の状態に有効である。

ABSTRACT

We study the motion of an ideal incompressible fluid in a perforated domain. The porous medium is composed of inclusions of size $a$ separated by distances $ ilde d$ and the fluid fills the exterior. We analyse the asymptotic behavior of the fluid when $(a, ilde d) o (0,0)$. If the inclusions are distributed on the unit square, this issue is studied recently when $\frac{ ilde d}a$ tends to zero or infinity, leaving aside the critical case where the volume fraction of the porous medium is below its possible maximal value but non-zero. In this paper, we provide the first result in this regime. In contrast with former results, we obtain an Euler type equation where a homogenized term appears in the elliptic problem relating the velocity and the vorticity. Our analysis is based on the so-called method of reflections whose convergence provides novel estimates on the solutions to the div-curl problem which is involved in the 2D-Euler equations.

研究の動機と目的

• 穴の大きさ a と穴間隔 ˜d が両方 0 に近づくとき、2次元非圧縮性Euler流れが穴あき領域で示す漸近的挙動を解析すること。このとき ˜d/a → ¯k > 0 である。
• 従来の研究が希薄（˜d/a → ∞）および稠密（˜d/a → 0）な状態にのみ焦点を当てていたため、その理論的ギャップを埋めること。
• 多孔質媒体が小さな非ゼロ体積分率を有する臨界的状態におけるEuler方程式の均質化系を導出すること。
• 反射法を用いて、穴あき領域における発散・回転問題に対する新たな推定を確立し、これにより非ゼロで小さな多孔質体積分率の範囲で均質化極限への収束を可能にする。

提案手法

• 解析は2段階に分けられる。まず、速度と回転度の関係をストリーム関数 ψN を介して記述する穴あき領域 FN における発散・回転問題に対する精密な一様推定を導出する。
• 反射法を用いて近似解を構成し、真の解 ψN と均質化解 ψc 間の誤差を評価する。
• 均質化問題は R² における div[(I₂ + kMK)∇ψc] = f として定式化される。ここで MK は穴の形状 K に依存する行列であり、k は体積分率の極限値である。
• 主な技術的革新は、重み付きノルムと回転度勾配の対数Sobolev型評価を用いて、真の速度 uN と均質化速度 uc 間の差に対する非等方的 L∞ 評価を導出することにある。
• エネルギー推定とvorticity輸送方程式にグロンウォール型の議論を適用することで、vorticityおよび速度場の収束を確立する。
• 解の存在時間の均一な延長を保証するため、ブートストラップ法を用い、収束が固定時間 T まで成り立つようにする。

実験結果

リサーチクエスチョン

• RQ1 多孔質媒体が小さな非ゼロ体積分率を有する状態、すなわち臨界的状態 ˜d/a → ¯k > 0 における2次元Euler方程式の均質化極限は何か？
• RQ2 非ゼロの多孔質媒体が、特に有効粘性係数や非等方性の観点から、Euler設定における巨視的流体力学にどのように影響を与えるか？
• RQ3 反射法が、臨界的状態における穴あき領域内の真の解と均質化解の間の誤差を十分に制御できるか？
• RQ4 均質化極限における流体運動を記述する有効方程式の構造は何か？また、古典的Euler方程式やDarcy型法則とはどのように異なるか？
• RQ5 穴あき領域におけるvorticityおよび速度場は、それぞれの均質化対応物へどのように収束するか？また、その収束に対する一様な境界は何か？

主な発見

• 臨界的状態における2次元Euler方程式の均質化極限は、R² における修正された楕円型問題で記述される：div[(I₂ + kMK)∇ψc] = f。ここで MK は穴の形状 K に依存する行列である。
• 単位円板の場合、MK = 2I₂ であるため、有効作用素は div[3I₂∇ψc] = f となり、古典的ラプラシアンに対して顕著な非等方的補正が生じる。
• 速度場 uN から均質化場 uc への収束は、L∞(OT) において、F(N,k) = (a/d)^{3−η} + ||µ−k||_{W^{−1,p}}^{(1−η)(p+2)/p} + ||µ−k||_{W^{−1,p}}^{1/2} + ||k||_{L∞}^2 で抑えられる収束速度を持つ。
• 対数グロンウォール推定により、vorticity勾配が時間に一様に有界であることが保証され、これはEuler方程式の非線形項を制御するために不可欠である。
• 穴あき領域における流体粒子の軌道 XN(t,x) は、[0,T] 上で均質化軌道 Xc(t,x) へ一様に収束し、誤差は O(F(N,k)) である。
• 初期vorticityがcompactlyな台を持ち、L∞ 範囲に有界である場合、収束は固定時間 T まで一様に成り立ち、体積分率 k が小さく、FV(ε₀) に属する限り成立する。