[論文レビュー] Superconvergence of Discontinuous Galerkin method for linear hyperbolic equations

Waixiang Cao, Zhimin Zhang+1 more|arXiv (Cornell University)|Nov 27, 2013

Advanced Numerical Methods in Computational Mathematics参考文献 11被引用数 8

ひとこと要約

本稿は、線形双曲型方程式に対する不連続ガレルキン（DG）法のダウンウィンド点および領域平均における長年の予想である $2k+1$ のスーパー収束率を、新規の補正関数アプローチを用いて証明する。また、初めてとして、内部左ラダウ点における導関数近似の $k+1$ のスーパー収束率を確立し、適切な初期離散化のもとで数値的検証により理論的予測が確認されている。

ABSTRACT

In this paper, we study superconvergence properties of the discontinuous Galerkin (DG) method for one-dimensional linear hyperbolic equation when upwind fluxes are used. We prove, for any polynomial degree $k$, the $2k+1$th (or $2k+1/2$th) superconvergence rate of the DG approximation at the downwind points and for the domain average under quasi-uniform meshes and some suitable initial discretization. Moreover, we prove that the derivative approximation of the DG solution is superconvergent with a rate $k+1$ at all interior left Radau points. All theoretical finding are confirmed by numerical experiments.

研究の動機と目的

• 線形双曲型方程式に対する不連続ガレルキン（DG）法において、ダウンウィンド点および領域平均での $2k+1$ のスーパー収束率が長年未解決であった問題を解明すること。
• 初めてとして、内部左ラダウ点における導関数近似の $k+1$ のスーパー収束率を確立すること。
• 修正されたラダウ射影との超接近性を達成するために、補正関数を用いた厳密な数学的枠組みを構築すること。
• 特に $k > 3$ の場合に、$2k+1$ のスーパー収束率を達成するためには、補正手順に基づいた適切な初期離散化が不可欠であることを示すこと。
• 1次元線形双曲型系におけるスーパー収束特性の包括的図を提示することで、従来のスーパー収束結果を統合・拡張すること。

提案手法

• 補正関数 $w$ がすべてのダウンウィンド点で消えるように構築され、DG解と $P_h^-u - w$ の間の超接近性が実現される。ここで $P_h^-u$ は、正確な解をダウンウィンド点で補間する切断ラダウ射影である。
• 解析はエネルギー推定と、ダウンウィンド点で正確な解を補間する特別な射影 $P_h^-$ に依存しており、誤差構造における高次近似を保証する。
• 時間積分は正確に行われると仮定し、従来の研究とは異なる独自の解析アプローチを採用。補正関数が双曲型問題における誤差伝搬を制御する役割を強調する。
• DG解と補正済みラダウ射影との間で、$2k+1$ 階の超接近性が確立され、これによりダウンウィンド点および領域平均でのスーパー収束が導かれる。
• 内部左ラダウ点における導関数スーパー収束は、同じ超接近性フレームワークから導出され、DG式の直交性および安定性特性を活用する。
• さまざまな初期離散化手法を用いた数値実験により理論的収束率を検証し、$k=4$ の場合にのみ提案された補正ベース手法が $2k+1$ の収束率を達成することを示す比較が行われた。

実験結果

リサーチクエスチョン

• RQ1 線形双曲型方程式に対するDG法において、ダウンウィンド点での $2k+1$ のスーパー収束率が厳密に証明可能か？
• RQ2 DG解の導関数近似は内部左ラダウ点でスーパー収束するか？ もしそうであれば、その収束率はどの程度か？
• RQ3 特に $k > 3$ の場合に、高次スーパー収束を達成するための初期離散化の役割は何か？
• RQ4 同じ枠組みのもとで、領域平均の $2k+1$ のスーパー収束率が証明可能か？
• RQ5 提案された補正関数アプローチは、既存手法と比較して、最適なスーパー収束率を達成する点でどのように異なるか？

主な発見

• 本稿は、線形双曲型方程式に対するDG法において、ダウンウィンド点での $2k+1$ のスーパー収束率が平均的意味で証明され、DG法における長年の予想が裏付けられた。
• 領域平均に対しても $2k+1$ のスーパー収束率が確立され、従来の結果よりも $1/2$ 階高い。
• 初めてとして、すべての内部左ラダウ点におけるDG解の導関数近似の $k+1$ のスーパー収束率が証明された。
• 数値実験により理論的収束率が確認され、$k=4$ の場合にのみ、提案された補正ベース初期離散化が $2k+1$ のスーパー収束率を達成することが示された。
• ダウンウィンド点における $2k+1/2$ のスーパー収束率は、理論的境界において最適でないことが判明した。数値的収束率は $2k+1$ であるため、理論的見積もりと実際の性能の間にギャップが存在する。
• セル平均における $2k+1$ のスーパー収束率の証明は、強力な数値的証拠があるものの、依然として未解決の問題のままである。