Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Validating methods for modeling composition gradients in planar shock experiments

An interface is Rayleigh–Taylor (RT) unstable when acceleration pushes a less dense material into a more dense one, and the growth of the instability is governed partly by the Atwood number gradient. Double-shell inertial confinement fusion capsules have a foam spacer layer pushing on an inner capsule composed of a beryllium tamper and high-Z inner shell, and so have RT unstable interfaces that require benchmarking. To this end, the results of a planar shock experiment with beryllium/tungsten targets are presented. One target had the normal bilayer construction of beryllium and tungsten in two distinct layers; the second target had the beryllium grading into tungsten with a quasi-exponential profile, motivated by the potential for reduced RT growth with the gradient profile. Simulations mimic the shock profiles for both targets and match the shock velocity to within 5%. These results validate the ability of our simulations to model double-shell capsules with bilayer or graded layer Be/W inner shells, which are needed to design future experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

The Bastille Day Magnetic Clouds and Upstream Shocks: Near Earth Interplanetary Observations

The energetic charged particle, interplanetary magnetic field, and plasma characteristics of the 'Bastille Day' shock and ejecta/magnetic cloud events at 1 AU occurring over the days 14-16 July 2000 are described. Profiles of MeV (WIND/LEMT) energetic ions help to organize the overall sequence of events from the solar source to 1 AU. Stressed are analyses of an outstanding magnetic cloud (MC2) starting late on 15 July and its upstream shock about 4 hours earlier in WIND magnetic field and plasma data. Also analyzed is a less certain, but likely, magnetic cloud (MC1) occurring early on 15 July; this was separated from MC2 by its upstream shock and many heliospheric current sheet (HCS) crossings. Other HCS crossings occurred throughout the 3-day period. Overall this dramatic series of interplanetary events caused a large multi-phase magnetic storm with min Dst lower than -300 nT. The very fast solar wind speed (greater than or equal to 1100 km/s) in and around the front of MC2 (for near average densities) was responsible for a very high solar wind ram pressure driving in the front of the magnetosphere to geocentric distances estimated to be as low as approx. 5 R(sub E), much lower than the geosynchronous orbit radius. This was consistent with magnetic field observations from two GOES satellites which indicated they were in the magnetosheath for extended times. A static force free field model is used to fit the two magnetic cloud profiles providing estimates of the clouds' physical and geometrical properties. MC2 was much larger than MCI, but their axes were nearly antiparallel, and their magnetic fields had the same left-handed helicity. MC2's axis and its upstream shock normal were very close to being perpendicular to each other, as might be expected if the cloud were driving the shock at the time of observation. The estimated axial magnetic flux carried by MC2 was 52 x 10(exp 20) Mx, which is about 5 times the typical magnetic flux estimated for other magnetic clouds in the WIND data over its first 4 years and is 17 times the flux of MC1. This large flux is due to both the strong axially-directed field of MC2 (46.8 nT on the axis) and the large radius (R(sub 0) = 0.189 AU) of the flux tube. MC2's average speed is consistent with the expected transit time from a halo-CME to which it is apparently related.

Lepping, R. P.↗
DOE OSTI2021

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↗
DOE OSTI2023

Isolation and phase-space energization analysis of the instabilities in collisionless shocks

We analyse the generation of kinetic instabilities and their effect on the energization of ions in non-relativistic, oblique collisionless shocks using a 3D-3V (three spatial with three velocity components) simulation by dHybridR, a hybrid particle-in-cell code. At sufficiently high Mach number, quasi-perpendicular and oblique shocks can experience rippling of the shock surface caused by kinetic instabilities arising from free energy in the ion velocity distribution due to the combination of the incoming ion beam and the population of ions reflected at the shock front. To understand the role of the ripple on particle energization, we devise a new instability isolation method to identify the unstable modes underlying the ripple and interpret the results in terms of the governing kinetic instability. We generate velocity-space signatures using the field–particle correlation technique to look at energy transfer in phase space from the isolated instability driving the shock ripple, providing a viewpoint on the different dynamics of distinct populations of ions in phase space. Together, the field–particle correlation technique and our new instability isolation method provide a unique viewpoint on the different dynamics of distinct populations of ions in phase space and allow us to completely characterize the energetics of the collisionless shock under investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Thermo-Mechanical Shock, X-Ray Transport and Radiation Response Analysis Experiments on NIF

The study of x-ray induced rapid heating and the resulting material responses - (thermomechanical shock (TMS) and thermo-structural response (TSR)) - is important to our understanding of material behavior in extreme environments. A series of high conversion efficiency (XRCE) x-ray sources have been developed for the NIF laser facility ranging in fluence and spectral content. In conjunction with source development, we have developed a new XTRRA (X-ray Transport and Radiation Response Analysis) test cassette to hold a set of samples, each at equal distance from target chamber center (TCC) where the x-ray source is placed. Design and performance details of the x-ray sources is discussed elsewhere. Here, we describe the development and use of this TMS XTRRA test cassette. The XTRRA test cassette is a snout mounted on a NIF Diagnostic Instrument Manipulator (DIM) designed to position six samples equidistant from the x-ray source to assure uniform irradiation of all samples. Photon Displacement interferometry (PDI) is used to measure the displacement of the rear surface of each sample resulting from x-ray exposure on the front of the sample. This will help determine the x-ray generated impulse on the sample. The samples are retrieved from the snout after the shot for inspection and post-shot analysis including the measurement of sample ablated mass and engineering tests examining changes in the material structure. The data from experiments are compared to predictions from LLNL’s ASC codes and used to inform models implemented in these codes. Thus, XRCE x-ray sources with the XTRRA snout at the National Ignition Facility (NIF) laser at the Lawrence Livermore National Laboratory provide the ability to expose a broad range of material samples directly to high x-ray fluence, measure shock transits, and recover the sample for further analysis post-shot.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Plasma processes in space

Elementary microscopic interactions in plasmas are described. The importance of plasma physics in space studies is illustrated by examining several phenomena which cannot be explained satisfactorily by MHD theory. These include kinetic instabilities, plasma turbulence in the bow shock, magnetic turbulence near the moon, VLF emissions in the magnetosphere, planetary and solar radio emissions, and interaction of planetary and cometary plasmas with the solar wind. Evidence for the existence of anomalous transport processes in terrestrial and planetary magnetospheres is presented.

Wu, C. S.↗
NASA NTRS1986

Nineteenth International Cosmic Ray Conference. Conference Papers: Invited Rapporteur, Highlight, Miscellaneous, Volume 9

Invited talks, rapporteur talks, and highlight talks are included. Topics of the invited and highlight talks include astrophysical jets, gamma-ray line astronomy, cosmic rays and gamma rays in astrophysics, the early universe, elementary particle physics, solar flares and acceleration of energetic particles, cosmogenic nuclei, extragalactic astronomy, composition of solar flare particles, very high energy gamma ray sources, gamma-ray bursts, shock acceleration in the solar wind, cosmic rays in deep underground detectors, spectrum of cosmic rays at 10 to the 19th power eV, and nucleus-nucleus interactions.

Jones, F. C.↗
NASA NTRS2010

Jet Noise Reduction by Microjets - A Parametric Study

The effect of injecting tiny secondary jets (microjets ) on the radiated noise from a subsonic primary jet is studied experimentally. The microjets are injected on to the primary jet near the nozzle exit with variable port geometry, working fluid and driving pressure. A clear noise reduction is observed that improves with increasing jet pressure. It is found that smaller diameter ports with higher driving pressure, but involving less thrust and mass fraction, can produce better noise reduction. A collection of data from the present as well as past experiments is examined in an attempt to correlate the noise reduction with the operating parameters. The results indicate that turbulent mixing noise reduction, as monitored by OASPL at a shallow angle, correlates with the ratio of jet to primary jet driving pressures normalized by the ratio of corresponding diameters (p d /pjD). With gaseous injection, the spectral amplitudes decrease at lower frequencies while an increase is noted at higher frequencies. It is apparent that this amplitude crossover is at least partly due to shock-associated noise from the underexpanded jets themselves. Such crossover is not seen with water injection since the flow in that case is incompressible and there is no shock-associated noise. Centerline velocity data show that larger noise reduction is accompanied by faster jet decay as well as significant reduction in turbulence intensities. While a physical understanding of the dependence of noise reduction on p d /pjD remains unclear, given this correlation, an analysis explains the observed dependence of the effect on various other parameters.

Zaman, K. B. M. Q.↗
NASA NTRS2011

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

Baurle, Robert A.↗
DOE OSTI2025

Investigation of warm dense matter by combining static and dynamic compression

Shock-compressing material samples that are precompressed with diamond-anvil cells (DACs) enables benchmarking equation of state measurements over a wide range of high-pressure conditions, in particular, warm dense matter relevant to planetary science, astrophysics, and inertial confinement fusion. These targets provide optical access to the shocked sample, so optical velocity interferometry and pyrometry are typical diagnostics. We present a scaling for the coupling of the laser drive to the shock pressure in DACs based on over 150 prior experiments that can be used to design future experiments. For initial static pressures ≤ 5 GPa, using a sapphire window mitigates the effect of photoionization compared to diamond windows and allows the measurement of shock pressures up to 1 TPa with optical diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Sound speed measurements in silver shock compressed to 300 GPa: Solid-state transition, melting, and liquid-state response

To gain insight into the thermodynamic response of shock compressed Ag states corresponding to the face-centered-cubic (fcc) to body-centered-cubic (bcc) transformation and melting observed using in-situ x-ray diffraction (XRD) measurements, we present longitudinal sound speed results and their analysis in Ag at peak stresses ranging from 60 to 300 GPa. The measured sound speeds increased linearly with density compression to 171 GPa, showing that the sound speeds (and longitudinal moduli) in the fcc and bcc phases are very similar. Between 171 and 218 GPa, the sound speed dropped significantly, consistent with the melting reported using XRD measurements in shock compressed Ag. From 218 to 300 GPa, the increasing sound speeds and Hugoniot states provide a determination of the liquid phase Ag response. In particular, determination of the Grüneisen parameter (Γ) showed that the density – Grüneisen parameter product (ρΓ) for liquid Ag is constant, but differs significantly from that for solid Ag at ambient conditions. Furthermore, the Mie- Grüneisen equation of state can be used to describe the Hugoniot and off-Hugoniot response of liquid Ag.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

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↗
DOE OSTI2025

Physics-based stabilized finite element approximations of the Poisson–Nernst–Planck equations

We present and analyze two stabilized finite element methods for solving numerically the Poisson–Nernst–Planck equations. The stabilization we consider is carried out by using a shock detector and a discrete graph Laplacian operator for the ion equations, whereas the discrete equation for the electric potential need not be stabilized. Discrete solutions stemmed from the first algorithm preserve both maximum and minimum discrete principles. For the second algorithm, its discrete solutions are conceived so that they hold discrete principles and obey an entropy law provided that an acuteness condition is imposed for meshes. Remarkably the latter is found to be unconditionally stable. We validate our methodology through transient numerical experiments that show convergence toward steady-state solutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coupling neutrino oscillations and simulations of core-collapse supernovae

At the present time even the most sophisticated, multidimensional simulations of core-collapse supernovae do not (self-consistently) include neutrino flavor transformation. This physics is missing despite the importance of neutrinos in the core-collapse explosion paradigm. Because of this dependence, any flavor transformation that occurs in the region between the protoneutron star and the shock could result in major effects upon the dynamics of the explosion. Here, we present the first hydrodynamic core-collapse supernova simulation which simultaneously includes flavor transformation of the free-streaming neutrinos in the neutrino transport. These oscillation calculations are dynamically updated and evolve self-consistently alongside the hydrodynamics. Using a M = 20 M ⊙ progenitor, we find that while the oscillations have an effect on the neutrino emission and the heating rates, flavor transformation alone does not lead to a successful explosion of this progenitor in spherical symmetry.

79 ASTRONOMY AND ASTROPHYSICS↗