• 한국연소학회:학술대회논문집
• /
• 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
• /
• 2008.03a
• /
• pp.364-370
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable and unstable pitch modes for the lower and higher activation energies, respectively. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of two modes. The maximum pressure history in the stable pitch remained nearly constant, and the single Mach leg existing on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the unstable pitch due to the generation and decay of complex Mach interaction on the shock front shape. The high frequency oscillation was self-induced because the intensity of the transverse wave was changed during propagation in one cycle. The high frequency behavior was not always the same for each cycle, and therefore the low frequency oscillation was also induced in the pressure history.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.367-373
• /
• 2008

Spinning detonations propagating in a circular tube were numerically investigated with a one-step irreversible reaction model governed by Arrhenius kinetics. Activation energy is used as parameter as 10, 20, 27 and 35, and the specific heat ratio and the heat release are fixed as 1.2 and 50. The time evolution of the simulation results was utilized to reveal the propagation mechanism of single-headed spinning detonation. The track angle of soot record on the tube wall was numerically reproduced with various levels of activation energy, and the simulated unique angle was the same as that of the previous reports. The maximum pressure histories of the shock front on the tube wall showed stable pitch at Ea=10, periodical unstable pitch at Ea=20 and 27 and unstable pitch consisting of stable, periodical unstable and weak modes at Ea=35, respectively. In the weak mode, there is no Mach leg on the shock front, where the pressure level is much lower than the other modes. The shock front shapes and the pressure profiles on the tube wall clarified the mechanisms of these stable and unstable modes. In the stable pitch at Ea=10, the maximum pressure history on the tube wall remained nearly constant, and the steady single Mach leg on the shock front rotated at a constant speed. The high and low frequency pressure oscillations appeared in the periodical unstable pitch at Ea=20 and 27 of the maximum pressure history. The high frequency was one cycle of a self-induced oscillation by generation and decay in complex Mach interaction due to the variation in intensity of the transverse wave behind the shock front. Eventually, sequential high frequency oscillations formed the low frequency behavior because the frequency behavior was not always the same for each cycle. In unstable pitch at Ea=35, there are stable, periodical unstable and weak modes in one cycle of the low frequency oscillation in the maximum pressure history, and the pressure amplitude of low frequency was much larger than the others. The pressure peak appeared after weak mode, and the stable, periodical unstable and weak modes were sequentially observed with pressure decay. A series of simulations of spinning detonations clarified that the unsteady mechanism behind the shock front depending on the activation energy.