[論文レビュー] Simulations of detonation waves with smoothed dissipative particle dynamics

Gérôme Faure, Jean‐Bernard Maillet|arXiv (Cornell University)|Sep 11, 2017

Gas Dynamics and Kinetic Theory被引用数 3

ひとこと要約

本稿では、化学反応系をモデル化するために、メソスコープ粒子ごとに反応進行度変数 λ を導入し、化学動力学に従って進化するように、Smoothed Dissipative Particle Dynamics (SDPD) を拡張した。この手法は、ニトロメタンにおける衝撃波から爆轟波への遷移および定常的爆轟波を成功裏にシミュレートし、爆轟速度や反応領域の幅といった主要な物理的性質を回復した。解像度が十分な場合、反応速度の前指数因子の変化に対しても結果が一貫しており、良好な性能を示した。

ABSTRACT

Smoothed Dissipative Particle Dynamics (SDPD) is a mesoscopic method which allows to select the level of resolution at which a fluid is simulated. The aim of this work is to extend SDPD to chemically reactive systems.To this end, an additional progress variable is attached to each mesoparticle and evolves according to chemical kinetics. This reactive SDPD model is illustrated with numerical studies of the shock-to-detonation transition in nitromethane as well as the stationary behavior of the reactive wave.

研究の動機と目的

• メソスコープスケールでの化学反応系をシミュレートできるように、SDPD手法を拡張すること。
• ニトロメタンのような反応性物質における爆轟波および衝撃波から爆轟波への遷移（STD）をシミュレート可能にする。
• 反応速度の前指数因子が爆轟波の構造およびダイナミクスに与える影響を評価すること。
• 粒子解像度が反応的SDPDシミュレーションの精度に与える影響を評価すること。
• 従来のMDやDPDEと比較して、粗い解像度でも爆轟波の主要な物理的性質を捉えられるかどうかを示すこと。

提案手法

• 反応系における化学的変換を追跡するために、メソスコープ粒子ごとに進行度変数 λ を導入する。
• DPDE（Maillet et al., 2007）の反応メカニズムをSDPDフレームワークに適応する。
• 温度依存の反応速度を有するアレニウス則を用いて、λ の時間的進化を支配する。
• 爆轟波の長時間ダイナミクスをシミュレートするために、移動ウィンドウ法を適用する。
• エネルギー保存則と流体力学的力とを満たす熱力学的に整合性のあるSDPD定式化を実装する。
• 滑らかな粒子補間を用いてナビエ＝ストークス方程式を解く有限個のメソスコープ粒子を用いてシミュレーションを実行する。

実験結果

リサーチクエスチョン

• RQ1 SDPDは、メソスケールにおける化学動力学を扱う反応系に効果的に拡張可能か？
• RQ2 アレニウスの前指数因子は、ニトロメタンにおける爆轟波の構造およびダイナミクスにどのように影響するか？
• RQ3 粒子解像度が、高速反応における反応的SDPDシミュレーションの精度に及ぼす影響はどの程度か？
• RQ4 反応的SDPDモデルは、衝撃波から爆轟波への遷移および定常的爆轟波の性質を正しく再現できるか？
• RQ5 解像度が十分な場合、反応速度の変化に対してスケール不変性を示せるか？

主な発見

• 解像度が十分な場合、全テスト済みの前指数因子に対して爆轟速度が理論的予測値から1%以内に収まる。
• 反応領域の幅は前指数因子の逆数に比例してスケーリングしており、時間および長さスケールの再スケーリングを示している。
• 前指数因子が 4×10¹⁵ s⁻¹ の場合、粗い解像度を用いると爆轟速度が15%ずれるため、高速動力学のシミュレーションには解像度が不十分であることが示された。
• 同じ高い前指数因子を用いて解像度 K=10 の高解像度シミュレーションを実行したところ、爆轟速度は6777 m·s⁻¹ となり、理論的期待値に近く、良好な一致を示した。
• 異なる前指数因子に対する圧力プロファイルは、スケーリングを施すことで同一曲線に重畳され、系のスケール不変性を確認した。
• メソスコープ粒子が大きすぎる場合には、反応メカニズムが高速反応を正確にモデル化できず、解像度依存のチューニングの必要性が浮き彫りになった。