Engineering topics

Miller, Ronald H.

Publications and source records attributed to Miller, Ronald H..

8 records

NASA NTRS1994

A Coulomb collision algorithm for weighted particle simulations

A binary Coulomb collision algorithm is developed for weighted particle simulations employing Monte Carlo techniques. Charged particles within a given spatial grid cell are pair-wise scattered, explicitly conserving momentum and implicitly conserving energy. A similar algorithm developed by Takizuka and Abe (1977) conserves momentum and energy provided the particles are unweighted (each particle representing equal fractions of the total particle density). If applied as is to simulations incorporating weighted particles, the plasma temperatures equilibrate to an incorrect temperature, as compared to theory. Using the appropriate pairing statistics, a Coulomb collision algorithm is developed for weighted particles. The algorithm conserves energy and momentum and produces the appropriate relaxation time scales as compared to theoretical predictions. Such an algorithm is necessary for future work studying self-consistent multi-species kinetic transport.

Miller, Ronald H.↗

NASA NTRS1993

Kinetic simulation of plasma flows in the inner magnetosphere

A one-dimensional hybrid particle code is used to study the interactions between upflowing thermal ions from conjugate ionospheres. The simulation model allows for multiple species, convection of plasmaspheric flux tubes, and Coulomb self-collisions which conserve momentum and energy locally. The model incorporates a variable-flux boundary condition where the flux, at the boundaries, approaches zero as the plasmasphere fills and equilibrium conditions are reached. The effects of two important processes on plasmaspheric refilling have been considered. The first includes convection of the plasmaspheric flux tube. The second is the interaction of ionospheric thermal plasma and particle injection from an external source. Particle injection seems to play an important role in the evolution of the total particle distribution on the early timescales (t less than 1 hour); however, for late timescales (t larger than 8 days) the thermal plasma from the ionosphere dominates the particle distribution.

Miller, Ronald H.↗

NASA NTRS1993

Self-consistent electrostatic potential due to trapped plasma in the magnetosphere

A steady state solution for the self-consistent electrostatic potential due to a plasma confined in a magnetic flux tube is considered. A steady state distribution function is constructed for the trapped particles from the constants of the motion, in the absence of waves and collisions. Using Liouville's theorem, the particle density along the geomagnetic field is determined and found to depend on the local magnetic field, self-consistent electric potential, and the equatorial plasma distribution function. A hot anisotropic magnetospheric plasma in steady state is modeled by a bi-Maxwellian at the equator. The self-consistent electric potential along the magnetic field is calculated assuming quasineutrality, and the potential drop is found to be approximately equal to the average kinetic energy of the equatorially trapped plasma. The potential is compared with that obtained by Alfven and Faelthammar (1963).

Miller, Ronald H.↗

NASA NTRS1993

Radiative properties of a plasma moving across a magnetic field. I - Theoretical analysis. II - Numerical results

The early-time evolution of plasmas moving across a background magnetic field is addressed with a 2D model in which a plasma cloud is assumed to have formed instantaneously with a velocity across a uniform background magnetic field and with a Gaussian density profile in the two dimensions perpendicular to the direction of motion. This model treats both the dynamics associated with the formation of a polarization field and the generation and propagation of electromagnetic waves. In general, the results indicate that, to zeroth order, the plasma cloud behaves like a large dipole antenna oriented in the direction of the polarization field which oscillates at frequencies defined by the normal mode of the system. Radiation damping is shown to play an important role in defining the plasma cloud evolution, causing a rapid decay of the polarizaiton field and a loss of plasma kinetic energy and momentum on time scales comprable to several ion gyroperiods. Scaling laws are derived for the plasma momentum and energy loss rates, and predictions for the braking time, the amplitude and spectrum of the radiation field, and the total radiated power are presented for conditions relevant to the recent Combined Release and Radiation Effects Satellite experiments.

Roussel-Dupre, Robert↗

NASA NTRS1992

A second-order theory for transverse ion heating and momentum coupling due to electrostatic ion cyclotron waves

A second-order theory for electrostatic instabilities driven by counterstreaming ion beams is developed which describes momentum coupling and heating of the plasma via wave-particle interactions. Exchange rates between the waves and particles are derived, which are suitable for the fluid equations simulating microscopic effects on macroscopic scales. Using a fully kinetic simulation, the electrostatic ion cyclotron instability due to counterstreaming H(+) beams has been simulated. A power spectrum from the kinetic simulation is used to evaluate second-order exchange rates. The calculated heating and momentum loss from second-order theory is compared to the numerical simulation.

Miller, Ronald H.↗

NASA NTRS1991

The directional dependence of magnetic fluctuations generated by cometary ion pickup

The properties of low frequency magnetic fluctuations generated by cometary ion pickup are examined by means of one-dimensional hybrid simulations, in which newborn ions are created at a constant rate. The helicity and direction of propagation of magnetic fluctuations are investigated for various cometary ion injection angles, alpha, relative to the solar wind magnetic field. The parameter eta represents the relative contribution of wave energy density propagating in the direction away from the comet, parallel to the beam. For small (quasi-parallel) injection angles, alpha is about 0 deg and eta is of order unity, while for larger (quasi-perpendicular) angles, alpha is about 90 deg and eta is about 0.5. At intermediate angles, alpha is about 60 deg, eta can vary between 0 and 1, depending on the wave number. The wave properties are consistent with the instabilities expected from linear theory.

Miller, Ronald H.↗

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

Pitch-angle scattering of cometary ions - Computer simulations

The relationships between the rate of pitch-angle scattering induced by wave-particle interactions and two parameters, alpha, and the cometary ion interjection rate, Lambda-j, are studied with the aid of one-dimensional electromagnetic hybrid computer simulations of homogeneous plasmas. The simulations are done with relatively weak ion injection rates commensurate with those in the distant environment of the Comet Halley. In response to the linear temporal growth of the fluctuating magnetic field energy, the injected ions pitch-angle scatter toward isotropy in both the quasi-parallel and quasi-perpendicular regimes. For the injection of cometary oxygen ions, the simulations show pitch-angle scattering rates that increase as the square root of the injection rate. The oxygen ion pitch-angle scattering rates are roughly twice as large at perpendicular injection, again in qualitative agreement with observations.

Gary, S. P.↗