Mixed Finite Element Methods for a Dirac Source: Divergence-Form Splitting and L^p Error Analysis

技術要約：ディラック源に対する混合有限要素法

問題設定
本論文は、有界で単連結なリプシッツ領域 $\Omega \subset \mathbb{R}^2$ における、ディラック源を持つ拡散方程式の数値解法を扱っている：
$$-\nabla \cdot (A\nabla u) = f_1 + \delta_{x_0} \quad \text{in } \Omega, \quad u = g \quad \text{on } \partial\Omega.$$
著者らは、混合有限要素法（MFEM）において、ディラック源による主要な困難は、単に解の正則性が低いこと（通常、$H^1$ ではなく $W^{1,p}$ ($p<2$) に属する）にあるのではなく、根本的な双対性の破綻にあることを指摘している。標準的な混合定式化では、ルベーグ空間（$L^q$）として設計されたスカラー空間に対して保存則をテストする。しかし、ディラック測度 $\delta_{x_0}$ は、一般的なテスト関数 $v \in L^q$ の双対空間には属さない。したがって、ペアリング $\langle \delta_{x_0}, v \rangle$ は定義されない。この破綻は連続レベルで発生するものであり、離散化に依存せず、メッシュの細分化や要素の選択のみでは修復できない。

手法：発散形式の分割（Divergence-Form Splitting）
メッシュや近似空間を変更することなくこの双対性の破綻を解決するために、著者らは、物理的なフラックスから明示的なベクトル場 $f_2$ を差し引く発散形式の分割を導入する。具体的には、$\nabla \cdot f_2 = \delta_{x_0}$ を満たす以下の $f_2$ を定義する：
$$f_2(x) := \frac{1}{2\pi} \frac{x - x_0}{|x - x_0|^2},$$
これは、分布の意味で $\nabla \cdot f_2 = \delta_{x_0}$ を満たし、$f_2 \in L^p(\Omega)^2$ （すべての $1 < p < 2$ に対して）に属する。

混合未知数は、修正フラックスとして再定義される：
$$\sigma := -A\nabla u - f_2.$$
この変換により、問題は以下の一次系へと書き換えられる：
$$A^{-1}\sigma + \nabla u = -A^{-1}f_2, \quad \nabla \cdot \sigma = f_1.$$
この新しい定式化において：

1. 保存則 $\nabla \cdot \sigma = f_1$ は、正則なデータ $f_1$ のみを扱うため、MFEMが要求する通り、完全なルベーグ空間を用いてテストすることが可能となる。
2. 特異性は構成方程式の中に閉じ込められ、フラックスのテスト関数とペアリングされる。$f_2 \in L^p$ であり、テスト関数 $\tau \in L^{p'}$ （ここで $1/p + 1/p' = 1$）であるため、ヘルダーの不等式により、このペアリングは定義可能であることが保証される。
3. ベクトル場 $f_2$ は、ソースの位置 $x_0$ にのみ依存し、係数 $A$ やメッシュには依存しない。したがって、係数依存の特異解や離散的なデルタ近似を必要としない。この手法は、$x_0$ が要素の内部にあるか頂点上にあるかにかかわらず機能する。

主要な貢献と解析

1. $L^p$ スケールにおける正則性と誤差解析：
   修正フラックス $\sigma$ は（係数 $A$ が $x_0$ 付近で単位行列のスカラ倍でない限り）一般に $|x-x_0|^{-1}$ 次の特異性を保持するため、$\sigma$ は $L^p(\Omega)^2$ ($1 < p < 2$) に属するが、$L^2(\Omega)^2$ には属さない。その結果、誤差解析はヒルベルトスケール ($p=2$) 未満で行われなければならない。

   • フラックス誤差： 著者らは、フラックスに対する**準最適近似界（quasi-best approximation bound）**を証明している。準一様メッシュにおいて、フラックス誤差は $O(h^{2/p - 1})$ の収束率を示す。このレートはシャープであり、極（pole）を含む要素上の特異性の局所的な近似によって決定される。
   • スカラー変数の誤差： 離散的な $L^{p'}$ インフ・スュップ条件と双対性議論を用いることで、$u$ の誤差は $u$ の最良の区分的定数近似によって漸近的に支配されることが示される。準一様メッシュにおいて、これは $L^2$ 誤差が $O(h |\log h|^{1/2})$ のオーダーであることを意味する。

2. 比較フラックスの構成：
   標準的なカノニカル補間（ラヴィアート・トーマス）は、特異性の角度成分によって $x_0$ を通るエッジ上の法線モーメントが非可積分になる場合、修正フラックスに対して失敗する可能性がある。著者らは、極の周囲のパッチにおける局所的なフラックス・バランスを通じて、**平衡化された比較フラックス（equilibrated comparison flux）**を構築する。これにより、カノニカルな自由度が定義できない場合でも、誤差解析のための有効な比較場の存在が保証される。

3. メッシュ・グレーディングと複雑性：
   準一様メッシュにおける劣最適なレート $h^{2/p - 1}$ は、手法の性質ではなく、特異性の性質である。**幾何学的に $x_0$ へ向かって細分化されるグレーデッドメッシュ（graded meshes）**を採用することで、著者らはフラックス誤差を、要素数 $N$ に関する一次の複雑性 $O(N^{-1/2})$ に回復できることを示している。

4. 事後誤差推定（A Posteriori Error Estimation）：
   本論文は、一次系最小二乗法と同様の、ルベーグスケールにおける残差ノルム等価性を確立している。これにより、拡張されたペア（修正フラックス、復元されたポテンシャル）に対する計算可能な $L^p$ 推定値が得られる。この推定値は、バブル関数やジャンプ項を必要とせず、信頼性と局所的な効率性を備えていることが示されている。

結果と数値実験
最低次ラヴィアート・トーマス（$RT_0$）および区分的定数（$P_0$）要素を用いた数値実験により、理論的知見が裏付けられた：

• 適応計算（Adaptive Runs）： $L^p$ 推定値によって駆動される適応的細分化は、フラックスに対して最適な収束率 $O(N^{-1/2})$ を達成し、特異性のペナルティを効果的に除去している。
• 一様計算（Uniform Runs）： 一様メッシュにおいて、フラックス誤差は予測されたレート $h^{2/p - 1}$ を示し、これは $p$（例：$p=1.2$ のとき $h^{0.667}$）によって変化する。
• スカラーの収束： スカラー変数は、グレーデッドメッシュ上で $O(h)$ のオーダー（対数因子は無視）で収束し、これは $p$ に依存しない。これは、双対性議論がポテンシャルに対して最適なレートを正常に回復したことを裏付けている。
• $L^2$ 発散： 実験により、離散フラックスの $L^2$ ノルムが対数的に発散することが確認された。これは、この問題が標準的なヒルベルトの枠組み内で解析できないことを裏付けている。

意義
本論文は、ディラック源を伴う混合法に内在する双対性の破綻に対する、ロバストでメッシュに依存しない解決策として、発散形式の分割を提示している。特異性を保存則から構成方程式へと移すことで、離散系のシステム行列を変更したり、人工的な離散デルタを導入したりすることなく、標準的な混合定式化をディラック問題に適用することを可能にしている。このアプローチの代償は $L^p$ 誤差解析（$p<2$）への移行であることを認めているが、適切なメッシュ・グレーディングと組み合わせることで、最適な収束率が得られることを示している。この手法は、ディラック源を $L^p$ 場の発散形式の負荷として扱うことで、広範な離散化手法（コンフォーマル、DG、有限体積法）に適用可能である。