• 한국연소학회:학술대회논문집
• /
• 2012.11a
• /
• pp.359-362
• /
• 2012

The results are shown highly unstable mode's detonation dynamics by compared with weakly unstable mode. And we investigate the difference and similar features of 2D and 3D results. By using PSD via FFT, the effects of pre-exponential factor difference and of unstable mode were investigated in this study. The result of PSD is shown pulsed features in weakly unstable mode, but noselike in highly unstable mode. By compared between Sheliren image and overlaid slice image, the irregular feature of detonation waves structure was discussed in highly unstable mode.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.15 no.4
• /
• pp.57-64
• /
• 2011

It has been suggested that turbulent effect should be considered for the study of highly unstable detonation of hydrocarbon fuels, as in the case of pulse detonation engine (PDE). A series of numerical study are carried out to understand the characteristics of the highly unstable detonation by considering viscosity, turbulence model and turbulence-combustion interaction model. Through studies of the different levels of modeling, it is understood that the viscosity and turbulence have negligible effects on low frequency characteristics, but tend to enhance the high frequency characteristics. It is also considered that the turbulence-chemistry interaction model should be taken the influence of the activation energy into account for detonation studies.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2007.11a
• /
• pp.367-370
• /
• 2007

Present study examines detonation wave propagation characteristics in annular channel. A normalized value of channel width to the annular radius was considered as a geometric parameter. A parametric study was carried out for a various regimes of detonation waves from weakly unstable to highly unstable detonation waves. Numerical approaches that used in the previous study of numerical requirements of the simulation of detonation wave propagations in 2D and 3D channel were used also for the present study with OpenMP parallization for multi-core SMP machines. The major effect of the curved geometry on the detonation wave propagation seems to be a flow compression effect, regardless of the detonation regimes. The flow compression behind the detonation wave by the curved geometry of the circular channel pushes the detonation wave front and results in the overdriven detonation waves with increased detonation speed beyond the Chapmann-Jouguet speed. This effect gets stronger as the normalized radius smaller, as expected. The effect seems to be negligible beyond the normalized radius of 10.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2006.05a
• /
• pp.314-318
• /
• 2006

The interest on the detonation in small tubes, which can be applied to the ignition devices of propulsion system, is increasing. However, the propagation dynamics of detonation waves in small tubes has not been investigated clearly yet. In the present experiments, propagations of detonation waves in stoichiometric propane-oxygen mixture through transparent tubes were recorded using a high speed camera and average velocities were measured as well. In terms of average velocity, there exists a transition regime where the waves show smooth transition from the normal Chapman-Jouguet(CJ) detonation to the low velocity detonation$(\sim0.5V_{CJ})$ along the decreasing initial pressure. In this transition regime, the detonation waves are highly unstable and show cyclic or intermittent longitudinal velocity fluctuation.

• 한국연소학회:학술대회논문집
• /
• 2014.11a
• /
• pp.169-174
• /
• 2014

Detonation propagation studies is recently getting more attention in these days for its feasibility in aerospace application. Another motivation for this study is the safety concern in industries, since the detonation can cause failure to the mechanical components particularly when the flame accelerates within a pipe or tubes. In this study we numerically simulated a Moderately unstable detonation case with various grid systems and fluid dynamic length scales and have compared in the contents. Moderately Unstable detonation case was selected for this study and detailed Hydrogen-Air Reaction Mechanisms proposed by Jachimowski was used in this study with N2 as inert species.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

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

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

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