[論文レビュー] Higher order derivative estimates for finite-difference schemes

István Gyöngy, Н. В. Крылов|ArXiv.org|May 20, 2008

Differential Equations and Numerical Methods参考文献 27被引用数 4

ひとこと要約

本稿では、2階の退化する放物型および楕円型PDEに対する有限差分スキームの解が、メッシュサイズ $ h $ に依存しない任意の順序の空間微分の推定値を有するための十分条件を確立する。主な貢献は、$ h $-依存のない高階微分の境界を導出するための枠組みを提供することであり、これにより、退化方程式に対してもリチャードソン補外および $ h $ におけるべき級数展開を用いた加速された有限差分スキームの構築が可能になる。

ABSTRACT

We give sufficient conditions under which solutions of finite-difference schemes in the space variable for second order possibly degenerate parabolic and elliptic equations admit estimates of spatial derivatives up to any given order independent of the mesh size.

研究の動機と目的

• 退化する放物型および楕円型PDEに対する有限差分スキームが、メッシュサイズ $ h $ に依存しない任意の順序の空間微分推定値を有するための十分条件を確立すること。
• 一次微分の境界に関する先行研究を高階微分に拡張し、離散解 $ u_h $ の空間的正則性をメッシュサイズ $ h $ にわたって一様に保証すること。
• 離散解 $ u_h $ の $ h $ におけるべき級数展開を活用することで、加速された有限差分スキームを構築する理論的基盤を構築することにより、高次の収束率を達成すること。
• 古典的PDE理論に依存せずに、退化PDEの解の存在と滑らかさが提示された推定値から直接導かれる可能性を示すこと。
• モデル例のクラスに対して主要定理の条件を検証し、実用的スキームへの適用可能性を確認すること。

提案手法

• 空間次元 $ \mathbb{R}^d $ における一般クラスの単調有限差分スキームを定式化し、メッシュサイズ $ h \in (0, h_0] $ でパラメータ化する。空間作用素 $ L_h $ は差分商 $ \delta_h $ と係数 $ p_\lambda, q_\lambda $ によって定義される。
• 係数 $ p_\lambda, q_\lambda $ 及びそれらの差分に対する構造的および有界性条件を課す仮定2.4および2.5を導入し、高階微分の制御を可能にする。
• 有限差分作用素の構造を活用して、第二階差分および高階差分商に対する事前推定値を用いたブートストラップ法を適用する。
• 第二階差分項を主項と余剰項に分解し、ヤングの不等式および重み付き和を用いて慎重に推定する。
• $ \varepsilon $-依存定数を用いて、高階差分項の和を低階項および二次形式で上から抑える重要な不等式 (4.8) を導出する。
• 得られた $ h $-依存のない推定値を用いて、$ u_h $ の $ h $ におけるべき級数展開の存在を正当化し、超収束を実現するリチャードソン補外のための基盤を提供する。

実験結果

リサーチクエスチョン

• RQ1 退化する放物型および楕円型PDEに対する有限差分スキームが、任意の順序の空間微分について $ h $-依存のない推定値を有するための条件は何か？
• RQ2 メッシュサイズ $ h \to 0 $ の全範囲にわたって、退化ケースにおいても離散解 $ u_h $ の滑らかさが一様に制御可能か？
• RQ3 微分推定値を用いて、$ u_h $ の $ h $ におけるべき級数展開を厳密に確立する方法は何か？ その結果、収束加速に与える影響は何か？
• RQ4 どの程度まで、退化PDEの解の存在と正則性を、連続PDE理論に依存せずに、離散的推定値から直接導けるか？
• RQ5 有限差分係数 $ p_\lambda, q_\lambda $ にどのような構造的条件を課すと、高階差分項が $ h $ にわたって一様に制御可能になるか？

主な発見

• 仮定2.4および2.5の下で、$ u_h $ の任意の指定順序 $ m $ までのすべての空間微分について、$ h $-依存のない推定値が確立される。
• 主結果である定理2.1は、再帰的推定スキームを用いて、$ h $ に依存しない高階差分商の事前推定値を提供する。
• 同一の仮定の下で、定理2.2は差分商の推定値から $ h $-依存のない高階微分の境界を導出する。
• 仮定2.4および2.5の条件は、モデル例のクラスに対して検証され、具体的なスキームへの適用可能性が確認された。
• $ u_h $ の $ h $ におけるべき級数展開の存在が確立され、リチャードソン補外を用いた加速スキームの構築が可能になる。
• 適切な異なる $ h $ の近似の線形結合を取ることで、退化方程式に対しても任意の収束率が得られる。