Detonation shock dynamics modeling and calibration of the HMX-based conventional high explosive PBX 9501 with application to the two-dimensional circular arc geometry

2020 ◽  
Vol 222 ◽  
pp. 213-232 ◽  
Author(s):  
Carlos Chiquete ◽  
Mark Short ◽  
Stephen J. Voelkel ◽  
Eric K. Anderson ◽  
Scott I. Jackson
2017 ◽  
Vol 828 ◽  
pp. 104-134 ◽  
Author(s):  
Brandon Lieberthal ◽  
D. Scott Stewart ◽  
Alberto Hernández

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.


2016 ◽  
Vol 807 ◽  
pp. 87-134 ◽  
Author(s):  
Mark Short ◽  
James J. Quirk ◽  
Chad D. Meyer ◽  
Carlos Chiquete

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.


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

Author(s):  
Asif Arefeen ◽  
Yujiang Xiang

Abstract A novel multibody dynamics modeling method is proposed for two-dimensional (2D) team lifting prediction. The box itself is modeled as a floating-base rigid body in Denavit-Hartenberg representation. The interactions between humans and box are modeled as a set of grasping forces which are treated as unknowns (design variables) in the optimization formulation. An inverse-dynamics-based optimization method is used to simulate the team lifting motion where the dynamic effort of two humans is minimized subjected to physical and task-based constraints. The design variables are control points of cubic B-splines of joint angle profiles of two humans and the box, and the grasping forces between humans and the box. Two numerical examples are successfully simulated with different box weights (20 Kg and 30 Kg, respectively). The humans’ joint angle, torque, ground reaction force, and grasping force profiles are reported. The joint angle profiles are validated with the experimental data.


Author(s):  
David L. Littlefield ◽  
Young-Ho Kim ◽  
David E. Lambert

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.


Sign in / Sign up

Export Citation Format

Share Document