scholarly journals Modeling two‐dimensional detonations with detonation shock dynamics

1989 ◽  
Vol 1 (7) ◽  
pp. 1261-1267 ◽  
Author(s):  
J. B. Bdzil ◽  
D. S. Stewart
Keyword(s):  
Two Dimensional ◽  
Shock Dynamics ◽  
Detonation Shock

Geometrical shock dynamics applied to condensed phase materials

2017 ◽  
Vol 828 ◽  
pp. 104-134 ◽  
Author(s):  
Brandon Lieberthal ◽  
D. Scott Stewart ◽  
Alberto Hernández
Keyword(s):  
Shock Front ◽  
Blast Wave ◽  
High Explosive ◽  
Spherical Particles ◽  
Shock Propagation ◽  
Inert Material ◽  
Two Dimensional ◽  
Shock Dynamics ◽  
Geometrical Shock Dynamics ◽  
Detonation Shock

Taylor blast wave (TBW) theory and geometrical shock dynamics (GSD) theory describe a radially expanding shock wave front through an inert material, typically an ideal gas, in the strong blast wave limit and weak acoustic limit respectively. We simulate a radially expanding blast shock in air using a hydrodynamic simulation code and numerically describe the intermediate region between these two limits. We test our description of the intermediate shock phase through a two-dimensional simulation of the Bryson and Gross experiment. We then apply the principles of GSD to materials that follow the Mie–Gruneisen equation of state, such as plastics and metals, and derive an equation that accurately relates the acceleration, velocity and curvature of the shock through these materials. Along with detonation shock dynamics (DSD), which describes detonation shock propagation through high explosive fluids, we develop a hybrid DSD/GSD model for the simulation of heterogeneous explosives. This model enables computationally efficient simulation of the shock front in high explosive/inert mixtures consisting of simple or complex geometric configurations. We simulate an infinite two-dimensional slab consisting of one half explosive, PBXN-9, and one half aluminium and model the boundary angle conditions using shock polar analysis. We also simulate a series of high explosive unit cells embedded with aluminium spherical particles, and we compare the propagation of the detonation shock front with a direct numerical simulation performed with the ALE3D code.

scholarly journals Shock dynamics in periodic two-dimensional composites

2020 ◽  
Author(s):  
Lauren L. Poole ◽  
Manny Gonzales ◽  
Abdel R. Moustafa ◽  
Austin R. C. Gerlt ◽  
Zachary C. Cordero
Keyword(s):  
Two Dimensional ◽  
Shock Dynamics

Steady detonation propagation in a circular arc: a Detonation Shock Dynamics model

2016 ◽  
Vol 807 ◽  
pp. 87-134 ◽  
Author(s):  
Mark Short ◽  
James J. Quirk ◽  
Chad D. Meyer ◽  
Carlos Chiquete
Keyword(s):  
Surface Wave ◽  
Angular Speed ◽  
Order Correction ◽  
Two Dimensional ◽  
Normal Surface ◽  
Circular Arc ◽  
First Order ◽  
Shock Dynamics ◽  
First Order Correction ◽  
Steady Detonation

We study the physics of steady detonation wave propagation in a two-dimensional circular arc via a Detonation Shock Dynamics (DSD) surface evolution model. The dependence of the surface angular speed and surface spatial structure on the inner arc radius ($R_{i}$), the arc thickness ($R_{e}-R_{i}$, where $R_{e}$ is the outer arc radius) and the degree of confinement on the inner and outer arc is examined. We first analyse the results for a linear $D_{n}$–$\unicode[STIX]{x1D705}$ model, in which the normal surface velocity $D_{n}=D_{CJ}(1-B\unicode[STIX]{x1D705})$, where $D_{CJ}$ is the planar Chapman–Jouguet velocity, $\unicode[STIX]{x1D705}$ is the total surface curvature and $B$ is a length scale representative of a reaction zone thickness. An asymptotic analysis assuming the ratio $B/R_{i}\ll 1$ is conducted for this model and reveals a complex surface structure as a function of the radial variation from the inner to the outer arc. For sufficiently thin arcs, where $(R_{e}-R_{i})/R_{i}=O(B/R_{i})$, the angular speed of the surface depends on the inner arc radius, the arc thickness and the inner and outer arc confinement. For thicker arcs, where $(R_{e}-R_{i})/R_{i}=O(1)$, the angular speed does not depend on the outer arc radius or the outer arc confinement to the order calculated. It is found that the leading-order angular speed depends only on $D_{CJ}$ and $R_{i}$, and corresponds to a Huygens limit (zero curvature) propagation model where $D_{n}=D_{CJ}$, assuming a constant angular speed and perfect confinement on the inner arc surface. Having the normal surface speed depend on curvature requires the insertion of a boundary layer structure near the inner arc surface. This is driven by an increase in the magnitude of the surface wave curvature as the inner arc surface is approached that is needed to meet the confinement condition on the inner arc surface. For weak inner arc confinement, the surface wave spatial variation with the radial coordinate is described by a triple-deck structure. The first-order correction to the angular speed brings in a dependence on the surface curvature through the parameter $B$, while the influence of the inner arc confinement on the angular velocity only appears in the second-order correction. For stronger inner arc confinement, the surface wave structure is described by a two-layer solution, where the effect of the confinement on the angular speed is promoted to the first-order correction. We also compare the steady-state arc solution for a PBX 9502 DSD model to an experimental two-dimensional arc geometry validation test.

Coupling a Detonation Shock Dynamics Code to an Eulerian Hydrocode

2009 ◽  
Author(s):  
David L. Littlefield ◽  
Young-Ho Kim ◽  
David E. Lambert
Keyword(s):  
Input File ◽  
User Input ◽  
Preliminary Work ◽  
Shock Dynamics ◽  
Boundary Information ◽  
Detonation Shock ◽  
Technical Effort

In this technical effort we have incorporated the Detonation Shock Dynamics (DSD) methodology into the Eulerian hydrocode CTH by coupling it to a DSD code called WAVETRACKER. Based on some preliminary work, we developed a strategy for coupling these two codes. Both CTH and WAVETRACKER are stand-alone codes. In the coupled code, CTH calls WAVETRACKER as a subroutine. Interface routines are developed to pass geometry, mesh and boundary information from CTH to WAVETRACKER. Additional parameters required by WAVETRACKER are read directly from a user input file. Results in the form of a burn table are interpolated and passed back from WAVETRACKER to CTH.

A numerical scheme for shock propagation in three dimensions

1988 ◽  
Vol 416 (1850) ◽  
pp. 179-198 ◽  
Keyword(s):  
Numerical Scheme ◽  
Three Dimensional ◽  
Experimental Results ◽  
Three Dimensions ◽  
Shock Propagation ◽  
Two Dimensional ◽  
Shock Surface ◽  
Shock Dynamics ◽  
Geometrical Shock Dynamics ◽  
Three Space

A numerical scheme for shock propagation in three space dimensions is presented. The motion of the leading shock surface is calculated by using Whitham’s theory of geometrical shock dynamics. The numerical scheme is used to examine the focusing of initially curved shock surfaces and the diffraction of shocks in a pipe with a 90° bend. Numerical and experimental results for the corresponding two-dimensional or axi-symmetrical cases are used to compare with the new and more complicated three-dimensional results.

Sign in / Sign up

Export Citation Format

Share Document