Search NASA⌕ Search

SEARCH · Search NASA

Results for “shock jumps”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quantum-cascade-laser-absorption-spectroscopy diagnostic for temperature, pressure, and NO X 2 Π 1/2 at 500 kHz in shock-heated air at elevated pressures

The design, validation, and application of a quantum-cascade-laser-absorption-spectroscopy diagnostic for measuring gas temperature, pressure, and nitric oxide (NO) in high-temperature air are presented. A distributed-feedback quantum-cascade laser (QCL) centered near 1976 c m − <#comment/> 1 was used to scan across two transitions of NO in its ground electronic state ( X 2 Π <#comment/> 1 / 2 ). A measurement rate of 500 kHz was achieved using a single QCL by: (1) performing current modulation through a bias-tee, and (2) targeting closely spaced transitions with a large difference in lower-state energy. The diagnostic was validated in a mixture of 95% argon and 5% NO, which was shock-heated to ≈ <#comment/> 2000 to 3700 K. The average mean percent differences between laser-absorption-spectroscopy (LAS) measurements and predictions from shock-jump relations for temperature, pressure, and NO mole fraction were 3.1%, 4.1%, and 6.5%, respectively. The diagnostic was then applied to characterize shock-heated air at high temperatures (up to ≈ <#comment/> 5500 K ) and high pressures (up to 12 atm) behind either incident or reflected shocks. The LAS measurements were compared to theoretical predictions from shock-jump relations, pressure sensors mounted in the wall of the shock tube, and equilibrium values of the NO mole fraction. The average mean percent differences between LAS measurements and their aforementioned reference values were 3.2%, 10.8%, and 10.4% for temperature, pressure, and NO mole fraction, respectively. Last, a comparison between a measured NO mole fraction time history and a time-stepped homogeneous reactor simulation performed using two different chemical kinetics mechanisms is presented.

Gilvey, Jonathan J. (ORCID:0000000317664383)↗
DOE OSTI2021

The antesonic condition for the explosion of core-collapse supernovae – II. Rotation and turbulence

In the problem of steady free fall on to a standing shockwave around a central mass, the ‘antesonic’ condition limits the regime of stable accretion to $c_T^2/v_\mathrm{esc}^2\le 3/16$, where cT is the isothermal sound speed in the subsonic post-shock flow, and vesc is the escape velocity at the shock radius. Above this limit, it is impossible to satisfy both the Euler equation and the shock jump conditions, and the system transitions to a wind. This physics explains the existence of a critical neutrino luminosity in steady-state models of accretion in the context of core-collapse supernovae. Here, we extend the antesonic condition to flows with rotation and turbulence using a simple one-dimensional formalism. Both effects decrease the critical post-shock sound speed required for explosion. While quite rapid rotation is required for a significant change to the critical condition, we show that the level of turbulence typically achieved in supernova simulations can greatly impact the critical value of $c_T^2/v_\mathrm{esc}^2$. A core angular velocity corresponding to a millisecond rotation period after contraction of the proto-neutron star results in only a ~5 percent reduction of the critical curve. In contrast, near-sonic turbulence with specific turbulent kinetic energy $K/c_T^2=0.5-1$, leads to a decrease in the critical value of $c_T^2/v_{\rm esc}^2$ by ~20 to 40 percent. This analysis provides a framework for understanding the role of post-shock turbulence in instigating explosions in models that would otherwise fail and helps explain why multidimensional simulations explode more easily than their one-dimensional counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic simulations verifying reconnection rates measured in the laboratory, spanning the ion-coupled to near electron-only regimes

The rate of reconnection characterizes how quickly flux and mass can move into and out of the reconnection region. In the Terrestrial Reconnection EXperiment (TREX), the rate at which the antiparallel asymmetric reconnection occurs is modulated by the presence of a shock and a region of flux pileup in the high-density inflow. Simulations utilizing a generalized Harris-sheet geometry have tentatively shown agreement with TREX's measured reconnection rate scaling relative to system size, which is indicative of the transition from ion-coupled toward electron-only reconnection. Here, in this study, we present simulations tailored to reproduce the specific TREX geometry, which confirm both the reconnection rate scale and the shock jump conditions previously characterized experimentally in TREX. The simulations also establish an interplay between the reconnection layer and the Alfvénic expansions of the background plasma associated with the energization of the TREX drive coils; this interplay has not yet been experimentally observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Temperature Measurements in Hypersonic Wind Tunnels via Femtosecond Coherent Anti-Stokes Raman Scattering

A femtosecond coherent anti-Stokes Raman scattering (fs CARS) instrument is developed to perform gas-phase thermometry in cold-flow hypersonic wind tunnels. Measurements are reported for Mach 8 and 14 pure-nitrogen flows. The fs CARS instrument includes a 100 fs pump/Stokes pulse and a spectrally narrow probe pulse from a second harmonic bandwidth compressor. Important experimental considerations such as limits on the pump/Stokes pulse energy are discussed. The fs CARS focusing and collimating optics are mounted on a two-axis translation stage system to scan the measurement location during a 30 second wind tunnel run. Single-laser-shot rotational CARS spectra are recorded at the laser repetition rate of 1 kHz in the wind tunnel freestream and near simple cone models. Spectral fitting is used to determine quantitative gas temperatures. Freestream temperatures at Mach 8 and 14 spanned ranges of 40–75 and 35–50 K, respectively, depending on tunnel operating conditions. Temperature variations across the central 100 mm span of the wind tunnel were quantified. Measured temperature jumps across conical bow shocks from various models varied by less than 1% from predicted values. Hypersonic boundary layer measurements were demonstrated. In conclusion, these measurements illustrate the utility and robustness of this instrument for the study of complex fluid flow phenomena in challenging ground test facilities.

Aerodynamics↗
DOE OSTI2022

Perpendicular Subcritical Shock Structure in a Collisional Plasma Experiment

We present a study of perpendicular subcritical shocks in a collisional laboratory plasma. Shocks are produced by placing obstacles into the supermagnetosonic outflow from an inverse wire array z pinch. We demonstrate the existence of subcritical shocks in this regime and find that secondary shocks form in the downstream. Detailed measurements of the subcritical shock structure confirm the absence of a hydrodynamic jump. We calculate the classical (Spitzer) resistive diffusion length and show that it is approximately equal to the shock width. As a result, we measure little heating across the shock (<10 % of the ion kinetic energy) which is consistent with an absence of viscous dissipation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Slow Shock Formation Upstream of Reconnecting Current Sheets

The formation, development, and impact of slow shocks in the upstream regions of reconnecting current layers are explored. Slow shocks have been documented in the upstream regions of magnetohydrodynamic (MHD) simulations of magnetic reconnection as well as in similar simulations with the kglobal kinetic macroscale simulation model. They are therefore a candidate mechanism for preheating the plasma that is injected into the current layers that facilitate magnetic energy release in solar flares. Of particular interest is their potential role in producing the hot thermal component of electrons in flares. During multi-island reconnection, the formation and merging of flux ropes in the reconnecting current layer drives plasma flows and pressure disturbances in the upstream region. These pressure disturbances steepen into slow shocks that propagate along the reconnecting component of the magnetic field and satisfy the expected Rankine–Hugoniot jump conditions. Plasma heating arises from both compression across the shock and the parallel electric field that develops to maintain charge neutrality in a kinetic system. Shocks are weaker at lower plasma β, where shock steepening is slow. While these upstream slow shocks are intrinsic to the dynamics of multi-island reconnection, their contribution to electron heating remains relatively minor compared with that from Fermi reflection and the parallel electric fields that bound the reconnection outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Tracker Fits for shot 1s1674

Gas gun shock-to-detonation experiments provide data for calibrating HE burn models. At LANL the principal diagnostic is a magnetic gauge package that provides data on the lead shock trajectory, which is used to determine the run-distance to detonation, and the lead shock speed with run distance, which together with jump-off velocity of the velocity gauges determines Hugoniot locus points for the reactants.

42 ENGINEERING↗
DOE OSTI2024

Electron Heating in 2D Particle-in-cell Simulations of Quasi-perpendicular Low-beta Shocks

Abstract We measure the thermal electron energization in 1D and 2D particle-in-cell simulations of quasi-perpendicular, low-beta ( β p = 0.25) collisionless ion–electron shocks with mass ratio m i / m e = 200, fast Mach number M ms = 1 –4, and upstream magnetic field angle θ Bn = 55°–85° from the shock normal n ˆ . It is known that shock electron heating is described by an ambipolar, B -parallel electric potential jump, Δ ϕ ∥ , that scales roughly linearly with the electron temperature jump. Our simulations have Δ ϕ ∥ / ( 0.5 m i u sh 2 ) ∼ 0.1 –0.2 in units of the pre-shock ions’ bulk kinetic energy, in agreement with prior measurements and simulations. Different ways to measure ϕ ∥ , including the use of de Hoffmann–Teller frame fields, agree to tens-of-percent accuracy. Neglecting off-diagonal electron pressure tensor terms can lead to a systematic underestimate of ϕ ∥ in our low- β p shocks. We further focus on two θ Bn = 65° shocks: a M s = 4 ( M A = 1.8 ) case with a long, 30 d i precursor of whistler waves along n ˆ , and a M s = 7 ( M A = 3.2 ) case with a shorter, 5 d i precursor of whistlers oblique to both n ˆ and B ; d i is the ion skin depth. Within the precursors, ϕ ∥ has a secular rise toward the shock along multiple whistler wavelengths and also has localized spikes within magnetic troughs. In a 1D simulation of the M s = 4 , θ Bn = 65° case, ϕ ∥ shows a weak dependence on the electron plasma-to-cyclotron frequency ratio ω pe /Ω ce , and ϕ ∥ decreases by a factor of 2 as m i / m e is raised to the true proton–electron value of 1836.

Tran, Aaron (ORCID:0000000334834890)↗
DOE OSTI2022

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗
DOE OSTI2023

Observation of subcritical shocks in a collisional laboratory plasma: scale dependence near the resistive length

We present a study of subcritical shocks in a highly collisional laboratory plasma with a dynamically significant magnetic field. Shocks were produced by placing cylindrical obstacles into the supermagnetosonic ( $M_{{\rm MS}} \sim 1.9$ ) outflow from an inverse wire array z-pinch at the MAGPIE pulsed power facility ( $n_{e} \sim 8.5 \times 10^{17}\,{\rm cm}^{-3}$ , $v \sim 45\,{\rm km}\,{\rm s}^{-1}$ ). We demonstrate the existence of subcritical shocks in this regime and find that secondary stagnation shocks form in the downstream which we infer from interferometry and optical Thomson scattering measurements are hydrodynamic in nature. The subcritical shock width is found to be approximately equal to the resistive diffusion length and we demonstrate the absence of a jump in hydrodynamic parameters. Temperature measurements by collective optical Thomson scattering showed little temperature change across the subcritical shock ( ${<}10\,\%$ of the ion kinetic energy) which is consistent with a balance between adiabatic and Ohmic heating and radiative cooling. We demonstrate the absence of subcritical shocks when the obstacle diameter is less than the resistive diffusion length due to decoupling of the magnetic field from the plasma. These findings are supported by magnetohydrodynamic simulations using the Gorgon and AstroBEAR codes and discrepancies between the simulations and experiment are discussed.

Physics↗
DOE OSTI2022

Effect of Topology on Transient Dynamic and Shock Response of Polymeric Lattice Structures

Architected cellular materials, such as lattice structures, offer potential for tunable mechanical properties for dynamic applications of energy absorption and impact mitigation. In this work, the static and dynamic behavior of polymeric lattice structures was investigated through experiments on octet-truss, Kelvin, and cubic topologies with relative densities around 8%. Here, dynamic testing was conducted via direct impact experiments (25–70 m/s) with high-speed imaging coupled with digital image correlation and a polycarbonate Hopkinson pressure bar. Mechanical properties such as elastic wave speed, deformation modes, failure properties, particle velocities, and stress histories were extracted from experimental results. At low impact velocities, a transient dynamic response was observed which was composed of a compaction front initiating at the impact surface and additional deformation bands whose characteristics matched low strain-rate behavior. For higher impact velocities, shock analysis was carried out using compaction wave velocity and Eulerian Rankine–Hugoniot jump conditions with parameters determined from full-field measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗
DOE OSTI2025

Biermann-Battery-Driven Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

This Letter reports the first complete observation of magnetized collisionless shock precursors formed through the compression of Biermann-battery magnetic fields in laser produced plasmas. At OMEGA, lasers produce a supersonic CH plasma flow which is magnetized with Biermann-battery magnetic fields. The plasma flow collides with an unmagnetized hydrogen gas jet plasma to create a magnetized shock precursor. The situation where the flowing plasma carries the magnetic field is similar to the Venusian bow shock. Imaging 2⁢𝜔 Thomson scattering confirms that the interaction is collisionless and shows density and temperature jumps. Proton radiographs have regions of strong deflections and FLASH magnetohydrodynamic (MHD) simulations show the presence of Biermann fields in the Thomson scattering region. Electrons are accelerated to energies of up to 100 keV in a power-law spectrum. Here, OSIRIS particle-in-cell (PIC) simulations, initialized with measured parameters, show the formation of a magnetized shock precursor and corroborate the experimental observables.

High-energy-density plasmas↗
DOE OSTI2024

A High-Order Eulerian–Lagrangian Runge–Kutta Finite Volume (EL–RK–FV) Method for Scalar Nonlinear Conservation Laws

Abstract We present a class of high-order Eulerian–Lagrangian Runge–Kutta finite volume methods that can numerically solve Burgers’ equation with shock formations, which could be extended to general scalar conservation laws. Eulerian–Lagrangian (EL) and semi-Lagrangian (SL) methods have recently seen increased development and have become a staple for allowing large time-stepping sizes. Yet, maintaining relatively large time-stepping sizes post shock formation remains quite challenging. Our proposed scheme integrates the partial differential equation on a space-time region partitioned by linear approximations to the characteristics determined by the Rankine–Hugoniot jump condition. We trace the characteristics forward in time and present a merging procedure for the mesh cells to handle intersecting characteristics due to shocks. Following this partitioning, we write the equation in a time-differential form and evolve with Runge–Kutta methods in a method-of-lines fashion. High-resolution methods such as ENO and WENO-AO schemes are used for spatial reconstruction. Extension to higher dimensions is done via dimensional splitting. Numerical experiments demonstrate our scheme’s high-order accuracy and ability to sharply capture post-shock solutions with large time-stepping sizes.

Chen, Jiajie↗
DOE OSTI2025

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING↗
DOE OSTI2022

Hybrid interferometric diagnostic for the refractive index measurement across nearly discrete shock waves

The measurement of high temperature gas properties is key for characterizing high speed flows. Nearly discrete changes in density across shock waves, in particular, are difficult to resolve through traditional fringe-counting interferometric methods. Existing techniques for estimating large fringe jumps are either resolution limited or require specialized window configurations. In this Letter, we describe a unique hybrid interferometric technique that combines narrowband fringes for high resolution and broadband fringes as an absolute reference to measure changes in refractive index with a resolution of up to 7 × 10 −8 across nearly discrete index changes of up to 1.5 × 10 −4 . By capturing fringes with an ultrahigh-speed camera, the refractive index changes across discrete shock fronts can be estimated inside a shock tube with high accuracy and time resolution. First, a novel hybrid calibration method for tracking finite fringes is discussed. Next, this technique is used to measure the post-initial-shock refractive indices for Mach 2.7 to 4.2 flows (pressures from 90.4 to 228.4 kPa). Results are then compared with theoretical values showing agreement within 2%.

Wang, Gwendolyn T. (ORCID:000000016271893X)↗