[論文レビュー] Numerical Solution of Cylindrically Converging Shock Waves

M. Shabouei, Reza Ebrahimi+1 more|arXiv (Cornell University)|Feb 8, 2016

Computational Fluid Dynamics and Aerodynamics参考文献 28被引用数 4

ひとこと要約

本稿では、TVD-MUSCL法とロウの近似リーマンソルバーを用いた有限体積法により、円柱座標系におけるEuler方程式を数値的に解き、円柱対称的に収束する衝撃波のシミュレーションを実施する。本研究では、軸近傍における衝撃波のダイナミクスを高精度に解像でき、パブリッシュされた解析的および数値的ベンチマークと照合することで、強い衝撃波収束を安定に扱う能力を示している。

ABSTRACT

The cylindrically converging shock wave was numerically simulated by solving the Euler equations in cylindrical coordinates with TVD scheme and MUSCL approach, using Roe's approximate Riemann solver and super-bee nonlinear limiter. The present study used the in house code developed for this purpose. The behavior of the solution in the vicinity of axis is investigated and the results of the numerical solution are compared with the computed data given by Payne, Lapidus, Abarbanel, and Goldberg, Sod, and Leutioff et al.

研究の動機と目的

• 円柱対称的に収束する衝撃波の挙動を、数値的に安定的かつ高精度にシミュレートすること。
• 特異点および急勾配が生じる軸近傍における数値解法の挙動を調査すること。
• ペイン、ソード、および他の研究者らが発表した解析的および数値的データと照合することで、自作計算コードの妥当性を検証すること。
• TVD-MUSCLスキームにスーパー・ビー制限子を適用した場合の、円柱幾何における強い衝撃波の解像度を評価すること。
• 慣性核融合および高エネルギー密度物理学への応用を想定した、衝撃波収束を研究するための信頼性の高い数値フレームワークを提供すること。

提案手法

• MUSCL再構成を用いた2次精度の有限体積法により、円柱座標系におけるEuler方程式を解く。
• 不連続性の近傍で不適切な振動を防ぐために、全変動減少（TVD）スキームを適用する。
• 数値フラックスを正確に計算し、衝撃波を捉えるために、ロウの近似リーマンソルバーを用いる。
• 数値的拡散を制御し、衝撃波フロントを鋭く保つために、スーパー・ビー非線形制限子を採用する。
• 円柱衝撃波シミュレーションに特化した自作計算コードを実装する。
• 特異性を扱い、保存則を維持するために、軸上で適切な境界条件を適用する。

実験結果

リサーチクエスチョン

• RQ1 TVD-MUSCL有限体積スキームは、軸近傍における円柱対称的収束衝撃波のダイナミクスをどの程度正確に捉えることができるか？
• RQ2 スーパー・ビー制限子を用いたロウのリーマンソルバーは、円柱幾何における強い衝撃波フロントの解像度において、どのような性能を示すか？
• RQ3 数値結果は、ペイン、ラピダス、アバーバネル、ゴールドバーグ、ソード、リューティオフらが発表した既存のベンチマーク解とどの程度一致するか？
• RQ4 軸での衝撃波収束をシミュレートする際に生じる数値的課題は何か？ 本研究のアプローチでは、それらはどのように緩和されているか？
• RQ5 自作コードは、既知の収束衝撃波の物理的挙動をどの程度正確に再現しているか？

主な発見

• 数値解法は、幾何的特異性が存在するにもかかわらず、衝撃波の収束を安定に捉え、軸近傍での安定性を維持している。
• ペイン、ラピダス、アバーバネル、ゴールドバーグ、ソード、リューティオフらのベンチマークデータと良好な一致を示し、シミュレーションの正確性が裏付けられた。
• スーパー・ビー制限子の適用により、数値的拡散が効果的に低減され、衝撃波フロントが鋭い形状を保っている。
• TVD-MUSCLアプローチにより、単調性が保たれ、強い勾配が存在する状況でも物理的に不適切な振動が発生しない。
• 自作コードは、円柱座標系における複雑な衝撃波ダイナミクスをシミュレートする上で、強固で信頼性が高いことを示している。
• シミュレーション結果は、圧力および密度の軸上での増幅という、期待される物理的挙動を確認している。