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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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 OSTI2021

The reflection of a blast wave by a very intense explosion

We demonstrate that the geometric similarity of Taylor’s blast wave persists beyond reflection from an ideal surface. Upon impacting the surface, the spherical symmetry of the blast wave is lost but its cylindrical symmetry endures. As the flow acquires dependence on a second spatial dimension, an analytic solution of the Euler equations becomes elusive. However, the preservation of axisymmetry, geometric similarity and planar symmetry in the presence of a mirror-like surface causes all flow solutions to collapse when scaled by the height of burst (HOB) and the shock arrival time at the surface. The scaled blast volume for any yield, HOB and ambient air density follows a single universal trajectory for all scaled time, both before and after reflection.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1986

Supershells and propagating star formation

Correlated supernovae from an OB association can carve large cavities (greater than 100 pc) in the interstellar medium (ISM), and can punch holes completely through the disk of a spiral galaxy. Supernova remnant energy within such a cavity is thermalized before the shock reaches the supershell. Thus stellar wind theory may be used to model these superbubbles. We describe how the evolution of the superbubble depends on the density distribution of the galactic disk gas and the rate of supernovae in the OB association. At a radius of 100 to 300 pc, the supershell becomes gravitationally unstable, forming giant molecular clouds which are the sites for new star formation. This gravitational instability of the supershells provides a physical mechanism for propagating star formation and may account for the observation of bursts of star formation in galaxies.

Maclow, M. M.↗

NASA NTRS1988

The stellar wind velocity function for red supergiants determined in eclipsing binaries

The potential for direct measurement of the acceleration of stellar winds from the supergiant component of Zeta Aurigae-type binary stars is discussed. The aberration angle of the interaction shock cone centered on the hot star provides a measure of the velocity of the cool star wind at the orbit of the secondary. This is confirmed by direct observations of stellar wind (P Cygni) line profile variations. This velocity is generally smaller than the final (terminal) velocity of the wind, deduced from the P Cygni line profiles. The contrast between these results and previously published supergiant wind models is discussed. The implication on the physics of energy source dissipation predicted in the theoretical models is considered.

Ahmad, Imad A.↗

NASA NTRS1994

Hydrogen and related materials at high density: Physics, chemistry and planetary implications

Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.

Hemley, R. J.↗

DOE OSTI2026

A stress-based fracture model for reacting metal ejecta

The evolution of reacting metal ejecta continues to be a topic of interest at the forefront of metals in reactive and extreme environments. Ejecta are small particles formed when the surface of a metal undergoes Richtmyer–Meshkov instability from a strong shock. Experiments have shown that in the case where ejecta are in ambient conditions that induce a reaction, the ejecta behave irregularly. The ejecta temperature rises and then plateaus, and the acceleration profile shows unexpected jumps. These variations are assumed to be related to the exothermic heat release and particle mass loss caused by the reaction. To explain this phenomenon, efforts to model this in simulations have increased. While current models can capture many of these physical processes, they currently assign a constant reaction shell thickness with little physical reasoning. This work remedies this problem by assigning a dynamic physically informed shell thickness to the reacting particles, using solid analysis. The shell thickness of the particles impacts the rate of change of reacted mass in the system, as well as the rate at which the particles react. The model is based on a simple stress–strain relationship and gives a dynamic assignment for when the reacting particle should begin to fracture. We compare our model to the previous computational and simulation data to analyze the effects of different model parameters.

42 ENGINEERING↗

NASA NTRS1981

On the role of vibrational excitation in dissociative recombination

An improved physical model of dissociative recombination is presented and applied to experimental data on the temperature dependence of rate coefficients for the rare-gas and atmospheric-gas ions. It is shown that in the charge neutralisation of the rare-gas dimer ions, autoionisation plays an important role (at least in comparison with the atmospheric-gas ions) and contributes to the fast fall-off in the rate coefficient with vibrational excitation observed in shock tube studies. Numerical estimates of the observed fall-off in rate coefficient with increasing vibrational excitation are also presented.

Cunningham, A. J.↗

DOE OSTI2022

Material effects on dynamics in triple-nozzle gas-puff Z pinches

Here, the gas-puff Z pinch has a long history with myriad applications as an efficient neutron or x-ray source. Its simplicity as a load configuration makes it suitable for studying fundamental plasma physics phenomena such as instabilities and energy transport. For example, the implosion of cylindrical shells onto a fusion fuel are inherently susceptible to instability growth on their external surfaces; if such instabilities are unmitigated, then the consequences in terms of degraded performance can be substantial. Similarly, mitigating heat transport from a hot fuel to its colder surrounding container can make fusion conditions more easily achievable. Here we have conducted a systematic study of triple-nozzle (outer liner, inner liner, fuel) gas puffs using two-dimensional (2D) magnetohydrodynamic simulations to investigate the effect of load material on the relevant dynamics. Analogous to past studies on spherical blast waves and converging shock waves, a trend emerges linking increased radiative cooling, lower adiabatic index, and increased magneto-Rayleigh-Taylor instability growth. Notably, our results suggest that, for the present configuration, Ar radiates less than both Ne and Kr during the early stages of the implosion while mass is being swept up and perturbations begin to seed instability growth. Consequently, pinches with Ar on the outer surface exhibit more stable 2D behavior. Here we also present a parameter scan of thermonuclear neutron yield, Y, as a function of peak current, I pk and dopant concentration with Ne or Ar, depending on the inner liner material. Above 6 MA, our results suggest Y ∝ $I$$^{5}_{pk}$ and even substantial mixing (10% by volume) of Ne into the fuel does not drastically reduce yield, suggesting an Ar/Ne/fuel configuration may reliably achieve DD thermonuclear yields of 10 13 – 10 14 /cm in the 10–20 MA range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Automated discovery of a robust interatomic potential for aluminum

Machine learning, trained on quantum mechanics (QM) calculations, is a powerful tool for modeling potential energy surfaces. A critical factor is the quality and diversity of the training dataset. Here we present a highly automated approach to dataset construction and demonstrate the method by building a potential for elemental aluminum (ANI-Al). In our active learning scheme, the ML potential under development is used to drive non-equilibrium molecular dynamics simulations with time-varying applied temperatures. Whenever a configuration is reached for which the ML uncertainty is large, new QM data is collected. The ML model is periodically retrained on all available QM data. The final ANI-Al potential makes very accurate predictions of radial distribution function in melt, liquid-solid coexistence curve, and crystal properties such as defect energies and barriers. We perform a 1.3M atom shock simulation and show that ANI-Al force predictions shine in their agreement with new reference DFT calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Simulated refraction-enhanced X-ray radiography of laser-driven shocks

Refraction-enhanced x-ray radiography (REXR) is used to infer shock-wave positions of more than one shock wave, launched by a multiple-picket pulse in a planar plastic foil. This includes locating shock waves before the shocks merge, during the early time and the main drive of the laser pulse that is not possible with the velocity interferometer system for any reflector. Simulations presented in this paper of REXR show that it is necessary to incorporate refraction and attenuation of x rays along with the appropriate opacity and refractive-index tables to interpret experimental images. Simulated REXR shows good agreement with an experiment done on the OMEGA laser facility to image a shock wave. REXR can be applied to design multiple-picket pulses with a better understanding of the shock locations. Finally, this will be beneficial to obtain the required adiabats for inertial confinement fusion implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Laser-plasma instabilities in long scale-length plasmas relevant to shock-ignition

We present 2D particle-in-cell simulations of laser-plasma instabilities (LPI) performed using conditions relevant to experiments on the OMEGA and NIF laser facilities. The laser intensity used is 2×10 15 W cm -2 , which allows comparison with previous experiments and simulations investigating shock ignition, though this is a relatively low intensity for the scheme. We find that the large convective gain of the instabilities leads to dynamics that are largely controlled by pump depletion and to scattering or absorption by LPIs predominantly occurring away from the quarter-critical density. The two-plasmon decay (TPD) instability is dominant for OMEGA-type conditions and remains important at the NIF-scale. In both cases, most absorption from TPD occurs near its Landau cutoff density, and we examine its nonlinear dynamics in this region, focusing on the ion-acoustic waves responsible for saturation. For NIF-type conditions, stimulated Raman scattering (SRS) and stimulated Brillouin scattering divert significant fractions of laser energy, and we find that SRS backscatter can occur with high gain at low densities due to kinetic effects. An in-depth analysis of hot-electron production is performed to identify the main sources and their characteristic electron temperature. For both simulations presented, the overall hot-electron temperature lies between 30 and 35 keV. This is relatively low and will likely increase shock pressure; however, the distributions also contain a significant number of high-energy (E k >100 keV) electrons that would likely cause unacceptable preheat. We, therefore, suggest some strategies that may be utilized to minimize this high-energy component. I. INTROD

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

GnarlyX: Eulerian multi-material hydrodynamics coupled to equation of state and hyperelastic, plastic constitutive models

GnarlyX is a new hydrocode for direct numerical simulations of the microstructural behavior of high explosives at the mesoscale. We summarize the computational framework for multidimensional, Eulerian multi-material hydrodynamics coupled to EOS and hyperelastic, plastic constitutive models. We present 1D verification tests of multi-material only and combined multimaterial and strength capabilities with comparisons to exact solutions of shock states resulting from an incident shock impacting the material interface of PBX9502 and aluminum. We show that GnarlyX performs well in capturing the resulting shock waves in examining numerical convergence with exact solutions. In later work, we will summarize the thermomechanics and multi-dimensional, parallel computing capabilities in GnarlyX with multi-dimensional verification tests.

42 ENGINEERING↗

NASA NTRS1986

The near-earth cross-tail current sheet - Detailed ISEE 1 and 2 case studies

Three near geomagnetic tail current sheet crossings of the ISEE 1 and 2 satellites, on April 5, 1979, are examined in detail. All are assoicated with the passage of an interplanetary shock and the region of variable solar wind pressure behind it. The general geometry of field reversing current sheets is discussed, and this geometry is examined for the cases studied, by using the ISEE 1 and 2 coorbiting satellite data sets. A new technique is employed which removes the effects of a variable sheet normal velocity for the first time. This allows us to calculate firm upper bounds on the current sheet thicknesses, and by utilizing certain physically motivated assumptions, determine the most probable actual sheet thicknesses, and inclinations of the field lines within these sheets. Current density profiles derived with this technique show the main cross-tail current sheet to be a structure that is many thermal ion-gyroradii thick and which is sometimes imbedded in a region that is three or more times thicker and contains much smaller current densities. These profiles also exhibit a considerable amount of fine structure in the sheet which appears as narrow peaks in the current density distributions. Possible explanations for these structures, and for the overall sheet structure itself, are examined.

Mccomas, D. J.↗

NASA NTRS1989

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗

NASA NTRS1993

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗

NASA NTRS1994

Lunar and martian meteorite delivery services

Launch mechanisms for lunar and martian meteorites have been investigated, by integrating physical modeling constraints, geochemical cosmic-ray exposure (CRE) constraints, and petrologic constraints. The potential source region for lunar meteorites is remarkably small compared to the final crater volume. CRE constraints indicate that most launches start at depths of less than or equal to 3.2 m, and cratering theory implies derivation of suitably accelerated objects from a subvolume with diameter only about 0.3 x the final crater diameter. The shallow depth provenance is probably related to shock-wave interference, enhanced by the lunar regolith's extremely low compressional wave velocity. CRE constraints alone imply that four to five separate launch events are represented among the eight well-studied lunar meteorites. Most of the lunar meteorites are regolith breccias, which tend to show only limited compositional diversity within any kilometer-scale region of the Moon. Several others are polymict breccias, which also show relatively subdued compositional diversity, compared to igneous rocks. The observed diversity among these samples in terms of abundances of mare basalt and KREEP, and in Mg/(Mg + Fe) ratio, implies that among eight well-studied lunar meteorites only two potential source craters pairings are plausible: between Asuka-881757 + Y-793169 (most probable) and between Y-793274 + EET875721. Altogether, these eight lunar meteorites apparently represent at least six separate source craters, including three in the past 10(exp 5) years and five in the past 10(exp 6) years. CRE constraints imply that SNC meteorites are launched from systematically greater than lunar meteorites. SNCs are also systematically bigger, and all nine well-studied SNCs are uncommonly young (by martian standards) mafic igneous rocks. Comparison between Viking and Apollo results reveals that rocks the size of common meteorites are remarkably scarce in the martian regolith, probably due to pervasive weathering. A plausible explanation for these trends is that most old, small, and shallow rocks on Mars have been weakened by a two-stage process of brecciation followed by pervasive weathering, to the point where they seldom survive the stresses of spallation off the planet. The scarcity of source-crater pairing among the lunar meteorites implies that these objects can be launched from craters much smaller than previously estimated and tends to support the suggestion of Rabinowitz (1993) that present-day cratering rates for the Earth-Moon region may be higher than previously estimated.

Warren, Paul H.↗