[论文解读] 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$ 超收斂率。
• 首次建立 DG 解在內部左勒達烏點處對導數逼近的 $k+1$ 超收斂率。
• 提出嚴謹的數學框架，利用校正函數實現 DG 解與修正勒達烏投影之間的超接近性。
• 展示基於校正步驟的適當初始離散化對於達成 $2k+1$ 超收斂率至關重要，特別是在 $k > 3$ 時。
• 統一並拓展先前的超收斂結果，提供一維線性雙曲系統中超收斂性質的完整圖像。

提出的方法

• 構造一個校正函數 $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$，表明理論界與實際結果之間存在差距。
• 儘管有強烈的數值證據，但 DG 解在單元平均值處的 $2k+1$ 超收斂率證明仍是開放問題。