[論文レビュー] The Consistent Multi-fluid Advection method

T. Plewa, Ewald Müller|ArXiv.org|Jul 23, 1998

Ionosphere and magnetosphere dynamics被引用数 13

ひとこと要約

本稿では、質量分率の正規化を強制し、修正された補間プロファイルにより数値拡散を低減することで、マルチフラッド流れを正確に運搬する保存的で高次精度のゴドノフ型スキームである一貫したマルチフラッド運搬（CMA）法を提案する。この手法は保存性と単調性を維持し、特に数値拡散が低減されることで、流体力学的シミュレーションにおける組成の分解能が顕著に向上する—とくに重元素の豊度が数倍変化するタイプII超新星モデルにおいて顕著である。

ABSTRACT

Simple modifications for higher-order Godunov-type difference schemes are presented which allow for accurate advection of multi-fluid flows in hydrodynamic simulations. The constraint that the sum of all mass fractions has to be equal to one in every computational zone throughout the simulation is fulfilled by renormalizing the mass fractions during the advection step. The proposed modification is appropriate for any difference scheme written in conservation form. Unlike other commonly used methods it does not violate the conservative character of the advection method. A new steepening mechanism, which is based on modification of interpolation profiles, is used to reduce numerical diffusion across composition discontinuities. Additional procedures are described, which are necessary to enforce monotonicity. Several numerical experiments are presented which demonstrate the capability of our Consistent Multi-fluid Advection (CMA) method in case of smooth and discontinuous distributions of fluid phases and under different hydrodynamic conditions. It is shown that due to the reduced diffusivity of the proposed scheme the abundance of some heavy elements obtained from hydrodynamic simulations of type II supernova explosions can change by a factor of a few in the most extreme cases.

研究の動機と目的

• マルチフラッド流体力学的シミュレーションにおける数値拡散と質量分率の不一致を解消すること。
• 複数の流体成分を運搬する際、高次精度ゴドノフスキームの保存的性質を維持すること。
• 保存性や単調性を損なわず、組成の不連続界面の分解能を向上させること。
• タイプII超新星のような天体物理学的環境における流体混合の正確なモデル化を可能とすること。
• 任意の保存形式の差分スキームに一般化可能な手法を開発すること。

提案手法

• 運搬ステップ中に再正規化を実行することで、すべての計算領域で質量分率の和が1に等しくなる制約を強制する。
• 補間プロファイルを修正し、組成界面を通過する際の数値拡散を低減する新たな鋭利化メカニズムを構築する。
• 保存形式で記述された任意の高次精度ゴドノフ型有限体積法と互換性がある。
• 解の物理的に不適切な振動を防ぐために、単調性を保証する追加手順を実装する。
• 元のスキームの保存的性質を維持し、人工的ソース項を回避する。
• 滑らかでない流れや不連続な流れの両方を想定した数値実験を通じて、性能を評価する。

実験結果

リサーチクエスチョン

• RQ1 組成の不連続界面を通過する際に、精度を保ちつつ保存的かつ正確にマルチフラッド運搬を実行する方法は何か？
• RQ2 高次精度ゴドノフスキームにどのような修正を加えると、保存性を損なわず数値拡散を低減できるか？
• RQ3 質量分率の再正規化がマルチフラッドシミュレーションの安定性と精度に与える影響は何か？
• RQ4 新しい補間ベースの鋭利化メカニズムは、流体界面の分解能をどの程度向上させるか？
• RQ5 低減された数値拡散は、超新星シミュレーションにおける重元素の予測豊度にどのような影響を与えるか？

主な発見

• CMA法は、シミュレーション領域全体で質量分率の正規化を機械精度まで維持している。
• 修正された補間プロファイルにより、急峻な組成不連続界面を通過する際の数値拡散が顕著に低減された。
• 元のゴドノフスキームの保存的性質を維持し、人工的ソース項を回避している。
• 追加のアルゴリズム的手順により単調性が保証され、物理的に不適切な振動が発生しなくなった。
• タイプII超新星のシミュレーションにおいて、低減された拡散性により極端な場合に重元素の豊度が数倍変化した。
• 数値実験により、滑らかでない流れや不連続な流れを含む多様な流体力学的条件下でも、安定した性能が確認された。