Embedded Trefftz DG method for reaction-diffusion problems on anisotropic meshes

技術要約：異方性メッシュ上での反応拡散問題に対する埋め込み型Trefftz DG法

問題設定
本論文は、有界領域 $\Omega \subset \mathbb{R}^2$ における、ディリクレ境界条件を伴う形式 $-\varepsilon^2 \Delta u + u = f$ の形で表される、特異摂動反応拡散問題の数値解法を扱う。焦点は $0 < \varepsilon \ll 1$ の領域にあり、これは解が鋭い境界層や強い異方性の特徴を示すシナリオである。標準的な等方性メッシュはこのような問題に対して非効率であることが多く、局所的な特徴をより少ない自由度で解像するために、異方性メッシュ細分化（特定の方向に引き伸ばされた四角形要素）が必要となる。不連続ガラーキン（DG）法は、一般的なメッシュ上で高次離散化を行うための柔軟性を備えているが、標準的なテンソル積DG定式化は、巨大な全体代数系をもたらす。Trefftz法は、微分作用子の核（kernel）に含まれる局所近似空間を用いることで自由度（DoF）の削減を実現するが、一般的な楕円型PDEに対して適切なTrefftz空間を構築することは、特に異方性のある曲面四角形要素においては容易ではない。

手法
著者らは、異方性四角形メッシュに合わせて設計された埋め込み型Trefftz不連続ガラーキン（DG）法を提案している。この手法は以下のステップを通じて構築される：

1. メッシュと空間： ドメインは、曲線的な四角形要素へと分割される。各要素 $K$ は、単位正方形からアフィン変換（異方性を処理するため）および $C^1$ 変換（曲率を処理するため）を通じて写像される。離散空間は、補助座標におけるテンソル積多項式で構成される。
2. 埋め込み型Trefftz構築： 明示的にグローバルなTrefftz空間を構築する代わりに、本手法は埋め込み戦略を採用している。各要素上に局所テスト空間 $L_h(K)$ が定義される。この空間は、プルバック要素の長い方の2辺上で消滅するテンソル積多項式からなる。
3. 局所作用子と制約： 局所的な反応拡散作用子 $A_K$ が定義される。埋め込み型Trefftz空間 $T_h(K)$ は、標準的なDG空間の部分空間として定義され、そこではすべての $q_h \in L_h(K)$ に対して、要素ごとの残差 $(A_K v_h, q_h)_K = 0$ が成立する。
4. グローバル定式化： グローバルな解 $u_h$ は、以下の2つの条件を満たすものとして求められる：
   • 局所制約： すべての局所テスト関数 $v_L \in L_h(K)$ に対して $(A_K u_h, v_L)_K = \ell_K(v_L)$ が成立する。
   • グローバルな整合性： すべてのグローバルテスト関数 $v_T$ に対して、標準的な内部ペナルティDG双線形形式 $a_h(u_h, v_T) = \ell_h(v_T)$ が成立する。
5. 実装： 実践においては、Trefftz空間は明示的に組み立てられない。代わりに、本手法は局所作用子行列の特異値分解（SVD）を利用して埋め込みを構築し、Trefftz空間に関連するグローバルな自由度のみを保持する。

主な貢献
本論文の具体的な技術的貢献は以下の通りである：

• 局所テスト空間の設計： 主要な新規性は、テンソル積多項式基底と互換性のある局所テスト空間の設計にある。単体基底を用いた従来の研究とは異なり、このアプローチにより、四角形および六面体メッシュ上での系統的な構築が可能となる。
• 異方性要素上の安定性： 著者らは、異方性のある（場合によっては曲面を持つ）四角形要素上での安定性と準最適性を証明している。重要な側面は、テスト空間の具体的な選択（長い方の辺で消滅すること）であり、これにより、安定性に不可欠な異方性逆不等式およびポアンカレ・ステックレフ不等式の導出が可能となる。
• 誤差解析： 著者らは、エネルギーノルムにおける異方性 a priori 誤差評価を導出している。これらの評価は、実際の異方性の方向における解の挙動を捉えており、等方的な評価よりも鋭い境界を提供する。
• 次元削減： 本手法は、局所近似空間の次元を、標準的なテンソル積の $(p+1)^2$ から $2p+2$ へと削減し、標準的なDG法やハイブリダイズドDG（HDG）法と比較して、グローバルに結合される自由度を大幅に減少させる（特に3次元において）。

結果
数値実験により、理論的知見が検証されている：

• $h$-収束： 境界層を持つ正方形および円形ドメインにおいて、埋め込み型Treffz法はエネルギーノルムにおいて最適な収束率を達成する。標準的なDG法と同じ精度に達するために、著しく少ない自由度を必要とする。
• 曲面ドメイン： 本手法は、混合された三角形および四角形要素（三角形要素は標準的な埋め込み技術を用いて処理される）を含む曲面メッシュ上でも堅牢に動作する。
• **$hp$-収束：** 境界層に適応した9要素メッシュにおいて、本手法は指数関数的な収束を示す。ただし、著者らは、高次の多項式次数（$p \gtrsim 5$）かつ大きな境界層サイズの場合、標準的な埋め込みが不安定になる可能性があると指摘している。これは、コーナー要素におけるテスト空間を修正してTrefftz条件を緩和することで軽減され、安定性と指数関数的収束が回復される。
• 異方性拡散： 本手法は、純粋な異方性拡散問題（反応項なし）に対してもテストされている。これは厳密な理論的枠組みの外にあるが、結果は、$h$-および$p$-細分化戦略の両方において、埋め込み型Trefftz法が一貫して標準的なDG法よりも低い誤差を与えることを示している。

意義と主張
本論文は、PDE作用子の核が明示的に知られていない、あるいは自明である場合でも効果的に機能する、Trefftz型DG法の一般的な枠組みを提供することを主張している。異方性四角形メッシュへの埋め込み型Trefftzアプローチの適用成功により、本研究は、特異摂動問題に対する高次離散化への道を開いた。著者らは、本手法がHDGに匹敵する競争力のある選択肢を提供し、3次元において、高次離散化の近似特性を維持しつつ、最小のグローバル系を実現することを強調している。解析に関しては、$p$-ロバスト性について控えめな姿勢をとっており、粗いメッシュ上で大きな境界層を持つ非常に高い多項式次数において安定性を維持するには、数値実験で示されたように、テスト空間の特定の修正が必要であることを認めている。