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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigation of novel twinning in additive body centered cubic titanium (Fiscal Year 2024 Final Report)

This report summarizes the result of a one year seedling project to investigate unusual twinning behavior in shock loaded additive Ti5552. The twinning behavior only occurs when the β phase of Ti5553 is metastable, and it appears to be a type of double twin involving two different twin variants, first a {332}⟨113⟩ twin forms before being consumed by a specific {112}⟨111⟩ twin variant to create a 20 o . This behavior has only been detected during shock loading around incipient spall damage. The twinning is investigated by performing postmortem EBSD and PED analysis of gas-gun loaded specimens and preliminary molecular dynamics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface Properties of Double-Fin Generated Shock-Wave/Boundary-Layer Interactions

Shock-Wave/Boundary-Layer Interaction (SBLI) is ubiquitous in high-speed air vehicles. SBLI surface flowfield properties of crossing shocks induced by symmetric Double-Fins (DFs) at fin angles of 8° and 10° at Mach 2 are experimentally investigated; a 10° Single-Fin (SF) is also explored as a comparison. Surface visualization is used to capture distinct flowfield structures such as separation and upstream influence. Steady and unsteady Pressure-Sensitive Paints (PSPs) are used to obtain global surface pressure fields. Characteristics of DF-SBLI surface topology are discussed by comparing the corresponding SF-SBLI. Here, the effect of fin angle on DF-SBLI is characterized, and unsteady dynamics of the pressure field are examined. Spectral Proper Orthogonal Decomposition (SPOD) is performed on unsteady PSP, and traveling surface pressure waves and their frequency-dependent behavior are discussed. Dispersion relation of surface pressure waves is revealed by performing two-dimensional space–time Fourier transform on unsteady PSP, and the results are compared with the phase velocity distributions obtained from SPOD. Accompanying group velocity behavior is examined for distinct SBLI regions wherein wave propagation is found dependent on its frequency. Steady surface pressures along the centerline are compared with previous experimental data. Fin-tip spacing is found to be an essential parameter for scaling of upstream influence length.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2026

3D probabilistic fracture mechanics / computational fluid dynamics simulation of a reactor pressure vessel under transient conditions

Reactor pressure vessels (RPVs) are safety-critical light-water-reactor components that, under irradiation, experience long-term material degradation in the form of embrittlement. This can increase their susceptibility to fracture under thermal-shock conditions, which could occur during off-normal transients such as loss-of-coolant accidents (LOCAs). During a LOCA, the most severe conditions for the RPV occur when emergency core cooling water is injected through the cold legs into the water-and-steam-filled RPV. The rapid cooling of the downcomer and internal RPV surface causes decreased temperature and elevated thermally driven tensile stresses in the RPV wall. This, combined with long-term material embrittlement, may cause fracture initiation at pre-existing flaws, challenging the integrity of the RPV. Assessing RPV integrity during transients with a large spatial variation in the coolant temperature requires a modeling approach that considers the effects of spatially varying coolant temperature on the fracture probability of a population of flaws distributed throughout the RPV, accounting for spatially varying embrittlement. Here, the present study addresses this need by demonstrating first-of-its-kind coupling of a high-fidelity 3D computational fluid dynamics code with 3D probabilistic fracture mechanics This was accomplished using representative models of a pressurized-water reactor subjected to small- and medium-break LOCA conditions, both of which can result in large spatial temperature variations. While the observed impact of accounting for 3D effects was minimal under the small-break LOCA this study indicates a significant increase in the probability of fracture initiation under the medium-break LOCA when 3D effects are considered, relative to a spatially uniform cooling scenario.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Probing the lattice structure of dynamically compressed and released single crystal iron through the alpha to epsilon phase transition

In this work, experiments using broadband Laue x-ray diffraction (XRD) were used to examine the lattice structure of dynamically compressed [100]-oriented single crystal iron samples at the Dynamic Compression Sector at the Advanced Photon Source. These experiments used 1 μm thick iron single crystal samples sandwiched between a polyimide ablator and a polycarbonate window. A 100 J, 10 ns duration laser pulse incident on the polyimide ablator was used to shock compress the iron samples to initial stresses greater than 25 GPa, exceeding the ~13 GPa alpha (body-centered-cubic or bcc structure) to epsilon (hexagonal-close-packed or hcp structure) phase transition stress. XRD measurements were performed at various times relative to the shock wave entering the iron sample: early times, <~150 ps while the initial shock waves propagated through the iron; intermediate times, after the iron equilibrated with the ablator and window reaching a plateau stress state (12 or 17 GPa) lasting several nanoseconds; and late times, during uniaxial strain release. The early time measurements show that in <~150 ps, the high-pressure hcp phase is relaxed with a c/a ratio of 1.61, contrary to previous laser shock experiments where a c/a ratio of 1.7 was inferred. In the plateau stress state and partially released states, XRD measurements showed that the hcp structure retained a c/a ratio of 1.61 with no observable changes in the microstructure. Upon stress release at ~1 GPa/ns release rate, the reverse phase transition (hcp to bcc) to the original single crystal orientation (implying a transformation memory effect) was observed to reach completion somewhere between 13 and 11 GPa, indicating little stress hysteresis under rapid uniaxial strain release. A similar memory effect for the reverse hcp to bcc transformation has been previously observed under hydrostatic compression. However, the bcc/hcp orientation relationships differ somewhat between dynamic and static compression experiments, implying that the transformation pathway under uniaxial dynamic strain differs from the Burgers mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Novel Hot-Spot Ignition Designs for Inertial Confinement Fusion with Liquid-Deuterium-Tritium Spheres

A new class of ignition designs is proposed for inertial confinement fusion (ICF) experiments. These designs are based on the hot-spot-ignition approach, but instead of conventional targets that comprise of spherical shells with thin frozen deuterium{tritium (DT) layers, homogeneous-density (solid or liquid) DT spheres are used where the lower-density central region and higher-density shell are created dynamically by appropriately shaping the laser pulse. These offer several advantages, including simplicity in target production and lower sensitivity to both laser imprint and physics uncertainty in shock interaction with the ice-vapor interface. The design evolution starts by launching an ~1-Mbar shock into a homogeneous DT sphere. After bouncing from the center, the reflected shock reaches the outer surface of the sphere and the shocked material starts to expand outward until its pressure drops below the ablation pressure. At this point, an adjustment shock is launched inward by supporting ablation pressure. This shock compresses the ablator and fuel, forming a shell. The shell is then accelerated and compressed by appropriately shaping the drive laser pulse, similar to the conventional thin-shell, hot-spot designs. Here, this paper demonstrates the feasibility of the new concept using hydrodynamic simulations and discusses the advantages and disadvantages of the concept compared with more-traditional ICF designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2026

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2022

Extemporaneous Mechanochemistry: Shock-Wave-Induced Ultrafast Chemical Reactions Due to Intramolecular Strain Energy

We report regions of energy localization referred to as hotspots are known to govern shock initiation and the run-to-detonation in energetic materials. Mounting computational evidence points to accelerated chemistry in hotspots from large intramolecular strains induced via the interactions between the shock wave and microstructure. However, definite evidence mapping intramolecular strain to accelerated or altered chemical reactions has so far been elusive. From a large-scale reactive molecular dynamics simulation of the energetic material 1,3,5-triamino-2,4,6-trinitrobenzene, we map decomposition kinetics to molecular temperature and intramolecular strain energy prior to reaction. Both temperature and intramolecular strain are shown to accelerate chemical kinetics. A detailed analysis of the atomistic trajectory shows that intramolecular strain can induce a mechanochemical alteration of decomposition mechanisms. The results in this paper could inform continuum-level chemistry models to account for a wide range of mechanochemical effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural evolution of liquid silicates under conditions in Super-Earth interiors

Molten silicates at depth are crucial for planetary evolution, yet their local structure and physical properties under extreme conditions remain elusive due to experimental challenges. In this study, we utilize in situ X-ray diffraction (XRD) at the Matter in Extreme Conditions (MEC) end-station of the Linear Coherent Linac Source (LCLS) at SLAC National Accelerator Laboratory to investigate liquid silicates. Using an ultrabright X-ray source and a high-power optical laser, we probed the local atomic arrangement of shock-compressed liquid (Mg,Fe)SiO 3 with varying Fe content, at pressures from 81(9) to 385(40) GPa. We compared these findings to ab initio molecular dynamics simulations under similar conditions. Results indicate continuous densification of theO-O and Mg-Si networks beyond Earth’s interior pressure range, potentially altering melt properties at extreme conditions. This could have significant implications for early planetary evolution, leading to notable differences in differentiation processes between smaller rocky planets, such as Earth and Venus, and super-Earths, which are exoplanets withmasses nearly three times that of Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Structural dynamics modeling of spent nuclear fuel during hypothetical package drop events

The response of spent nuclear fuel (SNF) to hypothetical package drop events is of particular interest in the scope of spent fuel storage and transportation because of the mechanical shock encountered in such scenarios. Previous testing and modeling by the U.S. Department of Energy has demonstrated that the shock and vibration environment of normal shipping and handling conditions (excluding package drop events) is relatively benign and does not challenge the integrity of spent nuclear fuel. Cask drop events are worth considering because SNF packages are required to withstand free drops onto unyielding surfaces as part of their licensing basis. The acceleration experienced during drop events can be orders of magnitude higher, and thus more advanced models are needed to encompass potential nonlinear behavior of the fuel, such as spacer grid buckling and rod-to-rod impact. This work describes a number of finite element models developed to calculate the response of spent nuclear fuel to various hypothetical drop events that have been validated by package and fuel assembly drop tests conducted in the last decade. Sensitivity of the model response to factors such as package drop orientation, secondary impacts, and irradiated material properties as well as their potential impacts to fuel cladding integrity, was also investigated. Cask drops are not expected as a regular occurrence during SNF transportation, but this work helps raise the understanding of SNF mechanical loads to the point of consistency with the package design requirements.

Kadooka, Kevin↗
DOE OSTI2024

Impacts of aluminum-cerium alloy deployment on different cerium commodities

An imbalance problem exists in the rare earth element (REE) market with cerium and lanthanum being more abundant and less valuable than REEs used in permanent magnets such as neodymium, praseodymium, and dysprosium. Research has been conducted on new applications for abundant REEs to improve mining economics. However, potential market impacts of these applications on individual REE commodities are unknown. A system dynamics model was built to investigate the impacts of a new aluminum-cerium (Al-Ce) alloy that can theoretically utilize three separate cerium forms: cerium carbonate, cerium oxide, or cerium metal. Cerium carbonate currently possesses the lowest price of the three available forms making it initially the most favorable choice. However, model results showed that cerium oxide could provide the best economics for future Al-Ce alloy deployment because its larger market size buffers major shocks to price.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Spontaneous shock waves in pulse-stimulated flocks of Quincke rollers

Active matter demonstrates complex spatiotemporal self-organization not accessible at equilibrium and the emergence of collective behavior. Fluids comprised of microscopic Quincke rollers represent a popular realization of synthetic active matter. Temporal activity modulations, realized by modulated external electric fields, represent an effective tool to expand the variety of accessible dynamic states in active ensembles. Here, we report on the emergence of shockwave patterns composed of coherently moving particles energized by a pulsed electric field. The shockwaves emerge spontaneously and move faster than the average particle speed. Combining experiments, theory, and simulations, we demonstrate that the shockwaves originate from intermittent spontaneous vortex cores due to a vortex meandering instability. They occur when the rollers’ translational and rotational decoherence times, regulated by the electric pulse durations, become comparable. The phenomenon does not rely on the presence of confinement, and multiple shock waves continuously arise and vanish in the system.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Nonlinear modeling and performance analysis of cracked beam microgyroscopes

Cracks constitute a common structural defect in microelectromechanical systems that may arise during manufacturing, mechanical fatigue, or shock loading. Areas weakened by cracks increase the flexibility of the microstructure. Furthermore, depending on the crack severity, the increased flexibility can cause dramatic changes in the static and dynamic behaviors of electrically actuated MEMS. In this study, a numerical investigation of the depth and location of multiple edge cracks on the performance of beam microgyroscopes is conducted. Applying the Griffith strain energy release theorem for an edge crack, the vibrational characteristics of microgyroscopes with multiple cracks including the static deflection, natural frequencies, and mode shapes are studied analytically. Then, the differential quadrature method is employed to model the nonlinear dynamic behaviors of the cracked gyroscope system. Cracks initiated along the driving or sensing directions of the microgyroscope are found to have differing effects on its nonlinear dynamic response. This numerical study reveals that the onset of severe single cracks or a series of multiple cracks can significantly affect the sensitivity of the microgyroscope to base rotations and therefore leads to a degradation in the performance of the damaged MEMS sensor. This study also shows that a potential design approach for broadband microgyroscope systems manufactured with cracks can be proposed.

42 ENGINEERING↗
DOE OSTI2021

Multidimensional Numerical Modeling of Combustion Dynamics in a Non-Premixed Rotating Detonation Engine With Adaptive Mesh Refinement

In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.

42 ENGINEERING↗