Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock compression”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Numerical solutions of Navier-Stokes equations for compressible turbulent two/three dimensional flows in terminal shock region of an inlet/diffuser

The multidimensional ensemble averaged compressible time dependent Navier Stokes equations in conjunction with mixing length turbulence model and shock capturing technique were used to study the terminal shock type of flows in various flight regimes occurring in a diffuser/inlet model. The numerical scheme for solving the governing equations is based on a linearized block implicit approach and the following high Reynolds number calculations were carried out: (1) 2 D, steady, subsonic; (2) 2 D, steady, transonic with normal shock; (3) 2 D, steady, supersonic with terminal shock; (4) 2 D, transient process of shock development and (5) 3 D, steady, transonic with normal shock. The numerical results obtained for the 2 D and 3 D transonic shocked flows were compared with corresponding experimental data; the calculated wall static pressure distributions agree well with the measured data.

Liu, N. S.↗
NASA NTRS1992

Studies in turbulence

Various papers on turbulence are presented. Individual topics addressed include: modeling the dissipation rate in rotating turbulent flows, mapping closures for turbulent mixing and reaction, understanding turbulence in vortex dynamics, models for the structure and dynamics of near-wall turbulence, complexity of turbulence near a wall, proper orthogonal decomposition, propagating structures in wall-bounded turbulence flows. Also discussed are: constitutive relation in compressible turbulence, compressible turbulence and shock waves, direct simulation of compressible turbulence in a shear flow, structural genesis in wall-bounded turbulence flows, vortex lattice structure of turbulent shear slows, etiology of shear layer vortices, trilinear coordinates in fluid mechanics.

Gatski, Thomas B.↗
NASA NTRS1983

Corotating pressure waves without streams in the solar wind

Voyager 1 and 2 magnetic field and plasma data are presented which demonstrate the existence of large scale, corotating, non-linear pressure waves between 2 AU and 4 AU that are not accompanied by fast streams. The pressure waves are presumed to be generated by corotating streams near the Sun. For two of the three pressure waves that are discussed, the absence of a stream is probably a real, physical effect, viz., a consequence of deceleration of the stream by the associated compression wave. For the third pressure wave, the apparent absence of a stream may be a geometrical effect; it is likely that the stream was at latitudes just above those of the spacecraft, while the associated shocks and compression wave extended over a broader range of latitudes so that they could be observed by the spacecraft. It is suggested that the development of large-scale non-linear pressure waves at the expense of the kinetic energy of streams produces a qualitative change in the solar wind in the outer heliosphere. Within a few AU the quasi-stationary solar wind structure is determined by corotating streams whose structure is determined by the boundary conditions near the Sun.

Burlaga, L. F.↗
NASA NTRS1983

Corotating pressure waves without fast streams in the solar wind

Voyager 1 and 2 magnetic field and plasma data are presented which demonstrate the existence of large scale, corotating, non-linear pressure waves between 2 AU and 4 AU that are not accompanied by fast streams. The pressure waves are presumed to be generated by corotating streams near the sun. For two of the three pressure waves that are discussed, the absence of a stream is probably a real, physical effect, viz., a consequence of deceleration of the stream by the associated compression wave. For the third pressure wave, the apparent absence of a stream may be a geometrical effect; it is likely that the stream was at latitudes just above those of the spacecraft, while the associated shocks and compression wave extended over a broader range of latitudes so that they could be observed by the spacecraft. It is suggested that the development of large-scale non-linear pressure waves at the expense of the kinetic energy of streams produces a qualitative change in the solar wind in the outer heliosphere. Within a few AU the quasi-stationary solar wind structure is determined by corotating streams whose structure is determined by the boundary conditions near the sun. Previously announced in STAR as N83-19694

Burlaga, L. F.↗
NASA NTRS1978

The artificial compression method for computation of shocks and contact discontinuities. III - Self-adjusting hybrid schemes

This paper presents a new computational method for the calculation of discontinuous solutions of hyperbolic systems of conservation laws, which deal effectively with both shock and contact discontinuities. The method consists of two stages: in the first stage a standard finite-difference scheme is hybridized with a nonoscillatory first order accurate method to provide for the monotonic variation of the solution near discontinuities, and in the second stage artificial compression is applied to sharpen transitions at discontinuities. This modification of a standard finite-difference method results in a scheme which preserves the order of truncation error of the original method and yet yields a sharp and oscillation free transition for both shocks and contact discontinuities. The modification can be easily implemented in existing computer codes.

Harten, A.↗
NASA NTRS1989

Separation shock motion in fin, cylinder, and compression ramp - Induced turbulent interactions

In conjunction with new experimental results at Mach 5, an examination has been made of published data on unsteadiness of shock-induced turbulent boundary-layer separation. The data are all wall pressure fluctuation measurements made under the unsteady separation shock and are from interactions induced by compression ramps, blunt and sharp fins, and circular cylinders. There is little evidence of a link between the separation shock zero-crossing frequency and characteristic frequency of the incoming boundary layer. The low shock frequencies and low shock speeds, and the trends with changes in model geometric parameters and incoming boundary layer, suggest that turbulent or global fluctuations at the upstream boundary of the separated flow drive the shock motion.

Dolling, D. S.↗
DOE OSTI2025

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗
NASA NTRS2004

Spectroscopy of the Perseus Cluster

We present preliminary results of a XMM-Newton 50 ks observation of the Perseus Cluster that provides an unprecedented view of the central 0.5 Mpc region. The projected gas temperature declines smoothly by a factor of 2 from a maximum value of approx. 7 keV in the outer regions to just above 3 keV at the cluster center. Over this same range, the heavy-element abundance rises slowly from 0.4 to 0.5 solar as the radius decreases from 14 ft. to 5 ft., and then it rises to a peak of almost 0.7 solar at 1&farcm;25 before declining to 0.4 at the center. Th global east-west asymmetry of the gas temperature and surface brightness distributions, approximately aligned with the chain of bright galaxies, suggests an ongoing merger, although the modest degree of the observed asymmetry certainly excludes a major merger interpretation. The chain of galaxies probably traces the filament along which accretion started some time ago and is continuing at the present time. A cold and dense (low-entropy) cluster core like Perseus is probably well "protected" against the penetration of the gas of infalling groups and poor clusters, whereas in non-cooling core clusters such as Coma and A1367, infalling subclusters can penetrate deeply into the core region. In Perseus, gas associated with infalling groups may be stripped completely at the outskirts of the main cluster and only compression waves (shocks) may reach the central regions. We argue, and show supporting simulations, that the passage of such a wave(s) can qualitatively explain the overall horseshoe shaped appearance of the gas temperature map (the hot horseshoe surrounds the colder, low-entropy core) as well as other features of the Perseus Cluster core. These simulations also show that as compression waves traverse the cluster core, they can induce oscillatory motion of the cluster gas that can generate multiple sharp "edges" on opposite sides of the central galaxy. Gas motions induced by mergers may be a natural way to explain the high frequency of "edges" seen in clusters with cooling cores.

Jones, Christine↗
NASA NTRS2014

MHD Modelling of Coronal Loops: Injection of High-Speed Chromospheric Flows

Context. Observations reveal a correspondence between chromospheric type II spicules and bright upward-moving fronts in the corona observed in the extreme-ultraviolet (EUV) band. However, theoretical considerations suggest that these flows are probably not the main source of heating in coronal magnetic loops. Aims. We investigate the propagation of high-speed chromospheric flows into coronal magnetic flux tubes and the possible production of emission in the EUV band. Methods. We simulated the propagation of a dense 104 K chromospheric jet upward along a coronal loop by means of a 2D cylindrical MHD model that includes gravity, radiative losses, thermal conduction, and magnetic induction. The jet propagates in a complete atmosphere including the chromosphere and a tenuous cool (approximately 0.8 MK) corona, linked through a steep transition region. In our reference model, the jet initial speed is 70 km per second, its initial density is 10(exp 11) per cubic centimeter, and the ambient uniform magnetic field is 10 G. We also explored other values of jet speed and density in 1D and different magnetic field values in 2D, as well as the jet propagation in a hotter (approximately 1.5 MK) background loop. Results. While the initial speed of the jet does not allow it to reach the loop apex, a hot shock-front develops ahead of it and travels to the other extreme of the loop. The shock front compresses the coronal plasma and heats it to about 10(exp 6) K. As a result, a bright moving front becomes visible in the 171 Angstrom channel of the SDO/AIA mission. This result generally applies to all the other explored cases, except for the propagation in the hotter loop. Conclusions. For a cool, low-density initial coronal loop, the post-shock plasma ahead of upward chromospheric flows might explain at least part of the observed correspondence between type II spicules and EUV emission excess.

Sun: chromosphere↗
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 OSTI2024

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Probabilistic flux limiters

The stable numerical integration of shocks in compressible flow simulations relies on the reduction or elimination of Gibbs phenomena (unstable, spurious oscillations). A popular method to virtually eliminate Gibbs oscillations caused by numerical discretization in under-resolved simulations is to use a flux limiter. A wide range of flux limiters have been studied in the literature, with recent interest in their optimization via machine learning methods trained on high-resolution datasets. The common use of flux limiters in numerical codes as plug-and-play blackbox components makes them key targets for design improvement. Even for deterministic dynamical models, numerical uncertainty is introduced via coarse-graining required by insufficient computational power to solve all scales of motion. Conventional flux limiters are deterministic and lack the capacity to address uncertainties, both aleatoric (inherent randomness) and epistemic (modeling uncertainty due to limited knowledge), which arise in coarse-grained numerical simulations. Here, we introduce a conceptually distinct type of flux limiter that is designed to handle the effects of randomness in the model and uncertainty in model parameters. Unlike traditional single-function flux limiters, these new probabilistic flux limiters incorporate multiple flux limiting functions, each applied with a learned probability drawn from high-resolution data to mitigate the effects of uncertainty in numerical simulations. This approach departs from traditional single-function limiters by explicitly modeling and incorporating uncertainty into the shock capturing process. Using the example of Burgers' equation as a testbed, we show that a machine learned, probabilistic flux limiter may be used in a shock capturing code to more accurately capture shock profiles. In particular, we show that our probabilistic flux limiter outperforms standard limiters and can be successively improved upon (up to a point) by expanding the set of probabilistically chosen flux limiting functions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Implementation of compound refractive lenses for large field-of-view x-ray phase-contrast imaging during hypervelocity impact experiments

Synchrotron x-ray phase-contrast imaging (XPCI) offers time-resolved visualization of dynamic compression phenomena, but its intrinsically small field-of-view (FOV) limits the time that key features remain in frame. A novel approach to enlarge the FOV is achieved by positioning a two-dimensional parabolic compound refractive lens (CRL) upstream of the sample to deliberately defocus the white beam. Ray-tracing simulations and XPCI measurements show that this CRL configuration can expand the beam by ∼50% vertically and ∼15% horizontally based on the full width at half-maximum of the beam. Implementing the CRL, however, attenuates the photon flux and lowers signal-to-noise ratio (SNR). Task-based analysis using a calibration grid (30 μm dots) showed that both setups fail to consistently meet the Rose criterion (SNR ≥ 5) for features of this size in single-bunch imaging. Extrapolating the measured SNR Rose values suggests that the minimum consistently detectable feature lies closer to 30–40 μm for the standard XPCI setup and above 40 μm for CRL-XPCI. Despite this limitation, the CRL configuration nearly doubles the illuminated area, enabling simultaneous tracking of front and rear observations of boron carbide targets subjected to rod and sphere impacts at 1.0–2.6 km/s. Image tracking algorithms and photonic Doppler velocimetry were used to measure penetration and rear-surface velocity histories. Together, these measurements capture crack fronts, penetration, and material breakout, offering new benchmark data for validating high-strain-rate constitutive models of ceramic materials.

Ceramic materials↗
DOE OSTI2025

A surprising proliferation of detwinning in β -tin at extreme loading rates

Integrating data from dynamic compression experiments of condensed matter across three national laboratories has led to insight and quantitative calibration of materials strength over decades of loading rate. For many materials, a single strength model (such as PTW) is sufficient to capture the flow-stress strain rate relationship which is monotonic. Here, we show here that β -tin, a tetragonal metal, exhibits dramatic deviations from this behavior. Naive fitting to a single PTW model is insufficient to capture the behavior; indeed, such resulting inferred flow stress versus strain exhibits a non-monotonic behavior. We suggest a resolution to this by proposing that in β -tin there are important Bauschinger effects arising from favorable conditions for twinning and detwinning. A simple yield surface model when paired with PTW hardening captures the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic strength of iron under pressure-temperature conditions of Earth’s inner core

Iron (Fe) is a primary constituent of terrestrial planetary cores, yet its rheological properties under extreme conditions remain uncertain. Here we present direct measurements of Fe strength at 310-430 GPa pressures and 3700-5800 K temperatures, obtained using Rayleigh-Taylor (RT) instability experiments at the National Ignition Facility. Single-crystal α-Fe samples with [001] and [111] orientations are shock-ramp compressed past the α-ε transition along paths approaching Earth’s inner core conditions. We find that ε-Fe derived from [001] α-Fe is consistently stronger (11-20 GPa) than that from [111] α-Fe (8-18 GPa), contrary to the trend at ambient conditions. Large-scale molecular dynamics simulations reproduce this atypical strength anisotropy and attribute it to microstructural evolution during the phase transition and subsequent ε-phase plasticity. Ripple growth analysis further constrains viscosities of 100-170 Pa·s under the driven conditions. Our results provide experimental benchmarks for Fe rheology at inner-core conditions, with implications for seismic anisotropy and the geodynamo.

Condensed-matter physics↗
DOE OSTI2024

High pressure suppression of plasticity due to an overabundance of shear embryo formation

Abstract High pressure shear band formation is a critical phenomenon in energetic materials due to its influence on both mechanical strength and mechanochemical activation. While shear banding is known to occur in a variety of these materials, the governing dynamics of the mechanisms are not well defined for molecular crystals. We conduct molecular dynamics simulations of shock wave induced shear band formation in the energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) to assess shear band nucleation processes. We find, that at high pressures, the initial formation sites for shear bands, “embryos”, form in excess and rapidly lower deviatoric stresses prior to shear band formation and growth. This results in the suppression of plastic deformation. A local cluster analysis is used to quantify and contrast this mechanism with a more typical shear banding seen at lower pressures. These results demonstrate a mechanism that is reversible in nature and that supersedes shear band formation at increased pressures. We anticipate that these results will have a broad impact on the modeling and development of high-strain rate application materials such as those for high explosives and hypersonic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Validation of 266 nm UV Raman thermometry using the alpha-to-beta quartz transition

Ultraviolet Raman scattering is a promising diagnostic for temperature measurements relevant to dynamic compression experiments, where improved signal strength and reduced background are critical. In this work, we experimentally validate Raman thermometry at 266 nm using crystalline quartz heated through the α to β phase transition in a tube furnace. Quartz provides a useful benchmark because of its well-characterized transition temperature near 846 K at ambient pressure. Temperatures inferred from the Stokes to anti-Stokes intensity ratio are compared with thermocouple measurements and with the observed phase transition. Single-shot Raman temperatures track the thermocouple closely, with an empirical single-shot scatter of approximately ±15 K under the present conditions. Temporal averaging reduces shot-to-shot scatter and improves agreement between the Raman temperatures and thermocouple readings to within approximately 8 K, while the transition is identified within 5 K of the accepted value. Measurements are performed using the fourth harmonic of a nanosecond pulsed Nd:YAG laser at 266 nm. The shorter excitation wavelength provides enhanced Raman scattering efficiency and reduced visible background relative to conventional visible excitation. Comparison with 532 nm Raman measurements acquired in the same experimental geometry demonstrates substantially improved signal-to-noise ratio for the UV configuration. These results establish 266 nm Raman thermometry as a reliable and accurate temperature diagnostic under controlled heating conditions and provide a quantitative benchmark for future application in dynamic compression and shock physics experiments.

Chemistry - Chemical explosives↗