Two-Dimensional Numerical Simulation of Detonation Transition with Multi-Step Reaction Model: Effects of Obstacle Height

2018 ◽  
Vol 191 (4) ◽  
pp. 659-675 ◽  
Author(s):  
Ayu Ago ◽  
Nobuyuki Tsuboi ◽  
Edyta Dzieminska ◽  
A. Koichi Hayashi
2014 ◽  
Vol 6 (06) ◽  
pp. 718-731
Author(s):  
Hua-Shu Dou ◽  
Zongmin Hu ◽  
Boo Cheong Khoo ◽  
Zonglin Jiang

AbstractNumerical simulation based on the Euler equation and one-step reaction model is carried out to investigate the process of deflagration to detonation transition (DDT) occurring in a straight duct. The numerical method used includes a high resolution fifth-order weighted essentially non-oscillatory (WENO) scheme for spatial discretization, coupled with a third order total variation diminishing Runge-Kutta time stepping method. In particular, effect of energy release on the DDT process is studied. The model parameters used are the heat release atq= 50,30,25,20,15,10 and 5, the specific heat ratio at 1.2, and the activation temperature atTi= 15, respectively. For all the cases, the initial energy in the spark is about the same compared to the detonation energy at the Chapman-Jouguet (CJ) state. It is found from the simulation that the DDT occurrence strongly depends on the magnitude of the energy release. The run-up distance of DDT occurrence decreases with the increase of the energy release forq= 50 ~ 20, and increases with the increase of the energy release forq= 20 ~ 5. This phenomenon is found to be in agreement with the analysis of mathematical stability theory. It is suggested that the factors to strengthen the DDT would make the detonation more stable, and vice versa. Finally, it is concluded from the simulations that the interaction of the shock wave and the flame front is the main reason for leading to DDT.


Materials ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3553
Author(s):  
Dengwang Wang ◽  
Yong Gao ◽  
Sheng Wang ◽  
Jie Wang ◽  
Haipeng Li

Carbon/Phenolic (C/P), a typical anisotropic material, is an important component of aerospace and often used to protect the thermodynamic effects of strong X-ray radiation. In this paper, we establish the anisotropic elastic-plastic constitutive model, which is embedded in the in-house code “RAMA” to simulate a two-dimensional thermal shock wave induced by X-ray. Then, we compare the numerical simulation results with the thermal shock wave stress generated by the same strong current electron beam via experiment to verify the correctness of the numerical simulation. Subsequently, we discuss and analyze the rules of thermal shock wave propagation in C/P material by further numerical simulation. The results reveal that the thermal shock wave represents different shapes and mechanisms by the radiation of 1 keV and 3 keV X-rays. The vaporization recoil phenomenon appears as a compression wave under 1 keV X-ray irradiation, and X-ray penetration is caused by thermal deformation under 3 keV X-ray irradiation. The thermal shock wave propagation exhibits two-dimensional characteristics, the energy deposition of 1 keV and 3 keV both decays exponentially, the energy deposition of 1 keV-peak soft X-ray is high, and the deposition depth is shallow, while the energy deposition of 3 keV-peak hard X-ray is low, and the deposition depth is deep. RAMA can successfully realize two-dimensional orthotropic elastoplastic constitutive relation, the corresponding program was designed and checked, and the calculation results for inspection are consistent with the theory. This study has great significance in the evaluation of anisotropic material protection under the radiation of intense X-rays.


Sign in / Sign up

Export Citation Format

Share Document