[論文レビュー] A parameter uniform fitted mesh method for a weakly coupled system of two singularly perturbed convection-diffusion equations

Saravana Sankar Kalaiselvan, John J. H. Miller+1 more|arXiv (Cornell University)|Sep 24, 2018

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

ひとこと要約

本稿では、異なる摂動パラメータを有する弱く結合された2つの特異摂動対流拡散方程式系に対して、パrameterに依存しない適合メッシュ法を提示する。一様なシャイシンメッシュ上での上流差分スキームを用いることで、最大ノルムにおいて両方の摂動パラメータに対して一様に1次収束を達成し、解析的誤差推定と数値例により、異なるパrameter値の下でも堅牢な収束が確認されている。

ABSTRACT

In this paper, a boundary value problem for a singularly perturbed linear system of two second order ordinary differential equations of convection- diffusion type is considered on the interval [0, 1]. The components of the solution of this system exhibit boundary layers at 0. A numerical method composed of an upwind finite difference scheme applied on a piecewise uniform Shishkin mesh is suggested to solve the problem. The method is proved to be first order convergent in the maximum norm uniformly in the perturbation parameters. Numerical examples are provided in support of the theory.

研究の動機と目的

• 異なる小さな摂動パラメータを有する2つの特異摂動対流拡散方程式系を解くための堅牢な数値法の開発。
• 弱い結合性と異なる摂動パラメータによる解成分における境界層の幅が異なるという課題の取り扱い。
• 特に層を多く含む解の下でも、両方の摂動パラメータに対して一様な収束を保証すること。
• 解の分解と微分係数の推定に関する厳密な解析的枠組みの構築により、先行研究を改善し、より明確な定義を提供すること。

提案手法

• 2階の特異摂動対流拡散方程式系を、区分的均等なシャイシンメッシュ上で上流差分スキームを適用して離散化する。
• 境界層を x=0 に有するシャイシンメッシュを構築し、より小さい摂動パラメータ ε₁ とより大きい ε₂ に適合したメッシュグレーディングを採用する。
• 解を正則部と層部に分解し、層関数は最大原理と安定性推定を用いて分析する。
• 離散問題に対して最大原理を確立し、数値スキームの正値性と安定性を保証する。
• 本手法が、ε₁ および ε₂ に対して最大ノルムにおいて一様に1次収束することを証明する。
• ベンチマーク問題に対する数値実験により、理論的収束性が検証される。

実験結果

リサーチクエスチョン

• RQ1 異なる摂動パラメータを有する弱く結合された2つの特異摂動対流拡散方程式系に対して、シャイシンメッシュ上での適合有限差分法がパラメータに依存しない収束を達成できるか。
• RQ2 結合性がある中で、解成分における境界層構造は個々の摂動パラメータ ε₁ と ε₂ にどのように依存するか。
• RQ3 先行研究における制限的仮定 ε₁ ≤ C N⁻¹ および ε₂ ≤ C N⁻¹ を除去することで、収束解析および手法設計にどのような影響が生じるか。
• RQ4 解の分解をどのように改善すれば、各摂動パラメータが層挙動に与える寄与をより正確に捉えることができるか。
• RQ5 実際の収束速度と誤差定数はどのように変化するか。また、ε₁ と ε₂ が0に近づく際にも安定しているか。

主な発見

• 提示されたシャイシンメッシュ上での上流差分スキームは、最大ノルムにおいて両摂動パラメータ ε₁ および ε₂ に対して一様に1次収束を達成する。
• 計算された収束次数は約 0.483 に安定し、すべてのテストパラメータ値において堅牢な収束が示された。
• 計算された (ε₁,ε₂)-一様誤差定数は N ≥ 128 の範囲で約 2.755 に安定し、パラメータに依存しない一様性が確認された。
• 数値結果から、誤差は N の増加とともに減少し、ε₁ と ε₂ の減少に伴い安定することが示され、理論的収束挙動が妥当であることが裏付けられた。
• 本手法は期待される層パターンを正しく捉えている：境界層条件に応じて、u₁ に ε₁-層、両成分に ε₂-層が現れる。
• 正則部と層部への解の分解が改善され、正しく定義されており、解構造に対するより明確な洞察が得られた。