Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JWL equation of state”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

17 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

JWL Equation of State

The JWL equation of state (EOS) is frequently used for the products (and sometimes reactants) of a high explosive (HE). Here we review and systematically derive important properties. The JWL EOS is of the Mie-Grueneisen form with a constant Grueneisen coefficient and a constants specific heat. It is thermodynamically consistent to specify the temperature at a reference state. However, increasing the reference state temperature restricts the EOS domain in the (V, e)-plane of phase space. The restrictions are due to the conditions that P ≥ 0, T ≥ 0, and the isothermal bulk modulus is positive. Typically, this limits the low temperature regime in expansion. The domain restrictions can result in the P-T equilibrium EOS of a partly burned HE failing to have a solution in some cases. For application to HE, the heat of detonation is discussed. Example JWL parameters for an HE, both products and reactions, are used to illustrate the restrictions on the domain of the EOS.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2014

JWL in a Nutshell

The Jones-Wilkins-Lee (JWL) equation of state is an analytical equation of state (EOS) used in the modeling of detonation products of high explosives. The JWL parameters for typical explosives are found in a series of reports by Brigitta Dobratz. The JWL parameters are usually derived from thermochemical codes or experimental cylinder tests using an elaborate and arcane procedure. It is often said that the JWL parameters need to be “self-consistent”, “thermodynamically consistent”, or “consistent with the conservation laws of mass, momentum, and energy.” The purpose of this document is to explain and illuminate what that means. Each parameter in the JWL equation of state is connected to the other parameters. An ideal simulation would produce detonation products at the Chapman-Jouguet state. I begin by describing the form of this equation of state.

Engineering(42)↗
DOE OSTI2015

An Exact, Compressible One-Dimensional Riemann Solver for General, Convex Equations of State

This note describes an algorithm with which to compute numerical solutions to the one- dimensional, Cartesian Riemann problem for compressible flow with general, convex equations of state. While high-level descriptions of this approach are to be found in the literature, this note contains most of the necessary details required to write software for this problem. This explanation corresponds to the approach used in the source code that evaluates solutions for the 1D, Cartesian Riemann problem with a JWL equation of state in the ExactPack package [16, 29]. Numerical examples are given with the proposed computational approach for a polytropic equation of state and for the JWL equation of state.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Optimizing Code Calibration of the JWL Explosive Equation-of-State to the Cylinder Test

Recent analysis of high explosives using the cylinder test has resulted in an increased understanding of the experimental uncertainty for the cylinder test. To fit the Jones-Wilkins-Lee equation-of-state model to less than experimental uncertainty, an optimization process is developed in this paper. This process, when applied to a large range of explosives, seems to show that the exponent R 1 increases with increasing energy density and is grouped based on ideality. Finally, the optimal parameters are shown to increase the accuracy of an analytic equation-of-state tool.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Improvements to the New CTH Code Verification & Validation Test Suite (FY2020)

The CTH multiphysics hydrocode, which is used for a wide range of important calculations, has undertaken in recent years to overhaul its software quality and testing processes. A key part of this effort entailed building a new, robust V&V test suite made up of traditional hydrocode verification problems, such as those listed in the ASC Tri-Lab Test Suite and the Enhanced Tri-Lab Test Suite, as well as validation problems for some of CTHs most frequently used equations of state, materials models, and other key capabilities. Substantial progress towards this goal was made in FY19. In FY20, this test suite has been expanded to include verification and validation tests of the Sesame and JWL equation of state models as well as the Mader verification problem from the Tri-Lab Test Suite and the Blake verification problem - a linear elastic analog to the Hunter problem from the Enhanced Tri-Lab Test Suite. This report documents CTH performance on the new test suite problems. Verification test results are compared to analytic solutions and, for most tests, convergence results are presented. Validation test results are compared to experimental data and mesh refinement studies are included. CTH performs well overall on the new test problems. Convergence rates for the Blake and Mader problems are comparable to those for similar ASC codes. The JWL and Sesame verification tests show good agreement with analytic solutions. Likewise, CTH simulation results show good agreement with experimental validation data for the Sesame and JWL equations of state for the materials tested. Future V&V work will focus on adding tests for other key capabilities like fracture and high explosive models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Validation of Reactive Flow Capabilities in PAGOSA

Several reactive flow modeling capabilities are studied and applied to physics problems of interest in the PAGOSA hydrocode with an emphasis on validation. The Davis equation of state and the Arrhenius shock temperature dependent WSD (AWSD) reactive flow model have been implemented by the author in PAGOSA over the last year. These high-explosive modeling options along with the previously existing HE-JWL equation of state and the SURF/SURFplus reactive flow models are discussed, exercised and compared. Validation of the models analyze numerical results against experimental data from cylinder tests, two-stage gas gun experiments and proton radiography (pRad) experiments. Applications of the modeling options are explored involving both HMX based PBX 9501 and TATB based PBX 9502. The effect of mesh size is considered for prediction of phenomena such as corner turning and dead zones, and the relative code performance of the modeling options is discussed.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2015

NSR&D FY15 Final Report. Modeling Mechanical, Thermal, and Chemical Effects of Impact

The main goal of this project is to develop a computer model that explains and predicts coupled mechanical, thermal and chemical responses of HE under impact and friction insults. The modeling effort is based on the LANL-developed CartaBlanca code, which is implemented with the dual domain material point (DDMP) method to calculate complex and coupled thermal, chemical and mechanical effects among fluids, solids and the transitions between the states. In FY 15, we have implemented the TEPLA material model for metal and performed preliminary can penetration simulation and begun to link with experiment. Currently, we are working on implementing a shock to detonation transition (SDT) model (SURF) and JWL equation of state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calibration of reactive burn and Jones-Wilkins-Lee parameters for simulations of a detonation-driven flow experiment with uncertainty quantification

Here, uncertainties in the explosive-specific parameters of the Jones-Wilkins-Lee (JWL) equation of state (EOS) are carefully considered in hydrodynamic simulations of an explosive experiment to minimize the error in the flow prediction. Experimental data of the leading shock position in the transverse direction over time serves as the prediction metric for quantifying simulation prediction error. The uncertainty quantification technique, global sensitivity analysis, is utilized to determine the JWL parameters to which the transverse shock propagation is most sensitive. A polynomial response surface (PRS) is constructed in the space of the most influential JWL parameters, and the point of minimum error between the experimental data and the PRS yields calibrated JWL parameters for the experimental flow. The simulation results following the parameter calibration show good agreement with the experimental data. It was found that two significant parameters, the heat release per unit mass of reactant Q and JWL model exponent R 1 are strongly related, which makes it difficult to identify accurate values.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Prediction of the Critical Curvature for LX-17 with the Time of Arrival Data from DNS

We extract the detonation shock front velocity, curvature and acceleration from time of arrival data measured at grid points from direct numerical simulations of a 50mm rate-stick lit by a disk-source, with the ignition and growth reaction model and a JWL equation of state calibrated for LX-17. We compute the quasi-steady (D, κ) relation based on the extracted properties and predicted the critical curvatures of LX-17. We also proposed an explicit formula that contains the failure turning point, obtained from optimization for the (D, κ) relation of LX-17.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2010

Ignition-and-Growth Modeling of NASA Standard Detonator and a Linear Shaped Charge

The main objective of this study is to quantitatively investigate the ignition and shock sensitivity of NASA Standard Detonator (NSD) and the shock wave propagation of a linear shaped charge (LSC) after being shocked by NSD flyer plate. This combined explosive train was modeled as a coupled Arbitrary Lagrangian-Eulerian (ALE) model with LS-DYNA hydro code. An ignition-and-growth (I&G) reactive model based on unreacted and reacted Jones-Wilkins-Lee (JWL) equations of state was used to simulate the shock initiation. Various NSD-to-LSC stand-off distances were analyzed to calculate the shock initiation (or failure to initiate) and detonation wave propagation along the shaped charge. Simulation results were verified by experimental data which included VISAR tests for NSD flyer plate velocity measurement and an aluminum target severance test for LSC performance verification. Parameters used for the analysis were obtained from various published data or by using CHEETAH thermo-chemical code.

Oguz, Sirri↗
DOE OSTI2014

Simulation Study of Near-Surface Coupling of Nuclear Devices vs. Equivalent High-Explosive Charges

A computational study was performed to examine the differences in near-surface ground-waves and air-blast waves generated by high-explosive energy sources and those generated by much higher energy - density low - yield nuclear sources. The study examined the effect of explosive-source emplacement (i.e., height-of-burst, HOB, or depth-of-burial, DOB) over a range from depths of -35m to heights of 20m, for explosions with an explosive yield of 1-kt . The chemical explosive was modeled by a JWL equation-of-state model for a ~14m diameter sphere of ANFO (~1,200,000kg – 1 k t equivalent yield ), and the high-energy-density source was modeled as a one tonne (1000 kg) plasma of ‘Iron-gas’ (utilizing LLNL’s tabular equation-of-state database, LEOS) in a 2m diameter sphere, with a total internal-energy content equivalent to 1 k t . A consistent equivalent-yield coupling-factor approach was developed to compare the behavior of the two sources. The results indicate that the equivalent-yield coupling-factor for air-blasts from 1 k t ANFO explosions varies monotonically and continuously from a nearly perfec t reflected wave off of the ground surface for a HOB ≈ 20m, to a coupling factor of nearly zero at DOB ≈ -25m. The nuclear air - blast coupling curve, on the other hand, remained nearly equal to a perfectly reflected wave all the way down to HOB’s very near zero, and then quickly dropped to a value near zero for explosions with a DOB ≈ -10m. The near - surface ground - wave traveling horizontally out from the explosive source region to distances of 100’s of meters exhibited equivalent - yield coupling - factors t hat varied nearly linearly with HOB/DOB for the simulated ANFO explosive source, going from a value near zero at HOB ≈ 5m to nearly one at DOB ≈ -25m. The nuclear-source generated near-surface ground wave coupling-factor remained near zero for almost all HOB’s greater than zero, and then appeared to vary nearly - linearly with depth-of-burial until it reached a value of one at a DOB between 15m and 20m. These simulations confirm the expected result that the variation of coupling to the ground, or the air, change s much more rapidly with emplacement location for a high-energy-density (i.e., nuclear-like) explosive source than it does for relatively low - energy - density chemical explosive sources. The Energy Partitioning, Energy Coupling (EPEC) platform at LLNL utilizes laser energy from one quad (i.e. 4-laser beams) of the 192 - beam NIF Laser bank to deliver ~10kJ of energy to 1mg of silver in a hohlraum creating an effective small-explosive ‘source’ with an energy density comparable to those in low-yield nuclear devices. Such experiments have the potential to provide direct experimental confirmation of the simulation results obtained in this study, at a physical scale (and time-scale) which is a factor of 1000 smaller than the spatial- or temporal-scales typically encountered when dealing with nuclear explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2015

Detonation Energy Densities from the Cylinder Test

The Cylinder test is a calibrated pipe-bomb used as the only means of obtaining detonation energy densities while the detonation is proceeding [1-2]. The test is close to being a precision test and is no longer the qualitative process found in much of the literature. Most copper tubes have a wall thickness that is 1/5 th the radius as the standard. The output is an annealed copper wall velocity as a function of time. The square of the velocity is proportional to the detonation energy density, E d , at a specific relative volume, v. We want to 1) convert copper wall velocity to detonation energy density, 2) have the process be analytical so it can be applied to hundreds of runs at once, 3) have the resulting JWL equation-of-state give a good result when run in a hydrocode, and 4) provide an accurate record for future use.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Proof-of-Concept Study for Uncertainty Quantification and Sensitivity Analysis using the BRL Shaped-Charge Example

These are the slides for a graduate presentation at Mississippi State University. It covers the following: the BRL Shaped-Charge Geometry in PAGOSA, mesh refinement study, surrogate modeling using a radial basis function network (RBFN), ruling out parameters using sensitivity analysis (equation of state study), uncertainty quantification (UQ) methodology, and sensitivity analysis (SA) methodology. In summary, a mesh convergence study was used to ensure that solutions were numerically stable by comparing PDV data between simulations. A Design of Experiments (DOE) method was used to reduce the simulation space to study the effects of the Jones-Wilkins-Lee (JWL) Parameters for the Composition B main charge. Uncertainty was quantified by computing the 95% data range about the median of simulation output using a brute force Monte Carlo (MC) random sampling method. Parameter sensitivities were quantified using the Fourier Amplitude Sensitivity Test (FAST) spectral analysis method where it was determined that detonation velocity, initial density, C1, and B1 controlled jet tip velocity.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2019

The diameter effect in Bullseye powder

Detonation velocity as a function of charge diameter is reported for Alliant Bullseye powder. Results are compared to those of mixtures of ammonium nitrate mixed with aluminum and ammonium nitrate mixed with fuel oil. In addition, measurements of free surface velocity of flyers in contact with detonating Bullseye are presented and results are compared to those of hydrocode calculations using a Jones–Wilkins–Lee equation of state generated in a thermochemical code. Comparison to the experimental results shows that both the free surface and terminal velocities were under-predicted.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2018

Combined Mini-Cylex & Disk Acceleration Tests in Type K Copper

In this work, we combined the miniature cylinder expansion test (Mini-Cylex), with the Disk Acceleration Test (DAX) using Type K copper, Picatinny Liquid Explosive, and photonic Doppler velocimetry. We estimated the CJ state using plate reverberation methods at the test cap. We extracted velocities at 2, 7, and 10 volume expansions to fit Jones-Wilkins-Lee Equation of State and estimated Gurney velocity at the tube wall. The test also provides an additional method to estimate reaction products Hugoniot through knowledge of the copper test cap. Our experiments and simulations are within expected uncertainty. Lastly, the test and the analysis effectively reduce costs while keeping or increasing fidelity.

42 ENGINEERING↗