Search NASA⌕ Search

SEARCH · Search NASA

Results for “Single 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Shock tube simulations for the three-layer Richtmyer–Meshkov instability with single-mode and multimode perturbations

While the canonical two-component, single-mode Richtmyer–Meshkov instability (RMI) has been extensively studied, relatively less work has focused on the effects of an additional intermediate-density middle layer. This work investigates such three-material RMI configurations at two Atwood number scenarios using the ARES hydrodynamics code. After validation against previous experimental and computational studies, setups corresponding to recent three-layer shock tube experiments are simulated. Cases with both single-mode and multimode perturbations are studied to quantify mixing across the interface between the materials with highest and intermediate density. In particular, this work is able to comprehensibly examine differences between two- and three-dimensional setups for the single-mode and multimode problems. Observations from previous two-layer investigations still apply in the three-layer setup, but over the time horizons considered, there appears to be insufficient nonlinear mode coupling to create significant differences between two- and three-dimensional simulations following the first passage of a shock. Finally, additional reshock simulations have additional nonlinear growth that does result in expected differences between two- and three-dimensional cases in this three-layer setup, but significant differences do not manifest during the time horizon studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1975

Stark-broadening measurement and calculations for a singly ionized aluminum line

Using an electromagnetically driven shock tube and rapid scanning Fabry-Perot spectrometer, the Lorentzian half-width of the 4663-A line of singly ionized aluminum was measured and found to be 0.45 A plus or minus 20% for an environment corresponding to an electron density of 10 to the 17th power and an electron temperature of 18,500 K. Semiclassical and semiempirical calculations yielded values of 0.56 A and 0.41 A, respectively, while the distorted-wave approximation gave 0.51 A.

Allen, A. W.↗
DOE OSTI2025

Thermal and kinematic properties of ejecta in SN1987A revealed by XRISM

We present an analysis of high-resolution spectra from the shock-heated plasmas in SN 1987A, based on an observation using the Resolve instrument onboard the X-Ray Imaging and Spectroscopy Mission (XRISM). The 1.7–10 keV Resolve spectra are accurately represented by a single-component, plane-parallel shock plasma model, with a temperature of $2.84_{-0.08}^{+0.09}$ keV and an ionization parameter of $2.64_{-0.45}^{+0.58}$ × $10^{11}\,\,{\rm s\,\, cm}^{-3}$. The Resolve spectra are also well reproduced by the 3D magneto-hydrodynamic simulation presented by Orlando et al. (2020, A&A, 636, A22) suggesting substantial contribution from the ejecta. The metal abundances obtained with Resolve align with the Large Magellanic Cloud value, indicating that the X-rays in 2024 originate from “non-metal-rich” shock-heated ejecta and the reverse shock has not reached the inner metal-rich region of ejecta. Doppler widths of the atomic lines from Si, S, and Fe correspond to velocities of 1500–1700 km s$^{-1}$, where the thermal broadening effects in this non-metal-rich plasma are negligible. Therefore, the line broadening seen in Resolve spectra is determined by the large bulk motion of ejecta. For reference, we determined a $90\%$ upper limit on non-thermal emission from a pulsar wind nebula at $4.3 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 2–10 keV range, aligning with NuSTAR findings by Greco et al. (2022, ApJ, 931, 132). Additionally, we searched for the $^{44}$Sc K line feature and found a $1\sigma$ upper limit of $1.0 \times 10^{-6}$ photons cm$^{-2}$ s$^{-1}$, which translates to an initial $^{44}$Ti mass of approximately $2 \times 10^{-4}\, M_{\odot }$, consistent with previous X-ray to soft gamma-ray observations (Boggs et al. 2015, Science, 348, 670; Grebenev et al. 2012, Nature, 490, 373; Leising 2006, ApJ, 651, 1019).

ISM: supernova remnants↗
NASA NTRS2000

Validation of a High-Order Prefactored Compact Scheme on Nonlinear Flows with Complex Geometries

Three benchmark problems are solved using a sixth-order prefactored compact scheme employing an explicit 10th-order filter with optimized fourth-order Runge-Kutta time stepping. The problems solved are the following: (1) propagation of sound waves through a transonic nozzle; (2) shock-sound interaction; and (3) single airfoil gust response. In the first two problems, the spatial accuracy of the scheme is tested on a stretched grid, and the effectiveness of boundary conditions is shown. The solution stability and accuracy near a shock discontinuity is shown as well. Also, 1-D nonlinear characteristic boundary conditions will be evaluated. In the third problem, a nonlinear Euler solver will be used that solves the equations in generalized curvilinear coordinates using the chain rule transformation. This work, continuing earlier work on flat-plate cascades and Joukowski airfoils, will focus mainly on the effect of the grid and boundary conditions on the accuracy of the solution. The grids were generated using a commercially available grid generator, GridPro/az3000.

Hixon, Ray↗
NASA NTRS1994

Explosions of infalling comets in Jupiter's atmosphere

In view of the expected collision of comet Shoemaker-Levy 9 (1993e) with Jupiter in 1994 July, we calculate basic properties of the initial interaction for a simplified Jovian atmosphere. The comet is expected to impact Jupiter at 60 km/sec and at an angle of 45 deg to the zenith. The shock wave generated by the bolide should be optically thick once it has penetrated to an atmospheric density approximately 10(exp -6) gr/cm(exp 3), and we calculate the post-shock conditions assuming local thermodynamic equilibrium (LTE) for shock velocities v(sub sh) in the range 10 to 60 km/sec and preshock densities rho(sub a) = 10(exp -6) to 10(exp -2) gr/cm(exp 3). Our shock calculations include molecular hydrogen, atomic hydrogen, ionized hydrogen, neutral helium, and singly ionized helium. Even at the highest shock velocity, the gas is only partially ionized and the postshock temperature rises with preshock density in order to maintain the ionization. The value of the effective shock adiabatic index gamma(sub sh) varies from 1.17 (at low v(sub sh) and rho(sub a) to 1.40 (at high v(sub sh) and rho(sub a). The ablation rate is limited by the radiative flux that reaches the bolide surface. We argue that the ablated gas does not efficiently transfer its kinetic energy to the atmosphere, and it ultimately slows in a similar fashion to the comet material. As the bolide initially falls through the atmosphere, the character of the shock emission changes. At rho(sub a) approximately 10(exp -8) gr/cm (exp 3), the gas is optically thin and we expect line emission; in the optical spectrum, Balmer emission is expected from the shocked atmosphere and low-ionization metal lines from ablated cometary material. At rho(sub a) approximately 10(exp -6) gr/m(exp 3), the shocked gas is optically thick and the shock front near the bolide produces a blackbody spectrum. The temperature is favorable for ultraviolet (1000 to 3000 A) emission and the luminosity may be approximately 5 x 10(exp 23) ergs/sec for approximately 0.6 sec for a bolide 1 km in radius. At rho(sub a) approximately 10(exp -4) gr/cm(exp 3), the bolide has passed below the ultraviolet photosphere. The shock front emits considerable ionizing radiation, but it is absorbed in a narrow preshock region. The bolometric correction for the optical luminosity is large and we expect a 3000 to 8000 A luminosity of approximately 3 x 10(exp 23) ergs/sec for approximately 1 sec. The optical emission is strongly peaked in the vicinity of the bolide. The bolide does have a somewhat less luminous, optically thick trail extending greater than or equal to 10 km, but the radiation is characterized by a temperature of 4000 to 5000 K. From the fragmentation model of Chyba, Thomas, & Zahnle (1993), the bolide deposits most of its kinetic energy at rho(sub a) approximately 10(exp -3) gr/cm(exp 3) and this is the effective explosion site. The shock wave from such an explosion can move up about one density scale height. We examine the breakout of the shock front from the Jovian atmosphere and find that the shock acceleration in the decreasing density region is slow, so that the energy flux in the shock front is small. Higher velocities might be generated by shock acceleration along the channel left by the bolide if the shock motion can occur before the channel closes off as a result of radiative cooling. Hot gas created by the explosion ultimately rises due to buoyancy on a timescale of a minute. The luminosity is highest when the bubble first rises into the optically thin part of the atmosphere and may be approximately 1 x 10(exp 25) ergs/sec in the near-infrared. Roughly 1% of the initial bolide energy may be radiated in this way; the rest of the energy is lost to sound waves from the initial explosion and to work done by the bubble on the surrounding atmosphere.

Chevalier, Roger A.↗
NASA NTRS1995

Theory and Observations of Non-Thermal Phenomena in Hot Massive Binaries

The shock between the colliding winds in binary systems containing two massive stars accelerates particles to relativistic energies. These energetic particles can produce observable non-thermal radiation from the radio to gamma-rays. The important physical processes in such systems are very similar to those we have proposed for non-thermal emissions from single hot stars, which have shocks generated by instabilities in the radiatively driven stellar winds. This paper discusses the theory and observations of non-thermal radiation in the radio, X-ray, and gamma-ray regions from both single stars and massive binaries. Similarities and differences between the two types of systems are outlined. We discuss two important physical effects that apparently have been neglected in previous theoretical work on colliding wind binaries.

White, Richard L.↗
NASA NTRS1976

Shock-induced fine-grained recrystallization of olivine - Evidence against subsolidus reduction of Fe/2+/

Electron spin resonance (ESR) studies have been carried out on three single grains of terrestrial olivine (Fo90) shock loaded along the 010 line to peak pressures of 280, 330, and 440 kbar. The results indicate that neither metallic Fe similar to that observed in returned lunar soils nor paramagnetic Fe(3+) caused by oxidation of Fe(2+) has been produced in these shock experiments. Trace amounts of Mn (2+) have been detected in both shocked and unshocked olivine. The ESR signals of Mn(2+) show spectral features which are found to correlate with the degree of shock-induced recrystallization observed petrographically. The increasing mass fraction of recrystallized olivine correlates with increasing shock pressures. This phenomenon is modelled assuming it results from the progressive effect of the shock-induced transformation of the olivine to a yet unknown high-pressure phase and its subsequent reversion to the low-pressure olivine phase. The mass fraction of recrystallized material is predicted to be nearly linear with shock pressure.

Ahrens, T. J.↗
NASA NTRS1987

Shock structure in jets of arbitrary exit geometry

The shock structure and hydrodynamic stability of single and twin jets with arbitrary exit geometry are studied. Both vortex sheets and realistic mean profiles are considered for the jet shear layer. The boundary element method is used to predict the shock spacing and screech tones in a vortex sheet model of a single jet. Results are given for circular, elliptic and rectangular jets and the results compared with analysis and experiment. A finite difference technique is used to examine the effect of finite mixing layer thickness for a single jet. The hydrodynamic stability of twin circular jets is also considered. The variation of the stability characteristics are presented for the vortex sheet case for several separation distances. The formulation of the problem to account for finite thickness effects is also given.

Morris, P. J.↗
NASA NTRS1973

A shock-induced phase change in orthoclase.

New shock compression data to 340 kb for single-crystal orthoclase (along (001)), demonstrate the onset of a shock-induced phase change at about 115 kb. Along the Hugoniot a mixed-phase region extends to about 300 kb, above which the data are believed to correspond to the properties of a high-pressure phase having the hollandite structure (zero pressure density of 3.84 g/cu cm) reported by Ringwood et al. If the hollandite value for the zero pressure density is used, the zero pressure bulk modulus of this phase is approximately 1.8 plus or minus 0.2 Mb.

Ahrens, T. J.↗
NASA NTRS1991

Periodic oscillations of the quasi-parallel bow shock as observed with the AMPTE/CCE spacecraft

The AMPTE/CCE spacecraft has, on occasion, encountered the earth's bow shock near apogee (about 9 earth radii). During one encounter, on November 1, 1984, the bow shock was observed to pass back and forth over the spacecraft several times in a periodic fashion. The period of the oscillation was about 20 seconds, and an upstream wave with a similar period was observed in the solar wind just prior to the multiple shock crossings. Minimum variance analysis of the shock crossings shows that the magnetic field carried out a single rotation through 360 deg for each pair of shock crossings. This implies that the shock was responding to the rotation of the field in the upstream wave, and it is shown how the plasma displacement associated with the upstream wave can produce the observed oscillations. These data may provide support for the wavelike structure of quasi-parallel shocks predicted in simulations.

Strangeway, R. J.↗
DOE OSTI2024

Shock compression of crystalline TeO 2 to the high-pressure fluid regime: Insights from ab initio molecular dynamics simulations

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ⁠) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α-TeO 2 paratellurite (⁠P4 1 2 1 2), TeO 2 cotunnite (⁠Pnma⁠), and TeO 2 post-cotunnite (⁠P2 1 /m⁠). Based on the discontinuity found in the calculated U s – u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of TeO 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α- TeO 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Full-field quantitative visualization of shock-driven pore collapse and failure modes in PMMA

The dynamic collapse of pores under shock loading is thought to be directly related to hot spot generation and material failure, which is critical to the performance of porous energetic and structural materials. However, the shock compression response of porous materials at the local, individual pore scale is not well understood. This study examines, quantitatively, the collapse phenomenon of a single spherical void in PMMA at shock stresses ranging from 0.4 to 1.0 GPa. Using a newly developed internal digital image correlation technique in conjunction with plate impact experiments, full-field quantitative deformation measurements are conducted in the material surrounding the collapsing pore for the first time. The experimental results reveal two failure mode transitions as shock stress is increased: (i) the first in situ evidence of shear localization via adiabatic shear banding and (ii) dynamic fracture initiation at the pore surface. Numerical simulations using thermo-viscoplastic dynamic finite element analysis provide insights into the formation of adiabatic shear bands (ASBs) and stresses at which failure mode transitions occur. Further numerical and theoretical modeling indicates the dynamic fracture to occur along the weakened material inside an adiabatic shear band. Finally, analysis of the evolution of pore asymmetry and models for ASB spacing elucidate the mechanisms for the shear band initiation sites, and elastostatic theory explains the experimentally observed ASB and fracture paths based on the directions of maximum shear.

42 ENGINEERING↗
NASA NTRS1995

Theoretical position-velocity diagrams of flux, electron density, and electron temperature in Herbig-Haro objects

Many features of Herbig-Haro objects can be reproduced using a kinematical bow shock model. We use the model to generate position-velocity (PV) diagrams of flux in H-alpha (O I) lambda lambda 6300+63, (S II) lambda lambda 6716+31, (O III) lambda lambda 4959+5007, and (C I) lambda lambda 9823+50, line ratios of (O I)/H-alpha, (O I)/(S II), (S II)/H-alpha, H-alpha(S II), and (O III)/H-alpha, electron density N(sub e), and electron temperature T(sub e). We show how position-velocity diagrams of N(sub e) and flux vary with shock velocity. By matching the diagrams from single lines, the line ratios, and N(sub e) with observations, we determined a narrow range of shock parameters for HH 1F, 2(A' + H), and 43 (B + C). We model the N(sub e) features of HH 2(A' + H) as a superposition of two bowshocks. We also show that the effects of slight misalignments of the two diagrams to be divided can produce artifacts in the line ratios and N(sub e) which obliterate the physical features. We show that N(sub e) in HH 1 can only be explained using the kinematical model by taking these misalignments into account.

Indebetouw, Remy↗
DOE OSTI2024

Robust 3D multi-material hydrodynamics using discontinuous Galerkin methods

A high-order discontinuous Galerkin (DG) method is presented for nonequilibrium multi-material (m ≥ 2) flow with sharp interfaces. Material interfaces are reconstructed using the algebraic THINC approach, resulting in a sharp interface resolution. The system assumes stiff velocity relaxation and pressure nonequilibrium. The presented DG method uses Dubiner's orthogonal basis functions on tetrahedral elements. This results in a unique combination of sharp multimaterial interfaces and high-order accurate solutions in smooth single-material regions. A novel shock indicator based on the interface conservation condition is introduced to mark regions with discontinuities. Slope limiting techniques are applied only in these regions so that nonphysical oscillations are eliminated while maintaining high-order accuracy in smooth regions. A local projection is applied on the limited solution to ensure discrete closure law preservation. The effectiveness of this novel limiting strategy is demonstrated for complex three-dimensional multi-material problems, where robustness of the method is critical. The presented numerical problems demonstrate that more accurate and efficient multi-material solutions can be obtained by the DG method, as compared to second-order finite volume methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1970

The aerodynamic characteristics of large angled cones with retrorockets

Analytical and experimental phases of the subject investigation are described. The analytical program for the single jet determines the terminal shock location, the jet boundary, the interface profile, the bow shock profile, the shear layer growth and the dead air region pressure. The experimental program described was conducted over the range from free stream Mach 0.4 to 2.0 at angles-of-attack up to 18 deg and at thrusting coefficients up to C sub T = T/q sub infinity A sub m = 30. Variables investigated included aeroshell angle, number of nozzles, engine thrust, size of nozzles, nozzle throttling and gas composition. The influence of these variables on the aeroshell stability, drag, and loads was determined by integrating pressure measurements on the aeroshell. The total system forces consist of components due to pure thrust and components due to pressure on the aeroshell arising from the jet-free stream interaction. Shadowgraphs provided flow field geometries which proved to be within 10% of those predicted analytically.

Jarvinen, P. O.↗
NASA NTRS1973

Shock-wave structure in a fully ionized gas.

The structure of a steady plane shock in a fully ionized gas has been investigated using the Boltzmann equation with the Gross-Krook type model as the governing equation and the discrete-ordinate method as a tool. The present results agree well with the results obtained by the continuum approach for weak shocks. For strong shocks the present approach gives results that are considerably different from those yielded by the continuum approach, particularly in the high pressure region. Also, the present method gives smooth and continuous shock profiles for high Mach numbers. On the other hand, the results from the continuum approach contain discontinuities in slope. The ion and electron distribution functions across the shock are found to be singly peaked for all cases considered and the effects of the induced electric field is small.

Lu, C. S.↗
NASA NTRS1979

Anorthite: Thermal equation of state to high pressures

The shock wave (Hugoniot) data on single crystal and porous anorthite (CaAl2Si208) to pressures of 120 GPa are presented. These data are inverted to yield high pressure values of the Grueneisen parameter, adiabatic bulk modulus, and coefficient of thermal expansion over a broad range of pressures and temperatures which in turn are used to reduce the raw Hugoniot data and construct an experimentally based, high pressure thermal equation of state for anorthite. The hypothesis that higher order anharmonic contributions to the thermal properties decrease more rapidly upon compression than the lowest order anharmonicities is supported. The properties of anorthite corrected to lower mantle conditions show that although the density of anorthite is comparable to that of the lower most mantle, its bulk modulus is considerably less, hence making enrichment in the mantle implausible except perhaps near its base.

Jeanloz, R.↗