[논문 리뷰] Detonation propagation in weakly confined gases

Youssef K. Wahba, XiaoCheng Mi외 2명|arXiv (Cornell University)|2026. 03. 25.

Combustion and Detonation Processes인용 수 0

한 줄 요약

본 논문은 음향 임피던스와 층 두께에 의해 약하게 제약된 구성에서 연소 폭발 행동을 분류하기 위해 두 층 가스 모델을 사용한 CFD를 이용하고, 예측 선구 충격, 저구동/과구동 상태, 그리고 위상 맵을 예측하기 위한 해석적 프레임워크를 개발한다.

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.

연구 동기 및 목표

• 두 층 채널에서 더 뜨거운 비활성 층에 의해 약하게 제약될 때 반응성 기체의 폭발이 어떻게 전파되는지 조사한다.
• 음향 임피던스 비와 층 두께에 의존하는 흐름장(flowfield) 영역을 식별한다.
• CFD로 관찰된 유동 거동을 해석하고 과구동(overdriven)과 저구동(underdriven) 폭발 사이의 전이를 예측하기 위한 해석 모델을 개발한다.

제안 방법

• 단일 스텝, 비 Arrhenius 모델을 사용하여 형성되는 세포 구조를 억제하기 위해 비활성 층에 의해 약하게 구속된 반응성 층의 2D CFD 시뮬레이션을 수행한다.
• 변수를 무차원화하고 조정 가능한 A2/A1 및 음향 임피던스 비 Z(모든 경우에서 Z<1)로 두 층 채널을 설정한다.
• 화학적 소스 항에 Strang 분할을 적용한 반응성 Euler 방정식을 풀기 위해 MUSCL-Hancock 스킴과 HLLC Riemann 솔버를 사용한다.
• 핫스팟을 통해 폭발을 시작하고 준정상 상태에 이를 때까지 실행하여 폭발 속도와 충격 구조를 측정한다.
• 관찰된 거동을 해석하기 위해 Kantrowitz-type 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 A2/A1 면적 비와 음향 임피던스 비 Z의 함수로서 비활성 층에서 선구 충격의 발현 임계조건은 무엇인가?
• RQ2 저구동 폭발과 과구동 폭발이 비활성 층의 흐름장과 반응성 층 상호작용에서 어떻게 나타나는가?
• RQ3 Z와 A2/A1의 함수로 서로 붙은(Attached), 선구(Precursor), 분리된(Detached) 충격의 거동 영역을 구분하는 해석적 위상 맵을 구성할 수 있는가?
• RQ4 전면 곡률과 충격 상호작용이 CJ 속도에 비해 속도 결손이나 과구동을 설명하는가?

주요 결과

• CFD 시뮬레이션은 Z와 A2/A1에 의존하는 다섯 가지 폭발 거동(Cases A–E)을 나타내며, 규칙적 또는 Mach 반사와 함께 저구동인 경우, 선구 충격이 반응성 층에서 MR/RR을 유발하는 경우, 그리고 분리된 충격이 있는 경우를 포함한다.
• Kantrowitz-type choking criterion이 면적 비와 임피던스를 M2 = M_CJ·Z로 연결하고 임계 A2/A1에 대한 유도 관계를 통해 선구충격 발현을 예측하기 위해 도출된다.
• 저구동 폭발(Z > Z_Kant)은 양의 전면 곡률과 반응층-비활성 층 경계에서 시작하는 음속 등위를 보이며, 전면 곡률은 Eyring 구성으로 분석된다.
• 과구동 경우(κ<0인 선구 충격)은 Mitrofanov 두 층 접근법과 극 분석을 사용하여 충격-폭발 상호작용을 예측한다.
• 위상 맵은 과구동/저구동 거동과 충격 상호작용을 Z 및 A2/A1에 연결하여 구성되고, CFD 결과와 대조하여 검증된다.

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