Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock dynamics”

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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Potential energy surfaces for high-energy N + O 2 collisions

Potential energy surfaces for high-energy collisions between an oxygen molecule and a nitrogen atom are useful for modeling chemical dynamics in shock waves. In the present work, we present doublet, quartet, and sextet potential energy surfaces that are suitable for studying collisions of O 2 ( 3 Σ$^{–}_{g}$) with N( 4 S) in the electronically adiabatic approximation. Two sets of surfaces are developed, one using neural networks (NNs) with permutationally invariant polynomials (PIPs) and one with the least-squares many-body (MB) method, where a two-body part is an accurate diatomic potential and the three-body part is expressed with connected PIPs in mixed-exponential-Gaussian bond order variables (MEGs). We find, using the same dataset for both fits, that the fitting performance of the PIP-NN method is significantly better than that of the MB-PIP-MEG method, even though the MB-PIP-MEG fit uses a higher-order PIP than those used in previous MB-PIP-MEG fits of related systems (such as N 4 and N 2 O 2 ). However, the evaluation of the PIP-NN fit in trajectory calculations requires about 5 times more computer time than is required for the MB-PIP-MEG fit.

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

The Modified Void Nucleation and Growth Model (MNAG) for Damage Evolution in BCC Ta

A void coalescence term was proposed as an addition to the original void nucleation and growth (NAG) model to accurately describe void evolution under dynamic loading. The new model, termed as modified void nucleation and growth model (MNAG model), incorporated analytic equations to explicitly account for the evolution of the void number density and the void volume fraction (damage) during void nucleation, growth, as well as the coalescence stage. The parameters in the MNAG model were fitted to molecular dynamics (MD) shock data for single-crystal and nanocrystalline Ta, and the corresponding nucleation, growth, and coalescence rates were extracted. The results suggested that void nucleation, growth, and coalescence rates were dependent on the orientation as well as grain size. Compared to other models, such as NAG, Cocks–Ashby, Tepla, and Tonks, which were only able to reproduce early or later stage damage evolution, the MNAG model was able to reproduce all stages associated with nucleation, growth, and coalescence. The MNAG model could provide the basis for hydrodynamic simulations to improve the fidelity of the damage nucleation and evolution in 3-D microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigating 3,4-bis(3-nitrofurazan-4-yl)furoxan detonation with a rapidly tuned density functional tight binding model

In this work, we describe a machine learning approach to rapidly tune density functional tight binding models for the description of detonation chemistry in organic molecular materials. Resulting models enable simulations on the several 10s of ps scales characteristic to these processes, with “quantum-accuracy.” We use this approach to investigate early shock chemistry in 3,4-bis(3-nitrofurazan-4-yl)furoxan, a hydrogen-free energetic material known to form onion-like nanocarbon particulates following detonation. We find that the ensuing chemistry is significantly characterized by the formation of large C x N y O z species, which are likely precursors to the experimentally observed carbon condensates. Beyond utility as a means of investigating detonation chemistry, the present approach can be used to generate quantum-based reference data for the development of full machine-learned interatomic potentials capable of simulation on even greater time and length scales, i.e., for applications where characteristic time scales exceed the reach of methods including Kohn–Sham density functional theory, which are commonly used for reference data generation.

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

Dynamic crushing of metal lattice metamaterials: Shock mode diagrams and transition to topology-independent compaction regime

Additively manufactured lattice metamaterials offer design versatility in strength and energy absorption and provide an additional degree of freedom through the selection of the lattice topology. Under quasistatic loading, the unit cell structure can strongly affect the stiffness, yield, and post-yield behavior, but whether and to what degree the effect of lattice topology persists into dynamic loading scenarios, up to the compaction shock regime, has not been established. LLNL’ s ALE3D hydrocode was used to perform a computational investigation of dynamic loading in multiple lattice types, including the gyroid, octet, Schwarz D, and rhombic dodecahedron, under impact velocities from 0.25 to 2.25 km/s. Shock Hugoniots for each lattice topology are generated and compared, suggesting that above a critical velocity, distinctions between architectures may not persevere and compacted lattices behave similarly. Here, to investigate the transition between topology-dependent quasistatic compression and the topology-independent regime above the critical velocity, a one-dimensional elastic-linear hardening plasticity-densified solid (E-LHP-DS) shock model for lattice materials was developed that relies upon confined compression to link the quasistatic and shock mechanics. Unlike similar works, the model does not assume rigid behavior prior to yield or locking behavior at densification, allowing a richer exploration of lattice mechanics. With only six parameters, the analytical model simultaneously fit quasistatic confined compression simulations for relative densities 0.1 $≤ \bar{ρ} ≤$ 0.9 and predicted dynamic compaction behavior to traverse several distinct shock modes, each defined by a critical impact speed (equivalently, critical stresses). Comparing the numerical results to the one-dimensional E-LHP-DS shock model predictions suggests that the topology-independence under strong shocks is linked to the onset of densification, which can be predicted based on quasistatic confined compression results.

Cellular material↗
NASA NTRS2015

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗
DOE OSTI2025

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2021

Biermann battery magnetic fields in ICF capsules: Total magnetic flux generation

This paper focuses on the process of magnetic flux generation in inertial confinement fusion implosions. Hot-spots are shown to be dominated by fields generated during stagnation when the temperature and density gradients are largest. A scaling of hot-spot magnetic flux is derived and compared with simulations, revealing that perturbations with both larger amplitudes and higher mode numbers generate more magnetic flux. Furthermore, the model allows for greater understanding of which target designs will be susceptible to magnetohydrodynamic effects. For example, the model can be used to ascertain the time when most magnetic flux is generated. If generation is weighted more toward early times, then more high-mode magnetic field loops will be present. A hot-spot with no high-mode perturbations at time of peak neutron production can still contain significant magnetic flux on those scales. By assuming that magnetic flux is deposited at the hot-spot edge by Nernst advection, the model can be used to post-process radiation-hydrodynamics data to estimate magnetic field strengths and magnetizations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inferring the temperature profile of the radiative shock in the COAX experiment with shock radiography, Dante, and spectral temperature diagnostics

Predicting and modeling the behavior of experiments with radiation waves propagating through low-density foams require a detailed quantification of the numerous uncertainties present. In regimes where a prominent radiative shock is produced, key dynamical features include the shock position, temperature, and curvature and the spatial distribution and temperature of the corresponding supersonic radiation wave. The COAX experimental platform is designed to constrain numerical models of such a radiative shock propagating through a low-density foam by employing radiography for spatial and shock information, Dante for characterizing the x-ray flux from the indirectly driven target, and a novel spectral diagnostic designed to probe the temperature profile of the wave. In this work, we model COAX with parameterized 2D simulations and a Hohlraum-laser modeling package to study uncertainties in diagnosing the experiment. The inferred temperature profile of the COAX radiation transport experiments has been shown to differ from simulations more than expected from drive uncertainties that have been constrained by simultaneous soft x-ray flux and radiography measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Enforcing detailed balance in the Borgnakke–Larsen redistribution method with temperature dependent relaxation models

For decades, it has been observed that the commonly used Borgnakke–Larsen method for energy redistribution in Direct Simulation Monte Carlo codes fails to satisfy the principle of detailed balance when coupled to a wide variety of temperature dependent relaxation models, while seemingly satisfying detailed balance when coupled to others. Many attempts have been made to remedy the issue, yet much ambiguity remains, and no consensus appears in the literature regarding the root cause of the intermittent compatibility of the Borgnakke–Larsen method with temperature dependent relaxation models. This paper alleviates that ambiguity by presenting a rigorous theoretical derivation of the Borgnakke–Larsen method's requirement for satisfying detailed balance. Specifically, it is shown that the Borgnakke–Larsen method maintains detailed balance if and only if the probability of internal-energy exchange during a collision depends only on collision invariants (e.g., total energy). Here, the consequences of this result are explored in the context of several published definitions of relaxation temperature, including translational, total, and cell-averaged temperatures. Of particular note, it is shown that cell-averaged temperatures, which have been widely discussed in the literature as a way to ensure equilibrium is reached, also fail in a similar, although less dramatic, fashion when the aforementioned relationship is not enforced. The developed theory can be used when implementing existing or new relaxation models and will ensure that detailed balance is satisfied.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Mathematical modelling of the concave front in the adjacent high explosive detonation problem

This study presents an analysis of the transition-zone in adjacent high explosive (HE) detonation problems which uses a $D, 𝜅, \dot{D}$ relationship, where $D$ is the detonation front-normal velocity, 𝜅 is the detonation front curvature and $\dot{D}$ is the time derivative of detonation front-normal velocity. Our approach extends the traditional $(D, 𝜅)$ model to accurately predict the behaviour of both diverging and converging detonation shock fronts. Our findings affirm that a hyperbolic type of front evolution equation, enhanced with wave acceleration, provides a robust framework for modelling complex shock front dynamics in HE materials. This approach not only captures the natural effects of straightness and boundary slope jumps in the transition-zone but also bridges the gap between mathematical predictions and experimental observations, offering insights into the behaviour of both diverging and converging detonation propagations in a homogeneous HE.

acceleration↗
DOE OSTI2024

Impact of mid- Z gas fill on dynamics and performance of shock-driven implosions at the OMEGA laser

Shock-driven implosions with 100% deuterium (D 2 ) gas fill compared to implosions with 50:50 nitrogen-deuterium (N 2 ⁢D 2 ) gas fill have been performed at the OMEGA laser facility to test the impact of the added mid-Ζ fill gas on implosion performance. Ion temperature (Τ ion ) as inferred from the width of measured DD-neutron spectra is seen to be 34%±6% higher for the N 2⁢ D 2 implosions than for the D 2 -only case, while the DD-neutron yield from the D 2 -only implosion is 7.2±0.5 times higher than from the N 2⁢ D 2 gas fill. The T ion enhancement for N 2 ⁢D 2 is observed in spite of the higher Z, which might be expected to lead to higher radiative loss, and higher shock strength for the D 2 -only versus N 2 ⁢D 2 implosions due to lower mass, and is understood in terms of increased shock heating of N compared to D, heat transfer from N to D prior to burn, and limited amount of ion-electron-equilibration-mediated additional radiative loss due to the added higher-Z material. Further, this picture is supported by interspecies equilibration timescales for these implosions, constrained by experimental observables. The one-dimensional (1D) kinetic Vlasov-Fokker-Planck code ifp and the radiation hydrodynamic simulation codes hyades (1D) and xrage [1D, two-dimensional (2D)] are brought to bear to understand the observed yield ratio. Comparing measurements and simulations, the yield loss in the N 2 ⁢D 2 implosions relative to the pure D 2 -fill implosion is determined to result from the reduced amount of D 2 in the fill (fourfold effect on yield) combined with a lower fraction of the D 2 fuel being hot enough to burn in the N 2 ⁢D 2 case. The experimental yield and T ion ratio observations are relatively well matched by the kinetic simulations, which suggest interspecies diffusion is responsible for the lower fraction of hot D 2 in the N 2 ⁢D 2 relative to the D 2 -only case. The simulated absolute yields are higher than measured; a comparison of 1D versus 2D XRAGE simulations suggest that this can be explained by dimensional effects. The hydrodynamic simulations suggest that radiative losses primarily impact the implosion edges, with ion-electron equilibration times being too long in the implosion cores. The observations of increased T ion and limited additional yield loss (on top of the fourfold expected from the difference in D content) for the N 2 ⁢D 2 versus D 2 -only fill suggest it is feasible to develop the platform for studying CNO-cycle-relevant nuclear reactions in a plasma environment.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Geophysical Observations of the 2023 September 24 OSIRIS-REx Sample Return Capsule Reentry

Sample return capsules (SRCs) entering Earth's atmosphere at hypervelocity from interplanetary space are a valuable resource for studying meteor phenomena. The 2023 September 24 arrival of the Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer SRC provided an unprecedented chance for geophysical observations of a well-characterized source with known parameters, including timing and trajectory. A collaborative effort involving researchers from 16 institutions executed a carefully planned geophysical observational campaign at strategically chosen locations, deploying over 400 ground-based sensors encompassing infrasound, seismic, distributed acoustic sensing, and Global Positioning System technologies. Additionally, balloons equipped with infrasound sensors were launched to capture signals at higher altitudes. This campaign (the largest of its kind so far) yielded a wealth of invaluable data anticipated to fuel scientific inquiry for years to come. The success of the observational campaign is evidenced by the near-universal detection of signals across instruments, both proximal and distal. This paper presents a comprehensive overview of the collective scientific effort, field deployment, and preliminary findings. The early findings have the potential to inform future space missions and terrestrial campaigns, contributing to our understanding of meteoroid interactions with planetary atmospheres. Furthermore, the data set collected during this campaign will improve entry and propagation models and augment the study of atmospheric dynamics and shock phenomena generated by meteoroids and similar sources.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗