[論文レビュー] Second order parameter-uniform convergence for a finite difference method for a singularly perturbed linear reaction-diffusion system

Joseph Paramasivam Mathiyazhagan, Valarmathi Sigamani+1 more|arXiv (Cornell University)|Apr 5, 2010

Differential Equations and Numerical Methods被引用数 6

ひとこと要約

本稿では、異なる小さなパrameterを有する線形反応拡散系に対して、シシュキンの区分的均等メッシュ上での有限差分法を提示する。パrameterに依存しない2次収束性を証明し、すべての摂動パrameterにわたって一様に $ C N^{-2} (\text{ln} N)^3 $ の誤差境界を達成する。グリーン関数手法に依存しない。

ABSTRACT

A singularly perturbed linear system of second order ordinary differential equations of reaction-diffusion type with given boundary conditions is considered. The leading term of each equation is multiplied by a small positive parameter. These singular perturbation parameters are assumed to be distinct. The components of the solution exhibit overlapping layers. Shishkin piecewise-uniform meshes are introduced, which are used in conjunction with a classical finite difference discretisation, to construct a numerical method for solving this problem. It is proved that the numerical approximations obtained with this method is essentially second order convergent uniformly with respect to all of the parameters.

研究の動機と目的

• 線形反応拡散系において、すべての特異的摂動パrameterにわたって一様に2次収束性を達成する数値解法を開発すること。
• 系に存在する異なる小さな摂動パrameterによって生じる重複する境界層を扱うこと。
• グリーン関数手法に依存せずにパrameterに依存しない誤差見積もりを確立することにより、高次元および非線形問題への応用可能性を高めること。
• 区分的均等なシシュキンメッシュと古典的有限差分法を用いて、安定性と精度を保証する強固な数値スキームを構築すること。

提案手法

• メッシュの遷移点は特異的摂動パrameter $ \varepsilon_i $ に基づいて定義され、境界層内での解像度が向上するように設計される。
• 2階常微分方程式系に古典的な中心差分離散化を適用し、離散線形系を形成する。
• シシュキンメッシュ上での離散作用素に対して、離散的最大原理および安定性特性を確立する。
• 局所切断誤差は2つの成分に分解して分析する：滑らか成分と特異成分。数学的帰納法とメッシュ依存推定を用いて鋭い境界を導出する。
• 誤差解析では成分ごとの見積もりを用い、遷移点における層の挙動を制御するための比 $ R(w_i(x_k)) $ を導入する。
• 証明はグリーン関数手法を避け、直接推定と帰納法に依存するため、多次元および非線形系への一般化の可能性が向上する。

実験結果

リサーチクエスチョン

• RQ1 異なる $ \varepsilon_i $ を有する線形反応拡散系において、シシュキンメッシュ上での有限差分法が、すべての特異的摂動パrameterにわたって一様に2次収束性を達成できるか。
• RQ2 異なる $ \varepsilon_i $ によって層が重複する場合、解の特異成分に対する鋭い境界をどのように導出できるか。
• RQ3 グリーン関数解析に依存せずにパrameterに依存しない誤差見積もりを確立できるか。これにより、複雑な問題への応用性が向上するか。
• RQ4 収束の正確な次数は何か。誤差境界に現れる対数因子 $ (\text{ln} N)^3 $ はどのようにして生じるか。

主な発見

• 離散解は、すべての $ \varepsilon_i \in (0,1] $ に対して、誤差境界 $ ||\vec{U} - \vec{u}|| \leq C N^{-2} (\text{ln} N)^3 $ を満たす2次パrameterに依存しない収束性を達成する。
• パrameter値が著しく異なる場合でも、境界層におけるメッシュの適応的細分化により、本手法はすべてのパrameter値に対して強固である。
• 特異成分の誤差は $ C N^{-2} (\text{ln} N)^3 $ で有界であり、対数因子は層の遷移およびメッシュ分割解析に起因する。
• 証明手法はグリーン関数を避けており、帰納法と成分ごとの推定に依存するため、多次元および非線形系への拡張の可能性が高まる。
• シシュキンメッシュ構造により、最小の $ \varepsilon_i $ に応じてメッシュ間隔が適応されるため、一様収束が保証される。すべての $ b_i = 0 $ の場合、メッシュは均等に還元される。
• 解析により、比 $ R(w_i(x_k)) $ が定数で有界であることが示され、これは遷移点における層の挙動を制御する上で重要である。