[論文レビュー] Detonation propagation in weakly confined gases

Youssef K. Wahba, XiaoCheng Mi+2 more|arXiv (Cornell University)|Mar 25, 2026

Combustion and Detonation Processes被引用数 0

ひとこと要約

この論文は、CFDと二層ガスモデルを用いて、音響インピーダンスと層厚さに基づく弱く閉塞された構成における着 detonations behaviorを分類し、 precursor shocks、underdriven/overdriven regimes、および位相図を予測する解析フレームワークを開発する。

ABSTRACT

This study investigates the propagation of detonations along a layered configuration where a reactive gas is weakly confined by a hotter inert layer. CFD simulations are performed using a single-step, non-Arrhenius reaction model designed to suppress cellular instabilities, enabling formulation of a theoretical framework directly compared with simulation results. The simulations reach a quasi-steady state, revealing distinct flowfield regimes that depend on the acoustic-impedance ratio and relative layer thicknesses, with some detonations exhibiting velocity deficits while others propagate above the ideal Chapman-Jouguet (CJ) speed. Analytical models are developed to interpret these regimes. When a precursor shock is observed in the inert layer, the detonation is overdriven; this is modeled using shock-polar analysis and velocity estimates based on the approach of Mitrofanov (Acta Astronaut. 3:995-1004, 1976). An analytical criterion for precursor shock onset is proposed. In underdriven scenarios, the detonation front exhibits positive curvature, analyzed using a geometric construction wherein the relationship between wave speed and front curvature is evaluated a priori. A simplified characteristic-based model captures the decay of the shock wave in the inert layer, after which shock-polar analysis determines the resulting wave interaction. Predictions from these models are assembled into a phase map delineating regions of overdriven and underdriven behavior, along with corresponding shock interactions, in the space of acoustic impedance and area ratios. This map is compared directly with CFD results. The combined numerical-theoretical framework clarifies transition mechanisms governing layered detonations and provides insights into detonation dynamics relevant to rotating detonation engines in which the detonation is bounded by hotter combustion products from a previous cycle.

研究の動機と目的

• hotter inert layer で弱く閉塞された二層チャネル内で反応性ガスの着 detonations propagates する様子を調べる。
• 音響インピーダンス比と層厚さに依存する流れ場の regime を特定する。
• CFD で観察された regime を解釈し、過駆動と過小駆動 detonations の遷移を予測する解析モデルを開発する。

提案手法

• 非 Arrhenius 型の単一ステップモデルを用いてセルラー不安定性を抑制し、弱い閉塞下の反応性層の 2D CFD を実行する。
• 変数を無次元化し、 adjustable A2/A1 および音響インピーダンス比 Z を持つ二層チャネルを設定する（いずれの場合も Z<1）。
• Strang 分割による化学反応源項と共に、MUSCL-Hancock スキームと HLLC リーマンソルバーを用いて反応性オイラー方程式を解く。
• hotspot により detonations を開始し、準定常状態まで実行して detonation 速度とショック構造を測定する。
• 観察された regime を解釈するための解析構成（Kantrowitz 型 choking criterion、Eyring 幾何構成、極解析、および Mitrofanov 二層法）を開発する。

Figure 1 : Schematic representation of the problem, showing terminal shock structures relative to the detonation front: (a) attached shock (behind the front); (b) precursor shock (ahead of the front).

実験結果

リサーチクエスチョン

• RQ1 不活性層における precursor shock の onset criterion は area ratio A2/A1 および音響- impedan ce ratio Z の関数としてどうなるか？
• RQ2 underdriven と overdriven detonations は不活性層の流れ場と反応性層の相互作用においてどう現れるか？
• RQ3 Z および A2/A1 の関数として、付着・ precursor・ detached shocks の regime を区分する解析的 phase map を構築できるか？
• RQ4 前方曲率とショックの相互作用は CJ 速度に対する速度 deficit や overdriving をどのように説明するか？

主な発見

• CFD シミュレーションは Z と A2/A1 に依存する 5 種の detonations behavior（Case A–E）を明らかにし、 regular または Mach 反射を伴う underdriven、反応性層における precursor shocks が MR/RR を駆動、および detachment された shocks を含む。
• Kantrowitz 型 choking criterion を導出し、 area ratio と impedance を M2 = M_CJ·Z に関連付け、臨界的な A2/A1 の関係を導出して precursor-shock onset を予測する。
• Underdriven detonations（Z > Z_Kant）は反応性層と不活性界面から発する正の前方曲率と sonic locus を示し、 front curvature は Eyring construction を用いて解析。
• Overdriven ケース（κ<0 の precursor shocks）は Mitrofanov 二層アプローチと極解析を用いてショック– detonations 相互作用を予測。
• 位相図を構築し、overdriven/underdriven の挙動とショック相互作用を Z および A2/A1 に結びつけ、CFD 結果で検証する。

Figure 2 : Initial conditions showing abrupt initiation via a high-pressure, high-temperature region.