Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shock initiation”

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

Tables and charts of equilibrium normal shock and shock-tube solutions for helium-hydrogen mixtures with velocities to 70 km/sec

Equilibrium thermodynamic and flow properties are presented in tabulated and graphical form for moving, standing, and reflected normal shock waves into helium-hydrogen mixtures representative of proposed outer planet atmospheres. The volumetric compositions of these mixtures are 0.35He-0.65H2, 0.20He-0.80H2, and 0.05He-0.95H2. Properties include pressure, temperature, density, enthalpy, speed of sound, entropy, molecular-weight ratio, isentropic exponent, velocity, and species mole fractions. Incident (moving) shock velocities are varied from 4 to 70 km/sec for a range of initial pressure of 5 N/sq m to 100 kN/sq m. The present results are applicable to shock-tube flows and to free-flight conditions for a blunt body at high velocities. A working chart illustrating idealized shock-tube performance with a 0.20He-0.80H2 test gas and heated helium driver gas is also presented.

Miller, C. G., III↗

Hypersonic shock structure with Burnett terms in the viscous stress and heat flux

The continuum Navier-Stokes and Burnett equations are solved for one-dimensional shock structure in various monatomic gases. A new numerical method is employed which utilizes the complete time-dependent continuum equations and obtains the steady-state shock structure by allowing the system to relax from arbitrary initial conditions. Included is discussion of numerical difficulties encountered when solving the Burnett equations. Continuum solutions are compared to those obtained utilizing the Direct Simulation Monte Carlo method. Shock solutions are obtained for a hard sphere gas and for argon from Mach 1.3 to Mach 50. Solutions for a Maxwellian gas are obtained from Mach 1.3 to Mach 3.8. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Shock density thickness, density asymmetry, and density-temperature separation are all more accurately predicted by the Burnett equations than by the Navier-Stokes equations.

Chapman, Dean R.↗

Scaling craters in carbonates: Electron paramagnetic resonance analysis of shock damage

Carbonate samples from the 8.9-Mt nuclear (near-surface explosion) crater, OAK, and a terrestrial impact crater, Meteor Crater, were analyzed for shock damage using electron paramagnetic resonance (EPR). Samples from below the OAK apparent crater floor were obtained from six boreholes, as well as ejecta recovered from the crater floor. The degree of shock damage in the carbonate material was assessed by comparing the sample spectra to the spectra of Solenhofen and Kaibab limestone, which had been skocked to known pressures. Analysis of the OAK Crater borehole samples has identified a thin zone of allocthonous highly shocked (10-13 GPa) carbonate material underneath the apparent crater floor. This approx. 5- to 15-m-thick zone occurs at a maximum depth of approx. 125 m below current seafloor at the borehole, sited at the initial position of the OAK explosive, and decreases in depth towards the apparent crater edge. Because this zone of allocthonous shocked rock delineates deformed rock below, and a breccia of mobilized sand and collapse debris above, it appears to outline the transient crater. The transient crater volume inferred in this way is found to by 3.2 +/- 0.2 times 10(exp 6)cu m, which is in good agreement with a volume of 5.3 times 10(exp 6)cu m inferred from gravity scaling of laboratory experiments. A layer of highly shocked material is also found near the surface outside the crater. The latter material could represent a fallout ejecta layer. The ejecta boulders recovered from the present crater floor experienced a range of shock pressures from approx. 0 to 15 GPa with the more heavily shocked samples all occurring between radii of 360 and approx. 600 m. Moreover, the fossil content, lithology and Sr isotopic composition all demonstrate that the initial position of the bulk of the heavily shocked rock ejecta sampled was originally near surface rock at initial depths in the 32 to 45-m depth (below sea level) range. The EPR technique is also sensitive to prehistoric shock damage. This is demonstrated by our study of shocked Kaibab limestone from the 49,000-year-old Meteor (Barringer) Crater Arizona.

Polanskey, Carol A.↗

A General Relativistic Magnetohydrodynamics Simulation of Jet Formation with a State Transition

We have performed the first fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity sim 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving shock wave. This evolution of jet-disk coupling suggests that the low/hard state of the jet system may switch to the high/soft state with a wind, as the accretion rate diminishes.

Nishikawa, K. I.↗

Effects of interplanetary shock waves on energetic charged particles

Experimental data on the influence of interplanetary perpendicular and oblique shock waves on the ambient energetic protons are presented along with a theoretical analysis of the acceleration of particles in almost perpendicular shock waves. It was found that low-energy protons can be accelerated in perpendicular shock waves by repeated crossings of the shock front up to a maximum energy given by the product of their initial energy times the ratio of the magnetic fields. High-energy protons need to stay at the shock front for longer times than low-energy protons in order to reach the same relative energy gain. In the theoretical study of proton acceleration at almost perpendicular shock waves, it was found that protons reflected at shock waves with the angle between the upstream magnetic field and the shock normal greater than about 80 deg achieve large energy gains at the shock front. The larger this angle, the higher the energy gain. However, the reflection and energization of protons at these shock waves is not 'instantaneous', neither is it a one-step process: it is performed through repeated crossings of the shock front.

Sarris, E. T.↗

Observational discrimination between modes of shock propagation in interstellar clouds: Predictions of CH+ and SH+ column densities in diffuse clouds

Considerable effort in recent years has been devoted to the study of shocks in the diffuse interstellar medium. This work has been motivated partly by the observations of rotationally excited states of H2, and partly by the realization that species such as CH(+), OH and H2O might be formed preferentially in hot, post-shock gas. The problem of CH(+) and the difficulties encountered when trying to explain the high column densities, observed along lines of sight to certain hot stars, have been reviewed earlier. The importance of a transverse magnetic field on the structure of an interstellar shock was also demonstrated earlier. Transverse magnetic fields above a critical strength give rise to an acceleration zone or precursor, in which the parameters on the flow vary continuously. Chemical reactions, which change the degree of ionization of the gas, also modify the structure of the shock considerably. Recent work has shown that large column densities of CH(+) can be produced in magnetohydrodynamic shock models. Shock speeds U sub s approx. = 10 km/s and initial magnetic field strengths of a few micro G are sufficient to produce ion-neutral drift velocities which can drive the endothermic C(+)(H2,H)CH(+) reaction. It was also shown that single-fluid hydrodynamic models do not generate sufficiently large column densities of CH(+) unless unacceptably high shock velocities (u sub s approx. 20 km/s) are assumed in the models. Thus, the observed column densities of CH(+) provide a constraint on the mode of shock propagation in diffuse clouds. More precisely, they determine a lower limit to the ion-neutral drift velocity.

Flower, D. R.↗

Acceleration of charged particles in oblique MHD shocks

The motion of high energy charged particles in ideal oblique MHD shocks has been studied extensively. The orbits of charged particles can be solved exactly from Lorentz force equation. It is found that a charged particle may cross the shock front many times before it leaves the shock and it is energized during crossings. The most energetic particles are those with initial pitch and phase near the boundary between crossing and reflected particles. The reflected particles cross the shock about twice as many times as the passing particles do. The energy gain increases for larger values of shock strength and speed. If reflection occurs, the energy gain is larger for smaller angles between the field and shock surface. The introduction of scattering results in a few more higher energy particles.

Chen, G.↗

Monte Carlo simulation of particle acceleration at astrophysical shocks

A Monte Carlo code was developed for the simulation of particle acceleration at astrophysical shocks. The code is implemented in Turbo Pascal on a PC. It is modularized and structured in such a way that modification and maintenance are relatively painless. Monte Carlo simulations of particle acceleration at shocks follow the trajectories of individual particles as they scatter repeatedly across the shock front, gaining energy with each crossing. The particles are assumed to scatter from magnetohydrodynamic (MHD) turbulence on both sides of the shock. A scattering law is used which is related to the assumed form of the turbulence, and the particle and shock parameters. High energy cosmic ray spectra derived from Monte Carlo simulations have observed power law behavior just as the spectra derived from analytic calculations based on a diffusion equation. This high energy behavior is not sensitive to the scattering law used. In contrast with Monte Carlo calculations diffusive calculations rely on the initial injection of supra-thermal particles into the shock environment. Monte Carlo simulations are the only known way to describe the extraction of particles directly from the thermal pool. This was the triumph of the Monte Carlo approach. The question of acceleration efficiency is an important one in the shock acceleration game. The efficiency of shock waves efficient to account for the observed flux of high energy galactic cosmic rays was examined. The efficiency of the acceleration process depends on the thermal particle pick-up and hence the low energy scattering in detail. One of the goals is the self-consistent derivation of the accelerated particle spectra and the MHD turbulence spectra. Presumably the upstream turbulence, which scatters the particles so they can be accelerated, is excited by the streaming accelerated particles and the needed downstream turbulence is convected from the upstream region. The present code is to be modified to include a better description of particle scattering (pitch-angle instead of hard-sphere) and as iterative procedure for treating the self-excitation of the MHD turbulence.

Campbell, Roy K.↗

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gas dynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel. Initial conditions are prescribed by assuming homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. The results suggest possible criteria for assessing the consequences of passage rotation on the wave processes, and they may therefore be applicable to pressure-exchange wave rotors. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface as well as the hub-to-tip radius ratio. Rarefaction fronts, shocks, and contact interfaces are observed to propagate faster with increasing wheel Mach number.

Larosiliere, L. M.↗

Analysis of unsteady wave processes in a rotating channel

The impact of passage rotation on the gas dynamic wave processes is analyzed through a numerical simulation of ideal shock-tube flow in a closed rotating-channel. Initial conditions are prescribed by assuming homentropic solid-body rotation. Relevant parameters of the problem such as wheel Mach number, hub-to-tip radius ratio, length-to-tip radius ratio, diaphragm temperature ratio, and diaphragm pressure ratio are varied. The results suggest possible criteria for assessing the consequences of passage rotation on the wave processes, and they may therefore be applicable to pressure-exchange wave rotors. It is shown that for a fixed geometry and initial conditions, the contact interface acquires a distorted three-dimensional time-dependent orientation at non-zero wheel Mach numbers. At a fixed wheel Mach number, the level of distortion depends primarily on the density ratio across the interface as well as the hub-to-tip radius ratio. Rarefaction fronts, shocks, and contact interfaces are observed to propagate faster with increasing wheel Mach number.

Larosiliere, Louis M.↗

Ionization fronts and shocked flows - The structure of the Orion Nebula at 0.1 arcsec

Initial results are presented from images of the Orion Nebula obtained with the HST Wide-Field Camera with a resolution of better than about 10 to the 15th cm. The morphology of the forbidden S II emission is markedly different from that seen in H-beta and forbidden O II. The forbidden S II is dominated by filamentary features with widths between 0.1 and 1 arcsesc which sharply highlight ionization fronts moving into dense neutral material. While the surfaces into which these fronts are moving are textured, subarcsec clumps with high density contrast are uncommon. The fields contains a number of HH objects whose structure is revealed by the HST, and which are described.

Hester, J. Jeff↗

Mechanical Properties of Shock-Damaged Rocks

Stress-strain tests were performed both on shock-damaged gabbro and limestone. The effective Young's modulus decreases with increasing initial damage parameter value, and an apparent work-softening process occurs prior to failure. To further characterize shock-induced microcracks, the longitudinal elastic wave velocity behavior of shock-damaged gabbro in the direction of compression up to failure was measured using an acoustic transmission technique under uniaxial loading. A dramatic increase in velocity was observed for the static compressive stress range of 0-50 MPa. Above that stress range, the velocity behavior of lightly damaged (D(sub 0) less than 0.1) gabbro is almost equal to unshocked gabbro. The failure strength of heavily-damaged (D(sub 0) greater than 0.1) gabbro is approx. 100-150 MPa, much lower than that of lightly damaged and unshocked gabbros (approx. 230-260 MPa). Following Nur's theory, the crack shape distribution was analyzed. The shock-induced cracks in gabbro appear to be largely thin penny-shaped cracks with c/a values below 5 x 10(exp -4). Moreover, the applicability of Ashby and Sammis's theory relating failure strength and damage parameter of shock-damaged rocks was examined and was found to yield a good estimate of the relation of shock-induced deficit in elastic modulus with the deficit in compressive strength.

He, Hongliang↗

On the fragmentation of protogalactic clouds

After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature of about 1 million K. This temperature lies near a very unstable region of the cooling curve, so that density enhancements will be expected to grow rapidly as the protogalaxy cools to lower temperatures. In this contribution, the dynamical effects on the growth of perturbations are examined. It is shown that perturbations with large initial amplitudes lead to the formation of H2, which results in rapid cooling to 100 K, and the formation of population III stars. Perturbations with small initial amplitudes, however, do not cool significantly ahead of the background, do not develop strong contrasts with the background, and so cannot induce the formtion of H2. They thus remain near 10,000 K, forming a population of protoclusters with masses of about 10 to the 7th solar masses, somewhat larger than typical globular cluster masses. Contamination by metals from a primordial generation of stars may lead to cooling and a significant reduction of the Jeans mass. Photodissociation induced by a background source of hard photons can also prevent cooling below 10,000 K, enhancing the formation of protoclusters relative to stars.

Murray, Stephen D.↗

Effects of Coupled Rolling and Pitching Oscillations on Transonic Shock-Induced Vortex-Breakdown Flow of a Delta Wing

Unsteady, transonic vortex-breakdown flow over a 65 deg. sharp edged, cropped-delta wing of zero thickness undergoing forced coupled pitching and rolling oscillations is investigated computationally. The initial condition of the flow is characterized by a transverse terminating shock which induces of the leading edge vortex cores to breakdown. The computational investigation uses the time-accurate solution of the laminar, unsteady, compressible, full Navier-Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme. The main focus is to analyze the effects of coupled motion on the wing response and vortex-breakdown flow by varying oscillation frequency and phase angle while keeping the maximum pitch and roll amplitude equal.

Kandil, Osama A.↗

Oxygen chemistry of shocked interstellar clouds. III - Sulfur and oxygen species in dense clouds

The chemical evolution of oxygen and sulfur species in shocked dense clouds is studied. Reaction rate constants for several important neutral reactions are examined, and revised values are suggested. The one-fluid magnetohydrodynamic shock structure and postshock chemical evolution are calculated for shocks of velocity v(s) = 10 km/s through clouds of initial number density n(0) = 100,000/cu cm and of molecule/atom ratios H2/H = 10, 1000, and 100,000 with most sulfur contained initially in molecules SO2 and SO. Abundances of SO2, SO, CS, and OCS remain near their preshock values, except in clouds containing substantial amounts of atomic hydrogen, where significant destruction of sulfur-oxygen species occurs. Abundances of shock-enhanced molecules HS and H2O are sensitive to the molecule/atom ratio. Nonthermal oxygen-hydrogen chemistry has a minor effect on oxygen-sulfur molecules in the case H2/H = 10.

Leen, T. M.↗

Stellar Collisions

Head-on collision between stars at high initial relative velocity, considering one-dimensional driven shock motion in inhomogeneous medium

SHOCK WAVE↗