SEARCH · Search NASA
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.
The motion of suprathermal electrons through the quasi-perpendicular, curved bow shock is considered, using a combined test-particle and hybrid-code approach. Whistler-wave generation and pitch-angle scattering of the derived distributions is studied with an implicit full-particle code. The results are compared to published ISEE observations at the Earth's bow shock.
Stationary transonic alpha-viscosity models of accretion disks surrounding nonrotating black holes have been investigated. The viscosity is modified such that it vanishes in the supersonic region to ensure its effect does not violate the causality condition. In contrast to previous studies, the viscous stress is taken to be explicitly proportional to the angular velocity gradient and is not assumed to depend solely on the local pressure in the disk. The numerical results reveal that the structure of the innermost regions of the disk are more sensitive to the modified form of the viscosity than to the form of the viscous stress. The critical sonic point is located inside the innermost stable circular orbit of a test particle at 3 Schwarzschild radii. In these solutions, the transition from subsonic to supersonic flow results from pressure effects and not viscous effects. The linear stability of these disks has been examined in the local approximation. It is found that radiative energy transport and viscous stresses in the radial direction can have important effects. As a result, it is shown that the growth rate of the inertial-acoustic mode reaches a maximum at a critical wavelength.
Electron transport and acceleration processes in the earth's magnetosphere have correspondences to analogous processes affecting electrons in the solar magnetosphere (i.e., heliosphere). Energetic electrons in planetary magnetospheres and the heliosphere are test particles probing transport and acceleration dynamics with minimal effects on dominant magnetic field configurations. Parallels are discussed relating to electron entry into the magnetospheres from interplanetary and interstellar space, circulatory transport processes, and acceleration by electric fields in boundary regions including shocks and magnetotails.
Measurements from the space shuttle flights have revealed that a large spacecraft in a low earth orbit is accompanied by an extensive gas cloud which is primarily made up of water. The charge exchange between the water molecule and the ionospheric O(+) ions produces a water ion beam traversing downstream of the spacecraft. In this report we present results from a study on the generation of plasma waves by the interaction of the water ion beams with the ionospheric plasma. Since velocity distribution function is key to the understanding of the wave generation process, we have performed a test particle simulation to determine the nature of H2O(+) ions velocity distribution function. The simulations show that at the time scales shorter than the ion cyclotron period tau(sub c), the distribution function can be described by a beam. On the other hand, when the time scales are larger than tau(sub c), a ring distribution forms. A brief description of the linear instabilities driven by an ion beam streaming across a magnetic field in a plasma is presented. We have identified two types of instabilities occurring in low and high frequency bands; the low-frequency instability occurs over the frequency band from zero to about the lower hybrid frequency for a sufficiently low beam density. As the beam density increases, the linear instability occurs at decreasing frequencies below the lower-hybrid frequency. The high frequency instability occurs near the electron cyclotron frequency and its harmonics.
Light gas gun studies have shown that 6 km/s solid mineral and glass test particles can be successively captured in 0.05 g cm(exp -3) aerogel without severe heating or fragmentation. In spite of this work, there has been uncertainty in the performance of aerogel for hypervelocity capture of real meteoroids. Natural impacts differ from simulations in that the particles are likely to be structurally weak and they typically impact at higher velocity that can be simulated in the laboratory. We are fortunate now to have had two successful capture experiments using aerogel exposed in space. These experiments provide fundamental data for the assessment of the value of silica aerogel for capture of hypervelocity meteoroids from spacecraft. The first experiment used 0.02 g cm(exp -3) aerogel flown on the lid of a Shuttle Get Away Special canister. During its 9 day exposure, the 0.165 m(exp 2) of aerogel in this Sample Return Experiment (SRE) captured two long 'carrot-shaped' tracks and one highly fractured bowl shaped 'crater'. The second collection was with 0.04 m(exp 2) of 0.05 g cm(exp -3) aerogel exposed on ESA's Eureca freeflying spacecraft that was exposed for 11 months before recovery by the Shuttle. The Eureca aerogel exposure consisted of four 10x10 cm module trays that were part of the TiCCE meteoroid collector built by the University of Kent at Canterbury. To date we have found ten 'carrot-shaped' tracks and two 'craters' on this experiment. The longest tracks in both exposures are over 2 mm long. Two of the TiCCE modules had a 0.1 micron Al film suspended a millimeter above the aerogel. On these modules several of the projectiles fragmented during passage through the film producing fields of carrot shaped tracks from the resulting miniature 'meteor' shower. Most of the tracks in these showers have observable particles at their ends. We have extracted one of the carrot track meteoroids and mounted it in epoxy for sectioning. So far the examination of these 14 impacts suggests that low density aerogel is a magic and highly effective media for intact capture of hypervelocity particles in space.
The structure of the quiet-time magnetotail is studied using a test particle simulation. Vlasov equilibria are obtained in the regime where v(D) = E(y) c/B(z) is much less than the ion thermal velocity and are self-consistent in that the current and magnetic field satisfy Ampere's law. Force balance between the plasma and magnetic field is satisfied everywhere. The global structure of the current sheet is found to be critically dependent on the source distribution function. The pressure tensor is nondiagonal in the current sheet with anisotropic temperature. A kinetic mechanism is proposed whereby changes in the source distribution results in a thinning of the current sheet.
The theory of cosmic-ray acceleration by subluminal magnetohydrodynamics (MHD) astrophysical shocks according to the test particle approximation is extended to highly oblique shocks where flow speeds which appear nonrelativistic in the shock rest frame appear on transformation to the de Hoffmann-Teller or E identically equal to 0 frame to have upstream flow speeds approaching c. Monte Carlo simulation shows that relative to the predictions of diffusion theory, as the upstream E identically equal to 0 frame velocity approaches c, flatter spectra and faster acceleration rates occur. These spectral and acceleration time changes are similar to those found for relativistic, parallel MHD shocks and may affect all nonthermal active galactic nuclei emission from relativistic electrons. Also there is an additional means of increasing the upper limit of the cosmic-ray spectrum expected from active galactic nuclei, although the approximation used may not accurately reproduce the spectral shape.
An identifying signature is proposed for the existence and location of the neutral line in the magnetotail. The signature, abrupt density and temperature changes in the Earth-tail direction, was first discovered in test particle simulations. Such temperature variations have been observed in International Sun Earth Explorer (ISEE) data (Huang et al., 1992), but their connection to the possible existence of a neutral line in the tail has not yet been established. The proposed signature develops earlier than the ion velocity space ridge of Martin and Speiser (1988) but can only be seen by spacecraft in the vicinity of the neutral line, while the latter can locate a neutral line remotely.
Spacecraft observations at the Earth's bow shock and at interplanetary shocks have established that the largest fluxes of accelerated suprathermal electrons occur in so-called shock spike events immediately downstream of the shock ramp. Previous theoretical efforts have mainly focused on explaining upstream energetic electron beams. Here we investigate the general motion and acceleration of energetic electrons in a curved, nearly perpendicular shock by numerically integrating the orbits of solar wind halo electrons in shock fields generated by a hybrid simulation (core electron fluid and kinetic ions). Close to the angle Theta(sub Bn) = 90 degs between the upstream magnetic field and shock normal, the calculations result in a (perpendicular) temperature increase proportional to the magnetic field ratio and give the highest phase space densities in the overshoot. For a steep distribution, the temperature change can correspond to an enhancement of the distribution by several orders of magnitude. These results are in agreement with predictions from adiabatic mapping. With smaller angles Theta(sub Bn), the overshoot and downstream densities fall off quickly, because the adiabatic energy gain is less and fewer electrons transmit. The shock curvature also leads to an accumulation of electrons close to 90 degs. Without pitch angle scattering, energization is only significant within a few (approximately 5 to 10 degs) degrees of the point of tangency. However, shock spike events appear to be observed more easily and farther away from 90 degs. Given that over a region of several degrees around 90 degs the theory gives enhancements of up to approximately 4 orders of magnitude, such electrons could in principle account for the typically observed enhancements of 1 to 2 orders of magnitude, if they were distributed over Theta(sub Bn). To test the idea that scattering could efficiently redistribute the energetic electrons, we have conducted test particle simulations in which artificial pitch angle scattering is included. We find that such a process can indeed be very efficient and can explain observations of shock spike events far away from Theta(sub Bn) is approximately 90 degs. The scattering naturally leads to much higher phase space densities at smaller Theta(sub Bn) than what a local one-dimensional mapping would predict and thus can account for an observed discrepancy with adiabatic theory stated in the literature.
Collection of electrons by a satellite at a relatively high potential in low Earth orbit, including the effects of the satellite's orbital motion, remains a poorly understood phenomenon. Using a test particle simulation in which charged particle motion is tracked in prescribed electric fields, we calculate here the current collection and the nature of the energy distribution function of the electrons collected by the satellite, including the effects of the satellite orbital motion. Calculations of the collected current without the orbital motion show an excellent agreement with the current from the Parker-Murphy model, but with the orbital motion of the satellite the current collection is enhanced, the degree of enhancement depending on the size of the satellite sheath extending along the magnetic field line. In the latter case, the flow of electrons around the satellite shows some interesting behaviors including azimuthal E x B drift around the satellite, axial trapping along the magnetic field, and formation of field-aligned flow of electrons in the wake region. The total energy of the collected electrons is approximately equal e phi(sub 0), where phi(sub 0) is the satellite potential, but the partition of the energy into components parallel (W(sub parallel)) and perpendicular (W(sub perpendicular)) to the magnetic field shows interesting features. For the magnetic field along the polar axis Z, the energy distributions near the poles are perfectly field aligned, that is, the parallel energy (W(sub parallel) = e phi(sub 0). The perpendicular component (W(sub perpendicular)) progressively increases toward the 'equator' of the spherical satellite. For a sufficiently large bias potential phi(sub 0), the equatorial energy distributions f(W(sub parallel)) and f(W(sub perpendicular)) show fine structures having multiple peaks. These features of the distribution functions result from the chaotic motion of electrons and the associated coupling between the electron motions parallel and perpendicular to the magnetic field in strongly nonuniform electric fields near the satellite.
Destruction of interstellar dust occurs predominantly in supernova shock waves in the warm neutral/ionized medium (density approximately = 0.25/cu cm, temperature approximately = 10(exp 4) K). Recent theoretical developments and laboratory data for sputtering processes and grain-grain collisional vaporization allows us to better evaluate the grain destruction rate in interstellar shocks in the warm medium. We find that, independent of composition, grain denstruction in supernova blast waves is dominated by nonthermal sputtering for shock velocities greater than 50 km/s and less than or equal to 150 km/s and thermal sputtering at higher shock velocities. We use a detailed scheme for the vaporization of grains colliding at high velocities (v(sub s) greater than or equal to 20 km/s) and show that the grain-grain collision destruction process is only dominant for shock velocities of less than or equal to 50-80 km/s and is less important than previously assumed. Nevertheless, the grain-grain destruction rates are of order 30%-90% of the sputtering rates at v(sub s) greater than 100 km/s and less than 200 km/s and are important in vaporizing the cores of grains. Detailed results for grain destruction as a function of grain size and composition are presented. We also present results for silicon carbide, iron, ice, and porous test particles. For carbonaceous grains we find that the fractional destruction is less than or equal to 0.29, and for silicate it is less than or equal to 0.45, for v(sub s) less than or equal to 200 km/s. We have calculated grain lifetimes, using the three-phase model of the interstellar medium, and find lifetimes of 4 x 10(exp 8) yr for carbonaceous grains and 2.2 x 10(exp 8) yr for silicate grains. Given that the typical stardust injection timescale of 2.5 x 10(exp 9) yr, we conclude that efficient mechanisms for grain growth in the interstellar medium must exist in order that a significant fraction of the refractory elements be incorporated in dust, as observed. Therefore, although our improved model has less vaporization of dust due to grain-grain collisions, sputtering still destroys dust efficiently and grain mantle growth in the interstellar medium is required, a conclusion reached in previous models of grain destruction in the interstellar medium. Carbonaceous mantles on silicate grains can protect the silicate cores from sputtering destruction in interstellar shock waves, provided that the protective mantles can efficiently reform in the interstellar medium. Also, if the grains are porous the postshock grain velocities are lower than for solid particles, and grain destruction is reduced. Porosity and mantling may increase the grain lifetime by factors of approximately 3 and approximately 3-4, respectively. The fraction of interstellar silicon in silicate stardust is therefore less than or equal to 0.25, but more silicon might be depleted in the form of a grain mantle.
The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.
The observations made during the encounter with comet Giacobini Zinner show that the character of MHD turbulence is governed by the magnetosonic (MS) waves generated by the pickup ions via a resonant cyclotron instability. The interaction of cometary ions in the mass loaded solar wind with MS waves propagating away from the comet and oblique to the interplanetary magnetic field (IMF) is investigated using the test particle approach. Ion trajectories, distribution functions, widths of pitch angle scattering and energy diffusion are obtained. We have also investigated the relationship between pitch angle scattering and the three parameters, namely, alpha, the angle between the solar wind flow direction and the ambient magnetic field, theta(sub BK), the angle between the ambient magnetic field and the wave propagation vector, and the ion injection velocity.
The effect of parallel ion acceleration associated with convection was first applied to energization of test particle polar ions by Cladis (1986). However, this effect is typically neglected in 'self-consistent' models of polar plasma outflow, apart from the fluid simulation by Swift (1990). Here we include approximations for this acceleration, which we broadly characterize as centrifugal in nature, in our time-dependent, semikinetic model of polar plasma outflow and describe the effects on the bulk parameter profiles and distribution functions of H+ and O+. For meridional convection across the pole the approximate parallel force along a polar magnetic field line may be written as F(sub cent, pole) = 1.5m(E(sub i))/B(sub i))squared (r(squared)/r(sup 3)(sub i)) where m is ion mass, r is geometric distance; and E(sub i), B(sub i) and r(sub i) refer to the electric and magnetic field magnitudes and geocentric distance at the ionosphere, respectively. For purely longitudinal convection along a constant L shell the parallel force is F(cent. long) = F(sub cent, pole)(1 - (r/(r(sub i)L))(sup 3/2)/(1 - 3r/(4 r(sub i)L))(sup 5/2). For high latitudes the difference between these two cases is relatively unimportant below approximately 5 R(sub E). We find that the steady state O+ bulk velocities and parallel temperatures strongly increase and decrease, respectively, with convection strength. In particular, the bulk velocities increase from near 0 km/s at 4000 km altitude to approximately 10 km/s at 5 R(sub E) geocentric distance for 50-mV/m ionospheric convection electric field. However, the centrifugal effect on the steady O+ density profiles depends on the exobase ion and electron temperatures: for low-base temperatures (T(sub i) = T(sub e) = 3000 K) the O+ density at high altitudes increases greatly with convection, while for higher base temperatures (T(sub i) = 5000 K, T(sub e) = 9000 K), the high-altitude O+ density decreases somewhat as convection is enhanced. The centrifugal force further has a pronounced effect on the escaping O+ flux, especially for cool exobase conditions; as referenced to the 4000-km altitude, the steady state O+ flux increases from 10(exp 5) ions/sq cm/s when the ionospheric convection field E(sub i) = 0 mV/m to approximately 10(exp 7) ions/sq cm/s when E(sub i) = 100 mV/m. The centrifugal effect also decreases the time scale for approach to steady-state. For example, in the plasma expansion for T(sub i) = T(sub e) = 3000 K, the O+ density at 7 R(sub E) reaches only 10(exp -7) of it final value approximately 1.5 hours after expansion onset for E(sub i) = 0. For meridional convection driven by E(sub i) = 50 mV/m, the density at the same time after initial injection is 30-50% of its asymptotic level. The centrifugal acceleration described here is a possible explanation for the large (up to approximately 10 km/s or more) o+ outflow velocities observed in the midlatitude polar magnetosphere with the Dynamics Explorer 1 and Akebono spacecraft.
Test particle calculations are used to compute the effects of gravity and ponderomotive acceleration by shear Alfven wave oscillations on the distribution function of O(+) ions along auroral field lines, assuming an ionospheric Maxwellian source of the ions at 2000 km altitude with approximately 0.5 eV of thermal energy in the parallel component of velocity. The electric field model corresponds to a standing wave oscillation with a frequency approximately 1 Hz in the azimuthal direction superimposed on the background dipole field, in which the wave amplitude is either increasing or decreasing in time. The electric field is taken to be primarily in the perpendicular direction. The time varying wave produces broad distributions with widths of 2 to 10 times the initial 0.5-eV thermal energy of the Maxwellian source, and the density and flux of upward going O(+) ions at one Earth radius are both enhanced in this model. The oxygen ion distribution functions at 1 R(sub E) altitude resulting from interaction with waves whose amplitudes are increasing in time have a more gradual lower energy cutoff than do the distribution functions resulting from decaying waves. The high-energy part of the distribution functions in growing waves reflects the temperature of the Maxwellian source, while the high-energy part of the distributions resulting from decaying waves steepens with time, independent of the source temperature.
The solar wind, a highly tenuous plasma streaming from the Sun into interplanetary space at supersonic speed, is roughly composed of 95% hydrogen and 5% helium by number. All other, heavy elements contribute less than 0.1% by number and thus are truly test particles Nevertheless, these particles provide valuable information not present in the main components. We first discuss the importance of the heavy ions as tracers for processes in the solar atmosphere. Specifically, their relative abundances are found to be different in the solar wind as compared to the photosphere. This fractionation, which is best organized as a function of the first ionization time (FIT) of the elements under solar surface conditions, provides information on the structure of the chromosphere. where it is imparted on the partially ionized material by an atom-ion separation mechanism. Moreover, the charge states of the heavy ions can be used to infer the coronal temperature, since they are frozen-in near the altitude where the expansion time scale overcomes the ionization/recombination time scales. Next, we review the published values of ion abundances in the solar wind, concentrating on the recent results of the SWICS instrument on Ulysses. About 8 elements and more than 20 charge states can be routinely analyzed by this sensor. There is clear evidence that both the composition and the charge state distribution is significantly different in the fast solar wind from the south polar coronal hole, traversed by Ulysses in 1993/94, as compared to the solar wind normally encountered near the ecliptic plane. The fractionation between low- and high-FIT elements is reduced, and the charge states indicate a lower, more uniform coronal temperature in the hole. Finally, we discuss these results in the framework of existing theoretical models of the chromosphere and corona, attempting to identify differences between the low- and high-latitude regions of the solar atmosphere.
Pickup ions in the solar wind are initially are born in ring-beam distributions, i.e. f(v) varies as delta(v(sub perpendicular) - V(sub sw)sin(Theta)) delta(v(sub parallel) - V(sub sw)cos(Theta)), where Theta is the angle between the solar wind velocity and the IMF(Interplanetary Magnetic Field), and V(sub sw) is the solar wind speed. Often the distribution has been presumed to relax to a distribution that is isotropic in Theta and essentially mono-energetic, a shell or a 'bi-spherical distribution.' However solar wind turbulence is capable of heating the ring distribution on the timescale of a few tens of gyroperiods, a timescale not greatly distinct from that required for pitch angle scattering to a shell. To describe this effect, we have performed test-particle studies of the heating/scattering of the ring beam distribution by MHD turbulence, adopting various models for the MHD fluctuations, including slab and fully dynamic 2D and 3D incompressible turbulence. Furthermore, a system composed of a cold ion ring and a background plasma is unstable to several kinetic plasma instabilities. We carried out kinetic simulations of the ring beam distribution, showing that plasma instabilities also rapidly energize and scatter particles. Results will be presented comparing relaxation and heating rates of the ring-beam distribution by the various mechanisms.
In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.