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 1,351 records · Page 75
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal plasticity finite element simulation of lattice rotation and x-ray diffraction during laser shock compression of tantalum

Here we present a crystal plasticity model tailored for high-pressure, high-strain-rate conditions that uses a multiscale treatment of dislocation-based slip kinetics. We use this model to analyze the pronounced plasticity-induced lattice rotations observed in shock-compressed polycrystalline tantalum via in situ x-ray diffraction. By making direct comparisons between experimentally measured and simulated texture evolution, we can explain how the details of the underlying slip kinetics control the degree of lattice rotation that ensues. Specifically, we show that only the highly nonlinear kinetics caused by dislocation nucleation can explain the magnitude of the rotation observed under shock compression. We demonstrate a good fit between our crystal plasticity model and x-ray diffraction data and exploit the data to quantify the dislocation nucleation rates that are otherwise poorly constrained by experiment in the dynamic compression regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

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

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

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

Measuring and simulating ice–ablator mix in inertial confinement fusion

Fuel–ablator mix has been established as a major performance degrading effect in the burning plasma regime of recent inertial confinement fusion (ICF) experiments. As such, the study of fuel–ablator mix with experiments and simulations can provide valuable insight for our understanding of these experiments and establish a path for even higher yields and increased robustness. Here we present a novel high-yield experimental ICF design that is motivated by recent experiments measuring ice–ablator mix with a CH ablator instead of a high-density carbon (HDC) ablator. Here we review these experiments in more detail and describe the modeling assumptions and parameters used to obtain agreement with the data from implosion and burn simulations with mix. Using this mix model calibrated a posteriori to the experimental data, we design an implosion that uses a CH ablator that is predicted to achieve better performance than a recent experiment that achieved net target gain of 1.5 in HDC. Because hydrodynamic instabilities are greatly reduced with this new design, we also expect a high reproducibility at the same implosion adiabat as current record yield experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Improvised Electrsopinning Set Up for Thicker Ceramic Nanofiber Mat for High Power Targets

In high energy particle physics there is a demand for multi-MW high performance particle production targets. Nanofiber microstructure will have better performance than current solid targets in mitigating increased thermal tress waves, radiation damage. The objective is to fabricate thicker ceramic/metallic nano-fiber with high strength, thermal shock resistance using low cost electrospinning process.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

The gasdynamics of compact relativistic jets

The gasdynamics of compact relativistic jets is explored by analyzing a specific idealized flow problem using the method of characteristics. The basic flow pattern of the gas and pressure waves within a jet experiencing a drop in external pressure is calculated, with analytic expressions given for many of the important parameters. Scaling laws which relate the intrinsic properties of the jet to the pressure of the surrounding medium are obtained and discussed. The physical properties of the jet depend critically on the value and abruptness of the decrease in external pressure, as well as on the initial Lorentz factor of the flow. A variety of the flow patterns can result, including jets of oscillating cross section, jets with standing shocks, and broad, nearly hollow beams which can break up into multiple jets. These results are discussed in relation to the observed characteristics of superluminal radio sources in general and the superluminal quasar 4C 39.25 in particular.

Daly, Ruth A.↗
NASA NTRS1990

Inlets, ducts, and nozzles

The internal fluid mechanics research program in inlets, ducts, and nozzles consists of a balanced effort between the development of computational tools (both parabolized Navier-Stokes and full Navier-Stokes) and the conduct of experimental research. The experiments are designed to better understand the fluid flow physics, to develop new or improved flow models, and to provide benchmark quality data sets for validation of the computational methods. The inlet, duct, and nozzle research program is described according to three major classifications of flow phenomena: (1) highly 3-D flow fields; (2) shock-boundary-layer interactions; and (3) shear layer control. Specific examples of current and future elements of the research program are described for each of these phenomenon. In particular, the highly 3-D flow field phenomenon is highlighted by describing the computational and experimental research program in transition ducts having a round-to-rectangular area variation. In the case of shock-boundary-layer interactions, the specific details of research for normal shock-boundary-layer interactions are described. For shear layer control, research in vortex generators and the use of aerodynamic excitation for enhancement of the jet mixing process are described.

Abbott, John M.↗
NASA NTRS1997

The Ambient and Perturbed Solar Wind: From the Sun to 1 AU

The overall objective of the proposed research was to use numerical solutions of the magnetohydrodynamic (MHD) equations along with comparisons of the computed results with observations to study the following topics: (1) ambient solar wind solutions that extend from the solar surface to 1 astronomical unit (AU), contain closed magnetic structures near the Sun, and are consistent with observed values; (2) magnetic and plasma structures in coronal mass ejections (CMES) as they propagate to the interplanetary medium; (3) relation of MHD shocks to CMEs in the interplanetary medium; (4) interaction of MHD shocks with structures (such as other shocks, corotating interaction regions, current sheets) in the interplanetary plasma; and (5) simulations of observed interplanetary structures. A planned close collaboration with data analysts served to make the model more relevant to the data. The outcome of this research program is an improved understanding of the physical processes occurring in solar-generated disturbances in the interplanetary medium between the Sun and 1 AU.

Steinolfson, R. S.↗
DOE OSTI2024

Pressure-induced transformation of Nb 2 O 5 under shock compression from first principles

Ab initio molecular dynamics (AIMD) simulations were carried out to investigate the equation of state of Nb 2 O 5 and its pressure-density relationship under shock conditions. Here, the focus of this study is on the monoclinic B-Nb 2 O 5 (C2/c) polymorph. Enthalpy calculations from AIMD trajectories at 300 K show that the pressure-induced transformation between the thermody-namically most stable crystalline monoclinic parent phase H-Nb 2 O 5 (P2/m) and B-Nb 2 O 5 occurs at ~1.9 GPa. This H→B transition is energetically more favorable than the H→L(Pmm2) pressure-induced transition recently observed at ~5.9-9.0 GPa. The predicted shock properties of Nb 2 O 5 polymorphs are also compared to their Nb and NbO 2 counterparts to assess the impact of niobium oxidation on shock response.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

A study of velocity, temperature, and density in the plasma generated by laser-induced breakdowns

The paper presents velocity measurements of shock-induced flow field, leading to vortex generation and plasma deformation in air. Femtosecond-laser electronic excitation tagging (FLEET) velocimetry was performed at the times Δt = 3, 20, 50, and 100 µs post laser-induced breakdown. Emissions over Δt = 3–10 µs showed the propagations of the initially elliptic shock, transitioning into a spherical front. The shock emanated along the laser axis causes the flow outward, and then the pressure gradient generated by the rarefaction wave drives the inward flow at later moments, with the velocity magnitude approaching a steady-state value of 40 m s -1 . Temporal velocity evolution was compared with non-self-similar solutions behind the propagating shock, which are sensitive to the size of the energy deposition, and the use of the measured initial plasma diameter reproduced the experiment. There establishes a region of uniform velocity around 35 m s -1 in the air-flow running through the plasma, which triggers the roll-up of the plasma surface by a large-scale vortex, providing the detail of flow field evolution from the shock propagation to the plasma deformation. A collective Thomson scattering and hybrid fs/ps pure rotational coherent anti-Stokes Raman scattering (CARS) were also performed to gain insight into the high-temperature plasma. An effective electron–ion recombination rate of 2 x 10 -12 cm 3 s -1 was measured at Δt = 0.5–10 µs, during a dynamic plasma expansion and compression. When the shock resides close to the plasma at Δt = 0.5–1 µs, the temperature distributions were found to follow the similarity law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

Molecular Cooling as a Probe of Star Formation: Spitzer Looking Forward to Herschel

We explore here the question of how cloud physics can be more directly probed when one observes the majority of cooling emissions from molecular gas. For this purpose we use results from a recent Spitzer Space Telescope study of the young cluster of embedded objects in NGC1333. For this study we mapped the emission from eight pure H2 rotational lines, from S(0) to S(7). The H2 emission appears to be associated with the warm gas shocked by the multiple outflows present in the region. The H2 lines are found to contribute to 25 - 50% of the total outflow luminosity, and can be used to more directly ascertain the importance of star formation feedback on the natal cloud. From these lines, we determine the outflow mass loss rate and, indirectly, the stellar infall rate, the outflow momentum and the kinetic energy injected into the cloud over the embedded phase. The latter is found to exceed the binding energy of individual cores, suggesting that outflows could be the main mechanism for cores disruption. Given the recent launch of Herschel and the upcoming operational lifetime of SOFIA we discuss how studies of molecular cooling can take a step beyond understanding thermal balance to exploring the origin, receipt, and transfer of energy in atomic and molecular gas in a wide range of physical situations.

Bergin, Edwin A.↗
DOE OSTI2020

Energy Flow in Thin Shell Implosions and Explosions

Energy flow and balance in convergent systems beyond petapascal energy densities controls the fate of late-stage stars and the potential for controlling thermonuclear inertial fusion ignition. Time-resolved x-ray self-emission imaging combined with a Bayesian inference analysis is used to describe the energy flow and the potential information stored in the rebounding spherical shock at 0.22 PPa (2.2 Gbar or billions of atmospheres pressure). This analysis, together with a simple mechanical model, describes the trajectory of the shell and the time history of the pressure at the fuel-shell interface, ablation pressure, and energy partitioning including kinetic energy of the shell and internal energy of the fuel. Here, the techniques used here provide a fully self-consistent uncertainty analysis of integrated implosion data, a thermodynamic-path independent measurement of pressure in the petapascal range, and can be used to deduce the energy ow in a wide variety of implosion systems to petapascal energy densities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

New Constraints on the Melting Temperature and Phase Stability of Shocked Iron up to 270 GPa Probed by Ultrafast X-Ray Absorption Spectroscopy

Studying the properties and phase diagram of iron at high-pressure and high-temperature conditions has relevant implications for Earth’s inner structure and dynamics and the temperature of the inner core boundary (ICB) at 330 GPa. Also, a hexagonal-closed packed to body-centered cubic (bcc) phase transition has been predicted by many theoretical works but observed only in a few experiments. The recent coupling of high-power laser with advanced x-ray sources from synchrotrons allows for novel approaches to address these issues. Here, we present a study on shock compressed iron up to 270 GPa and 5800 K probed by single-pulse (100 ps FWHM) x-ray absorption spectroscopy (XAS). Based on the analysis of the XAS spectra, we provide structural identification and bulk temperature measurements along the Hugoniot up to the melting. These results rule out the predicted transition to a high-temperature bcc phase and allow one to discriminate among existing equations of state models and melting curves. In particular, we report the first bulk temperature measurement in shock compressed iron on the melting plateau at 240(20) GPa and 5345(600) K. The melting curve resulting from our work extrapolates to a temperature of 6202(514) K at 330 GPa and represents a refined upper bound for the ICB temperature. Published by the American Physical Society 2024

Balugani, S.↗
NASA NTRS1990

On diamond, graphitic and amorphous carbons in primitive extraterrestrial solar system materials

Carbon is among the most abundant elements in the universe and carbon chemistry in meteorites and comets is an important key to understanding many Solar System and interstellar processes. Yet, the mineralogical properties and interrelations between various structural forms of elemental carbon remain ambiguous. Crystalline elemental carbons include rhombohedral graphite, hexagonal graphite, cubic diamond, hexagonal diamond (i.e., lonsdaleite or carbon-2H) and chaoite. Elemental carbon also occurs as amorphous carbon and poorly graphitized (or turbostratic) carbon but of all the forms of elemental carbon only graphite is stable under physical conditions that prevail in small Solar System bodies and in the interstellar medium. The recent discovery of cubic diamond in carbonaceous chondrites and hexagonal diamond in chondritic interplanetary dust particles (IDPs) have created a renewed interest in the crystalline elemental carbons that were not formed by shock processes on a parent body. Another technique, Raman spectroscopy, confirms a widespread occurrence of disordered graphite in the Allende carbonaceous chondrite and in chondritic IDPs. Elemental carbons have also been identified by their characteristic K-edge features in electron energy loss spectra (EELS). However, the spectroscopic data do not necessarily coincide with those obtained by selected area electron diffraction (SAED). In order to interpret these data in terms of rational crystalline structures, it may be useful to consider the principles underlying electron diffraction and spectroscopic analyses. Electron diffraction depends on electron scattering, on the type of atom and the distance between atoms in a crystal lattice. Spectroscopic data are a function of the type of atom and the energy of bonds between atoms. Also, SAED is a bulk sampling technique when compared to techniques such as Raman spectroscopy or EELS. Thus, it appears that combined analyses provide contradictory results and that amorphous, or short-range ordered, carbon identified by conventional TEM imaging and SAED may show evidence for sp(3) bonds in EELS spectra. It is suggested that complex, nanometer-scale, mineralogical interrelations are common to all elemental carbons irrespective of their origin. The subsequent thermal history, or energy balance, will determine the ultimate microstructure.

Rietmeijer, Frans J. M.↗
NASA NTRS1989

CFD for applications to aircraft aeroelasticity

Strong interactions of structures and fluids are common in many engineering environments. Such interactions can give rise to physically important phenomena such as those occurring for aircraft due to aeroelasticity. Aeroelasticity can significantly influence the safe performance of aircraft. At present exact methods are available for making aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for complex flows containing shock waves, vortices and flow separations, computational methods are still under development. Several phenomena that can be dangerous and limit the performance of an aircraft occur due to the interaction of these complex flows with flexible aircraft components such as wings. For example, aircraft with highly swept wings experience vortex induced aeroelastic oscillations. Correct understanding of these complex aeroelastic phenomena requires direct coupling of fluids and structural equations. Here, a summary is presented of the development of such coupled methods and applications to aeroelasticity since about 1978 to present. The successful use of the transonic small perturbation theory (TSP) coupled with structures is discussed. This served as a major stepping stone for the current stage of aeroelasticity using computational fluid dynamics. The need for the use of more exact Euler/Navier-Stokes (ENS) equations for aeroelastic problems is explained. The current development of unsteady aerodynamic and aeroelastic procedures based on the ENS equations are discussed. Aeroelastic results computed using both TSP and ENS equations are discussed.

Guruswamy, Guru P.↗
DOE OSTI2025

A Curated Dataset of Regional Meteor Events with Simultaneous Optical and Infrasound Observations (2006–2011)

We present a curated, openly accessible dataset of 71 regional meteor events simultaneously recorded by optical and infrasound instrumentation between 2006 and 2011. These events were captured during an observational campaign using the all-sky cameras of the Southern Ontario Meteor Network and the co-located Elginfield Infrasound Array. Each entry provides optical trajectory measurements, infrasound waveforms, and atmospheric specification profiles. The integration of optical and acoustic data enables robust linkage between observed acoustic signals and specific points along meteor trajectories, offering new opportunities to examine shock wave generation, propagation, and energy deposition processes. This release fills a critical observational gap by providing the first validated, openly accessible archive of simultaneous optical–infrasound meteor observations that supports trajectory reconstruction, acoustic propagation modeling, and energy deposition analyses. By making these data openly available in a structured format, this work establishes a durable reference resource that advances reproducibility, fosters cross-disciplinary research, and underpins future developments in meteor physics, atmospheric acoustics, and planetary defense.

astrometry↗