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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

From regolith to rock by shock

The purpose of this study is to relate the various mechanisms by which material is shock-lithified in terrestrial analogs of the lunar regolith to specific conditions of cratering, thereby making more specific the possible conditions of formation of rock from regolith by shock processes on the lunar surface. A model for shock-lithification of terrestrial and lunar regolith is proposed in which air or an air-water mixture initially in the pores of terrestrial soil affects the behavior of a soil-air-water system under shock-loading. Shock compression of porous terrestrial regolith by relatively small impact events give rise to three pressure regimes: (1) regime 1 - at pressures below 100 kb, material is compacted and weakly shock-lithified; (2) regime 2 - at pressures between 100 and 200 kb, material may be fragmented if the induced pore pressure exceeds the strength of the weak lithification mechanisms; and (3) at pressure above 200 kb, material is strongly lithified but may greatly expand in volume due to the pressure of pore gases. The three other regimes (below 50 kb, between 50 and 100 kb, and above 100 kb) associated with shock compression of lunar regolith are identified and discussed.

Kieffer, S. W.↗
NASA NTRS2022

An Investigation Into Some Important Aspects of Droplet Breakup/Vaporization Behavior Caused By a Gas Flow (Both Nonreacting & Detonative Combustion) in a Shock Tube

As a part of the rotating detonation engine (RDE) technology enablement project at NASA Glenn Research Center (GRC), an effort was undertaken to extend our current computational capabilities of OpenNCC in some important ways with the implementation of a modeling approach to account for the droplet breakup caused by a shock-induced gas motion, & a vaporization model valid over a wide range of pressure conditions encountered in multiphase detonation. With the modified code, a study was undertaken to investigate the individual droplet behavior followed by the passage of a shock front. The study is carried out by tracking a sparse group of droplets to gain some understanding of shock induced droplet behavior under various shock strengths & fuel injector conditions. The study also looks into the effect of randomization involved in determining the droplet breakup outcomes. Over a wide range of sub-critical conditions examined, larger droplets are observed to undergo significant changes in droplet behavior following their breakup. However, smaller droplets (10 µm or less ) remain unaffected by any shock induced breakup. In a follow-on work, we investigated the impact of shock and droplet interaction in a detonation study involving both gaseous as well as gas/liquid (droplet clouds) fuel/air stoichiometric mixtures in a simple 3D shock-tube configuration. The droplet clouds are made up of different initial droplet sizes of either 6, 10, or 30 µm. We also investigated the individual droplet behavior followed by the passage of a detonation front. The results represent conditions that lead to both overdriven and C-J (ChapmanJouguet) detonations. Under both test conditions, most of the droplet vaporization is completed within a short distance (duration) behind the detonation front & well within the region of complete combustion observed in a corresponding equivalent gas-phase fuel/air mixture. The impact of the shock-induced droplet breakup is found to be significant in the calculations involving the 30-µm droplets. Subsequent to the breakup, the drop sizes vary from 1 to 10 µm. Another factor that contributed to the observed rapid vaporization is the result of vaporization taking place under supercritical conditions. The overall detonation properties of various droplet clouds (made up of different initial sizes) are similar to those observed in a corresponding gaseous fuel/air mixture. In the calculation involving a gaseous fuel, the calculated C-J detonation velocity is 1822 m/s involving Jet-A/air and φ = 1. In the overdriven detonation, it is 2044 m/s. In the calculations involving droplet clouds, the corresponding detonation velocities are lower. The impact of increased droplet size is primarily seen in a higher reduction in the detonation velocity.

shock/droplet interaction↗
NASA NTRS2013

The First Pulse of the Extremely Bright GRB 130427A: A Test Lab for Synchrotron Shocks

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 s is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

Pulse↗
NASA NTRS1995

Ultra low frequency waves at the Earth's bow shock

The Earth's bow shock is a bountiful generator of waves. Some of these waves have group velocities that exceed the solar wind velocity directed into the shock and can propagate upstream against the flow. Upstream whistlers observed close to one Hertz in the spacecraft frame have been seen many Earth radii upstream. A second whistler mode wave, called the precursor, propagates upstream along the shock normal but is phase standing in the solar wind flow. The damping of both whistler mode waves is consistent with Landau damping. At low Mach numbers the precursor is connected to the non-coplanarity component in the shock ramp. At higher Mach numbers the upstream waves cannot propagate upstream and ion reflection becomes more important in providing free energy for wave particle interactions. The non-coplanarity component is still present but it now initiates a downstream wave train. Generally the waves just downstream from the bow shock are left hand circularly polarized ion cyclotron waves propagating along the magnetic field at the Alfven velocity. When the upstream Mach number is high and the helium content of the plasma is high, mirror mode waves are observed.

Russell, C. T.↗
NASA NTRS1995

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗
NASA NTRS1990

Two-state ion heating at quasi-parallel shocks

This paper documents the alternating occurrence of two different ion states of the ion distributions downstream from several quasi-parallel shocks: a cooler (and denser) core/shoulder type and a hotter (less dense) and more Maxwellian type. Three separate lines of evidence are presented to show that the two states are not related in an evolutionary sense, but that both are produced alternately at the shock: (1) the asymptotic downstream plasma parameters (density, ion temperautre, and flow speed) are intermediate between those characterizing the two different states closer to the shock, suggesting that the asymptotic state is produced by a mixing of the two initial states; (2) examples of apparently interpenetrating distributions are found during transitions from one state to the other; and (3) examples of both types of distributions are found at actual crossings of the shock ramp.

Thomsen, M. F.↗
NASA NTRS1976

Flux-corrected transport techniques for transient calculations of strongly shocked flows

New flux-corrected transport algorithms are described for solving generalized continuity equations. These techniques were developed by requiring that the finite difference formulas used ensure positivity for an initially positive convected quantity. Thus FCT is particularly valuable for fluid-like problems with strong gradients or shocks. Repeated application of the same subroutine to mass, momentum, and energy conservation equations gives a simple solution of the coupled time-dependent equations of ideal compressible fluid dynamics without introducing an artificial viscosity. FCT algorithms span Eulerian, sliding-rezone, and Lagrangian finite difference grids in several coordinate systems. The latest FCT techniques are fully vectorized for parallel/pipeline processing.

Boris, J. P.↗
NASA NTRS1972

Finite-amplitude pressure waves in the radial mode of a cylinder

A numerical study of finite-strength, isentropic pressure waves transverse to the axis of a circular cylinder was made for the radial resonant mode. The waves occur in a gas otherwise at rest, filling the cylinder. A method of characteristics was used for the numerical solution. For small but finite amplitudes, calculations indicate the existence of waves of permanent potential form. For larger amplitudes, a shock is indicated to occur. The critical value of the initial amplitude parameter in the power series is found to be 0.06 to 0.08, under various types of initial conditions.

Kubo, I.↗
NASA NTRS1980

A supersonic, three-dimensional code for flow over blunt bodies: User's manual

A computer code is described which may be used to calculate the steady, supersonic, three-dimensional, inviscid flow over blunt bodies. The theoretical and numerical formulation of the problem is given (shock-capturing, downstream marching), including exposition of the boundary and initial conditions. The overall flow logic of the program, its usage, accuracy, and limitations are discussed.

Chaussee, D. S.↗
NASA NTRS1980

Type 2 solar radio events observed in the interplanetary medium. Part 1: General characteristics

Twelve type 2 solar radio events were observed in the 2 MHz to 30 kHz frequency range by the radio astronomy experiment on the ISEE-3 satellite over the period from September 1978 to December 1979. These data provide the most comprehensive sample of type 2 radio bursts observed at kilometer wavelengths. Dynamic spectra of a number of events are presented. Where possible, the 12 events were associated with an initiating flare, ground based radio data, the passage of a shock at the spacecraft, and the sudden commencement of a geomagnetic storm. The general characteristics of kilometric type 2 bursts are discussed.

Cane, H. V.↗
NASA NTRS1980

Computerized three-dimensional aerodynamic design of a lifting rotor blade

A three-dimensional, inviscid, full-potential lifting rotor code was used to demonstrate that pressure distributions on both advancing and retreating blades could be significantly improved by perturbing local airfoil sections. The perturbations were described by simple geometric shape functions. To illustrate the procedure, an example calculation was made at a forward flight speed of 85 m/sec (165 knots) and an advance ratio of 0.385. It was found that a minimum of three shape functions was required to improve the pressures without producing undesirable secondary effects in high-speed forward flight on a hypothetical modern rotor blade initially having an NLR-1 supercritical airfoil. Reductions in the shock strength on the advancing blade could be achieved, while simultaneously lessening leading-edge pressure gradients on the retreating blade. The major blade section modifications required were blunting of the upper surface leading edge and some reshaping of the blade's upper surface resulting in moderately thicker airfoils.

Tauber, M. E.↗
NASA NTRS1986

Influence of numerical dissipation in computing supersonic vortex-dominated flows

Steady supersonic vortex-dominated flows are solved using the unsteady Euler equations for conical and three-dimensional flows around sharp- and round-edged delta wings. The computational method is a finite-volume scheme which uses a four-stage Runge-Kutta time stepping with explicit second- and fourth-order dissipation terms. The grid is generated by a modified Joukowski transformation. The steady flow solution is obtained through time-stepping with initial conditions corresponding to the freestream conditions, and the bow shock is captured as a part of the solution. The scheme is applied to flat-plate and elliptic-section wings with a leading edge sweep of 70 deg at an angle of attack of 10 deg and a freestream Mach number of 2.0. Three grid sizes of 29 x 39, 65 x 65 and 100 x 100 have been used. The results for sharp-edged wings show that they are consistent with all grid sizes and variation of the artificial viscosity coefficients. The results for round-edged wings show that separated and attached flow solutions can be obtained by varying the artificial viscosity coefficients. They also show that the solutions are independent of the way time stepping is done. Local time-stepping and global minimum time-steeping produce same solutions.

Kandil, O. A.↗
NASA NTRS1988

A solar-wind 'trigger' for the outer heliosphere radio emissions and the distance to the terminal shock

The solar wind data from the plasma science experiment on the Voyager spacecraft are examined to search for the source of the 2 to 3 kHz radio noise detected in the outer heliosphere. It is found that two anomalous high speed streams passed Voyager 2 before the noise was initially observed. It is suggested that the interaction of these streams with the terminal shock is responsible for the more intense emission. Using a time-of-flight argument, the distance to the shock is estimated at about 70 AU-140 AU. The larger value is consistent with an estimate using Voyager data and a standard pressure-balance argument.

Mcnutt, Ralph L., Jr.↗
NASA NTRS1994

Mineralogy and Ar-39 - Ar-40 of an old pristine basalt: Thermal history of the HED parent body

Previous investigations of mineral chemistry and Rb-Sr and Sm-Nd ages indicated that clast,84 from eucrite Yamato 75011 had preserved the pristine nature of its initial crystallization during an early stage of the HED parent body. Microscale mineralogy and Ar-39-Ar-40 ages of this clast, however, revealed local disturbance of microtextures and partially reset ages. This evidence suggests that, in addition to initial crystallization and rapid cooling, the Y75011,84 clast experienced shock deformation, reheating of short duration at higher temperature, and brecciation. These characteristics suggest two or more impact events. Fe-rich olivine filling fractures in pyroxene may have been introduced during the accompanying shock fracturing. The inferred Ar-39-Ar-40 degassing ages for Y75011 matrix and clast, 84 are 3.94 +/- 0.04 Ga and 3.98 +/- 0.03 Ga, respectively. The suggested degassing age for a clast from Y790020, believed to be paired with Y75011, is approximately 4.03 Ga, but could be younger. We consider it likely that all three samples experienced a common degassing event 3.95 +/- 0.05 Ga ago, but we cannot rule out two or more events spaced over a approximately 0.1 Ga interval. Higher temperature extractions of the two clast samples show significantly older apparent ages up to approximately 4.5 Ga and suggest that the time/temperature regime of this event was not sufficient to degas Ar totally. Most likely, the K-Ar ages were reset by thermal metamorphism associated with one or more impact events associated with shock fracturing, formation of Fe-rich olivine veins, and/or meteorite brecciation. The pyroxene annealing that commonly occurs in many eucrites is likely to be a much earlier process than the impact-produced textural changes and reset K-Ar ages observed in these meteorites. The existence of mineralogical and chronological evidence for metamorphism in an otherwise pristine eucrite suggests that the HED parent body experienced an extensive degree of early cratering.

Takeda, Hiroshi↗
NASA NTRS1994

Ion Chemistry in Atmospheric and Astrophysical Plasmas

There are many differences and also remarkable similarities between the ion chemistry and physics of planetary ionospheres and the ion chemistry and physics of astronomical environments beyond the solar system. In the early Universe, an expanded cooling gas of hydrogen and helium was embedded in the cosmic background radiation field and ionized by it. As the Universe cooled by adiabatic expansion, recombination occurred and molecular formation was driven by catalytic reactions involving the relict electrons and protons. Similar chemical processes are effective in the ionized zones of gaseous and planetary nebulae and in stellar winds where the ionization is due to radiation from the central stars, in the envelopes of supernovae where the ionization is initiated by the deposition of gamma-rays, in dissociative shocks where the ionization arises from electron impacts in a hot gas and in quasar broad-line region clouds where the quasar is responsible for the ionization. At high altitudes in the atmospheres of the Jovian planets, the main constituents are hydrogen and helium and the ion chemistry and physics is determined by the same processes, the source of the ionization being solar ultraviolet radiation and cosmic rays. After the collapse of the first distinct astronomical entities to emerge from the uniform flow, heavy elements were created by nuclear burning in the cores of the collapsed objects and distributed throughout the Universe by winds and explosions. The chemistry and physics became more complicated. Over 90 distinct molecular species have been identified in interstellar clouds where they are ionized globally by cosmic ray impacts and locally by radiation and shocks associated with star formation and evolution. Complex molecules have also been found in circumstellar shells of evolved stars. At intermediate and low altitudes in the Jovian atmospheres, the ion chemistry is complicated by the increasing abundance of heavy elements such as carbon, and an extensive array of complex molecules has been predicted. Reactions involving heavy elements dominate the structure of the ionspheres of the terrestrial planets and the satellites Titan and Triton.

Dalgarno, A.↗
NASA NTRS1994

Adiabatic invariants in stellar dynamics, 3: Application to globular cluster evolution

The previous two companion papers demonstrate that slowly varying perturbations may not result in adiabatic cutoffs and provide a formalism for computing the long-term effects of time-dependent perturbations on stellar systems. Here, the theory is implemented in a Fokker-Planck code and a suite of runs illustrating the effects of shock heating on globular cluster evolution are described. Shock heating alone results in considerable mass loss for clusters with R(sub g) less than or approximately 8 kpc: a concentration c = 1.5 cluster with R(sub g) kpc loses up to 95% of its initial mass in 15 Gyr. Only those with concentration c greater than or approximately 1.3 survive disk shocks inside of this radius. Other effects, such as mass loss by stellar evolution, will decrease this survival bound. Loss of the initial halo together with mass segregation leads to mass spectral indices, x, which may be considerably larger than their initial values.

Weinberg, Martin D.↗
NASA NTRS1981

Propagation of magnetically guided acoustic shocks in the solar chromosphere

The propagation of a train of acoustic shocks guided by diverging magnetic fields through a static model of the solar chromospheric network and transition region is investigated. For initial flux densities of 1,000,000 ergs/sq cm/sec in the lower chromosphere, the local efficiency of acoustic transmission into the corona can be much higher than that calculated for a plane parallel atmosphere. Acoustic energy will tend to be deposited at higher chromospheric levels in diverging magnetic fields, and magnetic guiding may influence the temperature profile of the network and plages. The total flux that can be transmitted into the corona along such diverging fields is severely limited since the magnetic elements occupy a small fractional area of the photosphere and the transmission efficiency is a rapidly decreasing function of initial flux density. Diverging magnetic fields and a varying ratio of specific heats are not likely to allow high frequency shocks to dissipate high enough in a static atmosphere to contribute significantly to the coronal energy balance.

Foukal, P.↗
NASA NTRS1991

Oblique MHD cosmic-ray modified shocks: Two-fluid numerical simulations

We present the first results of time dependent, two-fluid, cosmic-ray (CR) modified, MHD shock simulations. The calculations were carried out with a new numerical code for 1-D ideal MHD. By coupling this code with the CR energy transport equation we can simulate the time-dependent evolution of MHD shocks including the acceleration of the CR and their feedback on the shock structures. We report tests of the combined numerical method including comparisons with analytical steady state results published earlier by Webb, as well as internal consistency checks for more general MHD CR shock structures after they appear to have converged to dynamical steady states. We also present results from an initial time dependent simulation which extend the parameter space domain of previous analytical models. These new results support Webb's suggestion that equilibrium oblique shocks are less effective than parallel shocks in the acceleration of CR. However, for realistic models of anisotropic CR diffusion, oblique shocks may achieve dynamical equilibrium on shorter timescale than parallel shocks.

Frank, Adam↗