1 / 13

Detonation Failure during Diffraction

Detonation Failure during Diffraction. Marco Arienti, Joseph E. Shepherd CIT ASCI-ASAP FY01 Research Review October 22-23, 2001. Decoupled detonation in O 2 -H 2 mixture (E. Shultz). cross-sections. j = . 0. density. I. II. mass-fraction.

kaz
Download Presentation

Detonation Failure during Diffraction

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Detonation Failure during Diffraction Marco Arienti, Joseph E. Shepherd CIT ASCI-ASAP FY01 Research Review October 22-23, 2001

  2. Decoupled detonation in O2-H2 mixture(E. Shultz)

  3. cross-sections j = 0 density I II mass-fraction Decoupled Detonation in HMX(Eric Morano, FY00) pressure mass-fraction initial location ZND profile reflective VN II CJ Dirichlet reflective reflective I reflective Pressure-dependent reaction rate:

  4. Idea: examine the problem of detonation failure in terms of decoupling of the reaction zone from the shock front (in the mean flow). • Note: most reaction rate laws are strongly temperature dependent. When decoupling occurs, shocked particles are prevented from undergoing thermal runaway. • In a corner-turning geometry, the expansion propagating from the corner induces a curvature in the undisturbed detonation front and may trigger decoupling.

  5. Pressure field (z-axis) of failing detonation ideal gas and one-step Arrhenius chemistry • Examine the propagation of corner- generated expansionwaves. temperature isocontours corner

  6. Adapt Skews’ construction for diffracting shock(Shultz & Shepherd, 2001) for hydrocarbon and hydrogen mixtures

  7. Diagnostics & Modeling • Idea: follow a set of Lagrangian particles and perform “numerical dominant balance" of the terms contributing to the material derivative of temperature. • Cast Euler equation in reference frame attached to the shock (use Bertrand coordinates). • Track the detonation front and construct an intrinsic reference frame. • Express Lagrangian (material) derivative DT/Dt in this reference.

  8. Generalization to 2D settings • Bertrand (or intrinsic coordinates) x and h.The local velocity of the front is: • In cylindrical coordinates:

  9. Euler equations for reactive flow(Bdzil & Azlam, 2000) (1) (2) (3)

  10. Material derivative of Temperature (4)

  11. Special case: cylindrical symmetry(Eckett & Shepherd, JFM 2000) • When evaluated in cylindrical (or spherical) symmetry, Eq. (4) loses N1, N2, N3, the divergence, the centripetal and the transverse derivative term. (4)

  12. Analysis of detonation extinction in spherical symmetry Particle trajectory particle 10 fails to ignite time time time time

  13. Axial shock decay model • When evaluated along the axis of symmetry, Eq. (4) loses N1, N2, N3, the centripetal and the transverse derivative term. • In previous investigations, Taylor-Sedov strong spherical blast similarity solution (modified for post-shock energy release) has been used to model the evolution of the detonation front after the corner disturbance has reached the centerline.

More Related