• 한국연소학회:학술대회논문집
• /
• 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.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.11 no.1
• /
• pp.1-10
• /
• 2007

Three dimensional structures of detonation waves propagating in a square tube were investigated using a high resolution CFD code coupled with a conservation equation of reaction progress variable and an one-step irreversible reaction. The code were parallelized based on domain decomposition technique using MPI library. The computations were carried on an in-house Windows cluster with AMD processors. Three-dimensional unsteady analysis results in the smoked-foil records caused by the instabilities of the detonation waves, which showed the rectangular and diagonal modes of detonation instabilities depending on the initial condition of disturbances and the spinning detonation for case of small reaction constant.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.151-155
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied 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 a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• 한국연소학회:학술대회논문집
• /
• 2005.10a
• /
• pp.15-19
• /
• 2005

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied 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 a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

• 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.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2012.05a
• /
• pp.278-281
• /
• 2012

Numerical studies were performed to investigate the three-dimensional front structure and dynamics features of detonation wave propagating in a circular tube such as Pulse Detonation Engine (PDE). By carrying out a series of parametric study using one step irreversible Arrhenius kinetics model, mechanisms of the three-dimensional front structure were investigated for two-, three-, four and six-cell mode detonations. A comparison with two-dimensional results, the effects of slapping transverse waves in radial direction were confirmed. In the all muti-cell modes, the detonation front structures and smoked-records on the wall are formed by the propagation of transverse waves along the wall in clockwise and counter-clockwise while the slapping move in radial direction. And the strength of reflected waves on the curved wall is changed by the multi-dimensional confinement effect.

• 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.

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

Three-dimensional structures of detonation wave propagating in circular tube were investigated. Inviscid fluid dynamics equations coupled with a conservation equation of reaction progress variable were analyzed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Variable-$\gamma$ formulation was used to account for the variable properties between unburned and burned states and the chemical reaction was modeled by using a simplified one-step irreversible kinetics model. The computational code was parallelized based on domain decomposition technique using MPI-II message passing library. The computations were carried out using a home made Windows based PC cluster having 160 AMD AthloxXP and Athlon64 processor. The computational domain consisted of through a roundshaped tube with wall conditions. As an initial condition, analytical ZND solution was distributed over the computational domain with disturbances. The disturbances has circumferential large gradient. The unsteady computational results in three-dimension show the detailed mechanisms of multi-cell mode of detonation wave instabilities resulting diamond shape in smoked-foil record.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.341-345
• /
• 2006