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 199 records · Page 11

A transient MHD model applicable for the source of solar cosmic ray acceleration

A two-dimensional, time-dependent magnetohydrodynamic model is used to describe the possible mechanisms for the source of solar cosmic ray acceleration following a solar flare. The hypothesis is based on the propagation of fast mode MHD shocks following a sudden release of energy. In this presentation, the effects of initial magnetic topology and strength on the formation of MHD shocks have been studied. The plasma beta (thermal pressure/magnetic pressure) is considered as a measure of the initial, relative strength of the field. During dynamic mass motion, the Alfven Mach number is the more appropriate measure of the magnetic field's ability to control the outward motion. It is suggested that this model (computed self-consistently) provides the shock waves and the disturbed mass motion behind it as likely sources for solar cosmic ray acceleration.

Dryer, M.↗

Original size of the Vredefort structure, South Africa

The Vredefort structure is located approximately 120 km southwest of Johannesburg, South Africa, and is deeply eroded. Controversies remain on the origin of this structure with the most popular hypotheses being: (1) by impact cratering about 2.0 Ga; (2) as a cryptoexplosion structure about 2.0 Ga; and (3) by purely tectonic processes starting at about 3.0 Ga and ending with the Vredefort event at 2.0 Ga. In view of recent work in which the granophyre dikes are interpreted as the erosional remants of a more extensive impact melt sheet, injected downward into the underlying country rocks, the impact origin hypothesis for Vredefort is adopted. In order to estimate the original dimensions of the Vredefort impact structure, it is assumed that the structure was initially circular, that its predeformation center corresponds to the center of the granitic core, and that the pre-Vredefort geology of the area prior to approximately 2.0 Ga ago is as suggested by Fletcher and Reimold. The spatial relationship between shock metamorphic effects, the shock pressures they record, and the morphological features of the crater were established for a number of large terrestrial craters. The principles of crater formation at large complex impact structures comparable in size to Vredefort were also established, although many details remain unresolved. An important conclusion is that the transient crater, which is formed directly by excavation and displacement by the shock-induced cratering flow-field (i.e., the particle velocity flow field existing in the region of the transient crater but behind the initial outgoing shock front), is highly modified during the late stage processes. The original transient crater diameter lies well within the final rim of the crater, which is established by structural movements during late-stage cavity modification.

Therriault, A. M.↗

Initial results of a porous plug nozzle for supersonic jet noise suppression

As part of a continuing study of possible methods of jet noise reduction, some tests have been made on a porous plug type noise suppressor. Very little information exists on the aeroacoustic performance of jet nozzles having porous elements designed to eliminate the shock wave in the exhaust stream and by so doing to eliminate the shock associated noise and screech. Some initial results on the aeroacoustic performance of a model porous plug type jet noise suppressor are presented. Included are shadowgraph pictures of the flow exhausting from the porous plug nozzle with the comparable acoustic far-field spectra and cross correlations which illustrate the benefits of the test device.

Maestrello, L.↗

Collapse and fragmentation of molecular cloud cores. 2: Collapse induced by stellar shock waves

The standard scenario for low-mass star formation involves 'inside-out' collapse of a dense molecular cloud core following loss of magnetic field support through ambipolar diffusion. However, isotopic anomalies in presolar grains and meteoritical inclusions imply that the collapse of the presolar cloud may have been triggered by a stellar shock wave. This paper explores 'outside-in' collapse, that is, protostellar collapse initiated directly by the compression of quiescent dense cloud cores impacted by relatively slow stellar shock waves. A second-order accurate, gravitational hydrodynamics code has been used to study both the spherically symmetrical and three-dimensional evolution of initially centrally condensed, isothermal, self-gravitating, solar-mass cloud cores that are struck by stellar shock waves with velocities up to 25 km/s and postshock temperatures of 10 to 10,000 K. The models show that such mild shock waves do not completely shred and destroy the cloud, and that the dynamical ram pressure can compress the cloud to the verge of self-gravitational collapse. However, compression caused by a high postshock temperature is a considerably more effective means of inducing collapse. Shock-induced collapse produces high initial mass accretion rates (greater than 10(exp -4) solar mass/yr in a solar-mass cloud) that decline rapidly to much lower values, depending on the presence (approximately 10(exp -6) solar mass/yr) or absence (approximately 10(exp -8) to 10(exp -7) solar mass/yr) of an infinite reservoir of mass. Stellar mass accretion rates approximately 10(exp -7) solar mass/yr have been previously inferred from the luminosities of T Tauri stars; balanced mass accretion (stellar rate = envelope rate) at approximately 10(exp -7) solar mass/yr could then be possible if accretion occurs from a finite mass reservoir. Fluid tracers are used to determine what fraction of the stellar shock material is incorporated into the resulting protostellar object and disk; roughly half the impinging material is injected into the collapsing cloud core when there is a high postshock temperature. The models are consistent with a scenario where an AGB star wind triggered the collapse of the presolar cloud while injecting about 0.01 solar mass of matter derived from the AGB star envelope, as has been separately inferred on the basis of nucleosynthesis calculations.

Boss, Alan P.↗

A Comparison of EAST Shock-Tube Radiation Measurements with a New Air Radiation Model

This paper presents a comparison between the recent EAST shock tube radiation measurements (Grinstead et al., AIAA 2008-1244) and the HARA radiation model. The equilibrium and nonequilibrium radiation measurements are studied for conditions relevant to lunar-return shock-layers; specifically shock velocities ranging from 9 to 11 kilometers per second at initial pressures of 0.1 and 0.3 Torr. The simulated shock-tube flow is assumed one-dimensional and is calculated using the LAURA code, while a detailed nonequilibrium radiation prediction is obtained in an uncoupled manner from the HARA code. The measured and predicted intensities are separated into several spectral ranges to isolate significant spectral features, mainly strong atomic line multiplets. The equations and physical data required for the prediction of these strong atomic lines are reviewed and their uncertainties identified. The 700-1020 nm wavelength range, which accounts for roughly 30% of the radiative flux to a peak-heating lunar return shock-layer, is studied in detail and the measurements and predictions are shown to agree within 15% in equilibrium. The plus or minus 1.5% uncertainty on the measured shock velocity is shown to cause up to a plus or minus 30% difference in the predicted radiation. This band of predictions contains the measured values in almost all cases. For the highly nonequilibrium 0.1 Torr cases, the nonequilibrium radiation peaks are under-predicted by about half. This under-prediction is considered acceptable when compared to the order-of-magnitude over-prediction obtained using a Boltzmann population of electronic states. The reasonable comparison in the nonequilibrium regions provides validation for both the non-Boltzmann modeling in HARA and the thermochemical nonequilibrium modeling in LAURA. The N2 (+)(1-) and N2(2+) molecular band systems are studied in the 290 480 nm wavelength range for both equilibrium and nonequilibrium regimes. The non-Boltzmann rate models for these systems, which have significant uncertainties, are tuned to improve the comparison with measurements.

Johnston, Christopher O.↗

Determining particle density using known material Hugeniot curves

A method is detailed to determine the density of particles wherein the closing velocity is known between the impacting particles and a plate of known material. Either the shock wave velocity or the material velocity produced in the plate upon impact by an unknown material particle is determined and compared with the corresponding shock wave or material velocity that would by produced by different known material particles having the same closing velocity upon impact with the plate. The unknown material particle density is derived by obtaining a coincidence of the shock wave velocity or material velocity conditions initially produced upon impact between the known material plate and one of the different material particles and from the fact that shock wave velocity and material velocity are ordered on the impacting particle material density alone.

Dibattista, J. D.↗

Time-resolved spectroscopic measurements behind incident and reflected shock waves in air and xenon

Time-resolved spectra have been obtained behind incident and reflected shock waves in air and xenon at initial pressures of 0.1 and 1.0 torr using a rotating drum spectrograph and the OSU (The Ohio State University) arc-driven shock tube. These spectra were used to determine the qualitative nature of the flow as well as for making estimates of the available test time. The (n+1,n) and (n,n) band spectra of N2(+) (1st negative) were observed in the test gas behind incident shock waves in air at p1=1.0 torr and Us=9-10 km/sec. Behind reflected shock waves in air, the continuum of spectra appeared to cover almost the entire wavelength of 2,500-7,000 A for the shock-heated test gas. For xenon, the spectra for the incident shock wave cases for p1=0.1 torr show an interesting structure in which two intensely bright regions are witnessed in the time direction. The spectra obtained behind reflected shock waves in xenon were also dominated by continuum radiation but included strong absorption spectra due to FeI and FeII from the moment the reflected shock passed and on.

Yoshinaga, T.↗

A normal shock-wave turbulent boundary-layer interaction at transonic speeds

Experimental results, including surveys of the mean and fluctuating flow, and measurements of surface pressure, skin friction, and separation length, are compared with solutions to the Navier-Stokes equations utilizing various algebraic eddy viscosity models to describe the Reynolds shear stresses. The experimental data, obtained at a free-stream Mach number of 1.5 and Reynolds numbers between 10 million and 80 million, show that a separated zone forms near the foot of the shock and that its length is proportional to the initial boundary-layer thickness; that a supersonic region forms downstream of the shock; and that the shear stress increases significantly through the interaction and subsequently decays downstream. The computations adequately represent the qualitative features of the flow field throughout the interaction but quantitatively underpredict the extent of separation and the downstream level of skin friction.

Mateer, G. G.↗

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗

Heat shield erosion in a dusty Martian atmosphere

The study calculates the effects of dust particle impacts on the erosion of the forebody heat shield for a 26-m diameter aerobraking vehicle entering the Martian atmosphere at 8.6 km/s. An explicit thin-layer Navier-Stokes code is used to compute the dustless flowfield about the vehicle for the actual Martian atmospheric composition. The deceleration and melting of the dust particles within the forebody shock layer was computed for dust spherules having initial diameters from 1 to 19 microns. All the particles began vaporizing shortly after entering the shock layer, but most survived to hit the heat shield surface. The two different heat shield materials considered were Shuttle ceramic tiles with glassy surfaces and a moderate-density ablator. For a vehicle with a ballistic coefficient of 200 kg/sq m, a glassy heat shield surface experienced an average of 7 mm of surface erosion, or 20 times the thickness of the glassy layer on ceramic tiles.

Papadopoulos, Periklis↗

Weak shock propagation through a turbulent atmosphere

Consideration is given to the propagation through turbulence of transient pressure waveforms whose initial onset at any given point is an abrupt shock. The work is motivated by the desire to eventually develop a mathematical model for predicting statistical features, such as peak overpressures and spike widths, of sonic booms generated by supersonic aircraft. It is argued that the transient waveform received at points where x greater than 0 will begin with a pressure jump and a formulation is developed for predicting the amount of this jump and the time derivatives of the pressure waveform immediately following the jump.

Pierce, Allan D.↗

Solar Sentinels: Report of the Science and Technology Definition Team

The goal of NASA s Living With a Star (LWS) program is to develop the scientific understanding necessary to effectively address those aspects of the connected Sun Earth system that directly affect life and society. Along with the other elements of LWS, Solar Sentinels aims to discover, understand, and model the heliospheric initiation, propagation, and solar connection of those energetic phenomena that adversely affect space exploration and life and society here on Earth. The Solar Sentinels mission will address the following questions: (1) How, where, and under what circumstances are solar energetic particles (SEPs) accelerated to high energies and how do they propagate through the heliosphere? And (2) How are solar wind structures associated with these SEPs, like CMEs, shocks, and high-speed streams, initiated, propagate, evolve, and interact in the inner heliosphere? The Sentinels STDT recommends implementing this mission in two portions, one optimized for inner heliospheric in-situ measurements and the other for solar remote observations. Sentinels will greatly enhance the overall LWS science return.

Source record↗

High-Order Simulation of Non-Linear Oscillations and Shocks in the Solar Atmosphere

The solar atmosphere presents a rich source of highly non-linear magneto-hydrodynamic phenomena: strong gradients and forcing terms result in both large shocks and oscillations. The additional requirements of energy balance and initialization in hydrostatic equilibrium compound the challenge of this problem. A wealth of observational data allows us to check the results of our simulations. The problem of simulating the solar atmosphere provides, in addition to an interesting system in its own right, a challenging testbed for high-order shock-capturing methods. We discuss the challenge of simulating solar atmospheric phenomena, concentrating on various high-order central methods ranging from second to fourth order. Our method is based on the central-upwind scheme of Kurganov, Noelle and Petrova, which we extend to high order via various interpolants. We investigate various initial data for our simulations, corresponding to observed conditions in different regions of the solar surface: the normal quiet sun and sunspots. When non-oscillatory using second- and third-order methods, we are able to reproduce non-trivial observational results. In particular we find a correlation between initial data and both the shock speeds and particle oscillation spectra that match observations in the corresponding regions. When using fourth-order WENO interpolants, we find that while the individual shock profiles at any given time appear non-oscillatory, spurious oscillations appear in the fields after long time integrations. The issue of initialization in hydrostatic equilibrium raises difficult issues. Careful treatment of the gravitational source term can reduce violations of hydrostatic equilibrium, but difficulties remain, primarily due to discontinuities in the piecewise-polynomial reconstructions. In the case of the Euler equations in gravity, the use of high-order methods reduces the violation of hydrostatic equilibrium to a sufficiently low level for the above described results to be obtained. Including various energy flux terms in the equations, however, leads to severe instability when hydrostatic balance is violated. We discuss various strategies to address this problem, including adaptive grid methods.

Bryson, S.↗

Shock drift acceleration in the presence of waves

Attention is given to the initial results of a model designed to study the modification of the scatter-free, shock drift acceleration of energetic test particles by wave activity in the vicinity of a quasi-perpendicular, fast-mode MHD shock. It is emphasized that the concept of magnetic moment conservation is a valid approximation only in the perpendicular and nearly perpendicular regimes, when the angle theta-Bn between the shock normal and the upstream magnetic field vector is in the range from 70 deg to 90 deg. The present investigation is concerned with one step in a program which is being developed to combine the shock drift and diffusive processes at a shock of arbitrary theta-Bn.

Decker, R. B.↗

Ion acceleration at collisionless shock interactions

When two collisionless shocks collide, a population of high energy ions is produced. Using hybrid numerical simulations, it is shown that when the two shocks are quasi-perpendicular, ions with energies up to 15 E sub 0 (where E sub 0 is the initial kinetic energy of upstream ions in the shock frame) are produced by direct electric field acceleration. In the quasi-parallel regime, energies as high as 80 E sub 0 can be obtained. These ions are accelerated by both scattering off the approaching shocks and subsequently in the intense turbulence left behind by the shock collision.

Cargill, Peter J.↗

Distances to the termination shock and heliopause from a simulation analysis of the 1992-93 heliospheric radio emission event

A new heliospheric radio emission event observed by Voyagers 1 and 2 in mid-1992 is believed to have been produced by the interaction of an interplanetary shock with the heliopause. The shock is thought to have oriented near the Sun during a period of intense solar activity in late-May and early-June, 1991. The observed travel time of the shock to the heliopause is 408 days; the initial speed is estimated to be between 600 and 800 km/s. We use a numerical gasdynamic simulation of an interplanetary shock, propagating through an equilibrium solution of the solar wind/interstellar medium interaction, to compute the distances to the termination shock and the heliopause that are consistent with these observations. For a shock speed of 600 km/s, the termination shock is located at 92 AU, and the heliopause is located at 128 AU. These distances increase to 112 AU and 156 AU when the shock speed is increased to 800 km/s.

Steinolfson, R. S.↗

Focused interplanetary transport of approximately 1 MeV solar energetic protons through self-generated Alfven waves

We present a model of the focused transport of approximately 1 MeV solar energetic protons through interplanetary Alfven waves that the protons themselves amplify or damp. It is based on the quasi-linear theory but with a phenomenological pitch angle diffusion coefficient in the 'resonance gap.' For initial Alfven wave distributions that give mean free paths greater than approximately 0.5 AU for approximately 1 MeV protons in the inner heliosphere, the model predicts greater than roughly an order of magnitude amplification (damping) in the outward (inward) propagating resonant Alfven waves at less than or approximately equal to o.3 AU heliocentric distance. As the strength of proton source is increased, the peak differential proton intensity at approximately 1 MeV at 1 AU increases to a maximum of approximately 250 particles (/(sq cm)(s)(sr)(MeV)) and then decreases slowly. It may be attenuated by a factor of 5 or more relative to the case without wave evolution, provided that the proton source is sufficiently intense that the resulting peak differential intensity of approximately 1 MeV protons at 1 AU exceeds approximately 200 particles (/(sq cm)(s)(sr)(MeV)). Therefore, in large solar proton events, (1) one may have to take into account self-amplified waves in studying solar particle propagation, (2) the number of accelerated protons escaping from a flare or interplanetary shock may have been underestimated in past studies by a significant factor, and (3) accelerated protons escaping from a traveling interplanetary shock at r less than or approximately equal to 0.3 AU should amplify the ambient hydromagnetic waves siginificantly to make the shock an efficient accelerator, even if initially the mean free path is greater than or approximately equal to 1 AU.

Ng, C. K.↗