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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Calibration and characterization of the line-VISAR diagnostic at the HED-HIBEF instrument at the European XFEL

In dynamic-compression experiments, the line-imaging Velocity Interferometer System for Any Reflector (VISAR) is a well-established diagnostic used to probe the velocity history, including wave profiles derived from dynamically compressed interfaces and wavefronts, depending on material optical properties. Knowledge of the velocity history allows for the determination of the pressure achieved during compression. Such a VISAR analysis is often based on Fourier transform techniques and assumes that the recorded interferograms are free from image distortions. In this paper, we describe the VISAR diagnostic installed at the HED-HIBEF instrument located at the European XFEL along with its calibration and characterization. It comprises a two-color (532, 1064 nm), three-arm (with three velocity sensitivities) line imaging system. We provide a procedure to correct VISAR images for geometric distortions and evaluate the performance of the system using Fourier analysis. We finally discuss the spatial and temporal calibrations of the diagnostic. As an example, we compare the pressure extracted from the VISAR analysis of shock-compressed polyimide and silicon.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Compression-induced solidification of shock-melted cerium

Compression induced solidification has been observed in cerium on nanosecond timescales. For this study, a series of experiments were conducted on Sandia National Laboratories' Z Facility in which cerium was shock melted and subsequently shocklessly, or ramp, loaded across the melt line inducing solidification. The signature of solidification manifested in the recovery of material strength and the propagation of waves at the local elastic sound velocity. Density functional theory simulations of cerium along the experimental phase-space path exhibit spontaneous freezing to a tetragonal phase at the same pressure and closely predict the observed physical properties of solid and liquid cerium near melt.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

Here, we develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

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

Magnetopause Reconnection and Indents Induced by Foreshock Turbulence

Abstract Based on global hybrid simulation results, we predict that foreshock turbulence can reach the magnetopause and lead to reconnection as well as Earth‐sized indents. Both the interplanetary magnetic field (IMF) and solar wind are constant in our simulation, and hence, all dynamics are generated by foreshock instabilities. The IMF in the simulation is mostly Sun‐Earth aligned with a weak northward and zero dawn‐dusk component, such that subsolar magnetopause reconnection is not expected without foreshock turbulence modifying the magnetosheath fields. We show a reconnection example to illustrate that the turbulence can create large magnetic shear angles across the magnetopause to induce local bursty reconnection. Magnetopause reconnection and indents developed from the impact of foreshock turbulence can potentially contribute to dayside loss of planetary plasmas.

58 GEOSCIENCES↗
DOE OSTI2020

Dependence of Enstrophy Transport and Mixed Mass on Dimensionality and Initial Conditions in the Richtmyer–Meshkov Instability Induced Flows

Here we present a comparative study of the enstrophy budget and mixed mass between two- and three-dimensional flows induced by Richtmyer–Meshkov instability (RMI). Specifically, the individual contributions to the enstrophy budget due to the production from baroclinicity and from vortex stretching (which vanishes in two-dimensional (2D) flow) are delineated. This is enabled by a set of two- and three-dimensional computations at Atwood 0.5 having both narrow- and broad-band perturbations. A further three-dimensional (3D) computation is conducted at Atwood 0.9 using an identical narrowband perturbation to the Atwood 0.5 case to examine the sensitivity to density ratio. The mixed mass is also considered with the goal to obtain insight on how faithfully a simplified calculation performed in two dimensions can capture the mixed mass for an inertial confinement fusion (ICF) or other practical application. It is shown that the late time power law decay of variable density enstrophy is substantially different in two and three dimensions for the narrowband initial perturbation. The baroclinic production term is negligible in three dimensions (aside from the initial shock interaction), as vortex stretching is larger by two orders of magnitude. The lack of vortex stretching considerably reduces the decay rate in both narrowband and broadband perturbations in two dimensions. In terms of mixed mass, the lack of vortex stretching reduces the mixed mass in two dimensions compared to three in all cases. In the broadband cases, the spectral bandwidth in the 2D case is wider; hence, there is a longer time period of sustained linear growth which reduces the normalized mixed mass further.

42 ENGINEERING↗
DOE OSTI2020

Future impacts of ozone driven damages on agricultural systems

Current ozone (O 3 ) concentration levels entail significant damages in crop yields around the world. The reaction of the emitted precursors (mostly methane and nitrogen oxides) with solar radiation contribute to O 3 levels that exceed established thresholds for crop damage. This paper shows current and projected (up to 2080) relative yield losses (RYLs) driven by O 3 exposure for different crops and the associated economic damages applying dynamic crop production and prices that are calculated per region and period. We adjust future crop yields in the Global Change Assessment Model (GCAM) to reflect the RYLs and analyze the effects on agricultural markets. We find that the changes (generally reductions) in O 3 precursor emissions in a reference scenario would reduce the agricultural damages, compared to present, for most of the regions, with a few exceptions including India, where higher future O3 concentrations have large negative impacts on crop yields. The annual economic impact of O3 driven losses from 2010-2080 are, in billion US dollars at 2015 prices ($B), 5.0-6.0, 9.8-18.8, 6.7-10.6 and 10.4-12.5 for corn, soybeans, rice and wheat, respectively, with the large losses for wheat and soybeans driven by their comparatively high responses to O 3 . When O 3 effects are explicitly modelled as exogenous yield shocks in future periods, there is a direct impact in future agricultural markets. Therefore, the aggregated net present value (NPV) of crop production would be reduced around by $90.8B at a global level. However, these changes are not distributed evenly across regions, and the net present market value of the crops would increase by up to $118.2B (India) or decrease by up to $59.2B (China).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Solid-to-plasma transition of polystyrene induced by a nanosecond laser pulse within the context of inertial confinement fusion

Laser Direct Drive (LDD) Inertial Confinement Fusion (ICF) involves irradiating a spherical target of thermonuclear fuel coated with an ablator, usually made of polystyrene. Laser energy absorption near the target surface leads to matter ablation, hydrodynamic shocks, and ultimately capsule implosion. The conservation of spherical symmetry is crucial for implosion efficiency, yet spatial modulations in laser intensity can induce nonuniformities, causing the laser imprint phenomenon. Understanding laser imprint, especially considering the initial solid state, is essential for advancing LDD ICF. A first microscopic model of solid-to-plasma transition was built in 2019, accounting for laser absorption in the solid state with a band structure based ionization model. This model has been improved to include chemical fragmentation and a more accurate description of electron collision frequency in various matter states. The latest development involves assessing the model reliability by comparing theoretical predictions with experimental observations. Despite the success of this approach, questions remain, leading to further investigations and observations under different irradiation conditions. This work presents an experiment with a nanosecond pulse, taking into account hydrodynamic effects, and measures transmission dynamics over the entire laser beam area to observe two-dimensional effects. As a result, the objective is to adapt the theoretical model, couple it with a hydrodynamic code, and observe additional effects related to the initial solid state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Laser-driven flyer plate impact: Computational studies guided by experiments

We present a computational approach using a multimaterial, arbitrary Lagrangian–Eulerian code termed ALE3D to model the nanosecond/micrometer dynamics of the launch of 0.5–4.5 km/s laser-driven metal flyer plates and the impact with stationary targets of Pyrex and fused silica glasses, and Lexan and Plexiglas polymers, producing pressures in the target in the range of 5–20 GPa. The simulations are compared to experimental results where the flyer velocity profile and the velocity profile imparted to the target material were measured with high-speed velocimetry. The experimental flyer launch by a high-intensity pulsed laser is modeled by depositing heat into a thin vaporizable layer under the flyer plate. This model produces a flyer plate that has not been exposed to the laser pulse, allowing us to compare the properties of the real flyer to a simulated ideal flyer. The simulations of target impact are generally in good agreement with the experiment except at the highest impact velocities where the shock release process in the experiment is slower than that in the simulation. The cause of this disagreement is attributed to an inadequate description of the shock viscosity during the nanosecond unloading process.

Stekovic, Svjetlana (ORCID:0000000163878070)↗
DOE OSTI2019

A computational model for nanosecond pulse laser-plasma interactions

A multi-physics numerical model for laser-induced optical breakdown and laser-plasma interaction in a non-equilibrium gas is presented, accounting for: production of priming electrons via multi-photon ionization, energy absorption, cascade ionization, induced hydrodynamic response, and shock formation and propagation. The gas is governed by the Navier-Stokes equations, with non-equilibrium effects taken into account by means of a two-temperature model. The space-time dependence of the laser beam is modeled with a flux-tube formulation for the Radiative Transfer Equation. The flow governing equations are discretized in space using a second-order finite volume method. The semi-discrete equations are marched in time using an implicit-explicit (IMEX) dual time-stepping strategy, where diffusion and chemistry are solved implicitly, whereas convection is explicit. Application to a 20 ns long 50 mJ pulse laser-induced breakdown in quiescent O 2 shows the advantages of this temporal discretization during and just after the laser pulse, while a less-expensive symmetric Strang splitting (with implicit chemistry) is sufficient for the post-breakdown gas dynamics after ≃ 0.1 μs. The integrated model is shown to reproduce key features of corresponding experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Uncertainty Analysis in Multi‐Sector Systems: Considerations for Risk Analysis, Projection, and Planning for Complex Systems

Abstract Simulation models of multi‐sector systems are increasingly used to understand societal resilience to climate and economic shocks and change. However, multi‐sector systems are also subject to numerous uncertainties that prevent the direct application of simulation models for prediction and planning, particularly when extrapolating past behavior to a nonstationary future. Recent studies have developed a combination of methods to characterize, attribute, and quantify these uncertainties for both single‐ and multi‐sector systems. Here, we review challenges and complications to the idealized goal of fully quantifying all uncertainties in a multi‐sector model and their interactions with policy design as they emerge at different stages of analysis: (a) inference and model calibration; (b) projecting future outcomes; and (c) scenario discovery and identification of risk regimes. We also identify potential methods and research opportunities to help navigate the tradeoffs inherent in uncertainty analyses for complex systems. During this discussion, we provide a classification of uncertainty types and discuss model coupling frameworks to support interdisciplinary collaboration on multi‐sector dynamics (MSD) research. Finally, we conclude with recommendations for best practices to ensure that MSD research can be properly contextualized with respect to the underlying uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Post-explosion Evolution of Core-collapse Supernovae

Here we investigate the post-explosion phase in core-collapse supernovae with 2D hydrodynamical simulations and a simple neutrino treatment. The latter allows us to perform 46 simulations and follow the evolution of the 32 explosion models during several seconds. We present a broad study based on three progenitors (11.2, 15, and 27 M ⊙ ), different neutrino heating efficiencies, and various rotation rates. We show that the first seconds after shock revival determine the final explosion energy, remnant mass, and properties of ejected matter. Our results suggest that a continued mass accretion increases the explosion energy even at late times. We link the late-time mass accretion to initial conditions such as rotation strength and shock deformation at explosion time. Only some of our simulations develop a neutrino-driven wind (NDW) that survives for several seconds. This indicates that NDWs are not a standard feature expected after every successful explosion. Even if our neutrino treatment is simple, we estimate the nucleosynthesis of the exploding models for the 15 M ⊙ progenitor after correcting the neutrino energies and luminosities to get a more realistic electron fraction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Soliton Generation in Negative Thermal Expansion Materials

Strain solitons have been observed statically in several 2D materials and dynamically in substrate materials using ultrafast laser pulses. The latter case relies on lattice relaxation in response to ultrafast heating in a light-absorbing transducer material, a process which is sensitive to the thermal expansion coefficient. Here we consider an unusual case where the sign of the thermal expansion coefficient is negative, a scenario which is experimentally feasible in light of rapid and recent advances in the discovery of negative thermal expansion materials. We present numerical solutions to a nonlinear differential equation which has been repeatedly demonstrated to quantitatively model experimental data and discuss the salient results using realistic parameters for material linear and nonlinear elasticity. The solitons that emerge from the initial value problem with negative and positive thermal expansion are qualitatively different in several ways. The new case of negative thermal expansion gives rise to a nearly-periodic soliton train with chirped profile and free of an isolated shock front. We suggest this unanticipated result may be realized experimentally and assess the potential for certain applications of this generic effect.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Femtosecond temperature measurements of laser-shocked copper deduced from the intensity of the x-ray thermal diffuse scattering

We present 50-fs, single-shot measurements of the x-ray thermal diffuse scattering (TDS) from copper foils that have been shocked via nanosecond laser ablation up to pressures above ∼135 GPa. We hence deduce the x-ray Debye–Waller factor, providing a temperature measurement. The targets were laser-shocked with the DiPOLE 100-X laser at the High Energy Density endstation of the European X-ray Free-Electron Laser. Single x-ray pulses, with a photon energy of 18 keV, were scattered from the samples and recorded on Varex detectors. Despite the targets being highly textured (as evinced by large variations in the elastic scattering) and with such texture changing upon compression, the absolute intensity of the azimuthally averaged inelastic TDS between the Bragg peaks is largely insensitive to these changes, and allowing for both Compton scattering and the low-level scattering from a sacrificial ablator layer provides a reliable measurement of $T/Θ^2_D$, where Θ D is the Debye temperature. We compare our results with the predictions of the SESAME 3336 and LEOS 290 equations of state for copper and find good agreement within experimental errors. We, thus, demonstrate that single-shot temperature measurements of dynamically compressed materials can be made via thermal diffuse scattering of XFEL radiation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Broadband Rotational Spectroscopy in Uniform Supersonic Flows: Chirped Pulse/Uniform Flow for Reaction Dynamics and Low Temperature Kinetics

ConspectusThe study of gas-phase chemical reactions at very low temperatures first became possible with the development and implementation of the CRESU (French acronym for Reaction Kinetics in Uniform Supersonic Flows) technique. CRESU relies on a uniform supersonic flow produced by expansion of a gas through a Laval (convergent-divergent) nozzle to produce a wall-less reactor at temperatures from 10 to 200 K and densities of 1016-1018 cm-3 for the study of low temperature kinetics, with particular application to astrochemistry. In recent years, we have combined uniform flows with revolutionary advances in broadband rotational spectroscopy to yield an instrument that affords near-universal detection for novel applications in photodissociation, reaction dynamics, and kinetics. This combination of uniform supersonic flows with chirped-pulse Fourier-transform microwave spectroscopy (Chirped-Pulse/Uniform Flow, CPUF) permits detection of any species with a modest dipole moment, thermalized to the uniform temperature of the gas flow, with isomer, conformer, and vibrational state specificity. In addition, the use of broadband, high-resolution, and time-dependent (microsecond time scale) micro- and mm-wave spectroscopy makes it an ideal tool for characterizing both transient and stable molecules, as well as studying their spectroscopy and dynamics.In this Account, we review recent advances made using the CPUF technique, including studies of photodissociation, radical-radical reaction dynamics, and low temperature kinetics. These studies highlight both the strength of universal and multiplexed detection and the challenges of coupling it to a high-density collisional environment. Product branching and product evolution as a function of time have been measured for astrochemically relevant systems, relying on the detailed characterization of these flow conditions via experiments and fluid dynamics simulations. In the photodissociation of isoxazole, an unusual heterocyclic molecule with a very low-energy conical intersection, we have identified 7 products in 5 reaction channels and determined the product branching, pointing to both direct and indirect pathways. We have also approached the same system from separated NO and C3H3 reactants to explore a broader range of the potential energy surface, demonstrating the power of multichannel branching measurements for complex radical-radical reactions. We determined the product branching in the C3H2 isomers in the photodissociation of the propargyl radical and identified the importance of a hydrogen atom catalyzed isomerization to the lowest energy cyclic form. This then motivated a study of direct D-H exchange reaction in radicals, in which we demonstrate that it is an important and overlooked pathway for deuterium fractionation in astrochemical environments. Recently, we have shown the measurement of low temperature kinetics inside an extended Laval nozzle, after which a shock-free secondary expansion to low temperature and density affords an ideal environment for detection by rotational spectroscopy. These results highlight the power and potential of the CPUF approach, and future prospects will also be discussed in light of these developments.

Dias, Nureshan↗
DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗