[論文レビュー] Necessity of numerical smoothness

T. Sun|arXiv (Cornell University)|Jul 12, 2012

Probabilistic and Robust Engineering Design参考文献 4被引用数 3

ひとこと要約

本論文は、PDEの数値解法における最適収束率を達成するためには、局所多項式解の微分のジャンプを測る空間的滑らかさ指標の有界性が必須であることを確立している。著者らは、これらの指標が無限大に発散する場合には、期待される $\mathcal{O}(h^{p+1})$ の収束率を達成できないことを証明しており、これにより、事後誤差解析におけるこうした指標の使用が正当化され、数値的滑らかさの欠如に起因する収束喪失を実用的に診断するためのツールが提供される。

ABSTRACT

Numerical solutions of differential equations are usually not smooth functions. However, they should resemble the smoothness of the corresponding real solutions in one way or another. In two of our recent papers, a kind of spacial smoothness indicators was defined and subsequently applied on the {\it a posteriori} error analysis. Here we prove that the boundedness of those smoothness indicators is actually a necessary condition for a piecewise polynomial function to approximate a smooth function with optimal convergence rate. This should help in validating the error analysis in the related papers. Moreover, the result of this paper provides an efficient practical method to detect the loss of convergence rate due to the lack of numerical smoothness, hence it serves as a criterion for the qualities of many numerical schemes.

研究の動機と目的

• PDEの数値解法における最適収束率を達成するための数値的滑らかさの必要性を確立すること。
• 先行研究[2]および[3]で提唱された事後誤差解析における滑らかさ指標の使用を検証すること。
• DGやWENOのようなスキームにおける数値的滑らかさの不足に起因する収束喪失を検出するための実用的基準を提供すること。
• 従来の安定性概念が数値的滑らかさを強制しない場合には、それが最適収束を保証するのに不十分であることを示すこと。

提案手法

• セル界面における局所多項式解の微分のジャンプに基づいて、空間的滑らかさ指標 $ S^p = (\bar{M}, \bar{D}) $ を導入する。
• 滑らかさの3種類の定義を提示：$ W^{p+1}_\infty $、$ H^{p+1} $、$ W^{p+1}_1 $ は、それぞれ $ \bar{D} $、$ \|\bar{D}\|_{L^2} $、$ \|\bar{D}\|_{L^1} $ の有界性に基づく。
• スケーリング論法と近似理論を用いて、滑らかさ指標と正確解と再構成解の間の $ L^1 $、$ L^2 $、$ L^\infty $ 誤差の関係を確立する。
• 滑らかさ指標 $ Q(D^0_i, \dots, D^p_i) $ を用いて、誤差 $ \|u - u^R\| $ の下界を導出し、最適収束には $ Q \leq \mathcal{O}(1) $ が必須であることを示す。
• $ L^\infty $、$ L^2 $、$ L^1 $ ノルムにおける最適 $ \mathcal{O}(h^{p+1}) $ 収束を達成するためには、$ |D^k_i| \leq \mathcal{O}(1) $ が必須であることを証明する。
• Runge-Kutta DGおよびWENOなどの標準的スキームに結果を適用し、滑らかさ指標の制御が不十分な場合には収束率の喪失が生じることを示す。

実験結果

リサーチクエスチョン

• RQ1 数値解法がPDEの解法において最適収束率を達成するためには、滑らかさ指標の有界性が必須であるか？
• RQ2 [2]および[3]で用いられた滑らかさ指標は、収束劣化の信頼できる診断指標として機能するか？
• RQ3 微分のジャンプの有界性によって強制される数値的滑らかさは、最適収束のための必須なグローバル性質であるか？
• RQ4 従来の安定性概念は最適収束を保証するのに十分か、それとも数値的滑らかさがより根本的な要件であるか？
• RQ5 滑らかさ指標は、滑らかさの欠如に起因するスキームの故障を検出する実用的ツールとして使用可能か？

主な発見

• 最適収束 $ \|u - u^R\|_{L^\infty} \leq \mathcal{O}(h^{p+1}) $ を達成するためには、滑らかさ指標 $ Q(D^0_i, \dots, D^p_i) \leq \mathcal{O}(1) $ が必須であり、これはすべての $ k $ と $ i $ に対して $ |D^k_i| \leq \mathcal{O}(1) $ が成り立つことを意味する。
• $ L^1 $ 収束が $ \mathcal{O}(h^{p+1}) $ の割合で達成されるためには、$ \sum_{i=1}^{N-1} h \sqrt{Q(D^0_i, \dots, D^p_i)} \leq \mathcal{O}(1) $ が必須である。
• $ L^2 $ 収束が $ \mathcal{O}(h^{p+1}) $ の割合で達成されるためには、$ \sum_{i=1}^{N-1} h \, Q(D^0_i, \dots, D^p_i) \leq \mathcal{O}(1) $ が必須である。
• 最適収束のためには、$ J^k_i $（$ k $ 階微分のジャンプ）について $ |J^k_i| \leq \mathcal{O}(h^{p+1-k}) $ が必須である。
• 滑らかさ指標が $ \mathcal{O}(1) $ を超えると、スキームは最適収束率を達成できず、真の解の非滑らかさまたはスキームの欠陥を示している。
• 結果は、[2]および[3]で提唱された手法の妥当性を確認しており、誤差推定においてグローバル非線形解析を回避する目的で滑らかさ指標が使用されたことの正当性が裏付けられる。