Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusive shock acceleration”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Theoretical calculations of ion acceleration in the vicinity of comet Giacobini-Zinner

Ionization of cometary neutral molecules produces ions which are picked up by the solar wind. The cometary ion pickup process for comet Giacobini-Zinner is studied in two ways: (1) with a test particle method in which trajectories are numerically calculated for several thousand ions whose initial locations were chosen randomly with probability proportional to the neutral density and (2) with a quasi-linear diffusion model. The cometary ion distribution function was calculated with the test particle model at several locations upstream of the bow shock and for several types of magnetic fluctuations (or waves). These waves were allowed to propagate in both directions along the magnetic field at the Alfven speed. Both pitch angle scattering and energy diffusion are evident in the derived ion distributions. The monochromatic waves result in less ion acceleration than turbulent fluctuations with about the same amplitude. The calculated ion distribution functions are in reasonable agreement with the distributions measured in the vicinity of comet Giacobini-Zinner in 1985 by particle detectors on the ICE spacecraft when the ratio of power in sunward propagating Alfven waves to the power in antisunward propagating waves is assumed to lie between about 20% and 50%. However, the quasi-linear diffusion model results agree best with the measured distribution functions when the power ratio is only about 3 percent.

Mckenzie, M. L.↗
NASA NTRS1989

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.↗
NASA NTRS1989

The radiative signature of the shock acceleration process in active galactic nuclei

A time-dependent model in which plasma enters a shock region and undergoes first-order Fermi acceleration is considered. The time evolution of the particle distribution and the associated synchrotron emission has been calculated for an energy-dependent diffusion coefficient. When the emitting region is not spatially resolved, the spectrum is a power law with an exponential turnover in frequency. With increasing time, the power law and the turnover extend to higher frequencies. To test the model, observations at three or more frequencies should be made of time-resolved flux increases from the optically thin nonthermal emission region in variable active galactic nuclei. Existing data that can be used to test the model are scarce but offer qualitative support.

Bregman, Joel N.↗
NASA NTRS1989

Stochastic acceleration of cometary pickup ions - The classic leaky box model

The acceleration of cometary pickup ions by magnetohydrodynamic waves at P/Giacobini-Zinner is examined in a model where acceleration predominantly occurs downstream of the bow shock throughout the cometosheath where intense magnetic turbulence exists. The mean free path for scattering by the magnetic fluctuations in this region is less than the characteristic dimension of the cometosheath so that pickup ions are rendered isotropic and energized by a modest amount in the process. This principal loss mechanism for the ions is spatial diffusion out of the acceleration region moderated by the self-same accelerating waves. This particular feature constrains the model in a way that the predicted ion spectrum is uniquely determined by the power spectrum of the magnetic turbulence. At both P/Giacobini-Zinner and P/Halley, the turbulence is non-Kolmogoroff with a spectral index of 2 resulting in an ion spectral behavior that is approximately an exponential in ion speed, consistent with a recent analysis of the Giacobini-Zinner data (see Richardson et al.).

Barbosa, D. D.↗
NASA NTRS1981

Acceleration theory for 5-40 keV ions at interplanetary shocks

A theory is presented to explain the acceleration of suprathermal ions observed near propagating shocks in the solar wind. The hard power-law spectra of ions of energies 5 to 40 keV are shown to be accounted for by diffusive scattering across a plane shock with a small amount of adiabatic deceleration losses on both sides of the shock. The theory fits the observations of the event of November 25, 1977, which followed a large solar flare, and indicates a mean free path perpendicular to the shock of less than 0.0003 AU behind the shock and 0.01 AU in front of it. The theory also predicts a steepening ion energy spectrum at higher energies.

Forman, M. A.↗
NASA NTRS1994

Stream interfaces and energetic ions closer than expected: Analyses of Pioneers 10 and 11 observations

An empirical study of corotating interaction regions (CIRs) observed between 3.9 AU and 5.9 AU on Pioneers 10 and 11 shows that the main corotation energetic ion population (CEIP), which is associated with the trailing reverse shock, terminates within the CIR at a definite, structural boundary, which we show here is the stream interface. This new result has significant implications for solar wind and energetic particle modeling. In particular it implies either that the reverse shock forms closer to the stream interface than models suggest or that the theories that treat the generation and transport of these energetic ions, such as preshock Fermi acceleration and cross-field diffusion must be combined or extended. We test these scenarios by comparing the CEIP intensity profiles on the two sides of the stream interface. We find that while each automatically accounts for one or two aspects of the results none of them alone can account for all of our empirical results.

Intriligator, Devrie S.↗
NASA NTRS1992

Investigation of starting transients in the thermally choked ram accelerator

An experimental investigation of the starting transients of the thermally choked ram accelerator is presented in this paper. Construction of a highly instrumented tube section and instrumentation inserts provide high resolution experimental pressure, luminosity, and electromagnetic data of the starting transients. Data obtained prior to and following the entrance diaphragm show detailed development of shock systems in both combustible and inert mixtures. With an evacuated launch tube, starting the diffuser is possible at any Mach number above the Kantrowitz Mach number. The detrimental effects and possible solutions of higher launch tube pressures and excessive obturator leakage (blow-by) are discussed. Ignition of a combustible mixture is demonstrated with both perforated and solid obturators. The relative advantages and disadvantages of each are discussed. Data obtained from these starting experiments enhance the understanding of the ram accelerator, as well as assist in the validation of unsteady, chemically reacting CFD codes.

Burnham, E. A.↗
NASA NTRS1983

Active galaxies and the diffuse gamma-ray background

Active galaxies are shown to account for the observed gamma ray background radiation if a steepening of the spectra above about 100 keV is present. An analytical model is discussed in which protons undergo Fermi acceleration at a shock in a spherical accretion flow onto a massive black hole. Relativistic protons with power law spectra, nuclear interactions producing gamma rays from neutal pion decay and electrons from pion-mu meson-electron decay, with a power law spectrum above several hundred MeV, synchrotron and inverse Compton losses steepening the electron spectrum, a photon spectrum close to the pion gamma spectrum and a high-energy gamma ray spectrum steepened by photon-photon pair production interactions with X rays are covered in the model. Comparisons are made with HEAO 2 data on active galaxies, which have estimated luminosities and radii consistent with the compactness necessary for producing the steepening predicted by the model. The active galaxies spectra would be described by a spherical accretion-shock model.

Kazanas, D.↗
NASA NTRS1991

Finite element computation of compressible flows with the SUPG formulation

Finite element computation of compressible Euler equations is presented in the context of the streamline-upwind/Petrov-Galerkin (SUPG) formulation. The SUPG formulation, which is based on adding stabilizing terms to the Galerkin formulation, is further supplemented with a shock capturing operator which addresses the difficulty in maintaining a satisfactory solution near discontinuities in the solution field. The shock capturing operator, which has been derived from work done in entropy variables for a similar operator, is shown to lead to an appropriate level of additional stabilization near shocks, without resulting in excessive numerical diffusion. An implicit treatment of the impermeable wall boundary condition is also presented. This treatment of the no-penetration condition offers increased stability for large Courant numbers, and accelerated convergence of the computations for both implicit and explicit applications. Several examples are presented to demonstrate the ability of this method to solve the equations governing compressible fluid flow.

Le Beau, G. J.↗
NASA NTRS1984

Characteristics of the magnetospheric source of interplanetary energetic particles

The Earth's bow shock is frequently cited as an example of an astrophysical shock where particle acceleration is observed. However, because energetic particles observed upstream of the bow shock may be accelerated within the magnetosphere, it is important to understand the properties of the magnetospheric source. A first order picture of the spatial distribution of magnetospheric particles in the magnetosheath and upstream is obtained by mapping those magnetic field lines which drape over the magnetopause through the bow shock. Subsets of these field lines that connect to potential sites of magnetic merging on the magnetopause are also traced in the event that leakage occurs preferentially where normal components of the field are present across the boundary. The results can be used to determine whether the so-called diffuse particles observed upstream are accelerated locally or within the magnetosphere.

Luhmann, J. G.↗
NASA NTRS1988

Preshock region acceleration of implanted cometary H(+) and O(+)

A self-consistent, three-fluid model of plasma transport and implanted ion acceleration in the unshocked solar wind is presented. The solar wind plasma is depleted by charge exchange with the expanding cometary exosphere, while implanted protons and heavy ions are produced by photoionization and charge transfer and lost by charge exchange. A generalized transport equation describing convection, adiabatic and diffusive velocity change, and the appropriate production terms is used to describe the evolution of the two cometary ion components, while the moments of the Boltzmann equation are used to calculate the solar wind density and pressure. The flow velocity is obtained self-consistently by combining the conservation equations of the three ion species. The results imply that second-order Fermi acceleration can explain the implanted spectra observed in the unshocked solar wind. Comparison of measured and calculated distribution indicates that spatial diffusion of implanted ions probably plays an important role in forming the energetic particle environment in the shock vicinity.

Gombosi, Tamas I.↗
NASA NTRS1988

Observations related to the acceleration, injection, and interplanetary propagation of energetic protons during the solar cosmic ray event on February 16, 1984

This paper presents an analysis of data collected by the worldwide network of neutron monitors and from IMP-8 cosmic-ray telescopes, as well as by particle detectors on the GOES 5 and 6 and ICE satellites, on the solar cosmic ray event that took place on February 16, 1984. Using these data, the intensity-time (IT) profiles, the anisotropy-time profiles, the energy spectra, and the pitch angle distributions of the solar protons near earth were deduced. It was found that the solar protons propagated essentially scatter-free from the sun to the earth. The solar protons had easy access to the IMF lines to earth; the time from the onset to maximum intensity and the shape of the IT profiles at earth as a function of energy could be explained by the diffusion of the flare protons near the acceleration region. The energy spectrum of the solar flare protons injected into the undisturbed IMF at the sun was changing with time in both amplitude and shape. The observations suggest a shock acceleration process.

Debrunner, H.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1977

Nuclear burst plasma injection into the magnetosphere and resulting spacecraft charging

The passage of debris from a high altitude ( 400 km) nuclear burst over the ionospheric plasma is found to be capable of exciting large amplitude whistler waves which can act to structure a collisionless shock. This instability will occur in the loss cone exits of the nuclear debris bubble, and the accelerated ambient ions will freestream along the magnetic field lines into the magnetosphere. Using Starfish-like parameters and accounting for plasma diffusion and thermalization of the propagating plasma mass, it is found that synchronous orbit plasma fluxes of high temperature electrons (near 10 keV) will be significantly greater than those encountered during magnetospheric substorms. These fluxes will last for sufficiently long periods of time so as to charge immersed bodies to high potentials and arc discharges to take place.

Pavel, A. L.↗
NASA NTRS1990

First-order Fermi particle acceleration by relativistic shocks

Monte Carlo calculations of test particle spectra and acceleration times are presented from first-order Fermi particle acceleration for parallel shocks with arbitrary flow velocities and compression ratios r up to seven, shock velocities u1 up to 0.98c, and injection energies ranging from thermal to highly superthermal. Far above the injection energy, the spectra are well-approximated by a power law and the spectra are always harder than for nonrelativistic shocks. Approximate analytic expression are given for the spectral slope as a function of u1 and r. The acceleration time as a function of particle energy is less than for nonrelativistic shocks by a factor that increases with u1 and is about three for u1 = 0.98c. It is confirmed that the spectrum for pitch-angle diffusion is considerably steeper than for large-angle scattering for the same shock parameters.

Ellison, Donald C.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1993

Magnetic reconnection with large separatrix angles

The magnetic reconnection process is studied here using incompressible MHD simulations with different inflow boundary conditions and different magnetic Reynolds numbers R(m). The angle between the magnetic separatrices is in steady state reconnection depends mainly on the normal magnetic field on the inflow boundary. In steady state nonuniform reconnection with large separatrix angles, field-aligned plasma jets appear slightly downstream of the magnetic separatrices. The field-aligned plasma jet are stronger when R(m) is larger. Each field-aligned plasma jet consists of two parts: a slow shock and a fast-mode compressional wave. The slow shock converts the magnetic energy into plasma kinetic energy by acceleration and heating. The fast-mode compressional wave decelerates the plasma to a smaller outflow speed and heats it further. Nearly all the magnetic energy flowing into the diffusion region is converted into other forms. The length and width of the diffusion region depend on the values of the reconnection rate, R(m), and the normal magnetic field on the inflow boundary.

Yan, M.↗
NASA NTRS1987

Pick-up ions at Comet P/Halley's bow shock - Observations with the IIS spectrometer on Giotto

Gaseous material expanding from the nucleus of Comet Halley into space forms the neutral coma around the comet. Ionization in the solar UV radiation removes particles from the coma and injects them into the solar wind plasma. These freshly created ions are accelerated by the interplanetary electric field on cycloidal trajectories with gyrocenters moving with the speed of the magnetic field lines. In the solar wind frame of reference these particles move along the magnetic field lines with a fixed pitch angle. Pitch-angle scattering and energy diffusion reduce quickly the initial energy anisotropy which is associated with the narrow pick-up structures. First observations of heavy cometary pick-up ions (water group ions) at the bow shock are presented. The evolution of the distribution function in the vicinity of the shock and radial density profiles are discussed.

Wilken, B.↗