Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SHOCK”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Numerical investigation on target implosions driven by radiation ablation and shock compression in dynamic hohlraums

In a dynamic hohlraum driven inertial confinement fusion (ICF) configuration, the target may experience two different kinds of implosions. One is driven by hohlraum radiation ablation, which is approximately symmetric at the equator and poles. The second is caused by the radiating shock produced in Z-pinch dynamic hohlraums, only taking place at the equator. To gain a symmetrical target implosion driven by radiation ablation and avoid asymmetric shock compression is a crucial issue in driving ICF using dynamic hohlraums. It is known that when the target is heated by hohlraum radiation, the ablated plasma will expand outward. The pressure in the shocked converter plasma qualitatively varies linearly with the material temperature. However, the ablation pressure in the ablated plasma varies with 3.5 power of the hohlraum radiation temperature. Therefore, as the hohlraum temperature increases, the ablation pressure will eventually exceed the shock pressure, and the expansion of the ablated plasma will obviously weaken the shock propagation and decrease its velocity after propagating into the ablator plasma. Consequently, longer time duration is provided for the symmetrical target implosion driven by radiation ablation. In this paper these processes are numerically investigated by changing drive currents or varying load parameters. The simulation results show that a critical hohlraum radiation temperature is needed to provide a high enough ablation pressure to decelerate the shock, thus providing long enough time duration for the symmetric fuel compression driven by radiation ablation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effect of initial perturbation amplitude on Richtmyer-Meshkov flows induced by strong shocks

We systematically study the effect of the initial perturbation on Richtmyer-Meshkov (RM) flows induced by strong shocks in fluids with contrasting densities. Smooth Particle Hydrodynamics simulations are employed. A broad range of shock strengths and density ratios is considered. The amplitude of the initial single mode sinusoidal perturbation of the interface varies from 0% to 100% of its wavelength. The simulations results are compared, wherever possible, with four rigorous theories, and with other experiments and simulations, achieving good quantitative and qualitative agreement. Our study is focused on early time dynamics of the Richtmyer-Meshkov instability (RMI). We analyze the initial growth-rate of RMI immediately after the shock passage, when the perturbation amplitude increases linearly with time. For the first time, to the authors' knowledge, we find that the initial growth-rate of RMI is a non-monotone function of the initial perturbation amplitude, thus restraining the amount of energy that can be deposited by the shock at the interface. The maximum value of the initial growth-rate depends on the shock strength and the density ratio, whereas the corresponding value of the initial perturbation amplitude depends only slightly on the shock strength and density ratio.

42 ENGINEERING↗
DOE OSTI2015

Numerical study of heterogeneous mean temperature and shock wave in a resonator

When a frequency of gas oscillation in an acoustic resonator is sufficiently close to one of resonant frequencies of the resonator, the amplitude of gas oscillation becomes large and hence the nonlinear effect manifests itself. Then, if the dissipation effects due to viscosity and thermal conductivity of the gas are sufficiently small, the gas oscillation may evolve into the acoustic shock wave, in the so-called consonant resonators. At the shock front, the kinetic energy of gas oscillation is converted into heat by the dissipation process inside the shock layer, and therefore the temperature of the gas in the resonator rises. Since the acoustic shock wave travels in the resonator repeatedly over and over again, the temperature rise becomes noticeable in due course of time even if the shock wave is weak. We numerically study the gas oscillation with shock wave in a resonator of square cross section by solving the initial and boundary value problem of the system of three-dimensional Navier-Stokes equations with a finite difference method. In this case, the heat conduction across the boundary layer on the wall of resonator causes a spatially heterogeneous distribution of mean (time-averaged) gas temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sound speed measurements in lithium fluoride single crystals shock compressed to 168 GPa along [100]

The shock wave response of [100] lithium fluoride (LiF) single crystals at high stresses is of long-standing interest due to their extensive use as optical windows in dynamic compression experiments. The report of melting in shock compressed LiF single crystals between 134 GPa and 152 GPa – based on a single sound speed datum [J. Appl. Phys. 117, 045901 (2015)] – was surprising because good optical transmission was previously demonstrated in LiF shock compressed to ~200 GPa [J. Appl. Phys. 116, 033515 (2014)]. To address these apparent differences, we report on plate impact experiments on [100] LiF single crystals shock compressed to 168 GPa. Wave profiles were measured using laser interferometry to determine Hugoniot states and longitudinal sound speeds in shock compressed LiF. The measured Hugoniot states are in good agreement with those measured in previous studies. However, the measured sound speeds presented here show no evidence of melting up to 168 GPa. In particular, the abrupt drop reported previously in sound speed at 152 GPa was not observed in the present work. Furthermore, our results establish a lower bound of 168 GPa for the onset of melting in shock compressed LiF single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Experimental observations of detached bow shock formation in the interaction of a laser-produced plasma with a magnetized obstacle

The magnetic field produced by planets with active dynamos, like the Earth, can exert sufficient pressure to oppose supersonic stellar wind plasmas, leading to the formation of a standing bow shock upstream of the magnetopause, or pressure-balance surface. Scaled laboratory experiments studying the interaction of an inflowing solar wind analog with a strong, external magnetic field are a promising new way to study magnetospheric physics and to complement existing models, although reaching regimes favorable for magnetized shock formation is experimentally challenging. This paper presents experimental evidence of the formation of a magnetized bow shock in the interaction of a supersonic, super-Alfvenic plasma with a strongly magnetized obstacle at the OMEGA laser facility. The solar wind analog is generated by the collision and subsequent expansion of two counterpropagating, laser-driven plasma plumes. The magnetized obstacle is a thin wire, driven with strong electrical currents. Hydrodynamic simulations using the FLASH code predict that the colliding plasma source meets the criteria for bow shock formation. Spatially resolved, optical Thomson scattering measures the electron number density, and optical emission lines provide a measurement of the plasma temperature, from which we infer the presence of a fast magnetosonic shock far upstream of the obstacle. Proton images provide a measure of large-scale features in the magnetic field topology, and reconstructed path-integrated magnetic field maps from these images suggest the formation of a bow shock upstream of the wire and as a transient magnetopause. Finally, we compare features in the reconstructed fields to two-dimensional MHD simulations of the system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Strong reconnection electric fields in shock-driven turbulence

Turbulent magnetic reconnection in a quasi-parallel shock under parameters relevant to the Earth's bow shock is investigated by means of a two-dimensional particle-in-cell simulation. The addressed aspects include the reconnection electric field, the reconnection rate, and the electron and the ion outflow speeds. In the shock transition region, many current sheets are generated in shock-driven turbulence, and electron-only reconnection and reconnection where both ions and electrons are involved can occur in those current sheets. The electron outflow speed in electron-only reconnection shows a positive correlation with the theoretical speed, which is close to the local electron Alfvén speed, and a strong convection electric field is generated by the large electron outflow. As a result, the reconnection electric field becomes much larger than those in the standard magnetopause or magnetotail reconnection. In shock-driven reconnection that involves ion dynamics, both electron outflows and ion outflows can reach of the order of 10 times the Alfvén speed in the X-line rest frame, leading to a reconnection electric field the same order as that in electron-only reconnection. An electron-only reconnection event observed by the magnetospheric multiscale mission downstream of a quasi-parallel shock is qualitatively similar to those in the simulation and shows that the outflow speed reaches approximately half the local electron Alfvén speed, supporting the simulation prediction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An investigation of shock-induced phase transition in soda-lime glass

There exists a large body of evidence from experiments and molecular dynamics simulations to suggest the occurrence of phase transitions in soda-lime glass (SLG) and other silica glasses subject to shock compression to pressures above 3 GPa. In light of these findings, the current work investigated the existence of phase transition in SLG using shock and release experiments. The experiments employed symmetric SLG–SLG impact to achieve complete unloading to zero stress after shock compression to stresses in the range of 3–7 GPa. The stress–strain response and the Lagrangian release wave speed behavior of SLG obtained from these experiments are seen to reveal a mismatch between the loading and unloading paths of the pressure–strain curve for the material, which serves as compelling evidence for the occurrence of a shock-induced phase transition in the material at relatively low pressures. Furthermore, the release wave speed vs strain data obtained from experiments were used to construct a methodology for modeling the shock and release behavior of SLG. Lastly, this scheme implemented in numerical simulations was able to capture the release behavior of shock compressed SLG, for which a robust and satisfactory model was previously unavailable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

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

Shock wave bending around a dusty plasma void

We report on experimental observations of the bending of a dust acoustic shock wave around a dust void region. This phenomenon occurs as a planar shock wavefront encounters a compressible obstacle in the form of a void whose size is larger than the wavelength of the wave. As they collide, the central portion of the wavefront, that is the first to touch the void, is blocked while the rest of the front continues to propagate, resulting in an inward bending of the shock wave. The bent shock wave eventually collapses, leading to the transient trapping of dust particles in the void. Subsequently, a Coulomb explosion of the trapped particles generates a bow shock. The experiments have been carried out in a Direct Current (DC) glow discharge plasma, where the shock wave and the void are simultaneously created as self-excited modes of a three-dimensional dust cloud. The salient features of this phenomenon are reproduced in molecular dynamics simulations, which provide valuable insights into the underlying dynamics of this interaction.

Dust acoustic waves↗
DOE OSTI2011

INTERPRETING THE PROPERTIES OF SOLAR ENERGETIC PARTICLE EVENTS BY USING COMBINED IMAGING AND MODELING OF INTERPLANETARY SHOCKS

Images of the solar corona obtained by the Solar-Terrestrial Relations Observatory (STEREO) provide high-cadence, high-resolution observations of a compression wave forming ahead of a fast (940 km s{sup -1}) coronal mass ejection (CME) that erupted at {approx}9:00 UT on 2010 April 03. The passage of this wave at 1 AU is detected in situ by the Advanced Composition Explorer and Wind spacecraft at 08:00 UT on April 05 as a shock followed by a turbulent and heated sheath. These unprecedented and complementary observations of a shock-sheath region from the Sun to 1 AU are used to investigate the onset of a Solar Energetic Particle (SEP) event measured at the first Lagrange point (L1) and at STEREO-Behind (STB). The spatial extent, radial coordinates, and speed of the ejection are measured from STEREO observations and used as inputs to a numerical simulation of the CME propagation in the background solar wind. The simulated magnetic and plasma properties of the shock and sheath region at L1 agree very well with the in situ measurements. These simulation results reveal that L1 and STB are magnetically connected to the western and eastern edges of the driven shock, respectively. They also show that the 12 hr delay between the eruption time of the ejection and the SEP onset at L1 corresponds to the time required for the bow shock to reach the magnetic field lines connected with L1. The simulated shock compression ratio increases along these magnetic field lines until the maximum flux of high-energy particles is observed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

LONG-WAVELENGTH UNSTABLE MODES IN THE FAR UPSTREAM OF RELATIVISTIC COLLISIONLESS SHOCKS

The growth rate of long-wavelength kinetic instabilities arising due to the interaction of a collimated beam of relativistic particles and a cold unmagnetized plasma are calculated in the ultrarelativistic limit. For sufficiently culminated beams, all long-wavelength modes are shown to be Weibel-unstable, and a simple analytic expression for their growth rate is derived. For large transverse velocity spreads, these modes become stable. An analytic condition for stability is given. These analytic results, which generalize earlier ones given in the literature, are shown to be in agreement with numerical solutions of the dispersion equation and with the results of novel particle in cell simulations in which the electromagnetic fields are restricted to a given k-mode. The results may describe the interaction of energetic cosmic rays, propagating into the far upstream of a relativistic collisionless shock, with a cold unmagnetized upstream. The long-wavelength modes considered may be efficient in deflecting particles and could be important for diffusive shock acceleration. It is shown that while these modes grow in relativistic shocks propagating into electron-positron pair plasmas, they are damped in relativistic shocks propagating into electron-proton plasmas with moderate Lorenz factors {Gamma}{sub sh} {approx}< (m{sub p} /m{sub e}){sup 1/2}. If these modes dominate the deflection of energetic cosmic rays in electron-positron shocks, it is argued that particle acceleration is suppressed at shock frame energies that are larger than the downstream thermal energy by a factor of {approx}> {Gamma}{sub sh}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

MOLECULAR TRACERS OF TURBULENT SHOCKS IN GIANT MOLECULAR CLOUDS

Giant molecular clouds contain supersonic turbulence and simulations of magnetohydrodynamic turbulence show that these supersonic motions decay in roughly a crossing time, which is less than the estimated lifetimes of molecular clouds. Such a situation requires a significant release of energy. We run models of C-type shocks propagating into gas with densities around 10{sup 3} cm{sup -3} at velocities of a few km s{sup -1}, appropriate for the ambient conditions inside of a molecular cloud, to determine which species and transitions dominate the cooling and radiative energy release associated with shock cooling of turbulent molecular clouds. We find that these shocks dissipate their energy primarily through CO rotational transitions and by compressing pre-existing magnetic fields. We present model spectra for these shocks, and by combining these models with estimates for the rate of turbulent energy dissipation, we show that shock emission should dominate over emission from unshocked gas for mid to high rotational transitions (J > 5) of CO. We also find that the turbulent energy dissipation rate is roughly equivalent to the cosmic-ray heating rate and that the ambipolar diffusion heating rate may be significant, especially in shocked gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

ELECTRON ACCELERATIONS AT HIGH MACH NUMBER SHOCKS: TWO-DIMENSIONAL PARTICLE-IN-CELL SIMULATIONS IN VARIOUS PARAMETER REGIMES

Electron accelerations at high Mach number collisionless shocks are investigated by means of two-dimensional electromagnetic particle-in-cell simulations with various Alfven Mach numbers, ion-to-electron mass ratios, and the upstream electron {beta}{sub e} (the ratio of the thermal pressure to the magnetic pressure). We find electrons are effectively accelerated at a super-high Mach number shock (M{sub A} {approx} 30) with a mass ratio of M/m = 100 and {beta}{sub e} = 0.5. The electron shock surfing acceleration is an effective mechanism for accelerating the particles toward the relativistic regime even in two dimensions with a large mass ratio. Buneman instability excited at the leading edge of the foot in the super-high Mach number shock results in a coherent electrostatic potential structure. While multi-dimensionality allows the electrons to escape from the trapping region, they can interact with the strong electrostatic field several times. Simulation runs in various parameter regimes indicate that the electron shock surfing acceleration is an effective mechanism for producing relativistic particles in extremely high Mach number shocks in supernova remnants, provided that the upstream electron temperature is reasonably low.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

QUENCHED COLD ACCRETION OF A LARGE-SCALE METAL-POOR FILAMENT DUE TO VIRIAL SHOCKING IN THE HALO OF A MASSIVE z = 0.7 GALAXY

Using HST/COS/STIS and HIRES/Keck high-resolution spectra, we have studied a remarkable H I absorbing complex at z = 0.672 toward the quasar Q1317+277. The H I absorption has a velocity spread of {Delta}v = 1600 km s{sup -1}, comprises 21 Voigt profile components, and resides at an impact parameter of D = 58 kpc from a bright, high-mass (log M {sub vir}/M {sub Sun} {approx_equal} 13.7) elliptical galaxy that is deduced to have a 6 Gyr old, solar metallicity stellar population. Ionization models suggest the majority of the structure is cold gas surrounding a shock-heated cloud that is kinematically adjacent to a multi-phase group of clouds with detected C III, C IV, and O VI absorption, suggestive of a conductive interface near the shock. The deduced metallicities are consistent with the moderate in situ enrichment relative to the levels observed in the z {approx} 3 Ly{alpha} forest. We interpret the H I complex as a metal-poor filamentary structure being shock heated as it accretes into the halo of the galaxy. The data support the scenario of an early formation period (z > 4) in which the galaxy was presumably fed by cold-mode gas accretion that was later quenched via virial shocking by the hot halo such that, by intermediate redshift, the cold filamentary accreting gas is continuing to be disrupted by shock heating. Thus, continued filamentary accretion is being mixed into the hot halo, indicating that the star formation of the galaxy will likely remain quenched. To date, the galaxy and the H I absorption complex provide some of the most compelling observational data supporting the theoretical picture in which accretion is virial shocked in the hot coronal halos of high-mass galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

PHOSPHORUS CHEMISTRY IN THE SHOCKED REGION L1157 B1

We study the evolution of phosphorus-bearing species in one-dimensional C-shock models. We find that the abundances of P-bearing species depend sensitively on the elemental abundance of P in the gas phase and on the abundance of N atoms in the pre-shock gas. The observed abundance of PN and the non-detection of PO toward L1157 B1 are reproduced in C-shock models with shock velocity v = 20 km s{sup -1} and pre-shock density n(H{sub 2}) 10{sup 4}-10{sup 5} cm{sup -3}, if the elemental abundance of P in the gas phase is {approx}10{sup -9} and the N-atom abundance is n(N)/n{sub H} {approx}10{sup -5} in the pre-shock gas. We also find that P-chemistry is sensitive to O- and N-chemistry because N atoms are destroyed mainly by OH and NO. We identify the reactions of O-bearing and N-bearing species that significantly affect P-chemistry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

COLLISIONLESS SHOCKS IN A PARTIALLY IONIZED MEDIUM. III. EFFICIENT COSMIC RAY ACCELERATION

In this paper, we present the first formulation of the theory of nonlinear particle acceleration in collisionless shocks in the presence of neutral hydrogen in the acceleration region. The dynamical reaction of the accelerated particles, the magnetic field amplification, and the magnetic dynamical effects on the shock are also included. The main new aspect of this study, however, consists of accounting for charge exchange and the ionization of a neutral hydrogen, which profoundly change the structure of the shock, as discussed in our previous work. This important dynamical effect of neutrals is mainly associated with the so-called neutral return flux, namely the return of hot neutrals from the downstream region to upstream, where they deposit energy and momentum through charge exchange and ionization. We also present the self-consistent calculation of Balmer line emission from the shock region and discuss how to use measurements of the anomalous width of the different components of the Balmer line to infer cosmic ray acceleration efficiency in supernova remnants showing Balmer emission: the broad Balmer line, which is due to charge exchange of hydrogen atoms with hot ions downstream of the shock, is shown to become narrower as a result of the energy drainage into cosmic rays, while the narrow Balmer line, due to charge exchange in the cosmic-ray-induced precursor, is shown to become broader. In addition to these two well-known components, the neutral return flux leads to the formation of a third component with an intermediate width: this too contains information on ongoing processes at the shock.

79 ASTRONOMY AND ASTROPHYSICS↗