Using the Pagosa SURF Model to Simulate Fragment Impact on Energetic Materials for Safety Applications
Explore the source record for details and available documents.
Ma, Xia↗
SEARCH · Engineering Papers
Results for “SURF model”
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.
Re-calibration of PBX9501 SURF model
PBX 9501 is a plastic bonded explosive composed of 95 wt % HMX and a binder; see [Gibbs and Popolato, 1980, pp. 109–119]. SURF is a reactive burn model for shock initiation and propagation of detonation waves. It has previously been calibrated for PBX 9501. Here the SURF model is recalibrated for PBX 9501; specifically, lot 730-010 at ρ = 1.837 g/cc 3 . The new calibration uses the Davis reactants and products EOS calibrated for the AWSD model [Aslam et al., 2020]. The burn rate in the shock initiation regime is fit to the Pop plot from 5 embedded gauge shock-to-detonation transition (SDT) experiments from [Gustavsen et al., 1999, see fig 12 and table 5]. In the propagation regime, the burn rate is fit to curvature effect data (detonation speed as function of front curvature); see [Aslam, 2007]. Also the burn parameters are adjusted to fit the gap-stick experiment [Hill et al., 2018]. Simulating the detonation wave speed in this experiment requires a model that is accurate for initiation with complex shock loading; in particular, a pressure decreasing gradient behind a curved lead shock. This is more difficult than calibrating to the standard SDT experiments which are 1-D and driven by a sustained shock. Simulations of the gap-stick experiment for PBX 9501 with the SURF model will be discussed in a subsequent report.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Verification test of the SURF model in xRage: Part IV cylindrically converging detonation wave
Previous verification studies of the SURF reactive burn model in the xRage code used test problems for underdriven or self-supporting detonation waves. Here we consider overdriven detonation waves. As a test problem we use a cylindrical converging detonation wave. Since there is no exact solution, the numerical solutions at different resolutions are compared for 1-D cylindrical geometry and 2-D Cartesian geometry.
97 MATHEMATICS AND COMPUTING↗
SURF Model Resolution in xRage Code
A key issue for detonation wave simulations with a reactive burn model is the resolution needed for the reaction zone. The reaction-zone width is a physical length scale. For plastic-bonded explosives, it is small (on the order of 100 microns) but affects important detonation properties such as the detonation speed. Even with adaptive mesh refinement, a well resolved reaction zone is computationally expensive. Here we consider the trade-off between computational expense and accuracy when the resolution of the reaction zone is varied for simulations with the xRage code of PBX 9501 with the SURF burn model and of PBX 9502 with the SURFplus burn model. For assessing shock initiation, the numerical Pop plot is simulated, and for assessing propagating detonation waves, the cylinder test is used as a test problem. The reaction-zone resolution is varied from a cell size of 6.25 to 200 microns, and key properties of the solutions are compared. With 6 micron resolution, the numerical solution is near mesh converged and can be used to judge the accuracy of the coarser resolution simulations. For the propagating detonation wave in the cylinder test, at 200 microns the detonation speed is low by a couple of per cent. In terms of cpu-hours, the finest resolution is about a factor of 100 more expensive than the coarsest resolution.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Effects of Void Morphology on Detonation Initiation of PBX 9502
Accurate predictive modeling of high explosives in abnormal conditions, such as a fire, is critical to personnel safety. Modeling depends heavily on precise physical characterization of the high explosive in question. The aim of this thesis was to characterize the detonation response to shock and the microstructure at high-temperature (250°C) of PBX 9502 and use these data to attempt to model these results using a new reactive burn model called SURF (Scaled Uniform Reactive Front). In order to characterize the shock response of PBX 9502 at high-temperature a 1D gas gun experiment was designed and executed that provided shock to detonation transition data. These data were then used to do a preliminary calibration of SURF for PBX 9502 at 250°C. Small-angle neutron scattering (SANS) was employed to characterize the void morphology of PBX 9502 as a function of temperature. This was done in order to feed a modified version of SURF called physically informed - SURF (π-SURF) which uses the void size and quantity distribution data as measured by SANS along with the deflagration and detonation characteristics of a high explosive to inform the SURF model. The π-SURF model was then used to model both ambient and high-temperature PBX 9502 SDT response and the results were compared to the physically measured results. It was found that π-SURF works in ambient cases. At high-temperature, however, it was determined that not enough information is known about hot spot formation to use the π-SURF model for high-temperature scenarios.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
ZND Verification Tests for Reactive Burn Models in FLAG
The ZND theory, named after Zeldovich, von Neumann, and Doering, provides a simple model for one-dimensional ideal steady-state detonation. It assumes that the detonation wave front starts with a shock that is a discontinuous jump and is followed by a finite-length reaction zone. Reactive burn (also called reactive flow) models are based on ZND theory, as they model the shock initiation and detonation process with a finite reaction rate. The ZND wave propagation test is essentially the only available test case where an analytic solution exists for verification of reactive burn models in numerical codes. However, there are extensions and variants of the ZND test that have been devised for verification of multidimensional flows. The objective of this work is to provide verification of the reactive burn models currently implemented in the Lagrangian hydrocode FLAG and investigate the influence of mesh resolution, artificial viscosity models, and the Arbitrary Lagrangian-Eulerian (ALE) Euler relaxer on the simulation results. The burn models of interest are the Wescott-Stewart-Davis (WSD) model, the Scaled Uniform Reactive Flow (SURF) model (specifically with the SURFplus model extension), and the Arrhenius shock temperature state dependent WSD (AWSD) model. Previously, Ralph Menikoff has used ZND tests for verification of the SURF and SURFplus models in the Eulerian hydrocode xRAGE. The ZND tests here are somewhat different than the approach by Menikoff. In particular, we use a piston-driven ZND detonation wave (via a prescribed constant velocity boundary condition) in a Lagrangian framework whereas Menikoff had a ZND wave followed by a invariant rarefaction wave. The xRAGE simulations were carried out on uniform grids and adaptive mesh refinement (AMR) grids. Although AMR was recently implemented in FLAG and now fully functional for 2D simulations, it will not be evaluated with ZND tests at this time. Some work has been done previously for validating the reactive burn models in FLAG. For example, the validation studies of SURF with shock-to-detonation (SDT) tests, cylinder tests, and gap-stick tests. Further validation of the AWSD, WSD, and SURF models is described in for SDT, multi-shock, cylinder, and corner-turning tests. Recently, a large parameter study with approximately twenty SDT tests was performed to validate the AWSD, WSD, SURF, and SURFplus models while also investigating mesh resolution and artificial viscosity settings. To the best of our knowledge, the current work represents the first documented verification of these burn models in FLAG. The reader should note that many of the tables and figures in this report use units of cm/g/µs, which are the default units for FLAG. However, some lengths and velocities are expressed in µm and mm/µs, respectively, which are typical units for detonation analyses.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
A chemistry-informed hybrid machine learning approach to predict metal adsorption onto mineral surfaces
Historically, surface complexation model (SCM) constants and distribution coefficients (K d ) have been employed to quantify mineral-based retardation effects controlling the fate of metals in subsurface geologic systems. Our recent SCM development workflow, based on the Lawrence Livermore National Laboratory Surface Complexation/Ion Exchange (L-SCIE) database, illustrated a community FAIR data approach to SCM development by predicting uranium(VI)-quartz adsorption for a large number of literature-mined data. Here, we present an alternative hybrid machine learning (ML) approach that shows promise in achieving equivalent high-quality predictions compared to traditional surface complexation models. At its core, the hybrid random forest (RF) ML approach is motivated by the proliferation of incongruent SCMs in the literature that limit their applicability in reactive transport models. Our hybrid ML approach implements PHREEQC-based aqueous speciation calculations; values from these simulations are automatically used as input features for a random forest (RF) algorithm to quantify adsorption and avoid SCM modeling constraints entirely. Named the LLNL Speciation Updated Random Forest (L-SURF) model, this hybrid approach is shown to have applicability to U(VI) sorption cases driven by both ion-exchange and surface complexation, as is shown for quartz and montmorillonite cases. The approach can be applied to reactive transport modeling and may provide an alternative to the costly development of self-consistent SCM reaction databases.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Comments on the SURFplus reactive burn model
The SURFplus reactive burn model is intended for heterogeneous solid HE that produce an excess amount of carbon. It utilizes a fast hotspot reaction (SURF model) plus a second slow reaction for the energy release from carbon clustering. The fast reaction dominates shock initiation, though the fast rate parameters are slightly affected by the late energy release from the slow reaction. The slow reaction has a significant affect on propagating detonation waves. Propagating detonation waves are characterized by the curvature effect; i.e., detonation speed $D_n$ as function of the local front curvature κ. The $D_n(κ)$ curves have a qualitative different shape. This results from the sensitivity of the slope of the $D_n(κ)$ curve to the reaction-zone width up to the sonic point, and the sonic point moving as κ increases from the end of the slow reaction to the neighborhood of the end of the fast reaction. For PBX 9502, the steep slope for small curvature is determined by the slow reaction, while the lower slope for larger curvature is determined by the fast reaction rate in the high detonation pressure regime. The remainder of this paper is organized as follows. In section 2 we review the assumptions of the Shaw and Johnson carbon cluster model used by the SURFplus model for the slow reaction. Section 3 describes in detail the mixture model for a partly burned HE. The form of the slow rate used by the SURFplus model is described in section 4. We conclude with a brief discussion of applying a burn model with 2 rates to either a metalized or a composite PBX.
42 ENGINEERING↗
A Verification and Validation Study of the Cyclops I PBX 9502 Experiment: A validation of the SURFplus High Explosive Model [Slides]
SURF is an acronym for Scaled-Unified-Reactive-Front, invented by R. Menikoff and S. Shaw. What follows is a brief synopsis of this model as used for the validation study discussed here. For more details see the xrage user manual and the references cited there.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Simulating PBX 9501 Gap-stick Experiment
It has previously been shown that the rate parameters of the SURF reactive burn model can be adjusted to fit the gap-stick data for PBX 9501 [Johnson et al., 2018, fig 4]. A recent recalibration of the SURF model for PBX 9501 (lot 730-010 at ρ = 1.837 g/cc 3 ) used different reactants and products EOS and an updated fitting form for the burn rate [Menikoff, 2021]. The new calibration fit gap-stick data [Hill et al., 2021] together with data from SDT experiments and curvature effect data. Here the results of gap-stick simulations using the xRage code are shown. They provide insight into the reactive flow in the gap-stick experiment.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Approach of a planar detonation to a steady CJ detonation wave
Generating a planar steady detonation wave in a high explosive is important for accurately measuring key parameters such as the CJ detonation state and the reaction-zone width. A detonation wave can be promptly shock initiated in the TATB based explosive PBX 9502 due to a fast hotspot reaction. However, a slow carbon clustering reaction causes the detonation wave to have a slow approach to steady state. Here we study the approach of the reaction zone of a planar detonation to the steady ZND profile of an underdriven detonation wave. For PBX 9502 simulations were run with the SURFplus reactive burn model. For comparison, simulations were also run with the HMX based PBX 9501. The SURF reactive model is used for PBX 9501 since it has only a fast hotspot reaction.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SurFE-XD (Surface curvature-driven Finite Elements model for Diffusion under eXtreme conditions)
SurFE-XD is a mesoscale finite element framework to model surface diffusion under mutliphysics environments. The code uses legacy C++ library dolphin wrapped with python in a FEniCS driven unified form language and just-in-time (JIT) compilation setting. The purpose of the release is to attract wide-ranging usage of the code along with publication supplementation to support reproducibility of scientific data. SurFE-XD has been originally conceived under the LDRD-DR funding for “High-Gradient (C-BAND) Breakdown tolerant accelerator materials project. Currently SrFE-XD support electrostatics and Thermo-elasticity driven surface diffusion kernels. Releasing the code will also enable to include contributions from other physical regimes e.g., plasticity and electrodynamics etc as well portability to GPU-based platforms.
Bagchi, Soumendu↗
Synopsis of SURF and SURFplus
Developed by Menikoff and Shaw, the SURF reactive burn model builds on the Ignition and Growth concept by incorporating the lead shock pressure directly into the volumetric hot-spot burn rate. The plus-extension augments the model with a late time surface burn rate due to carbon clustering. Together, SURF and SURFplus allow for robust modeling of both conventional and insensitive high explosives. Practically, the SURF burn models require shock detection and the advection of the lead shock pressure as an additional material field. Calibrated parameters for given equations of state and thermodynamic closure are dependent on initial temperature and density. An implementation in Python is given.
42 ENGINEERING↗
Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502
High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Tools for 3D scientific visualization in computational aerodynamics at NASA Ames Research Center
Hardware, software, and techniques used by the Fluid Dynamics Division (NASA) for performing visualization of computational aerodynamics, which can be applied to the visualization of flow fields from computer simulations of fluid dynamics about the Space Shuttle, are discussed. Three visualization techniques applied, post-processing, tracking, and steering, are described, as well as the post-processing software packages used, PLOT3D, SURF (Surface Modeller), GAS (Graphical Animation System), and FAST (Flow Analysis software Toolkit). Using post-processing methods a flow simulation was executed on a supercomputer and, after the simulation was complete, the results were processed for viewing. It is shown that the high-resolution, high-performance three-dimensional workstation combined with specially developed display and animation software provides a good tool for analyzing flow field solutions obtained from supercomputers.
Bancroft, Gordon↗
U-Surf: a global 1 km spatially continuous urban surface property dataset for kilometer-scale urban-resolving Earth system modeling
High-resolution urban climate modeling has faced substantial challenges due to the absence of a globally consistent, spatially continuous, and accurate dataset to represent the spatial heterogeneity of urban surfaces and their biophysical properties. This deficiency has long obstructed the development of urban-resolving Earth system models (ESMs) and ultra-high-resolution urban climate modeling, over large domains. Here, we present U-Surf, a first-of-its-kind 1 km resolution present-day (circa 2020) global continuous urban surface parameter dataset. Using the urban canopy model (UCM) in the Community Earth System Model as a base model for satisfying dataset requirements, U-Surf leverages the latest advances in remote sensing, machine learning, and cloud computing to provide the most relevant urban surface biophysical parameters, including radiative, morphological, and thermal properties, for UCMs at the facet and canopy level. Generated using a systematically unified workflow, U-Surf ensures internal consistency among key parameters, making it the first globally coherent urban canopy surface dataset. U-Surf significantly improves the representation of the urban land heterogeneity both within and across cities globally; provides essential, high-fidelity surface biophysical constraints to urban-resolving ESMs; enables detailed city-to-city comparisons across the globe; and supports next-generation kilometer-resolution Earth system modeling across scales. U-Surf parameters can be easily converted or adapted to various types of UCMs, such as those embedded in weather and regional climate models, as well as air quality models. The fundamental urban surface constraints provided by U-Surf can also be used as features for machine learning models and can have other broad-scale applications for socioeconomic, public health, and urban planning contexts. We expect U-Surf to advance the research frontier of urban system science, climate-sensitive urban design, and coupled human–Earth systems in the future. The dataset is publicly available at https://doi.org/10.5281/zenodo.11247598 (Cheng et al., 2024).
Cheng, Yifan [Univ. of Illinois at Urbana-Champaig↗
Inhomogeneous Enrichment of Radioactive Nuclei in the Galaxy: Deposition of Live 53 Mn, 60 Fe, 182 Hf, and 244 Pu into Deep-sea Archives. Surfing the Wave?
Abstract While modeling the galactic chemical evolution (GCE) of stable elements provides insights to the formation history of the Galaxy and the relative contributions of nucleosynthesis sites, modeling the evolution of short-lived radioisotopes (SLRs) can provide supplementary timing information on recent nucleosynthesis. To study the evolution of SLRs, we need to understand their spatial distribution. Using a three-dimensional GCE model, we investigated the evolution of four SLRs: 53 Mn, 60 Fe, 182 Hf, and 244 Pu with the aim of explaining detections of recent (within the last ≈1–20 Myr) deposition of live 53 Mn, 60 Fe, and 244 Pu of extrasolar origin into deep-sea reservoirs. We find that core-collapse supernovae are the dominant propagation mechanism of SLRs in the Galaxy. This results in the simultaneous arrival of these four SLRs on Earth, although they could have been produced in different astrophysical sites, which can explain why live extrasolar 53 Mn, 60 Fe, and 244 Pu are found within the same, or similar, layers of deep-sea sediments. We predict that 182 Hf should also be found in such sediments at similar depths.
79 ASTRONOMY AND ASTROPHYSICS↗
Experimental and modelling study of the effect of pre-shock duration on the double shock initiation response of PBX 9501
Here, a series of double shock initiation gas-gun experiments have been performed on PBX 9501 (95% weight HMX, 5% binder) in which the input conditions were nominally the same, but the time delay between the two shocks was explicitly varied. The impact conditions generated pre- and main shock pressures of 2.9 and 6.2 GPa respectively, with time delays between them ranging from 0.56 to 1.10 µs. This resulted in different levels of reaction behind the pre-shock prior to the main shock traversing the explosive samples, with the growth of reaction monitored using embedded particle velocity and shock tracker gauges. Simulated results from three reactive burn models (CREST, AWSD and SURF) are compared with the measured data to assess their ability to predict the response of PBX 9501 under such double shock conditions. The described experiments add to the body of evidence to help improve our understanding of the double shock initiation response of explosives, and the data obtained provide an excellent test for reactive burn models.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗