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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

The central engine of quasars and active galactic nuclei Hadronic interactions of shock-accelerated relativistic protons

The detailed dynamics are presented of a model for AGNs and QSOs based on spherically symmetric accretion onto a black hole and dissipation of the accretion energy in a collisionless shock. The basic assumptions and equations relevant for the model are presented, and the relations between the physical parameters are given in terms of the shock parameters. The latter are calculated for varying Mach numbers, and the resulting scaling laws are discussed and compared to the data, in particular the recently discovered correlations between luminosity and mass. The model can naturally produce ultrarelativistic electrons with flat spectra, a feature which, when coupled with e(+)-e(-) pair production, can account for the observed continuum radiation from these objects over essentially the whole electromagnetic spectrum. The model also predicts the emission of high-energy neutrinos with a luminosity roughly equal to that of the photons.

Kazanas, D.↗
DOE OSTI2023

Drop Interactions with the Conical Shock Structure Generated by a Mach 4.5 Projectile

This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of $192$ $\mu$m. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier–Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant $C$ 2 = $2$ compared to the accepted value $C$ 2 = $2/3$. In conclusion, results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
NASA NTRS1973

Recent developments in shock tube research; Proceedings of the Ninth International Symposium, Stanford University, Stanford, Calif., July 16-19, 1973

Recent advances in shock tube research are described in papers dealing with the design and performance features of new devices as well as applications in aerodynamic, chemical, and physics experiments. Topics considered include a cryogenic shock tube for studying liquid helium fluid mechanics, studies of shock focusing and nonlinear resonance in shock tubes, applications in gas laser studies, very-low and very-high temperature chemical kinetic measurements, shock tube studies of ionization and recombination phenomena, applications in bioacoustic research, shock-tube simulation studies of sonic booms, and plasma research. Individual items are announced in this issue.

Bershader, D.↗
DOE OSTI2021

Geometry and dynamics of passive scalar structures in compressible turbulent mixing

We present a structure-based numerical analysis of passive scalar mixing in decaying homogeneous isotropic turbulence (DHIT) and shock-turbulence interaction canonical configurations. The analysis focuses on the temporal evolution of ensembles of passive scalar structures, initialized as spheres of different sizes relative to the Taylor microscale. An algorithm is introduced to track the evolution of each individual structure and the interactions with other structures in the ensemble, relating changes in the surface geometry and the underlying physical processes (turbulent transport, scalar dissipation, and shock compression). The tracking algorithm is applied to datasets from shock-capturing direct numerical simulations of DHIT, with Taylor microscale Reynolds number Re λ =40 and turbulence Mach number M t =0.2 ⁠, and STI cases in which the turbulence is processed by a shock wave at Mach numbers M = 1.5 and 3.0. Temporal surface convolution increases for initially larger structures, resulting in a higher probability of locally hyperbolic geometries where breakup into smaller structures occurs. Shock-induced deformation of the structures amplifies breakup processes, enhancing mixing, particularly for larger structures. Mixing enhancement by the shock is manifested as an amplification of the surface-averaged scalar gradient, which increases for initially larger structures. The alignment between the scalar gradient and the most extensional strain-rate eigendirection on the scalar isosurfaces also increases across the shock. Larger magnitudes of the scalar gradient and its alignment with the most compressive strain-rate eigendirection correlate with flatter surface regions. Shock-induced structure compression increases the area coverage of flat regions, where the amplification of scalar gradient is localized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Study of Shock Formation Parameters With Drive Conditions in Magnetically Accelerated Plasma Flows

We present experimental data regarding the formation of high-energy-density shocks in magnetically accelerated plasma flows using pulsed power drivers. We quantify the flow velocity and temperature of the ablated plasma using optical Thomson scattering and gated emission imaging across two different generators. We show that, regardless of the drive parameters, the plasma flows show continuous acceleration over centimeter spatial scales, in line with trends in published simulation work. When stationary targets are placed in these supersonic flows, bow-shock formation is observed at all drive parameters in a range of materials. In the higher density flow generated on the 1-MA COBRA generator at Cornell University, heating of the upstream flow ahead of the shock is observed and quantified, which is not observed at the lower density flow on the 0.2-MA Bertha driver at UC San Diego. Here, when combined with previous work on the XP generator at Cornell, we can show that these three experimental setups allow control of the effect of radiation loss and upstream absorption on the formation of the bow shock.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Direct simulation Monte Carlo with ionization and radiation

Improvements in the modeling of radiation in low density shock waves with Direct Simulation Monte Carlo (DSMC) are the subject of this study. The physical processes which determine the amount of radiation in a shock wave were investigated and the way in which they were modeled with DSMC was evaluated. Three physical processes were identified for which an improvement in the modeling technique could result in improved radiation predictions. New physical modeling schemes are introduced in this report for the three processes. First a method for determining the electric field and its effect on the flow is introduced. Second, a two step reaction process for electron impact ionization reactions is evaluated. Finally, a new scheme to determine the relaxation collision numbers for excitation of electronic states is proposed. Each new scheme attempts to move the DSMC method toward more accurate physics or more reliance on experimental data. The new schemes are all compared to the current modeling techniques and the differences in the results are evaluated. In all cases the results agree with the available data as well as, or better than the results from earlier schemes.

Carlson, Ann B.↗
DOE OSTI2022

Highly ordered graphite (HOPG) to hexagonal diamond (lonsdaleite) phase transition observed on picosecond time scales using ultrafast x-ray diffraction

Here, the response of rapidly compressed highly oriented pyrolytic graphite (HOPG) normal to its basal plane was investigated at a pressure of ~80 GPa. Ultrafast x-ray diffraction using ~100 fs pulses at the Materials Under Extreme Conditions sector of the Linac Coherent Light Source was used to probe the changes in crystal structure resulting from picosecond timescale compression at laser drive energies ranging from 2.5 to 250 mJ. A phase transformation from HOPG to a highly textured hexagonal diamond structure is observed at the highest energy, followed by relaxation to a still highly oriented, but distorted graphite structure following release. We observe the formation of a highly oriented lonsdaleite within 20 ps, subsequent to compression. This suggests that a diffusionless martensitic mechanism may play a fundamental role in phase transition, as speculated in an early work on this system, and more recent static studies of diamonds formed in impact events.

58 GEOSCIENCES↗
DOE OSTI2023

Same-sided successive-shock HED instability experiments

Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Large N algebras and generalized entropy

We construct a Type II ∞ von Neumann algebra that describes the large N physics of single-trace operators in AdS/CFT in the microcanonical ensemble, where there is no need to include perturbative 1/N corrections. Using only the extrapolate dictionary, we show that the entropy of semiclassical states on this algebra is holographically dual to the generalized entropy of the black hole bifurcation surface. From a boundary perspective, this constitutes a derivation of a special case of the QES prescription without any use of Euclidean gravity or replicas; from a purely bulk perspective, it is a derivation of the quantum-corrected Bekenstein-Hawking formula as the entropy of an explicit algebra in the G → 0 limit of Lorentzian effective field theory quantum gravity. In a limit where a black hole is first allowed to equilibrate and then is later potentially re-excited, we show that the generalized second law is a direct consequence of the monotonicity of the entropy of algebras under trace-preserving inclusions. Finally, by considering excitations that are separated by more than a scrambling time we construct a “free product” von Neumann algebra that describes the semiclassical physics of long wormholes supported by shocks. We compute Rényi entropies for this algebra and show that they are equal to a sum over saddles associated to quantum extremal surfaces in the wormhole. Surprisingly, however, the saddles associated to “bulge” quantum extremal surfaces contribute with a negative sign.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Ray-tracing image simulations of transparent objects with complex shape and inhomogeneous refractive index

Optical images of transparent three-dimensional objects can be different from a replica of the object’s cross section in the image plane, due to refraction at the surface or in the body of the object. Simulations of the object’s image are thus needed for the visualization and validation of physical models. We report ray tracing image simulations that achieved high physical fidelity, reproducing optical behaviors and image features not rendered in previous studies. We replicated brightfield microscopy images of drops with complex shapes, and images of pressure and shock waves traveling inside them. For high physical fidelity, the simulations must replicate the spatial and angular distribution of illumination rays, and both the experiment and the simulation must be designed for accurate optical modeling. The simulations are highly sensitive to the properties of the drops and can be used to diagnose and refine fluid dynamics models. The simulated images can also be optimized to extract multiple 3D properties from experimental images. Compared to specialized single-shot 3D imaging methods, this approach has the advantage that it preserves the experimental simplicity, the high resolution, and the visual interpretability characteristic to basic optical imaging. The techniques introduced here are directly applicable to optical microscopy, so they can be used in other fields, such as microfluidics and biology, to expand the type and the accuracy of three-dimensional information that can be extracted from basic optical images.

Cavitation↗
DOE OSTI2024

Diamond under extremes

Diamond is, by virtue of the covalent bonding between atoms and the very strong carbon to carbon bonds, the hardest natural material. It has been a fascinating material since its discovery, first as a decorative gem and more recently, for its numerous industrial uses because of its extreme hardness, elastic modulus, and optical transparency. In recent years, it has become a preferred ablator for laser shock experiments, and this has led to its choice as the capsule material for fusion experiments at the National Ignition Facility. Further, this review covers both experimental and computational (including machine learning) advancements in research on diamond subjected extreme conditions of temperature and pressure. The synergy between shock and ramp loading experiments and atomic level simulations is proving to be powerful in advancing our understanding of diamond under extremes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗
NASA NTRS1995

Low-energy ion acceleration at quasi-perpendicular shocks: Transverse diffusion

The problem of ion injection and acceleration at quasi perpendicular shocks has been the subject of some debate over the past two decades. It is widely known that these shocks efficiently accelerate particles that are well in the high-energy tail of the distribution. However, the issue of injection, or the acceleration of low-energy ions, has yet to reach a consensus. The fundamental issue is whether there is enough diffusion normal to the magnetic field for the particles to remain near the shock. Since transverse diffusion is a physical process that is not well understood in space plasmas, this is an important, and difficult issue to address. In this report, we will investigate the ion injection problem by performing test particle orbit integrations using synthesized turbulent fields. These fields are fully three-dimensional so that transverse diffusion is possible (cross-field diffusion is not possible in geometries where the electromagnetic fields are less than three dimensional). The synthesized fields are produced by superimposing a three-dimensional wave field on a background field. For completeness, we will compare the results from this model with the more well-established theories, such as the diffusive approximation and scatter-free shock drift acceleration. We will also compare these results with other numerical simulation techniques such as the well known hybrid simulation, and other test-particle calculations in which the shock fields are specified to have less than three dimensions. We will also discuss some recent relevant observations and how these compare with our results.

Giacalone, J.↗
DOE OSTI2026

Calculating shock Hugoniot and isentropes using multiphase equation of state tables and application to shock and release of diamond ablators in inertial confinement fusion implosions

Advances in shock and ramp compression techniques now allow experimental access to unprecedented extreme conditions of pressure and temperature, providing a means to test theoretical models. Here, we describe a simple methodology to compute multi-phase shock Hugoniot and isentropes using multiphase equation of state tables. We treat explicitly the phase coexistence along the phase boundary to reveal the evolution of the sample as it undergoes the phase transformation in adiabatic conditions. We illustrate the method by calculating the predicted shock and shock-and-release behavior of diamond at conditions relevant for the initial stage of inertial confinement fusion implosions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2017

THEMIS Satellite Observations of Hot Flow Anomalies at Earth's Bow Shock

Hot flow anomalies (HFAs) at Earth's bow shock were identified in Time History of Events and Macroscale Interactions During Substorms (THEMIS) satellite data from 2007 to 2009. The events were classified as young or mature and also as regular or spontaneous hot flow anomalies (SHFAs). The dataset has 17 young SHFAs, 49 mature SHFAs, 15 young HFAs, and 55 mature HFAs. They span a wide range of magnetic local times (MLTs) from approximately 7 to 16.5MLT. The largest ratio of solar wind to HFA core density occurred near dusk and at larger distances from the bow shock. In this study, HFAs and SHFAs were observed up to 6.3 RE and 6.1 RE (Earth radii), respectively, upstream from the model bow shock. HFA-SHFA occurrence decreases with distance upstream from the bow shock. HFAs of the highest event core ion temperatures were not seen at the flanks. The ratio of HFA ion temperature increase to HFA electron temperature increase is highest around 12MLT and slightly duskward. For SHFAs, (Tihfa=Tisw)/(Tehfa=Tesw) generally increased with distance from the bow shock. Both mature and young HFAs are more prevalent when there is an approximately radial interplanetary magnetic field. HFAs occur most preferentially for solar wind speeds from 550 to 600 km s-1. The correlation coefficient between the HFA increase in thermal energy density from solar wind values and the decrease in kinetic energy density from solar wind values is 0.62. SHFAs and HFAs do not show major differences in this study.

Interplanetary physics (planetary bow shocks; sola↗
DOE OSTI2024

Probing shocks interacting with radiation waves with the Radishock experiment

Both radiation flows and shocks have been extensively studied in the laboratory in the past few decades due to their critical roles in many astrophysical and high-energy density physics processes. In the Radishock experiment, a halfraum-powered radiation wave is driven into a low-density foam and interacts with an ablatively driven, counter-propagating shock. The interacting waves produce a spike in energy density with a temperature greater than the local temperature of the individual waves. As in the successful predecessor experiment, COAX, the primary diagnostic uses absorption spectroscopy at many locations down the cylindrical target, enabling a spatial temperature inference of the radiation wave and its interactions with the shock. Combined with a radiography diagnostic that is capable of imaging the shock and interaction features, we are able to study and inform model predictions of the interaction spike phenomenon. We describe the underlying physics behind the shock interactions with the radiation front and the implications of this experimental study for a broad range of astrophysical phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗