skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Hydrodynamic Modeling of Air Blast Propagation from the Humble Redwood Chemical High Explosive Detonations Using GEODYN

Abstract

Two-dimensional axisymmetric hydrodynamic models were developed using GEODYN to simulate the propagation of air blasts resulting from a series of high explosive detonations conducted at Kirtland Air Force Base in August and September of 2007. Dubbed Humble Redwood I (HR-1), these near-surface chemical high explosive detonations consisted of seven shots of varying height or depth of burst. Each shot was simulated numerically using GEODYN. An adaptive mesh refinement scheme based on air pressure gradients was employed such that the mesh refinement tracked the advancing shock front where sharp discontinuities existed in the state variables, but allowed the mesh to sufficiently relax behind the shock front for runtime efficiency. Comparisons of overpressure, sound speed, and positive phase impulse from the GEODYN simulations were made to the recorded data taken from each HR-1 shot. Where the detonations occurred above ground or were shallowly buried (no deeper than 1 m), the GEODYN model was able to simulate the sound speeds, peak overpressures, and positive phase impulses to within approximately 1%, 23%, and 6%, respectively, of the actual recorded data, supporting the use of numerical simulation of the air blast as a forensic tool in determining the yield of an otherwise unknown explosion.

Authors:
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1035964
Report Number(s):
LLNL-TR-501037
TRN: US1201308
DOE Contract Number:
W-7405-ENG-48
Resource Type:
Technical Report
Country of Publication:
United States
Language:
English
Subject:
42 ENGINEERING; 58 GEOSCIENCES; 73 NUCLEAR PHYSICS AND RADIATION PHYSICS; AIR; CHEMICAL EXPLOSIVES; EFFICIENCY; EXPLOSIONS; HYDRODYNAMIC MODEL; HYDRODYNAMICS; PRESSURE GRADIENTS; SIMULATION; VELOCITY

Citation Formats

Chipman, V D. Hydrodynamic Modeling of Air Blast Propagation from the Humble Redwood Chemical High Explosive Detonations Using GEODYN. United States: N. p., 2011. Web. doi:10.2172/1035964.
Copy to clipboard
Chipman, V D. Hydrodynamic Modeling of Air Blast Propagation from the Humble Redwood Chemical High Explosive Detonations Using GEODYN. United States. doi:10.2172/1035964.
Copy to clipboard
Chipman, V D. Tue . "Hydrodynamic Modeling of Air Blast Propagation from the Humble Redwood Chemical High Explosive Detonations Using GEODYN". United States. doi:10.2172/1035964. https://www.osti.gov/servlets/purl/1035964.
Copy to clipboard
@article{osti_1035964,
title = {Hydrodynamic Modeling of Air Blast Propagation from the Humble Redwood Chemical High Explosive Detonations Using GEODYN},
author = {Chipman, V D},
abstractNote = {Two-dimensional axisymmetric hydrodynamic models were developed using GEODYN to simulate the propagation of air blasts resulting from a series of high explosive detonations conducted at Kirtland Air Force Base in August and September of 2007. Dubbed Humble Redwood I (HR-1), these near-surface chemical high explosive detonations consisted of seven shots of varying height or depth of burst. Each shot was simulated numerically using GEODYN. An adaptive mesh refinement scheme based on air pressure gradients was employed such that the mesh refinement tracked the advancing shock front where sharp discontinuities existed in the state variables, but allowed the mesh to sufficiently relax behind the shock front for runtime efficiency. Comparisons of overpressure, sound speed, and positive phase impulse from the GEODYN simulations were made to the recorded data taken from each HR-1 shot. Where the detonations occurred above ground or were shallowly buried (no deeper than 1 m), the GEODYN model was able to simulate the sound speeds, peak overpressures, and positive phase impulses to within approximately 1%, 23%, and 6%, respectively, of the actual recorded data, supporting the use of numerical simulation of the air blast as a forensic tool in determining the yield of an otherwise unknown explosion.},
doi = {10.2172/1035964},
journal = {},
number = ,
volume = ,
place = {United States},
year = {Tue Sep 20 00:00:00 EDT 2011},
month = {Tue Sep 20 00:00:00 EDT 2011}
}
Copy to clipboard
Technical Report:
View Technical Report
DOI:  10.2172/1035964
Save / Share:
Export Metadata
Save to My Library
You must Sign In or Create an Account in order to save documents to your library.

  • ; ; ; ...

  • This report primarily concerns the use of two massively parallel finite element codes originally written and maintained at Lawrence Livermore National Laboratory. ALE3D is an explicit hydrodynamics code commonly employed to simulate wave propagation from high energy scenarios and the resulting interaction with nearby structures. This coupled response ensures that a structure is accurately applied with a blast loading varying both in space and time. Figure 1 illustrates the radial outward propagation of a pressure wave due to a center detonated spherical explosive originating from the lower left. The radial symmetry seen in this scenario is lost when instead amore » cylindrocal charge is detonated. Figure 2 indicates that a stronger, faster traveling pressure wave occurs in the direction of the normal axis to the cylinder. The ALE3D name is derived because of the use of arbitrary-Lagrange-Eulerian elements in which the mesh is allowed to advect; a process through which the mesh is modified to alleviate tanlging and general mesh distortion often cuased by high energy scenarios. The counterpart to an advecting element is a Lagrange element, whose mesh moves with the material. Ideally all structural components are kept Lagrange as long as possible to preserve accuracy of material variables and minimize advection related errors. Advection leads to mixed zoning, so using structural Lagrange elements also improves the visualization when post processing the results. A simplified representation of the advection process is shown in Figure 3. First the mesh is distorted due to material motion during the Lagrange step. The mesh is then shifted to an idealized and less distorted state to prevent irregular zones caused by the Lagrange motion. Lastly, the state variables are remapped to the elements of the newly constructed mesh. Note that Figure 3 represents a purely Eulerian mesh relaxation because the mesh is relocated back to the pre-Lagrange position. This is the case when the material flows through a still mesh. This is not typically done in an ALE3D analysis, especially if Lagrange elements exist. Deforming Lagrange elements would certainly tangle with a Eulerian mesh eventually. The best method in this case is to have an advecting mesh positioned as some relaxed version of the pre and post Lagrange step; this gives the best opportunity of modeling a high energy event with a combination of Lagrange and ALE elements. Dyne3D is another explicit dynamic analysis code, ParaDyn being the parallel version. ParaDyn is used for predicting the transient response of three dimensional structures using Lagrangian solid mechanics. Large deformation and mesh tangling is often resolved through the use of an element deletion scheme. This is useful to accommodate component failure, but if it is done purely as a means to preserve a useful mesh it can lead to problems because it does not maintain continuity of the material bulk response. Whatever medium exists between structural components is typically not modeled in ParaDyn. Instead, a structure either has a known loading profile applied or given initial conditions. The many included contact algorithms can calculate the loading response of materials if and when they collide. A recent implementation of an SPH module in which failed or deleted material nodes are converted to independent particles is currently being utilized for a variety of spall related problems and high velocity impact scenarios. Figure 4 shows an example of a projectile, given an initial velocity, and how it fails the first plate which generates SPH particles which then interact with and damage the second plate.« less