Search NASA⌕ Search

SEARCH · Search NASA

Results for “Shock instability”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Numerical simulations of time-dependent cosmic ray mediated shocks

Recent research in unsteady shocks that are influenced by nonlinear aspects of diffusive particle acceleration and that involve multidimensional flows is outlined, and some of the difficulties in making further progress are discussed. It is suggested that in many situations it is probably not appropriate to assume steady states for astrophysical CR shocks. This results both from the need to be concerned about the finite rate at which CR pressure builds and modifies a shock and from the existence of instabilities. Even so, it seems likely that CR-modified shocks can develop to a state with substantial CR pressure. The overall acceleration efficiency and the CR particle momentum distribution produced in time-dependent flows typically depend on several details such as the ratio of various time scales and the initial proportion of relativistic to nonrelativistic particles in the CR populations.

Jones, T. W.↗
NASA NTRS2011

Study of Unsteady, Sphere-Driven, Shock-Induced Combustion for Application to Hypervelocity Airbreathing Propulsion

A premixed, shock-induced combustion engine has been proposed in the past as a viable option for operating in the Mach 10 to 15 range in a single stage to orbit vehicle. In this approach, a shock is used to initiate combustion in a premixed fuel/air mixture. Apparent advantages over a conventional scramjet engine include a shorter combustor that, in turn, results in reduced weight and heating loads. There are a number of technical challenges that must be understood and resolved for a practical system: premixing of fuel and air upstream of the combustor without premature combustion, understanding and control of instabilities of the shock-induced combustion front, ability to produce sufficient thrust, and the ability to operate over a range of Mach numbers. This study evaluated the stability of the shock-induced combustion front in a model problem of a sphere traveling in a fuel/air mixture at high Mach numbers. A new, rapid analysis method was developed and applied to study such flows. In this method the axisymmetric, body-centric Navier-Stokes equations were expanded about the stagnation streamline of a sphere using the local similarity hypothesis in order to reduce the axisymmetric equations to a quasi-1D set of equations. These reduced sets of equations were solved in the stagnation region for a number of flow conditions in a premixed, hydrogen/air mixture. Predictions from the quasi-1D analysis showed very similar stable or unstable behavior of the shock-induced combustion front as compared to experimental studies and higher-fidelity computational results. This rapid analysis tool could be used in parametric studies to investigate effects of fuel rich/lean mixtures, non-uniformity in mixing, contaminants in the mixture, and different chemistry models.

Axdahl, Erik↗
NASA NTRS2005

Particle acceleration magnetic field generation, and emission in Relativistic pair jets

Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.

Nishikawa, K.-I.↗
NASA NTRS1991

A scale-invariant cellular-automata model for distributed seismicity

In the standard cellular-automata model for a fault an element of stress is randomly added to a grid of boxes until a box has four elements, these are then redistributed to the adjacent boxes on the grid. The redistribution can result in one or more of these boxes having four or more elements in which case further redistributions are required. On the average added elements are lost from the edges of the grid. The model is modified so that the boxes have a scale-invariant distribution of sizes. The objective is to model a scale-invariant distribution of fault sizes. When a redistribution from a box occurs it is equivalent to a characteristic earthquake on the fault. A redistribution from a small box (a foreshock) can trigger an instability in a large box (the main shock). A redistribution from a large box always triggers many instabilities in the smaller boxes (aftershocks). The frequency-size statistics for both main shocks and aftershocks satisfy the Gutenberg-Richter relation with b = 0.835 for main shocks and b = 0.635 for aftershocks. Model foreshocks occur 28 percent of the time.

Barriere, Benoit↗
NASA NTRS1985

Collisionless shocks in the heliosphere: Reviews of current research

The present conference on heliospheric collisionless shocks considers such macrostructure-, microstructure-, and particle acceleration-related topics as interplanetary shock phenomena near and within 1 AU, as well as beyond, planetary bow shocks, shock formation and evolution in the solar atmosphere, MHD and gasdynamic theories for planetary bow waves, and subcritical collisionless shock waves. Also discussed are ion reflection, gyration, and dissipation at supercritical shocks, the numerical simulation of quasi-perpendicular collisionless shocks, electron distributions near collisionless shocks, the microtheory of collisionless shock current layers, plasma waves and instabilities, the electron foreshock, upstream suprathermal ions, and both diffusive and shock drift acceleration.

Tsurutani, Bruce T.↗
NASA NTRS1973

TRW plasma wave experiment for the IMP-H mission

The IMP-H plasma wave experiment is designed to extend knowledge of wave-particle interactions in the disturbed cislunar region, the distant geomagnetic tail, the upstream solar wind, and the flanks of the magnetosheath-shock interface. It is expected to identify plasma instabilities, study particle acceleration and heating at collisionless shocks and other discontinuities, analyze turbulent conductivity and field line merging, and provide new information on dissipation processes for suprathermal particles. Instrumentation for the plasma wave experiment is designed to measure local electric and magnetic field oscillations over the frequency range 10 Hz to 100 kHz. A 24 inch electric dipole, a 7 inch diameter air core search coil, and the associated preamplifiers are mounted on a spacecraft counterweight boom. The frequency range of 10 Hz to 100 kHz for both E and B is processed using an eight-channel spectrum analyzer located in the instrument main-body package (a standard IMP trapezoidal module, 3 inches high). Electric fields as small as 10-100 microvolts/meter and magnetic signals as small as 1-3 milligamma will be detected.

Virobik, P. F.↗
NASA NTRS1983

Plasma and energetic particle structure of a collisionless quasi-parallel shock

The quasi-parallel interplanetary shock of November 11-12, 1978 from both the collisionless shock and energetic particle points of view were studied using measurements of the interplanetary magnetic and electric fields, solar wind electrons, plasma and MHD waves, and intermediate and high energy ions obtained on ISEE-1, -2, and -3. The interplanetary environment through which the shock was propagating when it encountered the three spacecraft was characterized; the observations of this shock are documented and current theories of quasi-parallel shock structure and particle acceleration are tested. These observations tend to confirm present self consistent theories of first order Fermi acceleration by shocks and of collisionless shock dissipation involving firehouse instability.

Kennel, C. F.↗
NASA NTRS1985

Plasma waves and instabilities

An account is given of the waves and instabilities occurring at collisionless shocks, with attention to the mechanisms responsible for the generation of these waves. The transition region of the shock usually involves an abrupt broadband burst of electrostatic noise that extends from below the lower hybrid resonance to near the electron plasma frequency, and by a broadband burst of whistler mode EM noise below the electron cyclotron frequency. Electrostatic lower hybrid waves are also noted. Upstream of the shock, electron plasma oscillations, ion acoustic waves, and intense ULF MHD waves are often observed. The region downstream of the shock is usually very chaotic; electrostatic waves often extend long distances into the downstream region, together with whistler mode emissions.

Gurnett, Donald A.↗
NASA NTRS1990

Lyman alpha radiation in external galaxies

The Ly alpha line of atomic hydrogen is often a luminous component of the radiation emitted by distant galaxies. Except for those galaxies which have a substantial central source of non-stellar ionizing radiation, most of the Ly alpha radiation emitted by galaxies is generated within regions of the interstellar medium which are photoionized by starlight. Conversely, much of the energy radiated by photoionized regions is carried by the Ly alpha line. Only hot, massive stars are capable of ionizing hydrogen in the interstellar medium which surrounds them, and because such stars are necessarily short-lived, Ly alpha emission traces regions of active star formation. Researchers argue that the strength of the Ly alpha emission observed from external galaxies may be used to estimate quantitatively the dust content of the emitting region, while the Ly alpha line profile is sensitive to the presence of shock waves. Interstellar dust particles and shock waves are intimately associated with the process of star formation in two senses. First, both dust particles and shock waves owe their existence to stellar activity; second, they may both serve as agents which facilitate the formation of stars, shocks by triggering gravitational instabilities in the interstellar gas that they compress, and dust by shielding star-forming molecular clouds from the ionizing and dissociative effects of external UV radiation. By using Ly alpha observations as a probe of the dust content in diffuse gas at high redshift, we might hope to learn about the earliest epochs of star formation.

Neufeld, David A.↗
NASA NTRS1993

Gasdynamics of relativistically expanding gamma-ray burst sources - Kinematics, energetics, magnetic fields, and efficiency

We calculate both analytically and numerically the evolution of highly relativistic fireballs through the stages of free expansion and coasting, and determine the dependence of the thermodynamic and radiation variables in the comoving and laboratory flames. The dynamics and the comoving geometry change at the (lab) expansion factors r/r(0) greater than eta and r/r(0) greater than eta-squared, respectively, where eta = E(0)/M(0)c-squared is the initial Lorentz factor. In the lab, the gas appears concentrated in a thin shell of width r(0) until r/r(0) of less than about eta-squared, and increases linearly after that. Magnetic fields may have been important in the original impulsive event. We discuss their effect on the fireball dynamics and also consider their effects on the radiation emitted when the fireball runs into an external medium and is decelerated. The inverse synchro-Compton mechanism can then yield high radiative efficiency in the reverse shock (and through turbulent instabilities and mixing also in the forward blast wave), producing a burst of nonthermal radiation mainly in the MeV to GeV range. The energy and duration depend on eta, the magnetic field strength, and the external density, and can match the range of properties observed in cosmic gamma-ray bursts.

Meszaros, P.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Kotov, Dmitry V.↗
NASA NTRS2013

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗
NASA NTRS2014

Comparitive Study of High-Order Positivity-Preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density ρ and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order Weighted Essentially Non-Oscillatory (WENO) schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative ρ may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more di cult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium ow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case as well as for the case consisting of two species and one reaction.. This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

Schemes↗
NASA NTRS1984

Effect of electron thermal anisotropy on the kinetic cross-field streaming instability

The investigation of the kinetic cross-field streaming instability, motivated by the research of collisionless shock waves and previously studied by Wu et al. (1983), is discussed more fully. Since in the ramp region of a quasi-perpendicular shock electrons can be preferentially heated in the direction transverse to the ambient magnetic field, it is both desirable and necessary to include the effect of the thermal anisotropy on the instability associated with a shock. It is found that Te-perpendicular greater than Te-parallel can significantly enhance the peak growth rate of the cross-field streaming instability when the electron beta is sufficiently high. Furthermore, the present analysis also improves the analytical and numerical solutions previously obtained.

Tsai, S. T.↗
NASA NTRS1981

Flow and acoustic properties of low Reynolds number supersonic underexpanded jets

Flow and acoustic measurements are made of cold model jets exhausting from a choked nozzle at pressure conditions corresponding to those of Mach 1.4 and 2.1 jets to investigate noise production properties of underexpanded supersonic jets. Mean flow measurements are made using pitot and static pressure probes, with flow fluctuation measurements made with a hot-wire probe and acoustic measurements made with a transversing microphone. Two convergent nozzles with exit diameters of 7.0 and 7.9 mm are used with an exciter consisting of a 0.8 mm tungsten electrode positioned 2 mm from the exit. Shock structure is observed as having a significant effect on the development of the flow field, while large-scale instabilities have higher growth rates in the shock containing underexpanded jets. The role of the asymmetric n = + or - 1 sinusoidal instability is clarified, and results suggest that the broadband shock associated noise of conventional high Reynolds number jets is not related to large-scale jet instability.

Hu, T. F.↗
NASA NTRS1982

Physical mechanisms for turbulent dissipation in collisionless shock waves

Dissipation processes in collisionless shock waves are discussed in terms of simple physical models. Two cases are considered: one is concerned with quasi-parallel shocks and the other deals with quasi-perpendicular shocks. In the former case it is explained why the magnetic field profile of a parallel shock usually appears turbulent. The dominant instability responsible for collisionless dissipation processes is discussed. The situation is quite different in a quasi-perpendicular shock in which cross-field currents are believed to be responsible for the collisionless dissipation. Important instabilities are reviewed and discussed. A brief comparison of the theory with recent observations will also be given.

Wu, C. S.↗
NASA NTRS1991

Alfven shock trains

The Cohen-Kulsrud-Burgers equation (CKB) is used to consider the nonlinear evolution of resistive, quasi-parallel Alfven waves subject to a long-wavelength, plane-polarized, monochromatic instability. The instability saturates by nonlinear steepening, which proceeds until the periodic waveform develops an interior scale length comparable to the dissipation length; a fast or an intermediate shock then forms. The result is a periodic train of Alfven shocks of one or the other type. For propagation strictly parallel to the magnetic field, there will be two shocks per instability wavelength. Numerical integration of the time-dependent CKB equation shows that an initial, small-amplitude growing wave asymptotes to a stable, periodic stationary wave whose analytic solution specifies how the type of shock embedded in the shock train, and the amplitude and speed of the shock train, depend on the strength and phase of the instability. Waveforms observed upstream of the earth's bowshock and cometary shocks resemble those calculated here.

Malkov, M. A.↗