[論文レビュー] Efficient Implementation of ADER Schemes for Euler and Magnetohydrodynamical Flows on Structured Meshes -- Comparison with Runge-Kutta Methods

Dinshaw S. Balsara, Chad Meyer+3 more|arXiv (Cornell University)|Jun 10, 2010

Computational Fluid Dynamics and Aerodynamics参考文献 47被引用数 5

ひとこと要約

本稿では、構造化メッシュ上でのEuler方程式および磁気流体力学（MHD）方程式を解くための、実用的で効率的なADER（任意の微分）スキームの実装を提示している。1段階の時間積分法を採用しており、従来の強安定性を保証するルンゲ＝クッタ法よりも優れた性能を発揮する。著者らはノード空間およびモーダル空間における詳細な定式化を提供し、発散自由MHDのための数値フラックスおよび電場計算を最適化しており、全テストされた次数においてADERスキームがルンゲ＝クッタ法のほぼ2倍の速度で動作することを示している。

ABSTRACT

ADER (Arbitrary DERivative in space and time) methods for the time-evolution of hyperbolic conservation laws have recently generated a fair bit of interest. The ADER time update can be carried out in a single step, which is desirable in many applications. However, prior papers have focused on the theory while downplaying implementation details. The purpose of the present paper is to make ADER schemes accessible by providing two useful formulations of the method as well as their implementation details on three-dimensional structured meshes. We therefore provide a detailed formulation of ADER schemes for conservation laws with non-stiff source terms in nodal as well as modal space along with useful implementation-related detail. We also provide details for the efficient use of ADER schemes in obtaining the numerical flux for conservation laws as well as electric fields for divergence-free magnetohydrodynamics. An efficient WENO-based strategy for obtaining zone-averaged magnetic fields from face-centered magnetic fields in MHD is also presented. The schemes catalogued here have been implemented in the first author's RIEMANN code. The speed of ADER schemes is shown to be almost twice as fast as that of strong stability preserving Runge-Kutta time stepping schemes for all the orders of accuracy that we tested.

研究の動機と目的

• 計算流体力学における理論的ADERSスキームと実装の間のギャップを埋めること。
• 構造化メッシュ上におけるノード空間およびモーダル空間におけるADERSメソッドの詳細で実用的な定式化を提供すること。
• 発散自由MHDシミュレーションのための数値フラックスおよび電場計算を最適化すること。
• 面中心データからゾーン平均磁束密度を求めるためのWENOに基づく戦略を提示すること。
• 強安定性を保証するルンゲ＝クッタ法と比較して、ADERSスキームの計算的優位性（速度および精度の面で）を実証すること。

提案手法

• ADERスキームは1段階の時間更新を用いて実装されており、ルンゲ＝クッタ法の段階的進行を回避している。
• 本手法は、3次元構造化メッシュ上での非剛性源項を伴う双曲型保存則を解くためにノードおよびモーダル定式化を採用している。
• 保存則に適した高次精度リーマンソルバを用いて、効率的な数値フラックスを計算している。
• MHDにおける電場は、発散自由な磁束密度を維持するためにフラックス積分から直接計算されている。
• ゾーン平均磁束密度を面中心の磁束密度データから求めるために、WENOに基づく再構成戦略が用いられている。
• すべてのスキームはRIEMANNコードに実装され、高次精度および性能評価が可能になっている。

実験結果

リサーチクエスチョン

• RQ1 計算流体力学における構造化メッシュ上でのADERSスキームを実用的に効率的に実装するにはどうすればよいか？
• RQ2 EulerおよびMHD方程式のADERSスキームで高い性能を発揮するための鍵となる実装詳細は何か？
• RQ3 異なる次数の精度において、ADERSスキームの計算コストは強安定性を保証するルンゲ＝クッタ法と比べてどの程度か？
• RQ4 ADERSに基づくMHDシミュレーションで発散自由な磁束密度を保証するための戦略は何か？
• RQ5 ADERSスキームは、高次精度を維持しつつ、ルンゲ＝クッタ法と比較して計算速度をほぼ2倍にできるか？

主な発見

• 全テストされた次数において、ADERSスキームは強安定性を保証するルンゲ＝クッタ法と比較して、ほぼ2倍の計算速度を達成している。
• ノード空間およびモーダル空間におけるADERSスキームの実装により、構造化メッシュ上での効率的で高次精度の解法が可能になっている。
• WENOに基づく再構成戦略の使用により、面中心データからのゾーン平均磁束密度が正確かつ安定に得られている。
• 数値フラックスおよび電場の効率的計算は、MHDにおける精度および発散自由制約の維持に不可欠である。
• ADERS法の1段階時間積分は、多段階のルンゲ＝クッタ法と比較して顕著な性能優位性を提供している。
• RIEMANNコードにおける実装により、提案されたADERS定式化の頑健性およびスケーラビリティが、複雑な流体力学問題に対して確認された。