Search NASA⌕ Search

SEARCH · Search NASA

Results for “SHOCK WAVE PROFILE”

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 55 records · Page 3

Fundamental problem in computing radiating flow fields with thick shock waves

Possible reasons for the failure of the Navier-Stokes equations to yield realistic profiles of the temperature and density through the structure of a hypersonic shock wave are investigated. Models for bulk viscosity in a monatomic gas are examined which yield a realistic thickness for the shock-wave density profile, but not the temperature profile, and hence are not satisfactory. A tentative computational model for nitrogen is explored which yields considerably more realistic results than the Navier-Stokes equations. This model involves a nonlinear stress-strain tensor, nonlinear heat flux vector, and nonequilibrium rotational energy.

Chapman, Dean R.↗

Laser Beam Propagation Through Inhomogeneous Media with Shock-Like Profiles: Modeling and Computing

Wave propagation in inhomogeneous media has been studied for such diverse applications as propagation of radiowaves in atmosphere, light propagation through thin films and in inhomogeneous waveguides, flow visualization, and others. In recent years an increased interest has been developed in wave propagation through shocks in supersonic flows. Results of experiments conducted in the past few years has shown such interesting phenomena as a laser beam splitting and spreading. The paper describes a model constructed to propagate a laser beam through shock-like inhomogeneous media. Numerical techniques are presented to compute the beam through such media. The results of computation are presented, discussed, and compared with experimental data.

Adamovsky, Grigory↗

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗

Shock waves in collisionless plasmas.

Theoretical models for collisionless plasma shock waves in terms of nonlinear, magnetosonic, constant profile waves, turbulent- and electrostatic-shock structures, etc

Tidman, D. A.↗

A method of atmospheric density measurements during Shuttle entry using UV laser Rayleigh scattering

A detailed study is described of the performance capabilities and the hardware requirements for a method in which ambient density is measured along the Space Shuttle flight path using on-board optical instrumentation. The technique relies on Rayleigh scattering of light from a pulsed, ultraviolet, ArF excimer laser operating at a wavelength of 193 nm. The method is shown to be capable of providing direct measurements of ambient density with an uncertainty of less than 1 percent and with a spatial resolution of 1 km, over an altitude range from 50 to 90 km. In addition, extensions of this concept are discussed that allow measurements of the shock wave location and the density profile within the shock layer. Two approaches are identified that appear to be feasible, in which the same laser system is used for the extended measurements as that required for the ambient density measurements.

Mckenzie, Robert L.↗

Generation, saturation, and convection of electrostatic waves in Jupiter's shock foot

In this paper, a model is developed for the analysis of the electrostatic waves produced in the shock foot at Jupiter. It is shown that an ion beam instability involving the ions reflected at the shock ramp and the incoming solar-wind electrons produces waves at the observed frequencies and that saturation via orbit diffusion limits the waves to amplitudes near to what is observed. Results from a two-dimensional model of the reflected beam in the foot indicate that the waves propagate against the solar wind away from the shock ramp and are amplified up to their saturation amplitudes. The saturation results, combined with the electron temperature profile due to wave-particle interactions predicted by quasi-linear theory, reproduce a wave amplitude profile for the shock foot that is in reasonable agreement with the observations.

Moses, S. L.↗

Shock-tube measurements of rotational relaxation in hydrogen.

Density profiles in hydrogen shock waves have been measured with an electron beam densitometer for Mach numbers from 1.2 to 3.0. The translational and rotational relaxation processes were found to be uncoupled in the shock wave allowing for a simple determination of the relatively large (from approximately 270 to approximately 30) number of collisions required for equilibration of the rotational degrees of freedom over the temperature range from 410-1240 K. An intermediate relaxation point and subsequent increase in the number of collisions required above 600 K was attributed to the participation of additional excited rotational quantum levels in the equilibration process.

Boitnott, C. A.↗

Molecular dynamics simulation of a piston driven shock wave in a hard sphere gas

Molecular dynamics simulation is used to study the piston driven shock wave at Mach 1.5, 3, and 10. A shock tube, whose shape is a circular cylinder, is filled with hard sphere molecules having a Maxwellian thermal velocity distribution and zero mean velocity. The piston moves and a shock wave is generated. All collisions are specular, including those between the molecules and the computational boundaries, so that the shock development is entirely causal, with no imposed statistics. The structure of the generated shock is examined in detail, and the wave speed; profiles of density, velocity, and temperature; and shock thickness are determined. The results are compared with published results of other methods, especially the direct simulation Monte-Carlo method. Property profiles are similar to those generated by direct simulation Monte-Carlo method. The shock wave thicknesses are smaller than the direct simulation Monte-Carlo results, but larger than those of the other methods. Simulation of a shock wave, which is one-dimensional, is a severe test of the molecular dynamics method, which is always three-dimensional. A major challenge of the thesis is to examine the capability of the molecular dynamics methods by choosing a difficult task.

Woo, Myeung-Jouh↗

A cost-effective strategy for nonoscillatory convection without clipping

Clipping of narrow extrema and distortion of smooth profiles is a well known problem associated with so-called high resolution nonoscillatory convection schemes. A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost effective because the wider stencil is used only where needed-in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗

Sound speed measurements in lithium fluoride single crystals shock compressed to 168 GPa along [100]

The shock wave response of [100] lithium fluoride (LiF) single crystals at high stresses is of long-standing interest due to their extensive use as optical windows in dynamic compression experiments. The report of melting in shock compressed LiF single crystals between 134 GPa and 152 GPa – based on a single sound speed datum [J. Appl. Phys. 117, 045901 (2015)] – was surprising because good optical transmission was previously demonstrated in LiF shock compressed to ~200 GPa [J. Appl. Phys. 116, 033515 (2014)]. To address these apparent differences, we report on plate impact experiments on [100] LiF single crystals shock compressed to 168 GPa. Wave profiles were measured using laser interferometry to determine Hugoniot states and longitudinal sound speeds in shock compressed LiF. The measured Hugoniot states are in good agreement with those measured in previous studies. However, the measured sound speeds presented here show no evidence of melting up to 168 GPa. In particular, the abrupt drop reported previously in sound speed at 152 GPa was not observed in the present work. Furthermore, our results establish a lower bound of 168 GPa for the onset of melting in shock compressed LiF single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shock Response of Polymethyl Methacrylate (PMMA) Under Explosive Loading

Abstract The dynamic response of polymethyl methacrylate (PMMA) is well understood for one-dimensional planar impact shocks, but limited research has been performed on the response of PMMA under spherical shock loading. In this work, the shock decay of an explosively-driven shock wave into PMMA was experimentally measured. PMMA cubes of various geometries were explosively loaded with an RP-80 detonator to produce the explosive shock wave. High-speed schlieren imaging was implemented to measure the explosively-driven shock wave velocity throughout the PMMA cubes. Photon Doppler velocimetry (PDV) was used to measure the particle velocity imparted by the shock wave at the surface of the cubes. The material shock response was studied at distances from 21.91 to 133.3 mm from the explosive source. The particle velocity history measured by PDV was compared to the wave profile visualized in the high-speed images. The shock wave pulse amplitude decreased with increased distance from the source. The conducted experiments extend the PMMA shock Hugoniot relating to the lower shock and particle velocity regime.

Torres, S. M. (ORCID:0009000287177206)↗

Acceleration of energetic particles by shock waves from large solar flares

Time-intensity profiles for solar proton events are analyzed with respect to the acceleration of energetic particles from the ambient solar wind by an interplanetary shock. The time-intensity profiles are derived from the Helio 1 particle data base for 1974 to 1984. It is found that many of the large proton events display a prolonged plateau. The Fe/O abundances in the large proton events are studied. It is observed that initially the Fe/O ratio is high (about 1.0), but after proton intensities increase to about 100 protons/sq cm sr sec MeV the Fe/O ratio decreases to a value near 0.1.

Reames, Donald V.↗

Sharp monotonic resolution of discontinuities without clipping of narrow extrema

A strategy is presented for accurately simulating highly convective flows containing discontinuities such as density fronts or shock waves, without distorting smooth profiles or clipping narrow local extrema. The convection algorithm is based on non-artificially-diffusive third-order upwinding in smooth regions, with automatic adaptive stencil expansion to (in principle, arbitrarily) higher order upwinding locally, in regions of rapidly changing gradients. This is highly cost-effective because the wider stencil is used only where needed - in isolated narrow regions. A recently developed universal limiter assures sharp monotonic resolution of discontinuities without introducing artificial diffusion or numerical compression. An adaptive discriminator is constructed to distinguish between spurious overshoots and physical peaks; this automatically relaxes the limiter near local turning points, thereby avoiding loss of resolution in narrow extrema. Examples are given for one-dimensional pure convection of scalar profiles at constant velocity.

Leonard, B. P.↗

A theoretical interpretation of Forbush decreases

An analysis is presented of the intensity variations of relativistic cosmic rays, at 1 AU, which are associated with the passage of a flare shock wave. The magnitude and time profile obtained is similar to that observed in flare-induced Forbush decreases, and the principal cause of the energy reduction is an increase in adiabatic cooling of the arriving particles due to prolonged containment behind the compressed field of the flare shock wave. The large decreases calculated for the smaller diffusion coefficients are inconsistent with observation, implying larger mean free paths than those traditionally assumed at these rigidities. The method presented may be used to study intensity variations associated with other large-scale structures of the interplanetary field, such as corotating interaction regions.

Gall, R.↗

Shock-induced chemical decomposition and overdriven detonation in hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) single crystals

Understanding the differences in the shock compression and detonation response of insensitive high explosives (IHEs) and conventional HEs (CHEs) is a long-standing need in HE science and technology. Having previously examined 1,1-diamino-2,2-dinitroethene (FOX-7) IHE single crystals, the shock and detonation response of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX)—a widely used CHE—was determined using wave profile measurements in ~250 μm thick single crystals shock compressed to 63 GPa. In marked contrast to FOX-7, RDX single crystals shocked along the [100] and [111] orientations showed wave profile features consistent with chemical decomposition onset at 15 GPa. These features were more pronounced for [100] RDX, suggesting a higher decomposition rate compared to [111] RDX. At 51 GPa and above, flat-topped single waves were observed for both orientations, showing the classic Chapman–Jouget (C–J) detonation response in which the decomposition is completed within the detonation front. The Hugoniot states and sound speeds determined for the detonation products were similar for both [100] and [111] orientations, showing that the overdriven detonation response for RDX does not depend on crystal orientation. The C–J pressure for RDX single crystals (35 GPa)—determined experimentally—is comparable to that of FOX-7. However, compared to FOX-7, chemical decomposition onset for RDX occurs at much lower pressures and the overdriven C–J detonation response occurs at higher pressures. Lastly, the present findings constitute the first experimental comparison of the shock and detonation response of conventional and insensitive HE single crystals over a broad pressure range below and above the C–J pressure.

36 MATERIALS SCIENCE↗

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗