Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK FRONT”

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 37 records · Page 2

Symmetry and scaling in one-dimensional compressible two-phase flow

Investigations of shock compression of heterogeneous materials often focus on the shock front width and overall profile. The number of experiments required to fully characterize the dynamic response of a material often belie the structure–property relationships governing these aspects of a shock wave. Recent observations measured a pronounced shock-front width on the order of 10 s of ns in particulate composites. We focus on particulate composites with disparate densities and investigate whether the mechanical interactions between the phases are adequate to describe this emergent behavior. The analysis proceeds with a general Mie–Grüneisen equation of state for the matrix material, a general drag force law with general power-law scaling for the particle-matrix coupling of the phases, and a volume fraction-dependent viscosity. Lie group analysis is applied to one-dimensional hydrodynamic flow equations for the self-consistent interaction of particles embedded in a matrix material. The particle phase is characterized by a particle size and volume fraction. The Lie group analysis results in self-similar solutions reflecting the symmetries of the flow. The symmetries lead to well-defined scaling laws, which may be used to characterize the propagation of shock waves in particle composites. An example of the derived scaling laws for shock attenuation and rise time is shown for experimental data on shock-driven tungsten-loaded polymers. A key result of the Lie analysis is that there is a relationship between the exponents characterizing the form of the drag force and the exponent characterizing the shock velocity and its attenuation in a particulate composite. Comparison to recent experiments results in a single exponent that corresponds to a conventional drag force.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The effects of mixture preburning on detonation wave propagation

Pressure gain combustion in the form of continuous detonations can provide a significant increase in the efficiency of a variety of propulsion and energy conversion devices. In this regard, rotating detonation engines (RDEs) that utilize an azimuthally-moving detonation wave in annular systems are increasingly seen as a viable approach to realizing pressure gain combustion. However, practical RDEs that employ non- premixed fuel and oxidizer injection need to minimize losses through a number of mechanisms, including turbulence-induced shock-front variations, incomplete fuel-air mixing, and premature deflagration. In this study, a canonical stratified detonation configuration is used to understand the impact of preburning on detonation efficiency. It was found that heat release ahead of the detonation wave leads to weaker shock fronts, delayed combustion of partially-oxidized fuel-air mixture, and non-compact heat release. Furthermore, large variations in wave speeds were observed, which is consistent with wave behavior in full-scale RDEs. Peak pressures in the compression region or near triple points were considerably lower than the theoretically-predicted values for ideal detonations. Analysis of the detonation structure indicates that this deflagration process is parasitic in nature, reducing the detonation efficiency but also leading to heat release far behind the wave that cannot directly strengthen the shock wave. As a result, this parasitic combustion leads to commensal combustion (heat release far downstream of the wave), indicating that it is the root cause of combustion efficiency losses.

42 ENGINEERING↗

Compensating slice emittance growth in high brightness photoinjectors using sacrificial charge

Achieving maximum electron beam brightness in photoinjectors requires detailed control of the 3D bunch shape and precise tuning of the beam focusing. Even in state-of-the-art designs, slice emittance growth due to nonlinear space charge forces and partial nonlaminarity often remains non-negligible. In this work, we introduce a new means to linearize the transverse slice phase space: a sacrificial portion of the bunch’s own charge distribution, formed into a wavebroken shock front by highly nonlinear space charge forces within the gun, whose downstream purpose is to dynamically linearize the desired bunch core. We show that linearization of an appropriately prepared bunch can be achieved via strongly nonlaminar focusing of the sacrificial shock front, while the inner core focuses laminarly. This leads to a natural spatial separation of the two distributions: a dense core surrounded by a diffuse halo of sacrificial charge that can be collimated. Multiobjective genetic algorithm optimizations of the ultracompact x-ray free electron laser injector employ this concept, and we interpret it with an analytic model that agrees well with the simulations. In simulation, we demonstrate a final bunch charge of 100 pC, peak current ∼ 30 A, and a sacrificial charge of 150 pC (250 pC total emitted from cathode) with normalized emittance growth of < 20 nm rad due to space charge. This implies a maximum achievable brightness approximately an order of magnitude greater than existing free electron laser injector designs. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

Role of opacity at the 9 keV back lighter energy used in measuring the equation of state of boron at pressures up to a Gbar

Many experiments have been conducted at the National Ignition Facility to measure the equation of state Hugoniot of plastic, boron, and diamond at extreme pressures up to a Gbar. The “Gbar” design employs a strong spherically converging shock launched through a solid ball of material using a hohlraum radiation drive. The shock front conditions are characterized using X-ray radiography, typically at energies near 9 keV. Here we examine how the opacity of boron at 9 keV changes at high pressures and temperatures. Understanding this is vital to unfolding the density in the shock front as pressures exceed 100 Mbar. We compare opacity calculations from a number of methods including the legacy XSN opacity tables, super transition array (STA), average atom (AVAT), and detailed configuration accounting (DCA) methods. We examine how the changing opacity is correlated with the K-shell occupation calculated using the opacity methods as well as other electronic structure methods such as MECCA, Purgatorio, and an all-electron Kohn-Sham density functional theory calculation that uses optimized all-electron norm-conserving Vanderbilt pseudopotentials (ONCV). We also examine the relative contribution of free-free opacity to the overall opacity, which is dominated by bound-free absorption.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Experimental characterization of hot-electron emission and shock dynamics in the context of the shock ignition approach to inertial confinement fusion

We report on planar target experiments conducted on the OMEGA-EP laser facility performed in the context of the shock ignition (SI) approach to inertial confinement fusion. The experiment aimed at characterizing the propagation of strong shock in matter and the generation of hot electrons (HEs), with laser parameters relevant to SI (1-ns UV laser beams with I ~10 16 W/cm 2 ). Time-resolved radiographs of the propagating shock front were performed in order to study the hydrodynamic evolution. The hot-electron source was characterized in terms of Maxwellian temperature, T h , and laser to hot-electron energy conversion efficiency η using data from different x-ray spectrometers. The post-processing of these data gives a range of the possible values for T h and η [i.e., T h [keV] ϵ (20, 50) and η ϵ (2%, 13%)]. These values are used as input in hydrodynamic simulations to reproduce the results obtained in radiographs, thus constraining the range for the HE measurements. According to this procedure, we found that the laser converts ~10% ± 4% of energy into hot electrons with T h = 27 ± 8 keV. The paper shows how the coupling of different diagnostics and numerical tools is required to sufficiently constrain the problem, solving the large ambiguity coming from the post-processing of spectrometers data. The effect of the hot electrons on the shock dynamics is then discussed, showing an increase in the pressure around the shock front. Furthermore, the low temperature found in this experiment without pre-compression laser pulses could be advantageous for the SI scheme, but the high conversion efficiency may lead to an increase in the shell adiabat, with detrimental effects on the implosion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The structure of 3-D collisional magnetized bow shocks in pulsed-power-driven plasma flows

Here we investigate three-dimensional (3-D) bow shocks in a highly collisional magnetized aluminium plasma, generated during the ablation phase of an exploding wire array on the MAGPIE facility (1.4 MA, 240 ns). Ablation of plasma from the wire array generates radially diverging, supersonic ($M_S$ $\sim 7$), super-Alfvénic ($M_A>1$) magnetized flows with frozen-in magnetic flux $(R_M\gg1)$. These flows collide with an inductive probe placed in the flow, which serves both as the obstacle that generates the magnetized bow shock, and as a diagnostic of the advected magnetic field. Laser interferometry along two orthogonal lines of sight is used to measure the line-integrated electron density. A detached bow shock forms ahead of the probe, with a larger opening angle in the plane parallel to the magnetic field than in the plane normal to it. Since the resistive diffusion length of the plasma is comparable to the probe size, the magnetic field decouples from the ion fluid at the shock front and generates a hydrodynamic shock, whose structure is determined by the sonic Mach number, rather than the magnetosonic Mach number of the flow. The 3-D simulations performed using the resistive magnetohydrodynamic (MHD) code Gorgon confirm this picture, but under-predict the anisotropy observed in the shape of the experimental bow shock, suggesting that non-MHD mechanisms may be important for modifying the shock structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron energization dynamics in interaction of self-generated magnetic vortices in upstream of collisionless electron/ion shocks

Abstract Relativistic collisionless shocks are considered responsible for particle energization mechanisms leading to particle acceleration. While electron energization in shock front region of electron/ion collisionless shocks are the most studied, the mechanism of electron energization in interaction with self-generated magnetic vortices (MVs) in the upstream region is still unclear. We investigate electron energization mechanism in the upstream region of electron/ion relativistic collisionless shocks, using two dimensional particle-in-cell (PIC) simulations. We discuss mechanism of electron energization which takes place in the upstream region of the shock, where the counter stream particles interact with incoming flow. The energy gain of electrons happens during their interaction with evolving fields of self-generated magnetic vortices in this region. Three Fermi-like electron energization scenarios are discussed. Stochastic acceleration of electrons in interaction with fields of MV leads to anisotropic heating of fast electrons due to diffusion in the momentum space of electrons and, finally, synergetic effect of evolving fields of MVs leads to the formation of a power-law tail of supra-thermal particles.

Naseri, N.↗

High-speed x-ray phase contrast imaging and digital image correlation analysis of microscale shock response of an additively manufactured energetic material simulant

The performance of energetic materials subjected to dynamic loading significantly depends on their micro- and meso-scale structural morphology. The geometric versatility offered by additive manufacturing opens new pathways to tailor the performance of these materials. Additively manufactured energetic materials (AMEMs) have a wide range of structural characteristics with a hierarchy of length scales and process-inherent heterogeneities, which are hitherto difficult to precisely control. It is important to understand how these features affect AMEMs’ response under dynamic/shock loading. Therefore, temporally and spatially resolved measurements of both macroscopic behavior and micro- and meso-level processes influencing macroscopic behavior are required. In this paper, we analyze the shock compression response of an AMEM simulant loaded under several impact conditions and orientations. Furthermore, x-ray phase contrast imaging (PCI) is used to track features across the observed shock front and determine the linear shock velocity vs particle velocity equation of state, as well as to quantify the interior deformation fields via digital image correlation (DIC) analyses. Photon Doppler velocimetry is simultaneously used to measure the particle velocities of the specimens, which are consistent with those obtained from x-ray PCI. The DIC analyses provide an assessment of the average strain fields inside the material, showing that the average axial strain depends on the loading intensity and reaches as high as 0.23 for impact velocities up to 1.5 km/s. The overall results demonstrate the utility of x-ray PCI for probing “in-material” equation of state and interior strains associated with dynamic shock compression behavior of the AMEM simulant.

3D printing↗

X-ray phase-contrast imaging of strong shocks on OMEGA EP

The ongoing improvement in laser technology and target fabrication is opening new possibilities for diagnostic development. An example is x-ray phase-contrast imaging (XPCI), which serves as an advanced x-ray imaging diagnostic in laser-driven experiments. In this work, we present the results of the XPCI platform that was developed at the OMEGA EP Laser-Facility to study multi-Mbar single and double shocks produced using a kilojoule laser driver. Two-dimensional radiation-hydrodynamic simulations agree well with the shock progression and the spherical curvature of the shock fronts. It is demonstrated that XPCI is an excellent method to determine with high accuracy the front position of a trailing shock wave propagating through an expanding CH plasma that was heated by a precursor Mbar shock wave. Furthermore, the interaction between the rarefaction wave and the shock wave results in a clear signature in the radiograph that is well reproduced by radiation-hydrodynamic simulations.

47 OTHER INSTRUMENTATION↗

Linear Grains

LIGR (Linear Grains) implements a linear (Fourier mode) analysis of a shock propagating through density perturbations. It has two primary capabilities. First, it calculates the shock front distortion and post-shock flow variables (density, velocity, pressure) for single Fourier mode pairs in two dimensions (a pair of x and y modes). Second, it calculates these same quantities for a grains-like pre-shock density perturbation, which consists of a sum of many Fourier modes. This grains-like perturbation is a periodic grid of rectangles at a specifiable density and size (the grains), separated by regions of a different specifiable size and density (interstitial spaces between the grains). The single-mode analysis is based on the work of Velikovich et al., Phys. Plasmas 14 072706 (2007).

Li, GraceJ↗

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)↗

Absolute Hugoniot measurements in low-density plastic foams on the NIKE shock compression platform

The propagation of Mbar-range shock waves in low-density foam materials is of interest to target design in all approaches to inertial confinement fusion, high energy density physics, and laboratory astrophysics. We report absolute Hugoniot measurements for CH 1.72 N 0.086 O 0.37 plastic foams with 73 and 94 mg/cm 3 densities in the 32–107 km/s shock velocity range. The experiments were performed on the shock compression platform developed on the NIKE KrF laser facility at the U.S. Naval Research Laboratory. NIKE's 4 or 8-ns long flat-top laser pulses drive steady shock waves into foam targets at the ablative pressures of 1–7 Mbar. The propagation of the ablation and shock fronts is tracked continuously in time using streaked side-on monochromatic x-ray imaging radiography. The straight x−t trajectories of the shock and ablation fronts in the recorded streak images confirmed their steadiness. The SESAME tabulated equation of state predictions generally agrees with our Hugoniot data within a 95% confidence band. The experimental uncertainty of the evaluated shock density compression ratios remains large, indicating the need for more experiments and improved theoretical understanding of the strong shock propagation mechanisms in dry foams.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron heat flux and propagating fronts in plasma thermal quench via ambipolar transport

The thermal collapse of a nearly collisionless plasma interacting with a cooling spot, in which the electron parallel heat flux plays an essential role, is both theoretically and numerically investigated. Here, we show that such thermal collapse, which is known as thermal quench in tokamaks, comes about in the form of propagating fronts, originating from the cooling spot, along magnetic field lines. The slow fronts, propagating with local ion sound speed, limit the aggressive cooling of plasma, which is accompanied by a plasma cooling flow toward the cooling spot. The extraordinary physics underlying such a cooling flow is that the fundamental constraint of ambipolar transport along the field line limits the spatial gradient of electron thermal conduction flux to the much weaker convective scaling, as opposed to the free-streaming scaling, so that a large electron temperature and, hence, pressure gradient can be sustained. The last ion front for a radiative cooling spot is a shock front where cold but flowing ions meet hot ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation to a Newborn Supernova Remnant from a Low-mass Iron Core Star

Supernova remnant observations show a high degree of asymmetry, mixing, and inhomogeneity. These asymmetries are seeded during the early seconds of the explosion and are further enhanced and modified as the shock and ejecta move through the stellar progenitor and into the circumstellar medium. We present simulations of a 9.6 M⊙ zero-metallicity progenitor initialized after shock revival and evolved for several years when the ejecta is in the circumstellar medium. A suite of 1D and 2D simulations examines the effects of neutron-star wind and radioactive decay heating. In 1D, decay heating forms a low-density bubble that suppresses the reverse shock. While in 2D, the heating is localized to metal-rich pockets, inflating them and compressing the surrounding material into dense shells. In 3D, the neutron-star wind and decay heating modify the plume morphology, producing more large-scale structures. The extended plume morphology leads to an asymmetrical shock breakout. After breakout, the leading plumes cannot keep up with the shock front, resulting in deceleration and fragmentation by the reverse shock while retaining the large-scale asymmetry. The projected ejecta morphology and velocities are strongly viewing angle dependent. The relatively uniform metal-rich distribution does not resemble the strongly inhomogeneous ejecta structure of Cas A. The 160-isotope decay network shows that 24.4% of the radioactive heating comes from decay chains other than the canonical 56Ni chain. The low explosion energy, low 56Ni yield, and Ni/Fe ratio greater than unity suggest an observational signature similar to an electron capture supernova.

Neopane, Sudarshan [University of Tennessee (UT)]↗

The Effect of Large-scale Magnetic Fluctuations on Energetic Particle Profiles across Interplanetary Shocks

In situ observations of energetic particles associated with interplanetary (IP) shocks, known as energetic storm particle (ESP) events, often show time-intensity profiles quite different from predictions of classical diffusive shock acceleration theory. We use numerical simulations, including test-particle simulations and hybrid simulations (with fluid electrons and kinetic protons) for shocks propagating into a turbulent magnetic field to study the ESP intensity-time profiles across a strong IP shock. We find that several types of energetic particle intensity profiles similar to in situ observations can be produced in our simulations. The peak of energetic particle count is often near, but can be shifted from, the locations of the shock front, with the peak usually being observed downstream of the shock. These findings may help understand particle acceleration at both traveling IP shocks and shocks in other heliospheric and astrophysical environments.

Lo, Ming-Yuan [New Mexico Consortium (NMC), Los Al↗

Development of an x-ray radiography platform to study laser-direct-drive energy coupling at the National Ignition Facility

A platform has been developed to study laser-direct-drive energy coupling at the National Ignition Facility (NIF) using a plastic sphere target irradiated in a polar-direct-drive geometry to launch a spherically converging shock wave. To diagnose this system evolution, eight NIF laser beams are directed onto a curved Cu foil to generate Heα line emission at a photon energy of 8.4 keV. These x rays are collected by a 100-ps gated x-ray imager in the opposing port to produce temporally gated radiographs. The platform is capable of acquiring images during and after the laser drive launches the shock wave. A backlighter profile is fit to the radiographs, and the resulting transmission images are Abel inverted to infer radial density profiles of the shock front and to track its temporal evolution. Here, the measurements provide experimental shock trajectories and radial density profiles that are compared to 2D radiation-hydrodynamic simulations using cross-beam energy transfer and nonlocal heat-transport models.

47 OTHER INSTRUMENTATION↗

Solar Energetic Particle Acceleration at a Spherical Shock with the Shock Normal Angle $θ$ $B_n$ Evolving in Space and Time

We present a 2D kinematic model to study the acceleration of solar energetic particles (SEPs) at a shock driven by a coronal mass ejection. The shock is assumed to be spherical about an origin that is offset from the center of the Sun. This leads to a spatial and temporal evolution of the angle between the magnetic field and the shock-normal direction (θ Bn ) as it propagates through the Parker spiral magnetic field from the lower corona to 1 au. We find that the high-energy SEP intensity varies significantly along the shock front due to the evolution of θ Bn . Generally, the west flank of the shock preferentially accelerates particles to high energies compared to the east flank and shock nose. This can be understood in terms of the rate of acceleration, which is higher at the west flank. Double power-law energy spectra are reproduced in our model as a consequence of the local acceleration and transport effects. These results will help us to better understand the evolution of SEP acceleration and provide new insights into large SEP events observed by multiple spacecraft, especially those close to the Sun, such as Parker Solar Probe and Solar Orbiter.

79 ASTRONOMY AND ASTROPHYSICS↗