[論文レビュー] A Parameter-Uniform Finite Difference Method for Multiscale Singularly Perturbed Linear Dynamical Systems

Valarmathi Sigamani, John J. H. Miller|ArXiv.org|Mar 10, 2009

Differential Equations and Numerical Methods参考文献 5被引用数 11

ひとこと要約

本稿では、異なる小さなパラメータが異なる解の層を支配するマルチスケール特異摂動線形動的システムに対するパrameterに依存しない有限差分法を提示する。シェィシュキンの非一様メッシュと古典的有限差分法を用いることで、すべての摂動パラメータに対して一様に1次収束が達成され、全誤差境界が $ O(N^{-1} ln N) $ となる。これは、すべてのスケールで頑健な精度を保証する。

ABSTRACT

A system of singularly perturbed ordinary differential equations of first order with given initial conditions is considered. The leading term of each equation is multiplied by a small positive parameter. These parameters are assumed to be distinct and they determine the different scales in the solution to this problem. A Shishkin piecewise--uniform mesh is constructed, which is used, in conjunction with a classical finite difference discretization, to form a new numerical method for solving this problem. It is proved that the numerical approximations obtained from this method are essentially first order convergent uniformly in all of the parameters.

研究の動機と目的

• 異なる小さなパラメータを有するマルチスケール特異摂動線形初期値問題に対する頑健な数値解法の開発。
• すべての摂動パラメータに対して収束が一様であることを保証し、パラメータの変動に伴う精度の低下を回避する。
• 異なるスケールに起因する複数の境界層を的確に扱える離散スキームの構築。
• 特異摂動パラメータの大きさに依存しない、厳密な誤差推定の確立。

提案手法

• 各異なる小さなパラメータ $ \varepsilon_i $ に起因する明確な境界層を解像できるように、シェィシュキンの非一様メッシュを構築する。
• シェィシュキンメッシュ上での古典的有限差分法を用いて系を離散化し、標準的手法で解ける線形方程式系を形成する。
• 離散的最大原理と安定性を確立することで誤差解析を支援し、離散解の有界性を保証する。
• 局所切り捨て誤差を滑らか成分と特異成分に分解し、数学的帰納法と層関数解析を用いてそれぞれに鋭い推定値を導出する。
• パラメータの順序付けに基づくすべての部分区間およびメッシュタイプにおいて、層の挙動を分析することで、局所切り捨て誤差に対するパrameterに依存しない境界を導出する。
• 滑らか成分と特異成分の境界を組み合わせることで、最終的な全誤差推定を得る。これにより、パrameterに依存しない一様収束結果が得られる。

実験結果

リサーチクエスチョン

• RQ1 複数の異なるスケールを有する系に対して、すべての特異摂動パラメータに対して一様収束を達成できる有限差分法を設計できるか？
• RQ2 異なる $ \varepsilon_i $ パラメータに起因するすべての層を同時に解像できるように、メッシュをどのように階層化すべきか？
• RQ3 このようなマルチスケール問題に対して、フィットメッシュ上での標準的有限差分スキームが達成可能な最良の収束次数は何か？
• RQ4 滑らか成分と特異成分の両方について、パrameterに依存しない鋭い誤差境界を導出できるか？

主な発見

• 提案手法は1次収束をパrameterに依存せず達成し、全誤差が $ C N^{-1} \nln N $ で抑えられる。これは摂動パラメータ $ \varepsilon_i $ の値に依存しない。
• 特異成分の局所切り捨て誤差は、すべての部分区間およびパラメータの順序付けに基づくすべての可能なメッシュタイプにおいて $ O(N^{-1} \nln N) $ で抑えられる。
• 解析により、最小のパラメータ $ \varepsilon_1 $ がゼロに近づいても、手法が安定かつ正確に保たれることを示した。これにより、頑健性が保証される。
• 適切に層をフィットさせたシェィシュキンメッシュの使用により、すべての層が解像され、遷移点 $ \sigma_i $ は $ \varepsilon_i / \alpha \nln N $ に基づいて選択される。
• バイナリーベクトル $ \vec{b} $ で定義される $ 2^n $ 通りのすべての可能なメッシュ構成に対して、誤差境界は鋭く一様に保たれる。
• 最大原理が安定性解析において中心的な役割を果たすとともに、帰納法と成分ごとの推定を用いた厳密な証明により理論的結果が裏付けられる。