Engineering topics

Reaugh, John E.

Publications and source records attributed to Reaugh, John E..

6 records

DOE OSTI2021

Observation of asymmetric explosive density evolution in the deflagration-to-detonation transition for porous explosives

We report deflagration-to-detonation transition phenomena in cylindrical columns of porous explosives based on octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX). The HMX powder and LX-14 (95.5 wt. % HMX and 4.5 wt. % estane polyurethane) prills were observed with high-speed cameras and flash x rays under conditions that allow for precise scaling and registration of the recorded images. These experimental results test the conventional assumption of axisymmetric, effective one-dimensional burning. Polycarbonate confinement tubes allowed this direct imaging with steel spheres employed as radio-opaque fiducials. Spheres embedded in the explosive column revealed the displacement of the explosive prior to deflagration. X-ray and fast camera images show that a dense plug of compact explosive develops ahead of the deflagration front. The internal fiducials register both internal changes in the plug and changes in the position between successive x-ray images. Tantalum witness strips placed on the inner wall of the explosive channel revealed radial expansion of the tube bore at successive x-ray exposures, indicating the internal pressure of about 100 MPa just before detonation. Images indicate that the deflagration is not only asymmetric but also involves newly observed transient phenomena just before detonation. Furthermore, we typically see a brief but continuous increase in deflagration speed and the formation of a deflagration channel that bypasses the plug just prior to detonation. We hypothesize that these phenomena play an important role in the deflagration-to-detonation transition itself.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Modeling The Effects of Shock Pressure and Pore Morphology on Hot Spot Mechanisms in HMX

Here, we investigate the effects of shock pressure and pore morphology on the formation and growth of hot spots in HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Both non-reactive and reactive ALE3D simulations are used in these studies. Our non-reactive simulations show a viscous-dominated pore collapse mode at lower shock pressures (2–10 GPa) with shear band formation and a hydrodynamic-dominated mode at higher shock pressures (20-40 GPa) due to bulk melting. When normalized by bulk shock heating, viscous-dominated pore collapse modes are more efficient at generating hot spots. Pore morphology influences the post-collapse temperature distributions and reaction rate for a fixed pore area and shock pressure. We find that multiple surface pores at the binder-grain interface tend to react the fastest. Due to their upstream location in the HMX grain, the surface pores collapse sooner than interior pores; thus, the extent of reaction will generally favor these morphologies because they have more time to grow. In general, multiple smaller hot spots tend to react faster than a single larger hot spot because they accelerate one another's burning. The rank order of morphology effects, however, is not the same for non-reactive and reactive simulations. For example, while multiple surface pores produce the highest reaction rates they do not produce the highest (non-reactive) hot spot temperatures. In conclusion, our numerical studies provide insights on hot spot mechanisms in lieu of direct measurements and can be used to develop advanced shock initiation models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

A Computer Model to Study the Response of Energetic Materials to a Range of Dynamic Loads

Over the past decade we developed a model to enable computer simulation of the mechanical and subsequent energetic response of explosives and propellants to mechanical insults such as impacts, perforations, drops, and falls. The model is embedded in computer simulation programs that solve the non-linear, large deformation equations of compressible solid and fluid flow in space and time. It is implemented as a user-defined model, which returns the updated stress tensor and composition that result from the simulation supplied strain tensor change. Although it is multi-phase, in that gas and solid species are present, it is single-velocity, in that the gas does not flow through the porous solid. More than 70 time-dependent variables are made available for additional analyses and plotting. The model encompasses a broad range of possible responses: mechanical damage with no energetic response, and a continuous spectrum of degrees of violence including delayed and prompt detonation. This work describes the basic workings of the model.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Theory, Solution Methods, and Implementation of the HERMES Model

The HERMES (high explosive response to mechanical stimulus) model was developed over the past decade to enable computer simulation of the mechanical and subsequent energetic response of explosives and propellants to mechanical insults such as impacts, perforations, drops, and falls. The model is embedded in computer simulation programs that solve the non-linear, large deformation equations of compressible solid and fluid flow in space and time. It is implemented as a user-defined model, which returns the updated stress tensor and composition that result from the simulation supplied strain tensor change. Although it is multi-phase, in that gas and solid species are present, it is single-velocity, in that the gas does not flow through the porous solid. More than 70 time-dependent variables are made available for additional analyses and plotting. The model encompasses a broad range of possible responses: mechanical damage with no energetic response, and a continuous spectrum of degrees of violence including delayed and prompt detonation. This paper describes the basic workings of the model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Direct numerical simulation of shear localization and decomposition reactions in shock-loaded HMX crystal

A numerical model is developed to study the shock wave ignition of HMX crystal. The model accounts for the coupling between crystal thermal/mechanical responses and chemical reactions that are driven by the temperature field. This allows for the direct numerical simulation of decomposition reactions in the hot spots formed by shock/impact loading. The model is used to simulate intragranular pore collapse under shock wave loading. In a reference case: (i) shear-enabled micro-jetting is responsible for a modest extent of reaction in the pore collapse region, and (ii) shear banding is found to be an important mode of localization. The shear bands, which are filled with molten HMX, grow out of the pore collapse region and serve as potential ignition sites. The model predictions of shear banding and reactivity are found to be quite sensitive to the respective flow strengths of the solid and liquid phases. In this regard, it is shown that reasonable assumptions of liquid-HMX viscosity can lead to chemical reactions within the shear bands on a nanosecond time scale.

36 MATERIALS SCIENCE↗

DOE OSTI2013

HERMES Model Modifications and Applications 2012

HERMES was developed to meet requirements for a model to describe explosive response to impacts that do not produce prompt detonation, but can include delayed detonations such as DDT and XDT. Such response is termed High Explosive Violent Response (HEVR). The HERMES model includes a constitutive description of strength as a function of pressure, strain-rate, strain, and ambient temperature; the development of damage including fragmentation and porosity; the conditions leading to ignition; the propagation of a burning reaction through the volume of damaged explosive; the burn-up of the damaged and undamaged explosive; and the development of detonation using the Computational Reaction Evolution dependent on Entropy (S) and Time (CREST) detonation model. CREST can distinguish gradual pressure increases from shock waves, and also predict the increase in shock sensitivity to detonation that results from increased porosity. This year HERMES has been implemented as a User Defined Function in ALE 3D. As a result, modifications to the model can be implemented in a timely fashion asynchronously from the ALE 3D distribution process. There are differences between the implementation of HERMES in ALE 3D and the implementation in LS DYNA. The major difference is the propagation of an ignition front from an element that has started to burn. In LS DYNA, the time at which an element begins to burn is based on the minimum line-of-sight distance between the element and any previously ignited element. In ALE 3D the time to burn is based on a level-set method that propagates the ignition front from all elements that were previously ignited. The ALE 3D implementation modifies the time of ignition when element deformation changes the distance from the element to the source. The LS DYNA implementation does not. There are also differences in the implementation of artificial viscosity, which is important for CREST since the artificial viscosity increases the element entropy in a shock. In general, however, the differences have proven minor. In collaboration with the Atomic Weapons Establishment, UK (AWE), we have characterized two additional explosives – PETN and Composition B, a mixture of 60% RDX and 40% TNT by weight. Applications of HERMES to these explosives are described below.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗