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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Spurious Numerical Oscillations in Numerical Simulation of Supersonic Flows Using Shock Capturing Schemes

The numerical simulation of transitional and turbulent processes in hypersonic boundary layers often involves a physical process of a shock-disturbance wave interaction in complex two-dimensional and three-dimensional flow fields. For such simulations, it is required that there be a high order of accuracy in capturing both the shock and the small disturbances. The purpose of this paper is to evaluate the viability of using high order shock capturing schemes to track small disturbances in a multi-dimensional steady hypersonic flow. The numerical methods that are to be studied are the Total Variation Diminishing (TVD) scheme, and Essentially Non-Oscillatory (ENO) scheme. This paper shows that the presence of numerical oscillations in the flow field solution may drastically hinder any attempt at tracking the propagation of any physical disturbances. It has been found that the numerical oscillations that exist for shock capturing methods may be significant enough to pollute a flow field containing small physical disturbances. The effects of the refinement of the grid do not reduce the oscillations, but rather they decrease the wavelength of the oscillations. It is shown that by aligning the shock with the grid, the amplitude of these spurious oscillations may be greatly reduced.

Lee, Theodore K.↗
DOE OSTI2020

A thermodynamic analysis of hypervelocity impacts on metals

Hypervelocity impacts on solid surfaces generate and propagate shock waves and create extreme density, pressure, and temperature conditions. Understanding hypervelocity impacts requires a multi-physics approach at the intersection of fluid dynamics, solid mechanics, and quantum theory, which has the capability of spanning multiple different dominating physics regimes. These strong shock conditions found in hypervelocity impacts are studied in the context of high-speed collisions between spacecraft and meteoroids and debris. We study the problem of hypervelocity impacts on metals with the lens of understanding the processes that lead to the creation of an impact generated plasma. In order to characterize the impacting particle as well as its threat to spacecraft from the subsequent radiating plasma an analysis was performed that highlights how the kinetic energy of the impactor is deposited into a target, how the shock wave and shocked material in the target evolve, and how these effects can lead towards the creation of an impact plasma through heating. The analysis suggests that at low velocities the majority of the energy from a hypervelocity impact is stored in reversible elastic energy and does not contribute to strong heating. At higher impact velocities, nucleus and electron thermal components dominate, increasing the energy available for ionization. Overall, this work provides the complete set of state variable initial conditions for any study of hypervelocity impact plasma expansion and its effects.

42 ENGINEERING↗
NASA NTRS1990

MHD intermediate shocks in coronal mass ejections

A simplified model of coronal mass ejections is considered in which at least a portion of the interaction with the background corona involves a shock wave, and the allowable shock solutions and their compressive signatures are examined. The MHD shock-jump equations have a maximum of three possible types of solutions with an entropy rise for fixed values of the physical variables (slow, intermediate, and fast shocks). However, one of the three solution classes (the intermediate shock) is widely believed to not occur in nature and is regarded as nonevolutionary or extraneous. Without the intermediate shock, there is no multiplicity of solutions in that only one shock (or none) can occur for given physical values. All three potential shock types are considered, and it is shown solely on the basis of the shock-jump equations, that intermediate shocks must exist along some segment of the shock front for certain parametric regimes and for conditions that probably occur in some coronal mass ejections.

Steinolfson, R. S.↗
NASA NTRS1974

Plasma instability modes related to the earth's bow shock

The present work examines the status of physical interpretations of some of the microscopic phenomena occurring in bow shock structures. A categorization of microscopic phenomena is given, and it is examined how various modes may or may not be invoked in explaining spacecraft measurements on bow structure. The macroscopic and observational context of the bow shock as presently understood is first defined, and then some of the microscopic plasma physical phenomena which might be expected to be found associated with certain macroscopic structures are outlined. Some problems in the use of the bow shock to test plasma shock theory are then discussed.

Greenstadt, W. W.↗
NASA NTRS1993

High-energy emission from the eclipsing millisecond pulsar PSR 1957+20

The properties of the high-energy emission expected from the eclipsing millisecond pulsar system PSR 1957+20 are investigated. Emission is considered by both the relativistic shock produced by the pulsar wind in the nebula surrounding the binary and by the shock constraining the mass outflow from the companion star of PSR 1957+20. On the basis of the results of microscopic plasma physical models of relativistic shocks it is suggested that the high-energy radiation is produced in the range from X-rays to MeV gamma rays in the binary and in the range from 0.01 eV to about 40 keV in the nebula. Doppler boost of the emission in the radiating wind suggests the flux should vary on the orbital time scale, with the largest flux observed roughly coincident with the pulsar's radio eclipse.

Arons, Jonathan↗
NASA NTRS2007

A Physical Relationship Between Electron-Proton Temperature Equilibration and Mach Number in Fast Collisionless Shocks

The analysis of Balmer-dominated optical spectra from nonradiative (adiabatic) SNRs has shown that the ratio of the electron to proton temperature at the blast wave is close to unity at v(s) approximate or smaller than 400 km s(-1) but declines sharply down to the minimum value of m(e) /m(p) dictated by the jump conditions at shock speeds exceeding 2000 km s(-1). We propose a physical model for the heating of electrons and ions in non-cosmic-ray-dominated, strong shocks (v(s) 400 km s(-1)) wherein the electrons are heated by lower hybrid waves immediately ahead of the shock front. These waves arise naturally from the cosmic ray pressure gradient upstream from the shock. Our model predicts a nearly constant level of electron heating over a wide range of shock speeds, producing a relationship (T(e)/T(p))sub 0 proportional to v(-2/s) (proportional to M(-2)) that is fully consistent with the observations.

SNR(ONRADIATIVE SUPERNOVA REMNANTS)↗
NASA NTRS1997

A Study of Fundamental Shock Noise Mechanisms

This paper investigates two mechanisms fundamental to sound generation in shocked flows: shock motion and shock deformation. Shock motion is modeled numerically by examining the interaction of a sound wave with a shock. This numerical approach is validated by comparison with results obtained by linear theory for a small-disturbance case. Analysis of the perturbation energy with Myers' energy corollary demonstrates that acoustic energy is generated by the interaction of acoustic disturbances with shocks. This analysis suggests that shock motion generates acoustic and entropy disturbance energy. Shock deformation is modeled numerically by examining the interaction of a vortex ring with a shock. These numerical simulations demonstrate the generation of both an acoustic wave and contact surfaces. The acoustic wave spreads cylindrically. The sound intensity is highly directional and the sound pressure increases with increasing shock strength. The numerically determined relationship between the sound pressure and the Mach number is found to be consistent with experimental observations of shock noise. This consistency implies that a dominant physical process in the generation of shock noise is modeled in this study.

Meadows, Kristine R.↗
DOE OSTI2025

A limit to strong shock behavior in the dynamic response of matter at pressure

Solids under high pressures experience a series of regimes, where their microstructure adapts to the applied compression and these key transitions are discussed in this paper. As strain increases, new forces emerge at extreme pressures. A previous study introduced the concept of the weak shock limit (WSL), at which the ambient theoretical shear strength is overcome. Above the WSL, further deformation under strong shock conditions results in electrons occupying higher energy levels as strain increases. As pressure rises further, shock melting occurs in the material and at around three times this melting pressure, the strong shock limit is reached where the driving physics under pressure switches, with electrons forced into higher energy states. This leads to significant reduction in their compressibility due to changes in electronic structure and developing electron degeneracy pressures. A derivation for conditions at this state is presented, which indicates that a dependence of the threshold pressure on the free electron number density defines the limit observed. This correlation suggests that ambient material moduli govern material compression up to nearly 50% strain. These observations show that models should account for different behaviors as dominant physics changes in each regime accessed as shock pressure increases.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Shock wave ignition of metal powders

The interaction between isolated metal particles and shock-induced flowfields is analyzed to characterize physical events and determine experimental feasibility for shock wave ignition studies. Emphasis is placed on incident shock wave phenomenon as opposed to reflected shock wave. A pre-ignition particle-flow interaction is hypothesized which presumes a vapor phase ignition reaction and an absence of significant vaporization prior to particle melting and erosion. Rapid vaporization of eroded mass is expected, with ignition reaction in the particle wakes. Experimental observation of the 6438-A Cd line and the 3721-A MgO band support the ignition mechanism hypothesized. Scattered cloud data are believed to substantiate the hypothesis that ignition delays could be correlated with predicted melting times.

Fox, T. W.↗
NASA NTRS2003

Mesoscale Computational Investigation of Shocked Heterogeneous Materials with Application to Large Impact Craters

The propagation of shock waves through target materials is strongly influenced by the presence of small-scale structure, fractures, physical and chemical heterogeneities. Pre-existing fractures often create craters that appear square in outline (e.g. Meteor Crater). Reverberations behind the shock from the presence of physical heterogeneity have been proposed as a mechanism for transient weakening of target materials. Pre-existing fractures can also affect melt generation. In this study, we are attempting to bridge the gap in numerical modeling between the micro-scale and the continuum, the so-called meso-scale. To accomplish this, we are developing a methodology to be used in the shock physics hydrocode (CTH) using Monte-Carlo-type methods to investigate the shock properties of heterogeneous materials. By comparing the results of numerical experiments at the micro-scale with experimental results and by using statistical techniques to evaluate the performance of simple constitutive models, we hope to embed the effect of physical heterogeneity into the field variables (pressure, stress, density, velocity) allowing us to directly imprint the effects of micro-scale heterogeneity at the continuum level without incurring high computational cost.

Crawford, D. A.↗
DOE OSTI2024

Temperature and density dependent pair potential for deuterium under shock

Large-scale classical molecular dynamics (CMD) simulations naturally include the microscopic physics necessary for atomistic modeling of shock release at the ablator-fuel interface in an inertial confinement fusion (ICF) capsule. Here, the multi-megabar shocks utilized in ICF experiments can drive the deuterium fuel from ambient to electron volt temperatures (T) and multi-fold compression. Modeling interatomic interactions over such an extreme range of conditions is challenging for empirical bond order potentials. We generate a pair potential for deuterium with explicit temperature and mass density dependence from ab initio density functional theory molecular dynamics using the iterative Boltzmann inversion method. This potential accurately reproduces the radial distribution functions and pressures from DFT in CMD equilibrium simulations across a wide range of thermodynamic conditions, yet fails to return the expected Hugoniot relations when used in direct CMD shock simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Is the common definition of the Mie-Grueneisen equation of state inconsistent

Anharmonic models are used to describe the high temperatures and the electronic effects characteristic of shock wave phenomena and planetary interiors. It is found that a generalized Grueneisen function dependent on temperature and volume is necessary for the self-consistent application of the Mie-Grueneisen equation at the physical conditions encountered in shock wave phenomena or in the interior of planets. It is noted that at temperatures of a few thousand degrees Kelvin, by deriving the function in terms of anharmonic lattice dynamics, the difference between gamma and gamma(th) is of the same order of a magnitude as gamma. This leads to the conclusions that (1) gamma is not a function of temperature and (2) gamma?gamma(th) does not explain either the equation of state of planetary interiors or high-temperature shock wave phenomena.

Mulargia, F.↗
NASA NTRS1999

HRI Observations of Balmer Dominated Filaments in the SNR RCW86

The objective of this project was to use the ROSAT HRI to study the recently discovered optical Balmer-dominated filaments in the young SNR R-CW86. The observations provide complete high-resolution X-ray coverage of the shell of the SNR. These X-ray observations are combined with new optical observations (both imaging and spectroscopic), and new high resolution radio observations to provide a better overall understanding of the state of the remnant (Sedov or reverse shock), its history (as the possible SNR of SN 185 AD), and the physics of non-radiative shocks.

Long, Knox S.↗
DOE OSTI2020

Thermal interfaces in dynamic compression experiments

Temperature is a difficult thermodynamic variable to measure in dynamic compression experiments. Optical pyrometry is a general-purpose technique for measuring temperature from a radiant surface, but that surface is often the interface between distinct materials with temperatures that vary spatially along the loading direction. This leads to a fundamental problem: how does the measured interface temperature relate to this temperature profile along the compression axis? Numerical analysis of loading history and thermal diffusion at these interfaces shows that seemingly subtle changes in experiment geometry can lead to very different temperature profiles. We compare these results to laboratory temperature measurements of shock-compressed tin.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1984

On the proper Mach number and ratio of specific heats for modeling the Venus bow shock

Observational data from the Pioneer Venus Orbiter are used to investigate the physical characteristics of the Venus bow shock, and to explore some general issues in the numerical simulation of collisionless shocks. It is found that since equations from gas-dynamic (GD) models of the Venus shock cannot in general replace MHD equations, it is not immediately obvious what the optimum way is to describe the desired MHD situation with a GD code. Test case analysis shows that for quasi-perpendicular shocks it is safest to use the magnetospheric Mach number as an input to the GD code. It is also shown that when comparing GD predicted temperatures with MHD predicted temperatures total energy should be compared since the magnetic energy density provides a significant fraction of the internal energy of the MHD fluid for typical solar wind parameters. Some conclusions are also offered on the properties of the terrestrial shock.

Tatrallyay, M.↗
DOE OSTI2026

Validation of SPH code Spheral to model interacting solid bodies in a supersonic flow

Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.

Airburst↗
NASA NTRS1994

Craters in aluminum 1100 targets using glass projectiles at 1-7 km/s

We report on impact experiments using soda-lime glass spheres of 3.2 mm diameter and aluminum targets (1100 series). The purpose is to assist in the interpretation of LDEF instruments and in the development of future cosmic-dust collectors in low-Earth orbit. Because such instruments demand understanding of both the cratering and penetration process, we typically employ targets with thicknesses that range from massive, infinite half-space targets, to ultrathin films. This report addresses a subset of cratering experiments that were conducted to fine-tune our understanding of crater morphology as a function of impact velocity. Also, little empirical insight exists about the physical distribution and shock-metamorphism of the impactor residues as a function of encounter speed, despite their recognized significance in the analysis of space-exposed surfaces. Soda-lime glass spheres were chosen as a reasonable analog to extraterrestrial silicates, and aluminum 1100 was chosen for targets, which among the common Al-alloys, best represents the physical properties of high-purity aluminum. These materials complement existing impact studies that typically employed metallic impactors and less ductile Al-alloys. We have completed dimensional analyses of the resulting craters and are in the process of investigating the detailed distribution of the unmelted and melted impactor residues via SEM methods, as well as potential compositional modifications of the projectile melts via electron microprobe.

Bernhard, R. P.↗