Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dense hydrogen layers for high performance MagLIF

We report that Magnetized Liner Inertial Fusion (MagLIF) experiments driven by the Z machine produce >10 13 deuterium-deuterium fusion reactions. Simulations indicate high yields and gains (1000) with increased current and deuterium-tritium layers for burn propagation. Such a coating also isolates the metal liner from the gaseous fuel, which should reduce mixing of liner material into the fuel. However, the vapor density at the triple point is only 0.3 kg/m3, which is not high enough for MagLIF operation. We present two solutions to this problem. First, a fuel wetted low-density plastic foam can be used to form a layer on the inside of the liner. The desired vapor density can be obtained by controlling the temperature. This does however introduce carbon into the layer which will enhance radiation losses. Simulations indicate that this wetted foam layer can significantly contribute to the fusion yield when the foam density is less than 35 kg/m 3 . Second, we show that a pure frozen fuel layer can first be formed on the inside of the liner and then low temperature gaseous fuel can be introduced just before the implosion without melting a significant amount of the ice layer. This approach is the most promising for MagLIF to produce high yield and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

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 OSTI2024

A comparison of smartphone and infrasound microphone data from a fuel air explosive and a high explosive

For prompt detection of large (>1 kt) above-ground explosions, infrasound microphone networks and arrays are deployed at surveyed locations across the world. Denser regional and local networks are deployed for smaller explosions, however, they are limited in number and are often deployed temporarily for experiments. With the expanded interest in smaller yield explosions targeted at vulnerable areas such as population centers and key infrastructures, the need for more dense microphone networks has increased. An “attritable” (affordable, reusable, and replaceable) and flexible alternative can be provided by smartphone networks. Explosion signals from a fuel air explosive (thermobaric bomb) and a high explosive with trinitrotoluene equivalent yields of 6.35 and 3.63 kg, respectively, were captured on both an infrasound microphone and a network of smartphones. The resulting waveforms were compared in time, frequency, and time-frequency domains. The acoustic waveforms collected on smartphones produced a filtered explosion pulse due to the smartphone's diminishing frequency response at infrasound frequencies (<20 Hz) and was found difficult to be used with explosion characterization methods utilizing waveform features (peak overpressure, impulse, etc.). However, the similarities in time frequency representations and additional sensor inputs are promising for other explosion signal identification and analysis. As an example, a method utilizing the relative acoustic amplitudes for source localization using the smartphone sensor network is presented.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

1-D coupled surface flow and transport equations revisited via the physics-informed neural network approach

The de Saint-Venant equation (SVE) and advection–diffusion equation (ADE) are commonly employed to solve solute transport problems in surface water. In this work, we propose a mesh-free method based on the physics-informed neural network (PINN) to solve the one dimensional (1-D) SVE, ADE, and the coupled SVE and ADE (SVE-ADE) under various initial and boundary conditions. The PINN model extends the architecture of deep neural networks (DNNs) with implementation of loss function, which are additionally subject to constraints imposed by the physical laws of SVE and ADE, along with their initial and boundary conditions. In such a manner, PINNs can be quickly steered to the true solution while obeying the physical laws. The results of PINN model are compared with the analytical and/or numerical solutions under various conditions to investigate its accuracy and efficiency in solving the SVE, ADE, and SVE-ADE. Our results indicate PINN can accurately simulate the shock wave morphology and avoid numerical dissipation in unsteady flow condition. The PINN method outweighed traditional numerical methods in several aspects, including its ability to function with small amounts of data, no grid discretization, and random selection of sampling points, etc. Additionally, the PINN method is also suitable for solving inverse problems with sparse and noisy data. With 1% noise and 2000 initial and boundary condition points (N u ), the errors of the estimated flow rate (v) and diffusion coefficient (D) are 0.003% and 0.105%, respectively, which indicate the accuracy and robustness of the proposed method. Finally, our results indicate the capability and robustness of the proposed PINN methodology for solving multi-physics problems, irrespective of the presence of sparse and noisy data in the training dataset.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Magnetic reconnection and kinetic waves generated in the Earth's quasi-parallel bow shock

Magnetic reconnection in quasi-parallel shocks, relevant to the Earth’s bow shock, is studied by means of two-dimensional full particle-in-cell simulations. As the Alfvénic Mach number increases, the propagation direction of the waves excited in the transition region changes, and the shock becomes more turbulent with more reconnection sites. In the higher Mach number shock, abundant electron-only reconnection sites are generated with scales on the order of the ion skin depth or less. Non-reconnecting current sheets can also generate electron jets and energy dissipation can occur there as well. However, non-reconnecting current sheets with the magnetic field reversal typically show a smaller energy dissipation rate than reconnecting current sheets. In the shock transition region, two types of waves are responsible for driving reconnection: one has a wavelength around three ion skin depths (d i ), and the other has a wavelength less than 1 d i . Here, electron and ion distribution functions show that in regions where the former type of waves is excited, there are two ion beams and a single-peaked electron distribution. In contrast, in regions where the latter type of waves is excited, there are multiple electron and ion beams. The waves propagating obliquely to the magnetic field bend the magnetic field lines, and magnetic reconnection occurs where oppositely directed field lines come into contact.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Observations of multi-ion physics and kinetic effects in a surrogate to the solar CNO reactions

The ‘CNO process’ occurs in heavier stars with finite metallicity in which hydrogen burning is catalyzed in the presence of 12 C. These reactions are more strongly dependent on temperature than the pp cycle reactions, and thus the CNO cycle dominates only in massive stars. For these types of reactions to be studied at ICF facilities such as OMEGA, an implosion platform using heavier nuclei in the fuel and capable of creating ion temperatures on the order of at least 20 keV is required. A potential route to reach these conditions is to take advantage of kinetic effects in low-convergence shock-driven ‘exploding pusher’ implosions. In this experiment, shots were conducted at the OMEGA laser facility using the surrogate reaction 13 C + D. Its cross section is substantially higher than the actual astrophysical CNO reactions. The yield of this reaction in these implosions was much lower than expected. Finally, physical explanations are discussed, with significant species stratification the likely explanation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Maximum Energy of Particles Accelerated in Gamma-Ray Burst Afterglow Shocks

Particle acceleration in relativistic collisionless shocks remains an open problem in high-energy astrophysics. Particle-in-cell (PIC) simulations predict that electron acceleration in weakly magnetized shocks proceeds via small-angle scattering, leading to a maximum electron energy significantly below the Bohm limit. This upper bound on electron energy manifests observationally as a characteristic synchrotron cutoff, providing a direct probe of the underlying acceleration physics. Gamma-ray burst (GRB) afterglows offer an exceptional laboratory for testing these predictions. Here, we model the spectral evolution of GRB afterglows during the relativistic deceleration phase, incorporating PIC-motivated acceleration prescriptions and self-consistently computing synchrotron and synchrotron self-Compton emission. We find that low-energy bursts in low-density environments, typical of short GRBs, exhibit a pronounced synchrotron cutoff in the GeV band within minutes to hours after the trigger. Applying our framework to GRB 190114C and GRB 130427A, we find that current observations are insufficient to discriminate between PIC-motivated acceleration and the Bohm limit, primarily due to poor photon statistics in the Fermi-Large Area Telescope band. Nevertheless, future MeV–TeV afterglow observations can break model degeneracies and place substantially tighter constraints on the mechanisms responsible for particle acceleration in relativistic shocks. To this end, we simulate a fiducial nearby short GRB as a promising probe of the cutoff location, for which the two acceleration scenarios are cleanly distinguishable and the detection of such an event in the near future remains feasible.

Wu, Zhao-Feng [Purdue University, West Lafayette, ↗
DOE OSTI2022

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Observations of twinning microstructure in iron ramp-compressed through the α – ε phase transition

Single crystal iron was ramp-compressed along the [100] axis at 4 GPa/ns through the α–ϵ phase transformation at 13 GPa. Unlike shock experiments, in these ramp compression experiments the bcc lattice is observed to be nearly isotropically relaxed before the phase transition. There was a mixed α/ϵ phase region starting at ~13 GPa which was observed for 0.75 ns until the peak stress of 18 GPa was reached. In situ x-ray diffraction measurements show the formation of new hcp orientations not reported in shock or quasistatic compression experiments. The new hcp orientations appear to be caused by {10¯11} <10¯1 ¯2> → {¯1011} <¯101¯2> sequential twins that occur during the phase transition. This twinning mechanism relieves the shear strain caused by the bcc-hcp phase transition from the isotropically relaxed bcc phase while subject to uniaxial compression. Furthermore, these results show that in iron the induced microstructure through a phase transition and the phase transition mechanism depend on the loading history.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Learning thermodynamically constrained equations of state with uncertainty

Numerical simulations of high energy-density experiments require equation of state (EOS) models that relate a material’s thermodynamic state variables—specifically pressure, volume/density, energy, and temperature. EOS models are typically constructed using a semi-empirical parametric methodology, which assumes a physics-informed functional form with many tunable parameters calibrated using experimental/simulation data. Since there are inherent uncertainties in the calibration data (parametric uncertainty) and the assumed functional EOS form (model uncertainty), it is essential to perform uncertainty quantification (UQ) to improve confidence in EOS predictions. Model uncertainty is challenging for UQ studies since it requires exploring the space of all possible physically consistent functional forms. Thus, it is often neglected in favor of parametric uncertainty, which is easier to quantify without violating thermodynamic laws. This work presents a data-driven machine learning approach to constructing EOS models that naturally captures model uncertainty while satisfying the necessary thermodynamic consistency and stability constraints. We propose a novel framework based on physics-informed Gaussian process regression (GPR) that automatically captures total uncertainty in the EOS and can be jointly trained on both simulation and experimental data sources. A GPR model for the shock Hugoniot is derived, and its uncertainties are quantified using the proposed framework. We apply the proposed model to learn the EOS for the diamond solid state of carbon using both density functional theory data and experimental shock Hugoniot data to train the model and show that the prediction uncertainty is reduced by considering thermodynamic constraints.

Sharma, Himanshu (ORCID:000900050235718X)↗
DOE OSTI2024

QCD-gravity double-copy in the Regge regime: Shock wave propagators

In a previous paper [H. Raj and R. Venugopalan, ], we demonstrated a double-copy relation between inclusive gluon radiation in shock wave collisions of ultrarelativistic nuclei and inclusive graviton radiation in trans-Planckian gravitational shock wave collisions. We compute here the corresponding gravitational shock wave propagators in general relativity and demonstrate that they too obey a double-copy relation to gluon shock wave propagators computed previously. These results provide key input in a renormalization group approach toward computing the high frequency radiation spectrum in close black hole encounters. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multimessenger measurements of the static structure of shock-compressed liquid silicon at 100 GPa

The ionic structure of high-pressure, high-temperature fluids is a challenging theoretical problem with applications to planetary interiors and fusion capsules. Here we report a multimessenger platform using velocimetry and angularly and spectrally resolved x-ray scattering to measure the thermodynamic conditions and ion structure factor of materials at extreme pressures. We document the pressure, density, and temperature of shocked silicon near 100 GPa with uncertainties of 6%, 2%, and 20%, respectively. The measurements are sufficient to distinguish between and rule out some ion screening models. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Rapid Construction of Validated Chemistry Models for Advanced Biofuels (Co-Optima Final Report)

Computer methods were developed and demonstrated that make it possible to predict the combustion and pyrolysis chemistry of proposed biofuels, even before the biofuels are ever manufactured or experimentally tested. The methods are based on first-principles calculations of reaction rates and equilibria. A large number of complicated calculations are required, so we developed a computer workflow that automates many of the steps. With the chemical structure of the biofuel and the system conditions (temperature, pressure, compositions, etc.) as the input, the computer methods can construct a detailed kinetic model with refined parameters without intense human intervention. The automation of model construction makes it practical for a single engineer to make the predictions in a reasonable time period. The accuracy of the computer predictions was tested by measuring the time-profiles of various reaction intermediates (carbon monoxide, carbon dioxide, water, formaldehyde, methane, and ethylene) formed as the biofuels pyrolyze or ignite using sophisticated laser shock tube methods. Several biofuels, including methyl propyl ether, cyclopentanone, butyl acetate isomers, etc., were carefully studied to demonstrate the validity of the proposed computational and experimental methodology. This research project opens up the possibility for assessing some proposed biofuels on the computer without expensive experimental campaigns. The new computer methods are open source and well-documented, allowing a knowledgeable modeler to easily access the developed methods and implement them in their work routine. Moreover, they are also convenient and powerful, promising for developing the fuel chemistry models needed for use in engine simulators, facilitating co-optimization of new engine designs with new fuels.

09 BIOMASS FUELS↗
DOE OSTI2024

Early photometric and spectroscopic observations of the extraordinarily bright INTEGRAL-detected GRB 221009A

Context. GRB 221009A, initially detected as an X-ray transient by Swift, was later revealed to have triggered the Fermi satellite about an hour earlier, marking it as a post-peak observation of the event’s emission. This GRB distinguished itself as the brightest ever recorded, presenting an unparalleled opportunity to probe the complexities of GRB physics. The unprecedented brightness, however, challenged observation efforts, as it led to the saturation of several high-energy instruments.Aims. Our study seeks to investigate the nature of the INTEGRAL-detected GRB 221009A and elucidate the environmental conditions conducive to these exceptionally powerful bursts. Moreover, we aim to understand the fundamental physics illuminated by the detection of teraelectronvolt (TeV) photons emitted by GRB 221009A.Methods. We conducted detailed analyses of early photometric and spectroscopic observations that span from the Fermi trigger through to the initial days following the prompt emission phase in order to characterize GRB 221009A’s afterglow, and we complemented these analyses with a comparative study.Results. Our findings from analyzing INTEGRAL data confirm GRB 221009A as the most energetic event observed to date. Early optical observations during the prompt phase negate the presence of bright optical emissions with internal or external shock origins. Spectroscopic analyses enabled us to measure GRB 221009A’s distance and line-of-sight properties. The afterglow’s temporal and spectral analysis suggests prolonged activity of the central engine and a transition in the circumburst medium’s density. Finally, we discuss the implications for fundamental physics of detecting photons as energetic as 18 TeV from GRB 221009A.Conclusions. Early optical observations have proven invaluable for distinguishing between the potential origins of optical emissions in GRB 221009A, underscoring their utility in GRB physics studies. However, the rarity of such data underscores the need for dedicated telescopes capable of synchronous multiwavelength observations. Additionally, our analysis suggests that the host galaxies of TeV GRBs share commonalities with those of long and short GRBs. Expanding the sample of TeV GRBs could further solidify these findings.Key words: techniques: photometric / techniques: spectroscopic / gamma-ray burst: general / gamma-ray burst: individual: GRB 221009A

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗