An Interface Green's Function Framework for Complete Discrete $W^{1,\infty}$ Analysis of Discontinuous Galerkin Methods

技術要約：不連続ガラーキン法における完全な離散 $W^{1,\infty}$ 解析のためのインターフェース・グリーン関数フレームワーク

問題提起
不連続ガラーキン（DG）法の点別誤差解析は重要な研究領域であり、既存の文献では $L^\infty$ ノルムおよび分解された（broken）$W^{1,\infty}$ 半ノルムにおける最適誤差評価が確立されている。しかし、決定的な理論的空白が残されている。すなわち、分解された $W^{1,\infty}$ 半ノルム $\|\nabla(u - u_h)\|_{L^\infty_h(\Omega)}$ は、メッシュ・インターフェースを横切る不連続性を制御できない。その結果、既存の解析では、不連続近似の本質であり、安定性、保存性、および後処理に不可欠なインターフェースのジャンプを直接制御することができない。さらに、凸多面体領域において、弱離散最大値原理が現在 DG 法では利用できない状況下で、多項式次数が 1 を超える場合に、劣最適な対数因子（例：$|\ln h|^\beta$）を伴わずに最適評価を導出することは、永続的な課題であった。

手法
著者らは、完全な離散 $W^{1,\infty}$ 解析を実現するために、新しいインターフェース・グリーン関数フレームワークを提案している。このフレーム手法の核心的な技術革新は、解またはその勾配を推定するために単一の離散デルタ関数を利用する古典的なグリーン関数解析から脱却している点にある。代わりに、このフレームワークは、隣接する離散デルタ関数の差によって生成されるインターフェース・グリーン関数を導入する。

1. インターフェース・グリーン関数の定義:

   • 内部メッシュ・インターフェースに対して、フレームワークは $-\Delta \Gamma_{if} = \delta_{h, K_+^*} - \delta_{h, K_-^*}$ を満たすグリーン関数 $\Gamma_{if}$ を定義する。ここで $\delta_{h, K}$ は要素 $K$ 上に支持される滑らかな離散デルタ関数を表す。このソース項は、自然に不連続近似のインターフェースを横切るジャンプを表している。
   • 境界面については、単一の離散デルタ関数を用いた同様のグリーン関数 $\Gamma_{bf}$ が定義される。

2. 表現公式:
   数値解（およびその勾配）のジャンプは、これらのインターフェース・グリーン関数とそれらの DG 近似を用いて表現される。これにより、インターフェース量の推定は、インターフェース・グリーン関数自体の近似誤差の推定へと還元される。

3. 主要な解析的要素：局所エネルギー評価:
   証明は、領域の**ダイアディック分解（dyadic decomposition）**を介して確立された、これらのインターフェース・グリーン関数の新しい局所エネルギー評価に基づいている。

   • 著者らは、ソース変数に関するテイラー展開の議論を用いている。
   • 決定的なのは、隣接する離散デルタ関数の差（$\delta_{h, K_+^*} - \delta_{h, K_-^*}$）が、特異なソースのレベルで**打ち消し効果（cancellation effect）**を生み出すことである。
   • この打ち消しにより、支配的な特異寄与が除去され、古典的な推定と比較してより鋭い局所的な推定（具体的には、$h$ の追加の因子）が得られる。
   • この改善により、既存の離散トレース不等式と最大値ノルム推定を直接適用した際に通常生じる対数因子を導入することなく、インターフェース推定を導出することが可能となる。

主要な貢献と結果
本論文は、以下のノルムによって定義される「完全な」離散 $W^{1,\infty}$ 理論を確立している：
$$|||v|||_{W^{1,\infty}_h(\Omega)} = \max_{K \in \mathcal{T}_h} \|\nabla v\|_{L^\infty(K)} + \max_{F \in \mathcal{E}_h} h_F^{-1} \|[v]\|_{L^\infty(F)}$$
このノルムは、分解された勾配とスケーリングされたインターフェース・ジャンプを同時に測定する。

このフレームワークは、ポアソン方程式に関する凸多面体領域上の以下の手法に適用されている：

• 対称内部ペナルティ DG (IPDG) 法
• ハイブリッド混合 DG 法
• 2つの代表的なハイブリッド化可能 DG (HDG) 法 (HDG-1 および HDG-2)

主な理論的結果:

• 最適なジャンプ推定: 著者らは、インターフェース・ジャンプ $[u - u_h]$ および数値フラックスの最適な最大値ノルム推定を証明し、凸多面体領域において以前に見られた対数因子 $\beta$ を除去した。
• 完全な誤差評価: 彼らは、分解された勾配とインターフェース・ジャンプを同時に制御する最適な誤差推定を導出し、これらの DG 法に関する点別解析を完結させた。
• 4次問題への拡張: フレームワークは、バイハーモニック方程式に対する $C^0$ 内部ペナルティ法へと拡張されている。著者らは、法線微分のジャンプに関する新しい最大値ノルム推定を導出した。これは、この手法の点別解析においてこれまで扱われていなかった量である。

意義と主張
本論文の主要な貢献は方法論的であると主張している。離散デルタ関数の差に基づくインターフェース・グリーン関数を導入することで、著者らは、凸多面体領域における DG 法における弱離散最大値原理の欠如という根本的な障壁を克服している。

• 対数因子の除去: フレームワークは、古典的なトレース不等式のアプローチから継承される劣最適な対数因子をうまく排除し、多項式次数が 1 より大きい場合でも最適な推定を提供している。
• 汎用性: この構成は、近似の具体的な数値定式化ではなく、局所的な不連続構造に依存している。著者らは、この手法が十分に一般的であり、広範なクラスの不連続および非適合有限要素法に適用可能であると主張している。
• 完全性: 本研究は、指定された DG 法に対して、分解された勾配とインターフェース・ジャンプを同時に制御する、最初の完全な離散 $W^{1,\infty}$ 解析を提供するものであり、不連続手法の点別近似に関する理論的理解における決定的な空白を埋めるものである。