[論文レビュー] Implicit-explicit, realizability-preserving first-order scheme for moment models with Lipschitz-continuous source terms

Florian Schneider|arXiv (Cornell University)|Nov 4, 2016

Gas Dynamics and Kinetic Theory参考文献 30被引用数 3

ひとこと要約

本稿では、Lipschitz連続な源項を伴うモーメントモデルに対して、実現可能性を保つ一次精度の陰-陽（IMEX）数値スキームを提示する。特に、最小エントロピー（M_N）モデルに応用する。剛性のある源項を陰的に、非剛性のフラックスを陽的に扱うことで、物理的パラメータに依存する時間刻みの制限を除去し、安定かつ効率的なシミュレーションを可能にするとともに、モーメントベクトルが実現可能（非負の運動論的分布と整合）であることを保証する。

ABSTRACT

We derive an implicit-explicit (IMEX), realizability-preserving first-order scheme for moment models with Lipschitz-continuous source terms. In contrast to fully-explicit schemes the time step does not depend on the physical parameters, removing the stiffness from the system. Furthermore, a wider class of collision operators (e.g. the Laplace-Beltrami operator) can be used. The derived scheme is applied to minimum-entropy models.

研究の動機と目的

• 剛性のある源項を伴う最小エントロピー（M_N）モデルにおいて、モーメントベクトルの実現可能性を維持する数値スキームの開発。
• 吸収係数や散乱係数などの物理的パラメータに依存する完全に陽的なスキームに見られる時間刻みの制限を克服すること。
• ラプラシアン＝ベルトラミ作用素を含む、より広いクラスの衝突作用素のモーメントモデルへの適用を可能にすること。
• 剛性のある源項を陰的に取り扱うことで、実現可能性を保ちつつ一次精度を維持するスキームの構築。
• 今後の研究において、完全に陽的な高次精度の運動論的および不連続ガレルキン法への高次IMEXスキームへの拡張の基盤を提供すること。

提案手法

• 剛性のある源項を陰的に、非剛性の双曲型フラックスを陽的に扱う陰-陽（IMEX）時間離散化を用いる。
• 空間均一なモーメント系に実現可能性を保つ陰離散化を適用し、再構築された運動論的密度の非負性を保証する。
• 各空間時間セル内で陰方程式系を局所的に解くことで、並列処理と計算効率を実現する。
• 衝突作用素が質量保存、局所的エントロピー散逸、時間ステップにおける非負性保存を満たすものと仮定する。
• 必要な仮定を満たし、より広範な適用性を可能にする代表的な衝突作用素としてラプラシアン＝ベルトラミ作用素を用いる。
• 理論的解析と数値的検証（製造された解とベンチマーク平面源問題を含む）を通じて、実現可能性の保存を検証する。

実験結果

リサーチクエスチョン

• RQ1 剛性があり、Lipschitz連続な源項を伴うモーメントモデルにおいて、実現可能性を保つIMEXスキームを構築可能か？
• RQ2 源項を陰的に取り扱うことで、完全に陽的なスキームに内在する時間刻みの制限が解消され、安定性が損なわれないか？
• RQ3 実現可能性を維持しつつ、ラプラシアン＝ベルトラミ作用素のようなより広いクラスの衝突作用素に適用可能か？
• RQ4 実際のシミュレーションにおいて、高次M_Nモデルにおいて、実現可能性境界にどれほど近づけるか？
• RQ5 理論的枠組みの仮定が実際には正確に満たされない場合でも、実現可能性を保つ性質は頑健か？

主な発見

• 標準的なCFL条件のもとで、剛性のある源項（強い吸収や散乱を含む）が存在しても、IMEXスキームは実現可能性を保つ。
• 時間刻みが物理的パラメータに依存しなくなったため、完全に陽的なスキームに比べてより大きな時間刻みが可能となり、計算がより効率的になる。
• 製造された解に対する収束テストにより、スキームが一次精度を達成していることが確認された。
• 平面源ベンチマークにおいて、解は実現可能性境界に近く保たれ、相対的距離はM₁で0.0039、M₂で2.2147×10⁻⁵、M₃で9.0981×10⁻⁷にとどまり、高次モデルではモーメント問題の解法が難しくなる傾向が示された。
• 位相空間におけるM₂およびM₃解のヒストグラムから、粒子が実現可能集合の下部に集中していることが確認され、境界から遠く離れた領域を避ける傾向が示され、スキームの頑健性が裏付けられた。
• 理論的仮定が完全に満たされない場合（例：ラプラシアン＝ベルトラミ作用素を用いた混合モーメントモデル）でも、本手法は良好に動作した。これは、仮定を弱めるか、より広いクラスのモデルに拡張する可能性を示唆している。