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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Numerical study on spray collapse process of ECN spray G injector under flash boiling conditions

This work features a numerical study of the spray collapse process that occurs under flash-boiling, gasoline-direct-injection conditions. The Engine Combustion Network's Spray G injector was modeled, and iso-octane, n-hexane, and propane were used as fuels with the goal to investigate fuel property effects. The development of individual jets and jet-to-jet interaction were analyzed to reveal the physical mechanisms of spray collapse. The results agreed qualitatively well with the experimental data of propane available from the literature and the predictions revealed the model's ability to capture the near-nozzle shock structures witnessed in the experiments. Simulations showed that the flow development of flashing sprays was strongly tied with the fuel properties. Specifically, under-expanded jets were not observed with iso-octane due to the fuel's low vapor pressure. In the case of n-hexane, the jets were under-expanded, and a low-pressure core was observed inside the jets due to the flow expansion. With propane, the very strong expansion caused the flow to become transonic and generate a unique shock structure near the nozzle exit. The shock structure revealed the formation of primary shock cells within the jets and secondary shock cells induced by jet-to-jet interaction. Based on the results, the spray collapse under flash boiling conditions was mainly attributed to two effects, i.e., the decrease in pressure at the spray center induced by the low-pressure cores inside and between the jets, and the isolation effect due to the jet-to-jet interaction which contributed to maintain the pressure differential between the spray center and the ambient.

42 ENGINEERING↗
DOE OSTI2020

Mixing in ICF implosions on the National Ignition Facility caused by the fill-tube

The micrometer-scale tube that fills capsules with thermonuclear fuel in inertial confinement fusion experiments at the National Ignition Facility is also one of the implosion's main degradation sources. It seeds a perturbation that injects the ablator material into the center, radiating away some of the hot-spot energy. This paper discusses how the perturbation arises in experiments using high-density carbon ablators and how the ablator mix interacts once it enters the hot-spot. Both modeling and experiments show an in-flight areal-density perturbation and localized x-ray emission at stagnation from the fill-tube. Simulations suggest that the fill-tube is degrading an otherwise 1D implosion by ~2×, but when other degradation sources are present, the yield reduction is closer to 20%. Characteristics of the fill-tube assembly, such as the through-hole size and the glue mass, alter the dynamics and magnitude of the degradation. Overall, these aspects point the way toward improvements in the design, some of which (smaller diameter fill-tube) have already shown improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Causes of fuel–ablator mix inferred from modeling of monochromatic time-gated radiography of OMEGA cryogenic implosions

Here, we present evidence, in the context of OMEGA cryogenic target implosions, that laser imprint, known to be capable of degrading laser-direct-drive target performance, plays a major role in generating fuel–ablator mix. OMEGA cryogenic target implosions show a performance boundary correlated with acceleration-phase shell stability; for sufficiently low adiabats (where the adiabat is the ratio of the pressure to the Fermi pressure) and high in-flight aspect ratios (IFAR's), the neutron-weighted shell areal density and neutron yield relative to the clean simulated values sharply decline. Direct evidence of Rayleigh–Taylor fuel–ablator mixing was previously obtained using a Si Heα backlighter driven by an ~20-ps short pulse generated by OMEGA EP. The shadow cast by the shell shortly prior to stagnation, as diagnosed using backlit radiographs, shows a softening near the limb, which is evidence of an ablator–fuel mix region for a low-adiabat implosion (α ~ 1.9, IFAR = 14) but not for a moderate adiabat implosion (α ~ 2.5, IFAR = 10). We find good agreement between experimental and synthetic radiographs in simulations that model laser imprint and account for uncertainty in the initial ablator thickness. Furthermore, we further explore the role of other mechanisms such as classical instability growth at the fuel–ablator interface, species concentration diffusion, and long-wavelength drive and target asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The effect of hardness on polymer-bonded pentaerythritol tetranitrate (PETN) explosive impact sensitivity

Handling sensitivity is an important property to assess when working with explosive samples and can be measured using a variety of tests, including drop-weight impact sensitivity. There exists a longstanding interest in the explosives community on the importance of measurable chemical, physical, and mechanical properties of explosives in impact sensitivity. However, most recent work in this area has explored chemical attributes rather than physical and mechanical properties of explosives. In this study, we measure hardness of explosive samples of pentaerythritol tetranitrate and Sylgard binder (XTX) during the curing process. The samples have been characterized for particle morphology through scanning electron microscopy and handling sensitivity through drop-weight impact testing. The relative importance of states of cure, methods of curing, morphology, and age of material are discussed. The data indicate that although there is a notable difference in morphology and mechanical properties for the samples as the polymer-bonded mixtures cure, the resulting changes to mechanical properties have a minimal effect on the sensitivity of the XTX.

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

Effect of laser preheat in magnetized liner inertial fusion at OMEGA

Integrated magnetized liner inertial fusion (MagLIF) experiments on the OMEGA laser system are modeled to investigate the effects of laser preheat on implosion performance. In both simulations and experiments, preheating the fuel with a laser enhances the neutron yield, with the maximum yield occurring at an optimal preheat laser energy. Increasing the preheat energy past the optimal value reduces the neutron yield. In simulations, the rate of the drop in neutron yield depends on the inclusion of the Nernst effect. The Nernst effect advects the magnetic field away from the center of the fuel region during the preheat stage of MagLIF on OMEGA and weakens magnetic-flux compression. If the Nernst effect is not included, the simulated drop in yield past the optimal preheat laser energy is more gradual as opposed to the sharper drop seen in experiments. Simulations are able to model trend in measured ion temperatures seen in experiments. Here, a mix model is used to show that including mix from the wall in the simulated fuel region further degrades the yield and lowers the optimal preheat laser energy. Simulations with mix predict that increasing the initial axial magnetic field may still enhance yield performance in integrated implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Viscosity Measurements in Extreme Conditions

Radiation-hydrodynamics simulations are critical to the NNSA complex for understanding high energy density physics and for addressing problems in national security science. However, current rad-hydro codes typically do not account for viscosity, which generally leads to large uncertainties in problems involving mixing of materials under shock loading. This report details progress made towards experimental measurements of fluid viscosity at high pressures and temperatures achieved by shock wave compression techniques, as well as complementary material model development and hydrodynamic simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulation of a strong steady-state plasma shock in a warm dense matter regime

The structure of collisional plasma shocks has been subject to an extensive, multi-decadal investigation—in the hydrodynamic, hybrid kinetic ion/electron fluid, and fully kinetic ion/electron limits. Despite this thoroughness, all of these studies apply exclusively to classical, weakly coupled plasmas. Here, we report the first results for a planar hydrodynamic simulation of a strong, steady-state shock in a subspace of the warm dense matter (WDM) regime. Specifically, we consider a plasma of fully degenerate electrons with moderate-to-strongly coupled ions. Since the WDM ion and electron transport coefficients and equation of state differ markedly from their non-degenerate, weak-coupling equivalents, we find that the structure of a WDM plasma shock notably deviates from the ideal plasma picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The effects of laser pulse length and collisional ionization on the acceleration of titanium ions

Abstract The interaction of a relativistic laser pulse ( > 10 18 W cm −2 µ m 2 ) with foil targets can accelerate ions to energies of tens of MeV u −1 with optimized laser and target parameters. We report the results on simulations of the interaction of a 6.0 × 10 20 W cm −2 laser pulse incident on ultrathin (10–500 nm) titanium foils to investigate the roles of laser pulse duration and ionization mechanisms in the acceleration of titanium ions. While holding peak intensity constant, two laser pulse durations were investigated, 140 and 650 fs. The optimum thickness is dependent on pulse duration, as it requires the concurrence of target transparency with the incidence of the peak laser intensity. The collisional processes do not play a significant role in Ti ion beam generation from the 140 fs laser pulse duration at the optimum thickness (30 nm). However, for the 650 fs laser optimum, collisions improve the conversion efficiency of highly energetic ( > 10 MeV u −1 ), high charge titanium (Ti 21 − 22 + ) by a factor of 20, and the titanium ion cutoff energy by ∼15%. This improvement is due to the fact that collisional ionization increases the electron density of the plasma, which delays the time of relativistic transparency, causing collisions to decrease the optimum foil thickness from 150 to 100 nm. Additionally, collisional ionization increases the charge-to-mass ratio of the titanium, and injects more electrons into the accelerating sheath field. At the optimum thickness, target normal sheath acceleration is the dominant mechanism of acceleration, with additional contributions from radiation pressure and shock wave acceleration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Using double pulse laser ablation in air to enhance the strength of laser-driven shocks

In the process of multi-pulse laser ablation, inter-pulse delay time, Δ t , is known to be an important parameter for maximizing ablation efficiency as well as impulse imparted to the target. In this work, using photon Doppler velocimetry, we show that for single pairs of colinear pulses (1064 nm, 8 ns, ∼ 60 J cm -2 per pulse) in air, the peak free surface velocity of the back surface of an aluminum target (125 µm thick) is increased, by a factor of nearly 3, when Δ t = 10 microseconds, compared with both pulses arriving simultaneously (Δ t = 0). Fast imaging of the ablation process suggests this enhancement is due to rarefaction of the contiguous air in the passage of the leading shock produced by ablation, which then in turn allows a larger fraction of the energy of the second pulse to reach the target surface. This interpretation is strengthened by additional experiments in which the two pulses do not overlap on the target surface, but the shock strength is nevertheless enhanced. Given a fixed energy budget this work suggests a prescription for maximizing laser-driven shock strength by judicious choice of inter-pulse delay.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fall SULI Report

This fall, I interned under Fan Guo and Xiaocan Li in T-2 at Los Alamos National Laboratory through the SULI program. I investigated how shock thickness is impacted by collisions between particles in plasmas and how particle acceleration in shocks is impacted by the mass ratio for protons and electrons used in simulations. Both investigations were carried out using the VPIC simulation code. I found that collisions do not meaningfully impact shock width and that higher mass ratios lead to somewhat more electron acceleration in shocks. This write-up is in three parts, the first discussing a description of what I was up to personally throughout the fall, the second discussing the results of the collisional shock study, and the third discussing the results of the high mass ratio shock study.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An efficient, conservative, time-implicit solver for the fully kinetic arbitrary-species 1D-2V Vlasov-Ampère system

In this paper, we consider the solution of the fully kinetic (including electrons) Vlasov-Ampère system in a one-dimensional physical space and two-dimensional velocity space (1D-2V) for an arbitrary number of species with a time-implicit Eulerian algorithm. The problem of velocity-space meshing for disparate thermal and bulk velocities is dealt with by an adaptive coordinate transformation of the Vlasov equation for each species, which is then discretized, including the resulting inertial terms. Mass, momentum, and energy are conserved, and Gauss's law is enforced to within the nonlinear convergence tolerance of the iterative solver through a set of nonlinear constraint functions while permitting significant flexibility in choosing discretizations in time, configuration, and velocity space. We mitigate the temporal stiffness introduced by, e.g., the plasma frequency through the use of high-order/low-order (HOLO) acceleration of the iterative implicit solver. We present several numerical results for canonical problems of varying degrees of complexity, including the multiscale ion-acoustic shock wave problem, which demonstrate the efficacy, accuracy, and efficiency of the scheme.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simulations of Precessing Jets and the Formation of X-shaped Radio Galaxies

Abstract Jet precession is sometimes invoked to explain asymmetries in radio galaxy (RG) jets and “X/S/Z-shaped” RGs, caused by the presence of a binary black hole companion to the source active galactic nucleus or by accretion instabilities. We present a series of simulations of RG jet precession to examine how these sources would evolve over time, including a passive distribution of cosmic-ray electrons so we can model radio synchrotron emissions and create synthetic radio maps of the sources. We find that a single source viewed from different angles can result in differing RG morphological classifications, confusing physical implications of these classifications. Additionally, the jet trajectories can become unstable due to their own self-interactions and lead to “reorientation events” that may look like the effects of external dynamics such as shocks, winds, or cold fronts in the medium. Finally, something akin to an “Odd Radio Circle” may be observed in the case of viewing the radio remnant of such a precessing source from a line of sight near the precession axis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

On the effects of reactant stratification and wall curvature in non-premixed rotating detonation combustors

The optimization of non-premixed rotating detonation combustors (RDCs) requires improved understanding of the coupled effects of reactant stratification, fluid property gradients, and complex shock-wave interactions on the detonation wave structure within annular geometries. In the current work, simultaneous orthogonal views of chemiluminescence and hydroxyl planar laser-induced fluorescence (PLIF) are utilized to establish the existence of a dual-wave system characterized by leading and trailing detonation waves that are closely coupled by the local flow physics. These features are persistent over a wide range of mass flow rates and are consistent with prior observations of non-premixed rotating detonations in annular geometries. The detailed instantaneous time sequences are compared with a 3D reactive unsteady Reynolds averaged Navier-Stokes (URANS) simulation to more clearly elucidate the in-situ combustion dynamics and the sensitivity to reactant inlet conditions. It is found that the dual-wave system results from unburned reactants that survive the leading detonation wave in the injector near field and are consumed within a trailing azimuthal reflected-shock combustion (ARSC) zone. By contrast, the injector far field is characterized by rapid mixing due to a sudden drop to subsonic conditions, and the bifurcated detonation wave structure collapses into a stronger, single-wave detonation front with higher overall pressure ratio as compared with the dual-wave system. While each RDC will have different inflow, mixing, and combustion characteristics, the underlying interactions between the stratified reactants and azimuthal wave dynamics identified through the combination of advanced MHz-rate diagnostics and 3D numerical simulations have important implications for the study of detonation wave stability, mode transition, and combustion efficiency in non-premixed annular RDCs.

3D URANS↗
DOE OSTI2021

A multishock model for the density variance of anisotropic, highly magnetized, supersonic turbulence

ABSTRACT Shocks form the basis of our understanding for the density and velocity statistics of supersonic turbulent flows, such as those found in the cool interstellar medium (ISM). The variance of the density field, $\sigma ^2_{\rho /\rho _0}$, is of particular interest for molecular clouds (MCs), the birthplaces of stars in the Universe. The density variance may be used to infer underlying physical processes in an MC, and parametrizes the star formation (SF) rate of a cloud. However, models for $\sigma ^2_{\rho /\rho _0}$ all share a common feature – the variance is assumed to be isotropic. This assumption does not hold when a trans-/sub-Alfvénic mean magnetic field, ${B}_0$, is present in the cloud, which observations suggest is relevant for some MCs. We develop an anisotropic model for $\sigma _{\rho /\rho _0}^2$, using contributions from hydrodynamical and fast magnetosonic shocks that propagate orthogonal to each other. Our model predicts an upper bound for $\sigma _{\rho /\rho _0}^2$ in the high Mach number $(\mathcal {M})$ limit as small-scale density fluctuations become suppressed by the strong ${B}_0$. The model reduces to the isotropic $\sigma _{\rho /\rho _0}^2\!-\!\mathcal {M}$ relation in the hydrodynamical limit. To validate our model, we calculate $\sigma _{\rho /\rho _0}^2$ from 12 high-resolution, three-dimensional, supersonic, sub-Alfvénic magnetohydrodynamical (MHD) turbulence simulations and find good agreement with our theory. We discuss how the two MHD shocks may be the bimodally oriented overdensities observed in some MCs and the implications for SF theory in the presence of a sub-Alfvénic ${B}_0$. By creating an anisotropic, supersonic density fluctuation model, this study paves the way for SF theory in the highly anisotropic regime of interstellar turbulence.

Beattie, James R.↗
DOE OSTI2024

Imaging the dynamics of initial laser-driven shocks and blowoff plasmas in polystyrene under laser-direct-drive fusion conditions

The dynamics of laser-driven shock propagation in solid ablator material, along with the concomitant processes of shinethrough and blowoff plasma plume expansion, at conditions relevant to laser-direct-drive fusion are measured using a tabletop system capable of supporting rapid dataset development with significantly higher resolution compared to existing methods. The experimental results are directly compared with current radiation-hydrodynamic simulations using the two-dimensional code . Discrepancies between experiments and simulations are evidenced and attributed to limitation of simulation to account for important processes involved in the solid-to-plasma transition. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Experimental Observations of Laser-Driven Tin Ejecta Microjet Interactions

The study of high-velocity particle-laden flow interactions is of importance for the understanding of a wide range of natural phenomena, ranging from planetary formation to cloud interactions. Experimental observations of particle dynamics are sparse given the difficulty of generating high-velocity flows of many particles. Ejecta microjets are micron-scale jets formed by strong shocks interacting with imprinted surfaces to generate particle plumes traveling at several kilometers per second. As such, the interaction of two ejecta microjets provides a novel experimental methodology to study interacting particle streams. In this Letter, we report the first time sequences of x-ray radiography images of two interacting tin ejecta microjets taken on a platform designed for the OMEGA Extended Performance (OMEGA EP) laser. We observe that the microjets pass through each other unattenuated for the case of 11.7±3.2 GPa shock pressures and jet velocities of 2.2±0.5 km/s but show strong interaction dynamics for 116.0±6.1 GPa shock pressures and jet velocities of 6.5±0.5 km/s. Furthermore, we find that radiation-hydrodynamic simulations of the experiments are able to capture many aspects of the collisional behavior, such as the attenuation of jet velocity in the direction of propagation, but are unable to match the full spread of the strongly interacting cloud.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗