[論文レビュー] Homogenization of reaction-diffusion equations in fractured porous media

Hermann Douanla, Jean Louis Woukeng|arXiv (Cornell University)|Jun 27, 2015

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

ひとこと要約

本稿では、周期的なマイクロ構造を有する破壊的多孔質媒体における反応項が大きい反応拡散方程式に対して、マルチスケール収束を用いた均質化フレームワークを構築する。解の収束をスケールパラメータ ε → 0 として、固定領域におけるマクロスケールの対流拡散反応方程式が得られ、孔と亀裂の周期的マイクロ構造が効果的輸送挙動を捉えている。

ABSTRACT

The paper deals with the homogenization of reaction-diffusion equations with large reaction terms in a multi-scale porous medium. We assume that the fractures and pores are equidistributed and that the coefficients of the equations are periodic. Using the multi-scale convergence method, we derive a homogenization result whose limit problem is defined on a fixed domain and is of convection-diffusion-reaction type.

研究の動機と目的

• 等分布する孔と亀裂を有するマルチスケールの破壊的多孔質媒体における拡散および反応プロセスのモデル化を目的とする。
• 複雑なマルチスケール幾何構造における数値シミュレーションの課題に対処するため、有効なマクロスケールモデルを導出することを目的とする。
• 微小スケールパラメータ ε → 0 における反応項が大きい反応拡散方程式の漸近的挙動を分析することを目的とする。
• スケールパラメータ ε → 0 における ε-解 u_ε が固定領域上に定義された極限問題に収束することを確立することを目的とする。

提案手法

• 周期的で破壊的多孔質媒体（マトリックス、孔、亀裂の相が明確に区別される）において、解の列を解析するため、マルチスケール収束法を用いる。
• 媒質を周期的複合体 Y = Y_m ∪ Y_c としてモデル化し、Y_m には孔（Z_s）と固体マトリックス（Z_p）が含まれ、Y_c は亀裂を表す。
• 拡散係数 A(y, τ) および反応項 g(y, τ, u) に対して周期性および一様強楕円性の仮定を適用し、有界性と正則性を保証する。
• 解 u_ε に対して二重スケール展開を導出し、y = x/ε および τ = t/ε² を用いて、u_0(x,t) + εu_1(x,t,y,τ) + ε²u_2(x,t,y,τ) + ... の形に分解する。
• セル Y_m × T 内で補正問題を解き、有効係数（マクロスケールの拡散テンソル Â(x,t) および対流項 L₁, L₂, L₃）を決定する。
• 極限問題として、マクロスケール領域 Ω × (0,T) における対流拡散反応方程式を確立し、セル問題から導かれた有効係数を用いる。

実験結果

リサーチクエスチョン

• RQ1 周期的マイクロ構造を有する破壊的多孔質媒体において、反応項が大きい場合、均質化プロセスにどのような影響を与えるか？
• RQ2 このようなマルチスケール媒体における反応拡散プロセスの長期挙動を捉える有効なマクロスケール方程式は何か？
• RQ3 孔および亀裂構造は、均質化モデルにおける有効拡散および反応係数にどのように寄与するか？
• RQ4 微小スケールパラメータ ε → 0 における ε-解 u_ε に対して、どのような収束性が成立するか？
• RQ5 均質化問題は、明示的に定義された有効係数を有する対流拡散反応方程式として定式化可能か？

主な発見

• 元の問題の解 u_ε は、ε → 0 のとき L²(Ω_T) 上で強い収束を示し、極限解 u₀ に収束する。
• 均質化問題は、固定領域 Ω × (0,T) 上の対流拡散反応方程式であり、時間微分項に |Z_s|∫_{Y_m} ρ(y) dy の係数がスケーリングされている。
• 有効拡散テンソル Â(x,t) は、マトリックス係数 Â(y,τ) および補正関数 θ と ω₁ を含むセル問題から導出される。
• 対流項 L₁, L₂, L₃ は、反応項 g とマイクロ構造の相互作用から生じ、L₂ および L₃ は反応ポテンシャルの勾配および補正関数に依存する。
• 有効係数は一様有界かつリプシッツ連続であり、これにより均質化問題の解の一意性および適切な定義が保証される。
• 極限問題は反応拡散方程式の構造を保ちつつ、マルチスケール効果を統合した修正された拡散項、対流項、反応項を有する。