[論文レビュー] Theoretical derivation of Darcy's law for fluid flow in thin porous media

Francisco J. Suárez‐Grau|arXiv (Cornell University)|Jan 1, 2020

Advanced Mathematical Modeling in Engineering参考文献 22被引用数 9

ひとこと要約

本稿は、膜厚さεと微細構造サイズaεの2つの小さなパラメータを伴う薄い多孔質媒体における非圧縮性ニュートン流体の流れに対して、Darcyの法則の厳密な理論的導出を提供する。条件aε ≪ εのもとで、新規の2パラメータアンフォールド法を用いて、3次元局所セル問題から導かれる透過率テンソルを有する2次元Darcyの法則を導出しており、従来の薄い媒体における均質化モデルとは異なる。

ABSTRACT

In this paper we study stationary incompressible Newtonian fluid flow in a thin porous media. The media under consideration is a bounded perforated $3D$ domain confined between two parallel plates. The description of the domain includes two small parameters: $\\varepsilon$ representing the distance between pates and $a_\\varepsilon$ connected to the microstructure of the domain such that $a_\\varepsilon\\ll \\varepsilon$. We consider the classical setting of perforated media, i.e. $a_\\varepsilon$-periodically distributed solid (not connected) obstacles of size $a_\\varepsilon$. The goal of this paper is to introduce a version of the unfolding method, depending on both parameters $\\varepsilon$ and $a_\\varepsilon$, and then to derive the corresponding $2D$ Darcy's law.

研究の動機と目的

• 穿孔多孔質媒体における古典的均質化理論を、膜厚さεと微細構造サイズaεという2つの小さなパラメータを有する薄い多孔質媒体に拡張すること。
• 膜厚さεと微細構造サイズaεの両方に依存する2パラメータアンフォールド法を構築し、薄い多孔質領域における流体の流れを解析すること。
• aε ≪ εの領域において、2次元Darcyの法則を厳密に導出すること。
• εとaεが同時に0に近づく際の均質化プロセスの収束を確立し、導出されたDarcy法則の数学的妥当性を保証すること。

提案手法

• 領域は穿孔された3次元シリンダQε = ω × (0, ε)としてモデル化され、周期的に配置された固体障害物（サイズaε）を有し、limε→0 aε/ε = 0 を満たす。
• 2パラメータアンフォールド作用素が導入され、周期的アンフォールド法を一般化して、膜厚さε（縦方向）と微細構造サイズaε（微細スケール）を同時に取り扱う。
• ストークス方程式が穿孔領域Ωε = Qε ∩ Tεに定式化され、障害物上ではスリップなし境界条件、上下面ではフラックスなし条件が適用される。
• 弱形式が導かれ、アンフォールド法を用いて極限に移行し、積空間ω × Yfにおける均質化問題が得られる。ここでYfは単位セルYの流体部である。
• 3次元流体領域Yfにおける局所セル問題を解き、透過率テンソルKを決定する。Kij = ∫Yf wⁱⱼ(y) dyであり、wⁱはYfにおけるストークス型問題の解である。
• 圧力pは2次元Darcy方程式：−∇x′·(K∇x′p) = f′ in ωを満たし、積分と極限の同定を通じて導出される。

実験結果

リサーチクエスチョン

• RQ1 薄い多孔質媒体における流体の流れの漸近的挙動は、膜厚さεと微細構造サイズaεの相対的スケーリングにどのように依存するか？
• RQ2 aε ≪ εの下で、特に微細構造が膜厚さ方向に整列していない場合に、Darcyの法則を導出する適切な数学的枠組みは何か？
• RQ3 横方向（ε）および微細スケール（aε）の両方の均質化を同時に処理できる2パラメータアンフォールド法を構築可能か？
• RQ4 薄い領域における有効Darcy法の形は何か？また、均質的に薄い媒体で導出された古典的2次元Darcy法とはどのように異なるか？
• RQ5 aε ≪ εの下で、透過率テンソルKは3次元局所問題からどのように導出されるか？また、その数学的性質は何か？

主な発見

• 本稿は、微細構造が膜厚さよりもはるかに小さいaε ≪ εの領域において、薄い多孔質媒体における2次元Darcyの法則を厳密に導出した。
• 透過率テンソルKは、単位セルYの流体部Yfにおける3次元局所セル問題から得られ、Kij = ∫Yf wⁱⱼ(y) dyであり、wⁱはYfにおけるストークス型問題の解である。
• 圧力pは2次元Darcy方程式：−∇x′·(K∇x′p) = f′ in ωを満たし、Kは対称かつ正定値であるため、解の存在と一意性が保証される。
• 速度場uはu(x′) = K(f′ − ∇x′p)として表され、アンフォールド法および弱収束の議論を通じて均質化極限から導出される。
• 列(aε⁻²ũε, pε)が均質化解(u, p)に収束することを証明し、導出されたDarcy法則の数学的妥当性を確立した。
• 導出されたDarcy法則は、従来の均質的薄い媒体モデル（例：Fabriciusら[9]）とは異なり、aε ≪ εのスケーリングの下で3次元局所問題を保持するため、他の薄い媒体領域とは異なる。