Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid simulations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Pitch angle scattering of cometary ions into monospherical and bispherical distributions

Low frequency magnetic fluctuations generated by solar wind/cometary ion distributions cause pitch-angle scattering of cometary ions. The nature of this pitch angle scattering is examined by means of one-dimensional hybrid simulations, in which newborn ions are created at a constant rate, and for various injection angles, alpha, relative to the magnetic field. Pitch angle scattering in the quasi-perpendicular regime (alpha above 60 deg) results in a bispherical velocity distribution while in the quasi-parallel regime (alpha between 0 and 60 deg), pitch angle scattering can result in either a monospherical or bispherical distribution, depending on the injection angle and initial beam velocity.

Miller, Ronald H.↗
NASA NTRS1991

Chaos and ion heating in a slow shock

An ion heating mechanism is proposed of slow shocks, which is associated with the chaotic motion of particles in the downstream wave field. For a coherent electromagnetic wave propagating along the downstream magnetic field, corresponding to switch-off shocks, the particle motions are not chaotic. For an oblique wave, the interaction between the particles and the wave field may lead to chaotic particle motions. Such particles may be greatly thermalized within one wavelength after they are incident into the downstream wave field. The results can be used to explain the existence of the critical intermediate Mach number observed in the hybrid simulations.

Lin, Y.↗
NASA NTRS1991

The directional dependence of cometary magnetic energy density in the quasi-parallel and quasi-perpendicular regimes

The direction of propagation of low frequency magnetic fluctuations generated by cometary ion pick-up is examined by means of 1D electromagnetic hybrid simulations. The newborn ions are injected at a constant rate, and the helicity and direction of propagation of magnetic fluctuations are explored for cometary ion injection angles of 0 and 90 deg relative to the solar wind magnetic field. The parameter eta represents the relative contribution of wave energy propagating in the direction away from the comet, parallel to the beam. For small (quasi-parallel) injection angles eta was found to be of order unity, while for larger (quasi-perpendicular) angles eta was found to be of order 0.5.

Miller, R. H.↗
NASA NTRS1991

Electrostatic potential jump across fast-mode collisionless shocks

The electrostatic potential jump across fast-mode collisionless shocks is examined by comparing published observations, hybrid simulations, and a simple model, in order to better characterize its dependence on the various shock parameters. In all three, it is assumed that the electrons can be described by an isotropic power-law equation of state. The observations show that the cross-shock potential jump correlates well with the shock strength but shows very little correlation with other shock parameters. Assuming that the electrons obey an isotropic power law equation of state, the correlation of the potential jump with the shock strength follows naturally from the increased shock compression and an apparent dependence of the power law exponent on the Mach number which the observations indicate. It is found that including a Mach number dependence for the power law exponent in the electron equation of state in the simple model produces a potential jump which better fits the observations. On the basis of the simulation results and theoretical estimates of the cross-shock potential, it is discussed how the cross-shock potential might be expected to depend on the other shock parameters.

Mandt, M. E.↗
NASA NTRS1992

Kinetic properties of mirror waves in magnetosheath plasmas

Linear and nonlinear properties of waves excited by the mirror instability in high beta, low anisotropy plasmas characteristic of the magnetosheath are investigated using linear theory and one-dimensional hybrid simulations. The mechanisms for wave growth and saturation at low amplitudes are discussed. A new method is considered for generating the large amplitude mirror waves observed in the magnetosheath based on external compression of magnetic flux tubes. Simulations in which the anisotropy is maintained by recycling the ions shows this process can inhibit the growth of ion cyclotron waves and enhances the growth of mirror waves.

Mckean, M. E.↗
NASA NTRS1992

Mirror and ion cyclotron anisotropy instabilities in the magnetosheath

The kinetic properties of the ion cyclotron anisotropy and the mirror instabilities are investigated using one-dimensional hybrid simulations. It is found that, for moderate values of the ion beta and the proton temperature anisotropy, the two instabilities produce similar levels of turbulence. For high values of beta or temperature anisotropy, the ion cyclotron instability produces higher fluctuation levels of turbulence than does the mirror instability.

Mckean, M. E.↗
NASA NTRS1993

Nonlinear evolution of electromagnetic ion beam instabilities

A comparative study of the ion/ion right-hand resonant instability and the ion/ion nonresonant instability is carried out to investigate and contrast their properties. Linear analysis demonstrates that the nonresonant instability becomes resonant if the density of the ion beam is sufficiently high. Hybrid simulations show that both the resonant and nonresonant instabilities result in the formation of nonlinear pulses, called pulsations, but with distinct features. For example, the pulsations generated by the resonant instability have a positive correlation with the ion density, while those generated by the nonresonant instability are likely to have a relatively weak negative correlation. The waves generated by the nonresonant instability are subject to a parametric decay instability and tend to form a state of condensate where the turbulence becomes nearly monochromatic.

Akimoto, K.↗
NASA NTRS1993

Near-specular reflection of ions at quasi-parallel shocks

One-dimensional hybrid simulations and a semianalytical model of the shock front are employed to investigate the source regions in the incident ion phase space of reflected and transmitted ions, the relative importance of the electric and magnetic forces in the reflection of incident ions, and how these characteristics change in time. The phase space origin of reflected particles and the fraction of incident ions that reflect are found to depend on the electromagnetic field structure of the shock front at the time the ions encounter it. The reflection fraction is maximized when the electric field along the shock normal (Ex) and the noncoplanar magnetic field (By) are at their maximum values. When Ex and By are large, the reflection process produces a beam that is cooler, more dense, and closer to specular than when they are small.

Mckean, M. E.↗
NASA NTRS1993

Propagation characteristics of waves upstream and downstream of quasi-parallel shocks

The propagation characteristics of waves upstream and downstream of quasi-parallel shocks are investigated by using 2D hybrid simulations. At low Alfven Mach numbers, M(A) below about 2, the shock is initially associated with upstream phase-standing whistlers. At later times, backstreaming ions excite longer-wavelength whistlers via the right-hand resonant ion/ion instability. These waves propagate along the magnetic field at a group velocity no smaller than the upstream flow speed, so that the waves remain in the upstream region. At higher MA (above about 3), these waves are convected back into the shock, causing its reformation and downstream perturbations. Shock transmitted waves mode-convert into Alfven/ion-cyclotron waves which have a wave vector along the shock normal (pointing upstream) and convect downstream. The 2D simulation results confirm our earlier suggestion that the upstream waves should be field aligned, and that their convection into the downstream is associated with linear mode conversion into the Alfven/ion-cyclotron branch.

Krauss-Varban, D.↗
NASA NTRS1994

The proton cyclotron instability and the anisotropy/B inverse correlation

Spacecraft observations in the strongly compressed subsolar magnetosheath show an inverse correlation between the proton temperature anistrophy (T(sub perpendicular to p)/T(sub parallel to p) greater than 1 where perpendicular to and parallel to denote directions perpendicular and parallel to the background magnetic field) and the parallel proton beta (beta(sub parallel to p). This manuscript uses one-dimensional hybrid simulations of the proton cyclotron anisotropy instability in homogeneous electron-proton plasmas to study this correlation which may represent a limited closure relation for fluid theories of anisotropic space plasmas. The emphasis is on driven simulation which increase the temperature anisotropy by periodically reducing the magnetic-field-aligned velocities of the protons. The late-time states from ensembles of both initial value and driven simulations yield very similar expressions for the proton anisotropy/beta(sub parallel to P) inverse correlation, and provide a basis for explaining differences between sheath observations from different spacecraft. The driven simulations also yield expressions for the maximum instability growth rate and the fluctuating field energy as functions of beta (sub parallel to p) and a parameter characterizing the anisotropy driver.

Gary, S. Peter↗
NASA NTRS1994

Reconnection layer at the flank magnetopause in the presence of shear flow

We present hybrid simulations of reconnection layer at the flank magnetopause, where a large plasma flow speed is present in the magnetosheath. It is found that there exists a threshold flow speed v(sub *) such that for the magnetosheath flow speed v(sub s) less than v(sub *) (v(sub s) greater than v(sub *)), the rotational discontinuity with a larger field rotation angle exists on the magnetosheath (magnetospheric) side of the reconnection layer. The threshold speed is found to be v(sub *) = v(sub Am) - v(sub As), where v(sub Am) (v(sub As)) is the Alfven speed on the magnetospheric (magnetosheath) side of the reconnection layer. Furthermore, for v(sub s) much less than v(sub *), the rotational discontinuity on the magnetosheath side is very thin, and an accelerated high-speed flow is located earthward of the rotational discontinuity, as observed at the dayside magnetopause. For v(sub s) approximately v(sub *), the magnetic field transition region is thick, and the accelerated flow is present in the entire field transition region, as observed at the flank magnetopause.

Lin, Y.↗
NASA NTRS1994

Low-freqency instabilities and the resulting velocity distributions of pickup ions at comet Halley

The interaction between the solar wind and newborn cometary ions is studied using a new analytical theory as well as one- and two-dimensional hybrid simulations. Using the observed parameters upstream of the comet Halley, a detailed study of wave excitation and the resulting particle distributions is presentd. Linear theory as well as simulations show that a variety of modes such as the fast magnetosonic mode, high frequency whistlers and obliquely propagating Alfven ion cyclotron waves can be excited. However, parallel propagating waves are found to be dominant in the wave spectrum and to control the scattering of the pickup ions. Several features of the observed distributions of pickup protons are explained. In particular, it is shown that the observed asymmetric pitch distribution for the pickup protons is due to the small saturation amplitude of the waves for the given parameters. Water group associated waves can lead to energy diffusion and further pitch angle scattering of protons. This effect is most likely to be important in the vicinity of the bow shock of comet Halley where the density of water group ions becomes comparable to that of protons. It is shown that the observed increase in the radius of the proton velocity shell just outside the bow shock can be due to water group waves. The nearby isotropic proton pitch angle distribution observed by Neugebauer et al. (1989) just outside the bow shock may, however, be related to the presence of a rotational discontinuity which has been identified in the magnetic field data. Just outside the bow shock, simulations show that parallel propagating water group waves can steepen with attached whistler wave packets. The steepening process at parallel propagation is a transient effect, in an important contrast to the case of steepening at oblique angles. The smaller beam densitites at comet Halley appears to be the main reason not only why waves at comet Halley have smaller amplitudes but also why oblique, steepening magnetosonic waves have not been detected at comet Halley, whereas they have been seen at comet Giacobini-Zinner. Halley, whereas they have been seen at comet Giacobini-Zinner.

Karimabadi, H.↗
NASA NTRS1994

Electron acceleration at nearly perpendicular collisionless shocks. 3: Downstream distributions

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.

Krauss-Varban, D.↗
NASA NTRS1994

Transition to whistler mediated magnetic reconnection

The transition in the magnetic reconnection rate from the resistive magnetohydrodynamic (MHD) regime where the Alfen wave controls reconnection to a regime in which the ions become unmagnetized and the whistler wave mediates reconnection is explored with 2-D hybrid simulations. In the whistler regime the electrons carry the currents while the ions provide a neutralizing background. A simple physical picture is presented illustrating the role of the whistler mediated reconnection is calculated analytically. The development of an out-of-plane component of the magnetic field is an observable signature of whistler driven reconnection.

Mandt, M. E.↗
NASA NTRS1994

Ion injection and acceleration at quasi-perpendicular shocks

We present and discuss results of a new model for ion injections and acceleration at quasi-perpendicular collisionless shocks. We use the one-dimensional hybrid simulation (kinetic ions/fluid electrons) and impose an assumption on the ion motion so that diffusion across the magnetic field (normal to the shock front) is possible. These motions are otherwise suppressed by both one- and two-dimensional simulations. We find that, even in strictly perpendicular shocks, when scattering normal to the field is included, a fraction of the incident ions are accelerated to suprathermal energies. When reasonable scattering times are considered, only pickup ions are injected, whereas thermal solar wind ions are not. The acceleration of these ions is very rapid. We have found that a few of the initially low-energy pickup ions can reach many tens to a few hundred times the plasma ramming energy in less than 100 gyroperiods. Furthermore, highly field aligned energetic ions are found to exist upstream of the slightly oblique shock. The most direct applications of this study is toward the interpretation of observations of both solar wind and interstellar pickup ion distributions in the vicinity of interplanetary shocks which are most often quasi-perpendicular. This work also directly addresses a fundamental issue with regard to our current understanding of the anomalous component of cosmic rays.

Giacalone, J.↗
NASA NTRS1993

Kinetic physics of the mirror instability

Kinetic mechanisms for the growth and saturation of the mirror instability are described using one-dimensional hybrid simulations. Two parameter regimes are considered. In the first regime, a relatively small ion anisotropy excites a slowly growing instability that produces small-amplitude waves; most ions respond to the waves as an adiabatic fluid. In the second regime a large anisotropy excites a rapidly growing instability that generates large-amplitude waves; the response of many ions in this case is nonadiabatic. The difference in ion response is due to the relative importance of two ion populations, resonant and nonresonant. Resonant ions, those ions with low velocities parallel to the background magnetic field, contribute to the growth of the instability as a result of their gyrointeractions with the noncoplanar component of the waves electric field and respond to the mirror waves nonadiabatically. Nonreasonant ions, those with large parallel velocities, respond as an adiabatic fluid. In both regimes, ion anisotropy is reduced by means of the magnetic mirror force; in the second regime, the anisotropy reduction reduces the free energy available for wave growth and leads to the saturation of the mirror instability.

Mckean, M. E.↗
NASA NTRS1995

Wave and ion evolution downstream of quasi-perpendicular bow shocks

Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Mckean, M. E.↗
NASA NTRS1994

Hot proton anisotropies and cool proton temperatures in the outer magnetosphere

The hot protons of the outer magnetosphere typically exhibit a temperature anistropy such that T(sub perp)/T(sub parallel) greater than 1, where perpendicular and parallel symbols denote directions relative to the background magnetic field. If this anisotrpy is sufficiently large, the electomagneitc proton cyclotron anistropy instability may be excited. This instability is studied using linear Vlasov theory and one-dimensional hybrid simulations for a homogeneous plasma model representative of conditions in the outer magnetosphere with a hot anisotropic proton component (denoted by subscript h) and a cool, initially isotropic proton component (subscript c). Linear theory yields an instability threshold condition on the hot proton temperature anistropy where as the simulations imply an upper bound on T(sub perp h)/T(sub parallel h); both the threshold and the upper bound have similar scaling with the maximum growth rate gamma (sub m), the parallel beta of the hot component, beta(sub parallel h), and the relative density of the hot component n(sub h)/n(sub e). An anlysis of plasma observations from the Los Alamos magnetospheric plasma analyzer (MPA) in geosynchronous orbits finds that the maximum value of the hot proton temperature anisotropy approximately satisfies the predicted scaling with beta(sub parallel h) and nu(sub h)/n(sub e) and yields the proportionality factor that quantifies this upper bound. The simulations are also used to examine the heating of the cool proton cyclotron instability. The simulations yield a scaling for the dimensionless late-time cool proton average temperature T(sub c)/T(sub parallel h) as (n(sub h)/n(sub e))/beta(sub parallel h exp 0.5). Analysis of MPA data shows that the observed values of T(sub c)/T(sub parallel h) have similar scaling and again yield the proportionality factor which quantifies this relationship.

Gary, S. Peter↗