A collocation scheme that is equivalent to discontinuous Galerkin discretizations

技術的概要：不連続ガラーキン離散化と同等の連合スキーム

問題提起
本論文は、非構造化グリッド上の偏微分方程式（PDE）に対する安定かつ高次の離散化の構築に取り組む。和による部分（SBP）演算子は安定性のための堅牢な枠組みを提供するが、単体要素（三角形など）上の多次元 SBP 演算子の構築は複雑であり、しばしば閉形式の表現を欠く。Chan による先行研究は、有限要素離散化の中間段階として、射影に基づく SBP 演算子（ここではモード連合、または MC 演算子と呼ばれる）を導入した。しかし、この演算子を連合スキームで直接使用することの具体的な特性、特に不連続ガラーキン（DG）法との関係については、十分に探求されていなかった。中心的な問いは、この特定の MC 演算子に基づく半離散化が、同じ基礎となる求積法を用いたモード DG 離散化と同等の解を生み出すかどうかである。

手法
著者は、正規直交多項式基底と正の内部求積則を用いて MC 演算子を定義する。核心的な手法は以下の通りである：

1. 演算子の定義：第一微分 MC 演算子 $D_{xd} = V_{xd}V^T W$ を定義する。ここで、$V$ は正規直交基底関数のヴァンデルモンド行列、$V_{xd}$ はそれらの微分を含む行列、$W$ は求積重みの対角行列である。
2. 理論的同等性の証明：論文は、MC 演算子を用いたスカラー非線形双曲型初期境界値問題（IBVP）の半離散化が、モード DG 解ベクトル $\tilde{u}(t)$ をノード空間への射影した解ベクトル $u(t)$（$u(t) = V\tilde{u}(t)$）を正確に生成することを示す証明を構築する。この証明は、解空間を多項式部分空間と演算子の零空間に分解することに基づいている。
3. スペクトル解析：著者は、結果として生じるシステム行列の固有値を解析する。MC 演算子のスペクトルは、対応する DG 演算子の固有値と、零空間に関連するゼロ固有値の集合から成ることが示される。
4. 零空間不整合の処理：MC 演算子は一般的に零空間不整合（その零空間は定数ベクトルだけでなく、より多くのものを含む）であることを論文は認めている。非定常問題では、これらの余分なモードはゼロに初期化され、ゼロのままとなる。定常問題については、著者は多項式の正確性を損なうことなく零空間の特異性を除去するための局所射影安定化（LPS）を提案する。
5. 数値的検証：理論は、三角形メッシュを用いた数値実験により検証される：
   • 定数係数の移流方程式。
   • 非線形バーガース方程式。
   • 「標準的」MC 離散化と「エントロピー保存的」（歪対称）定式化との比較。
   • 負の重みを含む求積法を用いたテスト。

主要な貢献

• 同等性定理：主要な貢献は、同じ求積法を共有する限り、MC 演算子を用いたモード連合離散化が数学的にモード DG 離散化と同等であるという証明である。これは、連合スキーム内の自由度（求積点）の数に関わらず成り立つ。
• スペクトル特性：論文は、MC 演算子のスペクトル半径が、ゼロ固有値の付加を除き、基礎となる DG 法によって実質的に決定されることを確立する。これは、ノードのクラスタリング（三角形上のスペクトル要素法でしばしば問題となる）が、従来のノード SBP 演算子の場合と同様に MC 演算子のスペクトル半径に悪影響を及ぼさないことを意味する。
• 零空間管理：MC 演算子は零空間不整合であるが、初期条件が正しく射影されていれば、時間依存シミュレーションにとって有害ではないことがこの研究で明確にされた。定常状態問題については、論文は、射影に基づく数値的減衰（LPS）を通じて整合性を回復する方法を実証する。
• エントロピー保存のニュアンス：この研究は、歪対称分割に基づくエントロピー保存的 MC 離散化は、残留項が多項式空間内に存在しない可能性があるため、本質的に DG と同等ではないことを明らかにする。しかし、残留項を多項式基底 onto 射影することで、同等性は回復される。

結果

• 定数係数移流：数値結果は、多項式次数や求積の正確性のレベルに関わらず、MC とモード DG 解間の $L_2$ 差が機械精度（$\approx 10^{-15}$）であることを示している。
• スペクトル半径：固有値分布は、三角形上のノードクラスタリングにより壊滅的な増大を被る最小ノルム SBP 演算子と比較して、MC 演算子が有利なスペクトル半径を示すことを確認する。
• 負の重み：同等性は、求積法が負の重みを含む場合でも成り立つ。これは、従来の SBP エネルギー安定性理論が安定性を予測することに失敗する可能性があるシナリオである。安定性は、エネルギー安定な DG 解との同等性によって保証される。
• バーガース方程式：非線形バーガース方程式において、「標準的」MC 離散化は DG と同等である。「エントロピー保存的」変種は、残留項が多項式基底によって張られないため、当初は不一致を示す。しかし、残留項を基底 onto 射影すると、同等性は機械精度まで回復する。

意義と主張
論文は、Chan によって導入されたモード連合演算子は単なる中間ツールではなく、実質的にモード DG 法である有効な離散化スキームであると主張する。その意義は、非構造化単体グリッド上で高次かつ安定なスキームを構築する際の簡素化にある。MC 演算子を使用することで、演算子の閉形式表現を活用しつつ、DG 法の安定性とスペクトル特性を維持することができる。

著者は、両手法は同等であるが、選択は応用次第であると謙虚に結論づける。質量行列と剛性行列が事前に計算される線形問題では DG が好まれる可能性がある。しかし、陽的時間進行または行列フリー陰的スキームで解かれる非線形問題については、MC 演算子が計算上の利点を提供しうる。論文は、多次元 SBP 演算子を構築するための好ましい手法として、モード連合の使用を強く推奨している。