[論文レビュー] Uniform convergence on a Bakhvalov-type mesh using the preconditioning approach: Technical report

Thái Anh Nhan, Relja Vulanović|arXiv (Cornell University)|Apr 16, 2015

Differential Equations and Numerical Methods参考文献 9被引用数 3

ひとこと要約

本稿では、バフカロフ型メッシュ上での特異摂動型対流拡散境界値問題に対して一様収束を達成するための前処理手法を提示する。離散化された系に前処理子を適用することにより、摂動パラメータ ε に依存しない条件数が得られ、境界層が存在する場合でも ε に対して一様に O(N⁻¹) の最適順序の点ごとの収束が可能になる。

ABSTRACT

The linear singularly perturbed convection-diffusion problem in one dimension is considered and its discretization on a Bakhvalov-type mesh is analyzed. The preconditioning technique is used to obtain the pointwise convergence uniform in the perturbation parameter.

研究の動機と目的

• 特異摂動型対流拡散問題における境界層を伴う一様収束性を有する数値解法の開発。
• 摂動パラメータ ε が小さい場合の精度と安定性を維持する課題に対処すること。
• バフカロフ型メッシュ上での離散化系に前処理を適用し、ε に対して一様収束を保証すること。
• 前処理系の条件数が ε に依存せずに有界であることを証明すること。
• 離散解が連続解に対して、ε に依存しない最適順序の点ごとの収束を確立すること。

提案手法

• 一次元の特異摂動型対流拡散問題をバフカロフ型メッシュ上で有限差分上流差分スキームで離散化する。
• 離散化系の条件を改善するために、スケーリングを行う前処理子行列を定義する。
• 前処理子を用いて元の系を、条件数が有界である新しい系に変換する。
• 一貫性・安定性の原則を適用し、一貫性誤差が O(N⁻¹) であり、前処理系が一様に有界な逆行列を持つことを証明する。
• 解を滑らかな部分と境界層部分に分解し、メッシュの異なる領域における切り捨て誤差を推定する。
• 離散最大原理を用いて、有界な条件数と一貫性誤差の推定値を組み合わせ、一様収束を確立する。

実験結果

リサーチクエスチョン

• RQ1 前処理を用いることで、バフカロフ型メッシュ上での対流拡散問題に対して一様収束を達成できるか？
• RQ2 前処理を適用した場合、離散化系の条件数が摂動パラメータ ε に依存せずに有界のままであるか？
• RQ3 バフカロフ型メッシュ上での有限差分スキームの一貫性誤差は、ε に対して一様に O(N⁻¹) であるか？
• RQ4 このアプローチを用いることで、ε に依存しない最適順序 O(N⁻¹) の精度で解を近似できるか？
• RQ5 前処理子は、境界層が存在する状況での離散化系の安定性および条件数にどのように影響を与えるか？

主な発見

• 前処理系の条件数は、摂動パラメータ ε に依存せず、O(N²) で有界である。
• バフカロフ型メッシュ上での有限差分スキームの一貫性誤差は、すべての ε に対して一様に O(N⁻¹) で有界である。
• a > 4/β を満たす限り、離散解は正確な解に対して点ごとの誤差 O(N⁻¹) で収束し、ε に依存しない。
• 前処理子は、ε が小さい場合に生じる悪条件化を効果的に軽減し、収束のロバスト性を保証する。
• 解析により、バフカロフ型メッシュと前処理の組み合わせが、急峻な境界層が存在する場合でも最適収束率を達成することが確認された。
• メッシュと前処理子の設計のおかげで、境界層領域に特別な処理を施さなくても一様収束が達成される。