[論文レビュー] Well-balanced and asymptotic preserving schemes for kinetic models

Casimir Emako, Min Tang|arXiv (Cornell University)|Mar 10, 2016

Mathematical Biology Tumor Growth参考文献 17被引用数 4

ひとこと要約

本稿では、ε → 0 のときの拡散極限と ε = O(1) のときの定常状態の両方のスケールで、良好にバランスと漸近保存（AP）の両方の性質を持つ数値スキームを設計する一般化されたフレームワークを提案する。この手法は、AP時間発展ステップと局所的定常状態の解法を組み合わせることで、WB特性を強制し、εに依存しない精度で両方の漸近的極限と定常状態を正しく捉えることに成功している。

ABSTRACT

In this paper, we propose a general framework for designing numerical schemes that have both well-balanced (WB) and asymptotic preserving (AP) properties, for various kinds of kinetic models. We are interested in two different parameter regimes, 1) When the ratio between the mean free path and the characteristic macroscopic length $ε$ tends to zero, the density can be described by (advection) diffusion type (linear or nonlinear) macroscopic models; 2) When $ε$ = O(1), the models behave like hyperbolic equations with source terms and we are interested in their steady states. We apply the framework to three different kinetic models: neutron transport equation and its diffusion limit, the transport equation for chemotaxis and its Keller-Segel limit, and grey radiative transfer equation and its nonlinear diffusion limit. Numerical examples are given to demonstrate the properties of the schemes.

研究の動機と目的

• 運動方程式に対して、良好にバランスと漸近保存の両方の性質を同時に満たす統合的数値フレームワークを開発すること。
• ε ≪ 1 のときの拡散極限と、ε = O(1) のときの強い源項の影響が顕著なスケールの両方において、運動方程式を正確にシミュレートする課題に対処すること。
• 特に、標準的なスキームが失敗する小ε領域においても、数値スキームが安定かつ正確に保たれるようにすること。
• ε → 0 のとき、離散化されたレベルで正しい巨視的極限（例：Keller-Segel方程式や非線形拡散方程式）を保持すること。
• 強い源項を有する双曲型系の長時間シミュレーションにおいて、定常状態を正確に捉えることで、高い精度を維持すること。

提案手法

• フレームワークは二段階の手続きを用いる：まず、APスキームを用いて時間依存運動方程式を解き、予測値を得る。
• 次に、予測された解を用いて、各計算セル内での局所的定常状態問題の流入境界条件を設定する。
• 定常状態問題は、巨視的極限と整合性を持つように調整されたAPスキームを用いて解き、その結果をもってして数値フラックスを補正することで、離散化されたレベルでの良好にバランスの性質を強制する。
• 得られたセル単位の流出量を用いて、数値フラックスを補正し、離散化されたレベルでの良好にバランスの性質を強制する。
• この手法は、中性子輸送、化学走化性（Othmer-Dunbar-Alt）、グレー放射線移動の3つのモデルに適用され、詳細な離散化スキームが提示されている。
• このアプローチは柔軟で、既存の任意のAPスキームを時間発展および定常状態問題に組み込むことが可能であり、効率的な実装を可能にする。

実験結果

リサーチクエスチョン

• RQ1 同一の数値フレームワークが、異なるパrameterスケールの運動方程式に対して、良好にバランスと漸近保存の両方の性質を同時に達成できるか？
• RQ2 ε → 0 のとき、離散スキームがεの大きさに依存せずに正しい巨視的極限（例：拡散または非線形拡散）を保持できるか？
• RQ3 ε = O(1) のとき、強い源項を有する双曲型系の定常状態解を、スキームが正確に捉えることができるか？
• RQ4 スキームの収束性は、メッシュサイズΔxおよびパrameter εに依存するか？また、εにわたって一様な収束性を示すか？
• RQ5 非線形系の繰り返し解法に依存せずに、離散化されたレベルでの良好にバランスの性質を強制する方法は何か？

主な発見

• 提案されたスキームは、すべてのεの値においてメッシュサイズΔxにわたる一様な収束性を示し、図7により収束次数が一貫していることから、漸近保存（AP）性質が裏付けられている。
• ε → 0 のときの正しい極限挙動が再現されている：放射強度Iεはψ⁰/2に収束し、ψ⁰ = ac(T⁰)⁴である。図8により、εにおける収束次数は約0.6に近く、良好な収束性が確認されている。
• ε = 1のとき、スキームは正確な定常状態解I(x,v) = |v|xおよびρ = xを正しく捉えており、T = 1における残差が減少し安定化していることから、良好にバランス（WB）の性質が確認されている（図9および10）。
• 数値結果により、ε ≪ 1 の拡散極限とε = O(1) の強い源項領域の両方において、性能劣化なしに精度と安定性を維持していることが確認された。
• このフレームワークは計算的に効率的であり、予測ステップに明示的または単純な陰解法APスキームを用いることができ、非線形系の高価なニュートン型ソルバーを回避できる。
• この手法は一般性に富み、他の運動方程式モデル、特に緩和項を有する双曲型系へも拡張可能であり、マルチスケールシミュレーションへの広範な応用が期待される。