Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock physics”

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

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

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Selected 3D Flash-X Checkpoints for the Long-Time Evolution of a 9.6 Solar-Mass Core-Collapse Supernova Model

This dataset contains selected 3D Flash-X checkpoint files from the long-time evolution of a low-energy core-collapse supernova explosion of a 9.6 Msun zero-metallicity, low-mass iron-core progenitor. The checkpoints span the shock-breakout phase through the young-remnant phase, ending at approximately 3 yr after core bounce. The dataset is intended for follow-up analysis and post-processing, especially radiation-transport calculations using the hydrodynamic and compositional structure of the ejecta. For a full description of the numerical setup, physical assumptions, limitations, and interpretation of the simulation, users should refer to the associated paper.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Long-term stability and properties of zirconia ceramics for heavy duty diesel engine components

Physical, mechanical, and thermal properties of commercially available transformation-toughened zirconia are measured. Behavior is related to the material microstructure and phase assemblage. The stability of the materials is assessed after long-term exposure appropriate for diesel engine application. Properties measured included flexure strength, elastic modulus, fracture toughness, creep, thermal shock, thermal expansion, internal friction, and thermal diffusivity. Stability is assessed by measuring the residual property after 1000 hr/1000C static exposure. Additionally static fatigue and thermal fatigue testing is performed. Both yttria-stabilized and magnesia-stabilized materials are compared and contrasted. The major limitations of these materials are short term loss of properties with increasing temperature as the metastable tetragonal phase becomes more stable. Fine grain yttria-stabilized material (TZP) is higher strength and has a more stable microstructure with respect to overaging phenomena. The long-term limitation of Y-TZP is excessive creep deformation. Magnesia-stabilized PSZ has relatively poor stability at elevated temperature. Overaging, decomposition, and/or destabilization effects are observed. The major limitation of Mg-PSZ is controlling unwanted phase changes at elevated temperature.

Larsen, D. C.↗
NASA NTRS1993

Comparisons between DSMC and the Navier-Stokes equations for reentry flows

A detailed comparison is made between Navier-Stokes and DSMC calculations for flows near the continuum limit to assess the accuracy of the continuum equations in this regime. Meaningful comparisons require the use of similar physical models. This necessitates the inclusion of a separate rotational energy equation and use of slip boundary conditions. Inclusion of slip boundary conditions resulted in improved agreement between surface properties. Moreover, good agreement was obtained for the various temperatures in the nonequilibrium portion of the flow field that does not contain the shock region. Departures are noted in the shock region and in regions where thermal diffusion effects are important.

Olynick, David R.↗
DOE OSTI2024

Simulations of n -dodecane/oxygen/nitrogen cellular detonations

In this work, two-dimensional n-dodecane/air/nitrogen cellular detonations are simulated with various equivalence ratios (ERs). A skeletal mechanism consisting of 54 species and 269 reactions is used. The lower and upper equivalence ratio boundaries for self-sustained detonation are 0.3 and 2.2, respectively. Detonation with different regimes characterized by the detonation cell patterns is observed, which aligns well with the category based on the stability parameter, i.e., weakly and highly unstable detonations, and extinction. Further, in terms of the frontal structure, non-negligible effect of diffusion on the cellular detonation is revealed, especially in the vicinity of the leading shock front. In highly unstable and quenched detonations, the alternation in reaction pathway within the induction zone accounts for the changes of detonation dynamics, such as the absence or extended sequence of important radical formation, e.g., OH. In addition, the composition of the unburned pockets depends on both pocket location from the leading shock front and the ER in the fresh mixture, because the former determines the residence time, whilst the latter affects the pocket reaction rate.

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

Experimental and modelling study of the effect of pre-shock duration on the double shock initiation response of PBX 9501

Here, a series of double shock initiation gas-gun experiments have been performed on PBX 9501 (95% weight HMX, 5% binder) in which the input conditions were nominally the same, but the time delay between the two shocks was explicitly varied. The impact conditions generated pre- and main shock pressures of 2.9 and 6.2 GPa respectively, with time delays between them ranging from 0.56 to 1.10 µs. This resulted in different levels of reaction behind the pre-shock prior to the main shock traversing the explosive samples, with the growth of reaction monitored using embedded particle velocity and shock tracker gauges. Simulated results from three reactive burn models (CREST, AWSD and SURF) are compared with the measured data to assess their ability to predict the response of PBX 9501 under such double shock conditions. The described experiments add to the body of evidence to help improve our understanding of the double shock initiation response of explosives, and the data obtained provide an excellent test for reactive burn models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Two self-similar Reynolds-stress transport models with anisotropic eddy viscosity

Two Reynolds-averaged Navier-Stokes models with full Reynolds-stress transport (RST) and tensor eddy viscosity are presented. These new models represent RST extensions of the $k−2L−a−\mathscr{C}$ and $k−ϕ−L−a−\mathscr{C}$ models by Morgan. Self-similarity analysis is used to derive constraints on model coefficients required to reproduce expected growth parameters for a variety of canonical flows, including Rayleigh-Taylor (RT) and Kelvin-Helmholtz (KH) mixing layers. Both models are then applied in one-dimensional simulation of RT and KH mixing layers, and the expected self-similar growth rates and anisotropy are obtained. Next, models are applied in two-dimensional simulation of the so-called “tilted rocket rig” inclined RT experiment and in simulation of a shock-accelerated localized patch of turbulence. Here it is found that RST is required to capture the qualitative growth of the shock-accelerated patch, and an anisotropic eddy viscosity provides substantial improvement over a Boussinesq treatment for the tilted rocket rig problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Computational shock formation & development: An arbitrary Lagrangian-Eulerian characteristics approach [Slides]

We propose a new computational shock formation-development algorithm. We use a “good” geometry adapted to the evolving solution, along with a good set of variables defined in this geometry. We are able to: (1) Accurately capture the pre-shock; (2) Track distinguished characteristics and capture weak discontinuities; (3) Approximate solutions to classical Riemann problems; and (4) Accurately solve challenging problems for which standard methods fail.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A Foreshock Model for Interstellar Shocks of Solar Origin: Voyager 1 and 2 Observations

The Voyager 1 (V1) and Voyager 2 (V2) spacecraft were launched in 1977 on a mission to explore the outer planets and reach the heliopause, the boundary between the hot solar plasma and the relatively cool interstellar plasma. V1 reached the heliopause on 2012 August 25, at 121.6 au, and V2 reached the heliopause on 2018 November 5, at 119.0 au. One of their remarkable discoveries was the detection of shocks propagating into the interstellar plasma from energetic solar events. These shocks are typically preceded by electron plasma oscillations excited by electron beams streaming along interstellar magnetic field lines ahead of the shocks. The frequencies of the plasma oscillations have now provided radial electron density profiles in the outer heliosphere and in the interstellar medium to radial distances of more than 145 au. The oscillations are typically preceded by bursts of high-energy ∼5–100 MeV electrons. These electron bursts are interpreted as being due to the reflection (and acceleration) of cosmic-ray electrons by the shock at the time the shock first contacts the magnetic field line that passes through the spacecraft. Relative timing between the cosmic rays reflected by the shock and the onset of the plasma oscillations allow us, for the first time, to estimate the energy, ∼20–100 eV, of the electron beams responsible for the plasma oscillations. These observations are combined into a self-consistent model called the foreshock model that describes the interaction of shocks of solar origin with the interstellar plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Physics of the interaction of the solar wind with the ionosphere of Venus - Flow/field models

Venus deflects incident solar wind by means of thermal and magnetic pressure, both of which are produced by ionospheric currents that are in turn driven by pressure gradients and electric fields induced in the ionosphere by the flowing magnetized solar wind plasma. The ionospheric currents connect across the ionopause into the ionosheath region, and close on, and behind, the bow shock. The computed ionospheric plasma velocity distribution is found to be spatially asymmetric, with the fastest flow velocities lying perpendicular to the ionospheric magnetic fields. Computations also indicate that the ionosphere is susceptible to an MHD shear instability whose extent and location is a function of external solar wind conditions. The present Venus ionosphere model compares well with in situ Pioneer Venus Orbiter observations.

Cloutier, P. A.↗
DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗
NASA NTRS1982

Trajectories of solid particles spalled from a carbonaceous heat shield

Trajectories are calculated of solid carbon particles that spall from a carbonaceous heat shield and travel through a given flow field. The mathematical model takes into account mass, momentum, and energy conservation during evaporation of the particles in an effort to understand spallation phenomena and their consequences in a physical way. The solution technique is applied to available Galileo Probe flow field solutions. Plots of trajectories and other particle parameters are presented for a range of initial particle sizes and velocities. It is shown that a significant amount of gaseous carbon is deposited in the inviscid region and ahead of the bow shock as well as in the ablation layer.

Davies, C. B.↗
NASA NTRS1983

Magnetospheric and interplanetary physics 1979-1982

Major trends in the study of magnetospheric and interplanetary physics during the 1979-1982 period are surveyed. Topics discussed include the exploration of the Saturnian and Jovian magnetospheres by Voyagers 1 and 2, the behavior of different ions in the earth magnetosphere, auroral kilometric radiation, computer modeling of global magnetospheric MHD flow, the magnetic substorm, the quiet state, the earth's bow shock, the heliospheric current sheet, and new techniques such as electron beam experiments, 'active' injection experiments, auroral radars, and observations of the earth's distant magnetic tail. The future of this area of research is seen in the combination of data from different spacecraft and ground observations in a single correlated data set, and in the consolidation of past gains by analysis of the large data backlog, while a small number of new missions goes forward.

Stern, D. P.↗
NASA NTRS1982

Jet noise suppression by porous plug nozzles

Jet noise suppression data presented earlier by Maestrello for porous plug nozzles were supplemented by the testing of a family of nozzles having an equivalent throat diameter of 11.77 cm. Two circular reference nozzles and eight plug nozzles having radius ratios of either 0.53 or 0.80 were tested at total pressure ratios of 1.60 to 4.00. Data were taken both with and without a forward motion or coannular flow jet, and some tests were made with a heated jet. Jet thrust was measured. The data were analyzed to show the effects of suppressor geometry on nozzle propulsive efficiency and jet noise. Aerodynamic testing of the nozzles was carried out in order to study the physical features that lead to the noise suppression. The aerodynamic flow phenomena were examined by the use of high speed shadowgraph cinematography, still shadowgraphs, extensive static pressure probe measurements, and two component laser Doppler velocimeter studies. The different measurement techniques correlated well with each other and demonstrated that the porous plug changes the shock cell structure of a standard nozzle into a series of smaller, periodic cell structures without strong shock waves. These structures become smaller in dimension and have reduced pressure variations as either the plug diameter or the porosity is increased, changes that also reduce the jet noise and decrease thrust efficiency.

Bauer, A. B.↗
NASA NTRS1984

Observations of Absorption Lines from Highly Ionized Atoms

In the ultraviolet spectra of hot stars, absorption lines can be seen from highly ionized species in the interstellar medium. Observations of these features which have been very influential in revising the perception of the medium's various physical states, are discussed. The pervasiveness of O 6 absorption lines, coupled with complementary observations of a diffuse background in soft X-rays and EUV radiation, shows that there is an extensive network of low density gas (n approx. fewX 0.001/cucm) existing at coronal temperatures, 5.3 or = log T or = 6.3. Shocks created by supernova explosions or mass loss from early-type stars can propagate freely through space and eventually transfer a large amount of energy to the medium. To create the coronal temperatures, the shocks must have velocities in excess of 150 km/sec; shocks at somewhat lower velocity 9v or = 100 km/sec) can be directly observed in the lines of Si3. Observations of other lines in the ultraviolet, such as Si 4V and C 5, may highlight the widespread presence of energetic uv radiation from very hot, dward stars. More advanced techniques in visible and X-ray astronomical spectroscopy may open up for inspection selected lines from atoms in much higher stages of ionization.

Jenkins, E. B.↗
NASA NTRS2016

Development of Multi-Physics Dynamics Models for High-Frequency Large-Amplitude Structural Response Simulation

An analytic approach is demonstrated to reveal potential pyroshock-driven dynamic effects causing power losses in the Thermo-Electric (TE) module bars of the Mars Science Laboratory (MSL) Multi-Mission Radioisotope Thermoelectric Generator (MMRTG). This study utilizes high-fidelity finite element analysis with SIERRA/PRESTO codes to estimate wave propagation effects due to large-amplitude suddenly-applied pyro shock loads in the MMRTG. A high fidelity model of the TE module bar was created with approximately 30 million degrees-of-freedom (DOF). First, a quasi-static preload was applied on top of the TE module bar, then transient tri-axial acceleration inputs were simultaneously applied on the preloaded module. The applied input acceleration signals were measured during MMRTG shock qualification tests performed at the Jet Propulsion Laboratory. An explicit finite element solver in the SIERRA/PRESTO computational environment, along with a 3000 processor parallel super -computing framework at NASA AMES, was used for the simulation. The simulation results were investigated both qualitatively and quantitatively. The predicted shock wave propagation results provide detailed structural responses throughout the TE module bar, and key insights into the dynamic response (i.e., loads, displacements, accelerations) of critical internal spring/piston compression systems, TE materials, and internal component interfaces in the MMRTG TE module bar. They also provide confidence on the viability of this high-fidelity modeling scheme to accurately predict shock wave propagation patterns within complex structures. This analytic approach is envisioned for modeling shock sensitive hardware susceptible to intense shock environments positioned near shock separation devices in modern space vehicles and systems.

Shock Waves↗
DOE OSTI2020

Insight into the Chemistry of PETN Under Shock Compression Through Ultrafast Broadband Mid-Infrared Absorption Spectroscopy

In this study, thin films of pentaerythritol tetranitrate (PETN) were shock compressed using the laser driven shock apparatus at Los Alamos National Laboratory (LANL). Two spectroscopic probes were available to this apparatus: visible white light transient absorption spectroscopy (VIS) from 400 to 700 nm and mid-infrared transient absorption spectroscopy (MIR) from 1150 to 3800 cm –1 . Important PETN vibrational modes are the symmetric and antisymmetric NO 2 stretches at 1280 and 1650 cm –1 , respectively, as well as CH stretches at ~2900 cm –1 . Shock strength was varied from approximately 3 to 55 GPa to span from the chemically unreactive regime to the regime in which fast chemical reaction took place on the 250 ps time scale of the measurements. VIS and MIR results suggest irreversible chemistry was induced in PETN at pressures above 30 GPa. At lower shock pressures, the spectroscopy showed minimal changes attributable to pressure induced effects. Under the higher-pressure reactive conditions, the frequency region at the antisymmetric NO 2 stretch mode had a significantly increased absorption while the region around the symmetric NO 2 stretch did not. No observable increased absorption occurred in the higher frequency regions where CH-, NH-, and OH- bond absorptions would be observed. A broad absorption appeared on the shoulder at the red-edge of the CO 2 vibrational band around 2200 cm –1 . In addition to the experiments, reactive molecular dynamics were carried out under equivalent shock conditions to correlate the evolution of the infrared spectrum to molecular processes. The simulations show results consistent to experiments up to 30 GPa but suggest that NO and NO 2 related features provided the strongest contributions to the shocked infrared changes. Proposed mechanisms for shocked PETN chemistry are analyzed as consistent or inconsistent with the data presented here. Our experimental data suggests C≡O or N 2 O bond formation, nitrite formation, and absence of significant hydroxyl or amine concentrations in the initial chemistry steps in PETN shocked above 30 GPa.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗