Effect of initial disturbance on the detonation front structure of a narrow duct

Dou, Hua-Shu; Khoo, Boo Cheong

doi:10.1007/s00193-009-0240-8

Cited by 17 publications

(11 citation statements)

References 43 publications

Supporting

Mentioning

Contrasting

Unclassified

Order By: Relevance

“…Figure 4 shows, in the duct with width 2×L 1/2 , the typical front feature and density gradient on the walls of the detonation in a cycle. The front structure is similar to that in [48,49], suggesting spinning detonations have formed. Due to the differences in the type of perturbation, spinning detonation rotates counter-clockwise, opposite to that in [48].…”

Section: Verification Of the Numerical Resolution Of Three-dimensionasupporting

confidence: 69%

“…The front structure is similar to that in [48,49], suggesting spinning detonations have formed. Due to the differences in the type of perturbation, spinning detonation rotates counter-clockwise, opposite to that in [48]. It can be seen that, the triple lines and transverse waves collide with the walls, and strong explosions take place near the walls.…”

Section: Verification Of the Numerical Resolution Of Three-dimensionasupporting

confidence: 69%

See 1 more Smart Citation

High Resolution WENO Simulation of 3D Detonation Waves

Wang¹,

Shu²,

Han³

et al. 2012

View full text Add to dashboard Cite

In this paper, we develop a three-dimensional parallel solver using the fifth order highresolution weighted essentially non-oscillatory (WENO) finite difference scheme to perform extensive simulation for three-dimensional gaseous detonations. A careful study is conducted for the propagation modes of the three-dimensional gaseous detonation wave-front structure in a long square duct with different widths under different initial perturbations. The numerical results indicate that, with a transverse sinusoidal perturbation of the initial ZND profile, when the width of the duct is less than the cellular width (4.5×L 1/2 ), unstable detonations can trigger a spinning motion in the duct. The detonation wave propagates in a single-headed spinning motion, with a distinctive "ribbon" displayed on the four walls. In this case, the measured pitch-to-diameter ratio is approximately 3.42, which is slightly larger than the theoretically predicted value 3.128 for a round duct. When the channel width is greater than the cellular width, detonation waves propagate in an out-of-phase rectangular mode. With a transverse cosine perturbation of the initial ZND profile, the front of the stable detonation has a rectangular structure, and regular cellular patterns and in-phase "slapping waves" can be observed clearly on the four walls. The width-to-length ratio of the cellular patterns is approximately 0.5. For a mildly unstable detonation, its front has an in-phase rectangular structure at the early stage, then the wave-front becomes flat. Over time, but it still maintains an in-phase rectangular structure after reigniting. For highly unstable detonations, the wave-front has a rectangular structure at the early stage. After a low pressure stage for a very long time, detonation occurs once again. At this time, the detonation front structure becomes very twisted, and the triple-lines become asymmetrical. Finally, a spinning detonation mode is triggered. With a symmetrical perturbation mode along the diagonals of the detonation front, for the stable detonation, an diagonal detonation is formed and the detonation front maintains a diagonal structure, but no "slapping waves" appears on the walls. The width-to-length ratio of the cellular structure is equal to that in the rectangular structure. For mildly unstable and highly unstable detonations, the front has a diagonal structure at the early stage. After a short period of time, the diagonal structure of the detonation front cannot be maintained, and it ultimately evolves into a spinning detonation.Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Oper...

show abstract

Section: Verification Of the Numerical Resolution Of Three-dimensionasupporting

confidence: 69%

Section: Verification Of the Numerical Resolution Of Three-dimensionasupporting

confidence: 69%

High Resolution WENO Simulation of 3D Detonation Waves

Wang¹,

Shu²,

Han³

et al. 2012

View full text Add to dashboard Cite

show abstract

“…The governing equations describing the fluid flow and the detonation propagation are the threedimensional Euler equations with a source term which represents the chemical reaction process. In conservative form, these are written as [29][30][31]…”

Section: Governing Equationsmentioning

confidence: 99%

Computational Study of Deflagration to Detonation Transition in a Straight Duct: Effect of Energy Release

Dou

Khoo

2012

Int. J. Mod. Phys. Conf. Ser.

Self Cite

View full text Add to dashboard Cite

Numerical 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 nonoscillatory 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 at q =50, 30, 25, 20, 15, 10 and 5, the specific heat ratio at 1.2, and the activation temperature at Ti =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 for q =50~20, and increases with the increase of the energy release for q =50~20. 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.

show abstract

“…High-order finite-difference methods based on essentially non-oscillatory (ENO) [1,2] and weighted ENO (WENO) [3][4][5][6] schemes for spatial discretisation coupled with explicit [7,8] as well as implicit [9,10] time integration techniques are very popular in computational fluid dynamics (CFD) of compressible flows. These schemes are very robust, even in the presence of strong discontinuities like shocks [11] or detonation waves [7]. This property can be traced back to the adaptivity of the stencil to the present flow conditions.…”

Section: Introductionmentioning

confidence: 99%

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

Hermes

Klioutchnikov

Olivier

2011

Numerical Methods in Fluids

View full text Add to dashboard Cite

SUMMARY This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.

show abstract

Effect of initial disturbance on the detonation front structure of a narrow duct

Cited by 17 publications

References 43 publications

High Resolution WENO Simulation of 3D Detonation Waves

High Resolution WENO Simulation of 3D Detonation Waves

Computational Study of Deflagration to Detonation Transition in a Straight Duct: Effect of Energy Release

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

Contact Info

Product

Resources

About