Code-to-code comparison of an explosively fired projectile
Explore the source record for details and available documents.
Subramanian, Gopinath Santhana Ganapathy↗
Engineering topics
Ticknor, Christopher
Publications and source records attributed to Ticknor, Christopher.
Uncertainty studies of diameter-effect models applied to HMX-based explosive
Explore the source record for details and available documents.
high explosive, HMX, diameter effect, uncertainty ↗
Atomic cluster expansion potential for large scale simulations of hydrocarbons under shock compression
We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene. A novel approach to selecting structures from DFT MD is also presented, which allows for the rapid selection of unique DFT MD frames from complex trajectories. Comparisons to DFT and experimental data demonstrate that the presented ACE potential accurately reproduces isotherms, carbon melting curves, radial distribution functions, and shock Hugoniots for carbon and hydrocarbon systems for pressures up to 100 GPa and temperatures up to 6000 K for hydrocarbon systems and up to 9000 K for pure carbon systems. This work delivers a potential that can be used for accurate, large-scale simulations of shocked hydrocarbons and demonstrates a methodology for fitting and validating machine learning interatomic potentials to complex molecular environments, which can be applied to energetic materials in future works.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Weak shock compaction on granular salt
This study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories. Mesoscale simulations further enabled advanced analysis of the shock wave’s substructure. In simulations, the shock front’s precursor was shown to have a heterogeneous nature, which is usually modeled as uniform in continuum analyses. The presence of force chains results in a spread out of the shock precursor over a greater ramp distance. With increasing impact velocity, the shock front thickness reduces, and the precursor of the shock front becomes less heterogeneous. Furthermore, mesoscale modeling suggests the formation of force chains behind the shock front, even under the conditions of weak shock. This study presents novel mesoscale simulation results on salt corroborated with data from experiments, thereby characterizing the compaction front speeds in the weak shock regime.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Calibration and uncertainty quantification for Davis Equation of State models for the High Explosive PBX 9501 products
Abstract This paper investigates the uncertainty in the parameters used in the calibration of an Davis Equation Of State (EOS) for the detonation products of the High Explosive PBX 9501. The procedure sought to make use of all available information about this HE to inform the best set of calibration parameters as well as the uncertainty in these parameters. The procedure made use of historical experimental data, the results from thermo‐chemical modeling as well as data on the best isentrope function fit to cylinder test experimental data. Combining all these heterogeneous data sources together in a Bayesian calibration, yielded a posterior mean and covariance. Sampling from the posterior distribution and evaluating an important Quantity Of Interest (QOI) in the EOS model, the detonation speed of a one‐inch rate stick, produced a distribution which showed variations which were in agreement with experiments. The uncertainty in the EOS was reported as eleven sets of model calibrations which spanned the range of this QOI.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models
GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.
42 ENGINEERING↗
On the determination of shocked surface velocities
Laser Doppler velocimetry is commonly fielded in shock physics experiments to diagnose the times of arrival, magnitudes, and durations of shock impulses as they arrive at shock loaded surfaces. However, in the presence of ejecta, the resulting frequency noise often obscures the surface and subsequent impulses at later times. In this context, we evaluate a dominant wave analysis technique applied to velocimetry data from explosively shocked Sn subjected to a second shockwave, on which ejecta are emitted subsequent to both the first and subsequent shockwave impulses. The analysis reveals surface velocities, the times of arrival of subsequent impulses, and their durations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Dominant wave analysis
Explore the source record for details and available documents.
PDV, ejecta physics, signal analysis↗
Instabilities of a vortex-ring-bright soliton in trapped binary three-dimensional Bose-Einstein condensates
Instabilities of vortex-ring-bright coherent structures in harmonically trapped two-component three-dimensional Bose-Einstein condensates are studied numerically within the coupled Gross-Pitaevskii equations and interpreted analytically. Interestingly, the filled vortex core with a sufficiently large amount of the bright component is observed to reduce the parametric interval of stability of the vortex ring. Here, we have identified the mechanisms of several linear instabilities and one nonlinear parametric instability in this connection. Two of the linear instabilities are qualitatively different from ones reported earlier and are associated with azimuthal modes of m = 0 and m = 1, i.e., deviations of the vortex from the stationary ring shape. Our nonlinear parametric resonance instability occurs between the m = 0 and m = 2 modes and signals the exchange of energy between them.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A reactive flow model for the 3,3'-diamino-4,4'-azoxyfurazan based plastic bonded explosive (PBX 9701)
Here, this paper describes the calibration of the Arrhenius Wescott–Stewart–Davis (AWSD) reactive flow model for the recently developed high explosive PBX 9701, which consists of 97% 3,3'-diamino-4,4'-azoxyfurazan (DAAF) and 3% FK-800 binder by weight. DAAF-based explosive formulations have several desirable qualities as they are relatively insensitive to non-shock insults but have higher performance than triaminotrinitrobenzene based formulations. Equations of state for the explosive reactants and products are calibrated using a combination of existing and new experimental data and theoretical calculations. The AWSD rate law calibration utilizes both one-dimensional shock-to-detonation and multidimensional rate stick experiments to capture the shock initiation and propagating detonation regimes. Validation of the calibrated model is demonstrated through comparison with recent gas-gun experiments.
45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
Shock to detonation transition of pentaerythritol tetranitrate (PETN) initially pressed to 1.65 g/cm 3
A novel set of experiments and reactive flow modeling of pentaerythritol tetranitrate (PETN) is presented. In this work, the specific phenomenon of shock to detonation transition is examined, where an initial, relatively weak shock is propagated into pressed PETN powder at 1.65 g/cm 3 and the subsequent buildup to detonation is observed experimentally. These experiments, in conjunction with reactant and products’ equations of state, are utilized for building reactive flow models.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
First principles reactive simulation for equation of state prediction
The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman–Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. Furthermore we demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Static and dynamic properties of multi-ionic plasma mixtures
Complex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multi-component systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion and viscosity for ternary mixtures at extreme conditions. These quantities can be produced from first principles using orbital free molecular dynamics at the computational expense of very intensive simulations to reach good statistics. Utilizing the first-principles results as reference data, we assess the merit of a global analytic model for transport coefficients, "Pseudo-Ions in Jellium" (PIJ), based on an iso-electronic assumption (iso-n e ). With a multi-component hypernetted-chain integral equation, we verify the quality of the iso-n e prescription for describing the static structure of the mixtures. This semi-analytical modeling compares well with the simulation results and allows one to consider plasma mixtures not accessible to simulations. Further, applications are given for the mix of materials in ICF experiments. A reduction of a multicomponent mixture to an effective binary mixture is also established in the hydrodynamic limit and compared with PIJ estimations for ICF relevant mixtures.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Model for the electrical conductivity in dense plasma mixtures
A new density functional theory, average atom based model for the electrical conductivity of dense plasmas with a mixture of ion species, containing no adjustable parameters, is presented herein. The model takes the temperature, mass density and relative abundances of the species as input. It takes into account partial ionization, ionic structure, and core-valence orthogonality, and uses quantum mechanical calculations of cross sections. Comparison to an existing high fidelity but computationally expensive method reveals good agreement. The new model is computationally efficient and can reach high temperatures. A new mixing rule is also presented that gives reasonably accurate conductivities for high temperature plasma mixtures.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Equations of State with Magpie [Slides]
Magpie is a thermochemical equilibrium code developed at LANL. It calculates equations of state for high explosive (HE) mixtures quickly using statistical mechanics.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Equations of state for polyethylene and its shock-driven decomposition products
We construct new equations of state (EOS) for high density and ultrahigh molecular weight polyethylene and their chemical decomposition products under shock loading. The former were built using the SESAME framework, based in part on new specific heat and thermal expansion data reported here. The products EOS was based on thermochemical modeling under the assumption of full thermodynamic and chemical equilibrium. The products are represented as the ideal mixture of bulk carbon in the form of diamond, H 2 , H, and CH 4 . In the process of building a new EOS for the products, we recalibrated our exponential-6 pair potential for methane in order to better agree with data that have appeared since its original parameterization. The polyethylene EOS were calibrated to thermal, thermomechanical, and shock data, and their performance was evaluated in hydrodynamic modeling of deep release experiments reported previously.
36 MATERIALS SCIENCE↗
Quasistable quantum vortex knots and links in anisotropic harmonically trapped Bose-Einstein condensates
The long-time existence of topologically nontrivial configurations of quantum vortices in the form of torus knots and links in trapped Bose-Einstein condensates is demonstrated numerically within the three-dimensional Gross-Pitaevskii equation with an external anisotropic parabolic potential. We identify here parametric domains of trap anisotropy, characterized by the axial over planar frequency ratio ${\lambda}{\approx}1.5-1.6$, where the lifetime of such quasistationary rotating vortex structures is many hundreds of typical rotation times. This suggests the potential experimental observability of the structures. We quantify the relevant lifetimes as a function of the model parameters (e.g., ${\lambda}$) and initial condition parameters of the knot profile.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Enhancement factor of nuclear reactions in partially ionized plasmas and asymmetric mixtures
When hydrogen is mixed with heavy (high‐ Z ) elements, even as traces, such as iron in the sun, the plasma mixture transits from a weakly coupled state to the intermediate coupling regime of simultaneous weak and strong couplings due to charge asymmetry. In those circumstances, partial ionization is an important aspect to consider. We present a global approach for evaluating enhancement factors for thermonuclear reactions in such mixtures. We grounded our approach on extensive orbital free molecular dynamics (OFMD) simulations of hydrogen‐silver mixtures, chosen as a prototype of hydrogen—high‐ Z mixtures, in a large range of temperature (100 < T < 1600 eV). By comparing a series of OFMD simulations to hyper‐netted chain calculations of effective binary ionic mixtures (eBIM), we set a general procedure for computing enhancement factors. The electron screening contribution is also included implicitly in the definition of the effective ionizations within the average atom framework. The agreement between the pair distribution functions of OFMD and eBIM validates a prescription in which each average atom representing a given species shares the same external electron density with the others (iso‐ n e rule). We checked that an interpolation formula between weak and strong couplings in BIM, due to Chugunov and DeWitt, is accurate enough to address the limit of extreme dilution of a heavy element.
Arnault, Philippe↗