• Journal of the Korean Society of Propulsion Engineers
• /
• v.12 no.3
• /
• pp.68-75
• /
• 2008

The three-dimensional structure of detonation wave propagating in a circular tube was investigated using a parallel computational code developed previously. A series of parametric study for a circular tube of a fixed diameter gave the formation mechanism of the detonation cell structures depending on pre-exponential factor, k. The unsteady results in three-dimension showed the mechanisms of two, three and four cell mode of detonation wave front structures. The detonation cell number was increased but cell width and length were decreased with increased pre-exponential factor k. In the all multi-cell mode, the detonation wave structure and smoked-foil records on the wall are made by the moving of transverse waves. The detonation wave front structures have the regular polygon and windmill shapes periodically.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• 한국연소학회:학술대회논문집
• /
• 2006.04a
• /
• pp.205-212
• /
• 2006

Three-dimensional numerical study was carried out for the investigation of the detonation wave structures propagating in tubes. Fluid dynamics equations and conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The variable gas properties between unburned and burned states were considered by using variable specific heat ratio formulation. The unsteady computational results in three-dimension show the detailed mechanisms of rectangular and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record. The results for the small reaction constant shows the spinning mode of three-dimensional detonation wave dynamics, which was rarely observed in the previous numerical simulation of the detonation waves.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.431-434
• /
• 2005

A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the cell structure of detonation wave by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for capturing cell structures of the different regimes of the detonation phenomena.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.232-235
• /
• 2007

A comprehensive numerical study was carried out to identify the on-set condition of the cell structures of oblique detonation waves (ODWs). Mach 7 incoming flow was considered with all other flow variables were fixed except the flow turning angles varying from 35 to 38. For a given flow conditions theoretical maximum turning angle is $38.2^{\circ}$ where the oblique detonation wave may be stabilized. The effects of grid resolution were tested using grids from $255{\times}100$ to $4,005{\times}1,600$. The numerical smoked foil records exhibits the detonation cell structures with dual triple points running opposite directions for the 36 to 38 turning angles. As the turning angle get closer to the maximum angle the cell structures gets finer and the oscillatory behavior of the primary triple point was observed. The thermal occlusion behind the oblique detonation wave was observed for the $38^{\circ}$ turning angle.

• Nuclear Engineering and Technology
• /
• v.49 no.7
• /
• pp.1423-1430
• /
• 2017

Detonation cell width is an important parameter in hydrogen explosion assessments. The experimental data on gas detonation are statistically analyzed to establish a universal method to numerically predict detonation cell widths. It is commonly understood that detonation cell width, ${\lambda}$, is highly correlated with the characteristic reaction zone width, ${\delta}$. Classical parametric regression methods were widely applied in earlier research to build an explicit semiempirical correlation for the ratio of ${\lambda}/{\delta}$. The obtained correlations formulate the dependency of the ratio ${\lambda}/{\delta}$ on a dimensionless effective chemical activation energy and a dimensionless temperature of the gas mixture. In this paper, support vector regression (SVR), which is based on nonparametric machine learning, is applied to achieve functions with better fitness to experimental data and more accurate predictions. Furthermore, a third parameter, dimensionless pressure, is considered as an additional independent variable. It is found that three-parameter SVR can significantly improve the performance of the fitting function. Meanwhile, SVR also provides better adaptability and the model functions can be easily renewed when experimental database is updated or new regression parameters are considered.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.10 no.2
• /
• pp.1-14
• /
• 2006

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma} $ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.45-52
• /
• 2008

In order to investigate the initiation and propagation processes of a spherical detonation wave induced by direct initiation, numerical simulations were carried out using two-dimensional compressible Euler equations with an axisymmetric assumption and a one-step reaction model based on Arrhenius kinetics with various levels of ignition energy. By varying the amount of ignition energy, three typical initiation behaviors, which were subcritical, supercritical and critical regimes, were observed. Then, the ignition energy of more than $137.5{\times}10^6$ in non-dimensional value was required for initiating a spherical detonation wave, and the minimum ignition energy(i.e., critical energy) was less than that of the one-dimensional simulation reported by a previous numerical work. When the ignition energy was less than the critical energy, the blast wave generated from an ignition source continued to attenuate due to the separation of the blast wave and a reaction front. Therefore, detonation was not initiated in the subcrtical regime. When the ignition energy was more than the minimum initiation energy, the blast wave developed into a multiheaded detonation wave propagating spherically at CJ velocity, and then a cellular pattern radiated regularly out from the ignition center in the supercritical regime. The influence on ignition energy was observed in the cell width near the ignition center, but the cell width on the fully developed detonation remained constant during the expanding of detonation wave due to the consecutive formation of new triple points, regardless of ignition energy. When the ignition energy was equal to the critical energy, the decoupling of the blast wave and a reaction front appeared, as occurred in the subcrtical regime. After that, the detonation bubble induced by the local explosion behind the blast wave expanded and developed into the multiheaded detonation wave in the critical regime. Although few triple points were observed in the vicinity of the ignition core, the regularly located cellular pattern was generated after the onset of the multiheaded detonation. Then, the average cell width on the fully developed detonation was almost to that in the supercritical regime. These numerical results qualitatively agreed with previous experimental works regarding the initiation and propagation processes.

• 한국연소학회:학술대회논문집
• /
• 2015.12a
• /
• pp.169-172
• /
• 2015

This paper presents a numerical investigation on kerosene-air mixture detonation and behaviors of thermal elasto-plstic thin metal tube under detonation loading based on multi-material analysis. The detonation loading is modeled by the kerosene-air mixture detonation which is compared with CJ condition and experimental cell size. And the thermal softening effect on elasto-plstic model of metal tube is indicated by different dynamic response of detonation loaded tube in various temperature and tube thickness.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• v.y2005m4
• /
• pp.343-346
• /
• 2005

Three dimensional structures of detonation wave propagating through a square-shaped duct were investigated using computational method and parallel processing. Inviscid fluid dynamics equations coupled with $variable-{\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by MUSCL-type TVD scheme and four stage Runge-Kutta time integration. The unsteady computational results in three dimension show the detailed mechanism of transverse mode and diagonal mode of detonation wave instabilities resulting same cell length but different cell width in smoked-foil record.