Evolving finite elements for advection diffusion with an evolving interface

技術的概要：進化する界面を伴う移流拡散問題のための進化有限要素法

問題定義
本論文は、移動する界面 $\Gamma(t)$ を含む静止領域 $\Omega$ 上で定義された放物型方程式の数値近似を取り扱う。この界面は領域を 2 つの部分領域 $\Omega_1(t)$ および $\Omega_2(t)$ に分割する。支配方程式は、以下の条件のもとで、これらの部分領域内のスカラー場 $u_1$ と $u_2$ の輸送および拡散をモデル化する：

• 時間依存係数（$A_i, B_i, C_i$）を有する各部分領域内での移流 - 拡散 - 反応ダイナミクス。
• 外境界 $\partial\Omega$ におけるディリクレ境界条件。
• 界面を跨ぐ解の連続性（$u_1 = u_2$）。
• 界面を跨ぐ prescribed なフラックスのジャンプ（$A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$）。

特定された主な課題は、時間経過に伴う部分領域の進化と、移動する界面を跨ぐフラックスのジャンプの存在である。流体力学、材料科学、細胞生物学に応用例はあるが、本論文は問題の数学的および数値的定式化に焦点を当てている。

手法
著者らは、[21] で確立された抽象的な枠組みに基づき、**進化有限要素空間（ESFEM）**を利用する任意ラグランジュ - オイラー（ALE）アプローチを提案する。

1. 関数解析的設定：連続的な問題は、進化ヒルベルト空間の文脈内で定式化される。著者らは、$\phi_t$ を参照領域から進化領域へ関数を輸送するフロー写像とする、互換性のある空間のペア $(X(t), \phi_t)$ を定義する。これにより、弱物質微分の定義と、適切に定式化された弱形式の導出（定理 2.16）が可能となる。
2. メッシュ構築：
   • 初期メッシュ：領域は、線形片多項式単体を用いて初期に三角分割される。極めて重要なのは、界面ファセットの頂点が初期界面 $\Gamma(0)$ 上に正確に存在するように制約されている点である（仮定 M5）。
   • 等パラメータ写像：高次精度を達成するために、線形メッシュは微分同相写像 $\Psi_h$ を用いて等パラメータメッシュへ写像される。この写像は、メッシュの連続性を保ちながら、界面ファセットを正確な曲線界面 $\Gamma(0)$ 上に押しやる。
   • 時間進化：メッシュは、既知の滑らかな大域的速度場 $w$ に従ってラグランジュ節点を移動させることで進化させる。離散フロー $\Phi^h$ は、連続フローの多項式補間として構成され、メッシュが常に適合する許容分割として維持されることを保証する。
3. リフト写像：離散的計算領域（曲線界面を近似するもの）から正確な連続領域へ関数を写像するリフト作用素 $\Lambda_h$ が導入される。これは誤差解析において離散解と厳密解を比較するために不可欠である。
4. 離散スキーム：この手法は、次数 $k$ の等パラメータ有限要素を採用する。離散変分問題は、離散物質微分とリフトされた係数を用いて定式化される。このスキームは、時間ステップ法（数値実験では具体的には後退差分公式、BDF）を用いて解かれる。

主な貢献

• 適切性：本論文は、データに関する特定の仮定のもとでの追加的正則性結果を含め、連続問題の解の存在と一意性を証明する厳密な関数解析的設定を提供する（定理 2.16）。
• 等パラメータ ESFEM 枠組み：著者らは、フラックスのジャンプを伴う両側界面問題を処理するために、抽象的な進化有限要素枠組みを拡張する。彼らは、曲線境界を正確に近似し、高次収束を可能にする等パラメータ要素の構築を詳述する。
• 最適誤差 bound：主要な理論的結果（定理 5.1）は、最適次数の誤差 bound を確立する。具体的には、$L^2$ ノルムにおける誤差は $O(h^{k+1})$ によって、$H^1$ 半ノルムにおける誤差は $O(h^k)$ によって抑えられる。ここで、$h$ はメッシュサイズ、$k$ は多項式の次数である。これらの bound は定数まで成り立ち、既知の滑らかな解の補間誤差と一致する。
• 実装と検証：本論文には、Firedrake パッケージと GMSH を用いた堅牢な実装が含まれている。2 次元および 3 次元の数値実験は、様々な多項式次数に対する理論的収束率を検証する。

結果
数値実験は、界面を跨いで連続的だが微分不可能な既知の厳密解を持つモデル問題において、手法の有効性を示す。

• 収束：次数 $k$ の空間離散化に対して、数値結果は $L^2$ ノルムにおいて次数 $k+1$ の収束と、$H^1$ 半ノルムにおいて次数 $k$ の収束を確認する。
• 堅牢性：結果は、メッシュが一様細分化ではなく独立して生成された場合であっても、2 次元および 3 次元の両方の領域で成り立ち、手法の安定性を示している。
• コードの入手可能性：著者らは再現性とコミュニティによる更なる利用を可能にするため、シミュレーションコードを提供する。

意義と範囲
本論文は、高次等パラメータ有限要素を用いて移動界面を伴う移流 - 拡散問題を解くための、完全な理論的および実践的枠組みを提供すると主張している。その意義は、フラックスのジャンプを伴う界面問題に対する抽象的な進化空間理論と、実用的な実装との間のギャップを埋める点にある。

著者らは速度場に関して控えめな範囲を維持している：彼らは、大域的で滑らかな速度場が事前に与えられると仮定している。彼らは明示的に、そのような速度場を計算する方法（例えば、メッシュの絡みつきを防ぐための再パラメータ化など）には取り組んでいないと述べており、そのような手法は文献 [8, 18, 19] に存在すると指摘している。さらに、理論的には領域の進化が解の一部である結合された PDE システムにも適用可能であるが、現在の解析はそうした結合問題にはカバーしていない。焦点は、既知の領域進化が与えられた場合の放物型伝達問題の近似に厳密に置かれている。