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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗
NASA NTRS1993

Prediction and control of vortex-dominated and vortex-wake flows

This progress report documents the accomplishments achieved in the period from December 1, 1992 until November 30, 1993. These accomplishments include publications, national and international presentations, NASA presentations, and the research group supported under this grant. Topics covered by documents incorporated into this progress report include: active control of asymmetric conical flow using spinning and rotary oscillation; supersonic vortex breakdown over a delta wing in transonic flow; shock-vortex interaction over a 65-degree delta wing in transonic flow; three dimensional supersonic vortex breakdown; numerical simulation and physical aspects of supersonic vortex breakdown; and prediction of asymmetric vortical flows around slender bodies using Navier-Stokes equations.

Kandil, Osama↗
NASA NTRS1994

Development of processing techniques for advanced thermal protection materials

The effort, which was focused on the research and development of advanced materials for use in Thermal Protection Systems (TPS), has involved chemical and physical testing of refractory ceramic tiles, fabrics, threads and fibers. This testing has included determination of the optical properties, thermal shock resistance, high temperature dimensional stability, and tolerance to environmental stresses. Materials have also been tested in the Arc Jet 2 x 9 Turbulent Duct Facility (TDF), the 1 atmosphere Radiant Heat Cycler, and the Mini-Wind Tunnel Facility (MWTF). A significant part of the effort hitherto has gone towards modifying and upgrading the test facilities so that meaningful tests can be carried out. Another important effort during this period has been the creation of a materials database. Computer systems administration and support have also been provided. These are described in greater detail below.

Selvaduray, Guna S.↗
NASA NTRS1994

Bolt cutter functional evaluation

The Aerospace Corporation has been implementing finite difference and finite element codes for the analysis of a variety of explosive ordnance devices. Both MESA-2D and DYNA3D have been used to evaluate the role of several design parameters on the performance of a satellite separation system bolt cutter. Due to a lack of high strain rate response data for the materials involved, the properties for the bolt cutter and the bolt were selected to achieve agreement between computer simulation and observed characteristics of the recovered test hardware. The calculations provided insight into design parameters such as the cutter blade kinetic energy, the preload on the bolt, the relative position of the anvil, and the anvil shape. Modeling of the cutting process clarifies metallographic observation of both cut and uncut bolts obtained from several tests. Understanding the physical processes involved in bolt cutter operation may suggest certain design modifications that could improve performance margin without increasing environmental shock response levels.

S. Goldstein↗
NASA NTRS1996

X-ray emission from superbubbles

The results of Rosat observations of superbubbles in the Large Magellanic Cloud provided X-ray spectra for bright objects and data on upper limits of faint objects. The X-ray emission from superbubbles can be used to determine the physical condition of the superbubble interior, place constraints on the ambient interstellar medium, and diagnose the existence of supernova remnant shocks.

Chu, You-Hua↗
NASA NTRS2005

Constraints on Grain Formation Around Carbon Stars from Laboratory Studies of Presolar Graphite

We report the results of an investigation into the physical conditions in the mass outflows of asymptotic giant branch (AGB) carbon stars that are required for the formation of micron-sized presolar graphite grains, either with or without internal crystals of titanium carbide (TiC). In addition to providing detailed information about stellar nucleosynthesis, the structure and composition of presolar grains give unique information about the conditions of grain formation. In the present work we use laboratory observations of presolar graphite to gain insight into the physical conditions in circumstellar outflows from carbon AGB stars. The periodic pulsation of AGB stars enhances the gas density through shocks in the stellar atmosphere above the photosphere, promoting the condensation of dust grains. Copious mass outflow occurs largely because grains are coupled to the radiation field of the star, which accelerates them by radiation pressure; momentum is in turn transferred to gas molecules by collisions with grains. The dust/gas mixture is effectively a two-component fluid whose motion depends on atmospheric structure and which, in turn, influences that structure. In particular, the radiation pressure on the grains determines the velocity field of the outflow and thus the density distribution, while the density distribution itself determines the conditions of radiative transfer within the outflow and thus the effective radiation pressure.

Bernatowicz, T. J.↗
NASA NTRS2012

Collisional Histories of Comets and Trojan Asteroids: Insights from Forsterite and Enstatite Impact Studies

Impacts into forsterite and orthoenstatite at speeds typically encountered by comets demonstrate that shock imparted by collisions is detectable in the infrared signatures of their dust. The spectral signatures can be traced to physical alterations in their crystalline structures, as observed in TEM imaging and modeled using a dipole approximation. These results yield tantalizing insights into the collisional history of our solar system, as well as the history of individual comets and Trojan asteroids.

Lederer. S. M.↗
NASA NTRS2017

The Mars Crustal Magnetic Field Control of Plasma Boundary Locations and Atmospheric Loss: MHD Prediction and Comparison with MAVEN

We present results from a global Mars time-dependent MHD simulation under constant solar wind and solar radiation impact considering inherent magnetic field variations due to continuous planetary rotation. We calculate the 3-D shapes and locations of the bow shock (BS) and the induced magnetospheric boundary (IMB) and then examine their dynamic changes with time. We develop a physics-based, empirical algorithm to effectively summarize the multidimensional crustal field distribution. It is found that by organizing the model results using this new approach, the Mars crustal field shows a clear, significant influence on both the IMB and the BS. Specifically, quantitative relationships have been established between the field distribution and the mean boundary distances and the cross-section areas in the terminator plane for both of the boundaries. The model-predicted relationships are further verified by the observations from the NASA Mars Atmosphere and Volatile EvolutioN (MAVEN) mission. Our analysis shows that the boundaries are collectively affected by the global crustal field distribution, which, however, cannot be simply parameterized by a local parameter like the widely used subsolar longitude. Our calculations show that the variability of the intrinsic crustal field distribution in Mars-centered Solar Orbital itself may account for approx.60% of the variation in total atmospheric loss, when external drivers are static. It is found that the crustal field has not only a shielding effect for atmospheric loss but also an escape-fostering effect by positively affecting the transterminator ion flow cross-section area.

Global↗
DOE OSTI2022

Dynamic compression of [100] MgF 2 single crystals: Shock-induced polymorphism to highly coordinated structures

AX 2 -type compounds exhibit a rich and complex series of phase transitions under compression. The high-pressure polymorphism of rutile-type AX 2 difluorides can be characterized by a typical sequence from rutile type → Ca Cl 2 type → HP – Pd F 2 type → cotunnite type . Relative to dioxides, the reduced valence and ionic radius of the F – anion in rutile-type MgF 2 (sellaite) results in lowered transition pressures, making it a useful analog for SiO 2 and other dioxides. In this work, MgF 2 single crystals were shock-compressed along the [100] direction to 24 to 120 GPa using plate impact techniques, and wave profiles were measured using laser interferometry. At low stresses (24 to 44 GPa), we observe features consistent with an elastic–inelastic response, followed by a phase transformation. Peak stress-density states in this stress range are consistent with those expected for the modified fluorite-type (HP- PdF 2 ) phase or a related structure. At higher stresses (69 to 91 GPa), we observe a two-wave structure with peak stress-density states consistent with transformation to a denser and likely more highly coordinated phase such as the cotunnite-type structure. Furthermore, at 120 GPa, only a single wave structure is observed, indicating that the wave profile features observed at lower stresses are overdriven.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2021

Temperature and Density on the Forsterite Liquid‐Vapor Phase Boundary

Abstract The physical processes during planet formation span a large range of pressures and temperatures. Giant impacts, such as the one that formed the Moon, achieve peak pressures of 100s of GPa. The peak shock states generate sufficient entropy such that subsequent decompression to low pressures intersects the liquid‐vapor phase boundary. The entire shock‐and‐release thermodynamic path must be calculated accurately in order to predict the post‐impact structures of planetary bodies. Forsterite (Mg 2 SiO 4 ) is a commonly used mineral to represent the mantles of differentiated bodies in hydrocode models of planetary collisions. Here, we performed shock experiments on the Sandia Z Machine to obtain the density and temperature of the liquid branch of the liquid‐vapor phase boundary of forsterite. This work is combined with previous work constraining pressure, density, temperature, and entropy of the forsterite principal Hugoniot. We find that the vapor curves in previous forsterite equation of state models used in giant impacts vary substantially from our experimental results, and we compare our results to a recently updated equation of state. We have also found that due to under‐predicted entropy production on the principal Hugoniot and elevated temperatures of the liquid vapor phase boundary of these past models, past impact studies may have underestimated vapor production. Furthermore, our results provide experimental support to the idea that giant impacts can transform much of the mantles of rocky planets into supercritical fluids.

58 GEOSCIENCES↗
DOE OSTI2021

First-principles equation of state database for warm dense matter computation

We put together a first-principles equation of state (FPEOS) database for matter at extreme conditions by combining results from path integral Monte Carlo and density functional molecular dynamics simulations of the elements H, He, B, C, N, O, Ne, Na, Mg, Al, and Si as well as the compounds LiF , B 4 C , BN , CH 4 , CH 2 , C 2 H 3 , CH , C 2 H , MgO , and MgSiO 3 . For all these materials, we provide the pressure and internal energy over a density-temperature range from ~ 0.5 to 50 g cm - 3 and from ~ 10 4 to 10 9 K, which are based on ~ 5000 different first-principles simulations. We compute isobars, adiabats, and shock Hugoniot curves in the regime of L - and K -shell ionization. Invoking the linear mixing approximation, we study the properties of mixtures at high density and temperature. Furthermore, we derive the Hugoniot curves for water and alumina as well as for carbon-oxygen, helium-neon, and CH-silicon mixtures. We predict the maximal shock compression ratios of H 2 O , H 2 O 2 , Al 2 O 3 , CO , and CO 2 to be 4.61, 4.64, 4.64, 4.89, and 4.83, respectively. Finally we use the FPEOS database to determine the points of maximum shock compression for all available binary mixtures. We identify mixtures that reach higher shock compression ratios than their end members. We discuss trends common to all mixtures in pressure-temperature and particle-shock velocity spaces. In the Supplemental Material, we provide all FPEOS tables as well as computer codes for interpolation, Hugoniot calculations, and plots of various thermodynamic functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Reduced instability growth and improved radiation trapping with optimized shock timing in double-shell inertial confinement fusion capsules

The double shell is a volume-burn inertial confinement fusion concept consisting of two concentric shells: a low-Z outer shell that collides with and transfers momentum to a high-Z inner shell which compresses and heats the thermonuclear fuel. The increased number of capsule interfaces and severe hydrodynamic instability of the high-density pusher during its acceleration phase provide challenges to the success of the double shell. Two-dimensional radiation-hydrodynamics simulations predict the hydrodynamic instability growth on the outer surface of the pusher can be greatly reduced through appropriate timing of two shocks that cross this interface. One of these shocks, unique to multi-shell designs, arises from radiation-driven ablation of the inner shell ahead of the main shock, the second shock of concern. The shock timing is optimized by increasing the thickness of a low-Z tamper layer exterior to the pusher, resulting in only minimal changes to the implosion timing. Reducing the instability growth on the outer surface of the high-Z pusher can dramatically decrease the modulations that feedthrough to the pusher inner surface, improving the efficacy of radiation trapping in the thermonuclear fuel and increasing the predicted thermonuclear yield by ≳20×.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗
DOE OSTI2022

Phase transitions in Zr at sub-nanosecond time scales

Here solid-solid phase transitions are investigated in Zr thin films using shock compression induced by a short laser pulse (<1 ns). Shock wave profiles are measured at free surfaces for films of different thicknesses (a few micrometers) using chirped-pulse line velocimetry with 10 ps time resolution. Experiments are performed at pressures up to ~50 GPa, which is sufficient to reach the ω and β phases under equilibrium conditions. The shock wave structures are analyzed using a general Lagrangian analysis method, which allows for the calculation of stress–strain paths and assessments of phase transition behavior. In agreement with recent short laser pulse experiments using ultra-fast x-ray diffraction, we do not find any clear evidence of the α–ω transition, though this would be expected from the phase diagram. Instead, we infer a direct transformation to a metastable β-phase at lower shock pressures (<20 GPa) and equilibrium β at higher pressures. Through the velocimetry analysis, we find α–β transformation onset times of less than ~100 ps and completion times of less than ~200 ps.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Double sided measurements and associated error bars for planar samples or subassemblies

Many of the high energy density (HED) experiments use targets with planar components usually in a multilayer stack as the physics package which is at the heart of the experiment. The physics package contains the sample part under study as well as other components as shown in Fig. 1, that when exposed to the laser drive applies pressure onto or shock the main component to take it to the desired state in the phase diagram. These other components may include the ablator, pusher, x-ray pre-heat shield, window or other special components. The components are usually glued together. Currently, experiments on NIF and Omega in the Dynamic Materials campaign use several platforms such as Diffraction, EOS and Strength (Ref 1-3) which all use such planar stack physics packages.

42 ENGINEERING↗