Search NASA⌕ Search

Engineering topics

Pontius, D. H., Jr.

Publications and source records attributed to Pontius, D. H., Jr..

Field lines and magnetic surfaces in a two-component slab/2D model of interplanetary magnetic fluctuations

A two-component model for the spectrum of interplanetary magnetic fluctuations was proposed on the basis of ISEE observations, and has found an intriguing level of application in other solar wind studies. The model fluctuations consist of a fraction of 'slab' fluctuations, varying only in the direction parallel to the locally uniform mean magnetic field B(0) and a complement of 2D (two-dimensional) fluctuations that vary in the directions transverse to B(0). We have developed an spectral method computational algorithm for computing the magnetic flux surfaces (flux tubes) associated with the composite model, based upon a precise analogy with equations for ideal transport of a passive scalar in planar two dimensional geometry. Visualization of various composite models will be presented, including the 80 percent 2D/ 20 percent slab model with delta B/B(0) approximately equals 1 and a minus 5/3 spectral law, that is thought to approximately represent a snapshot of solar wind turbulence. Characteristically, the visualizations show that flux tubes, even when defined as regular on some plane, shred and disperse rapidly as they are viewed along the parallel direction. This diffusive process, which generalizes the standard picture of field line random walk, will be discussed in detail. Evidently, the traditional picture that flux tubes randomize like strands of spaghetti with a uniform tangle along the axial direction is in need of modification.

Matthaeus, W. H.↗

Test-particle simulations in increasingly strong turbulence

Quasi-linear theory supposes that the energy in resonant fluctuations is small compared to that in the mean magnetic field. This is evident in the fact that the zeroth-order particle trajectories are helices about a mean field B(sub o) that is spatially uniform over many correlation lengths. However, in the solar wind it is often the case that the fluctuating part of the field is comparable in magnitude to the mean part. It is generally expected that quasi-linear theory remains viable for particles that are in resonance with a region of the fluctuation spectrum having only small energy density, but even so, care must be taken when comparing simulations to theoretical predictions. We have performed a series of test-particle simulations to explore the evolution of ion distributions in turbulent situations with varying levels of magnetic fluctuations. As delta-B/B(sub o) is increased the distinctions among absolute pitch angle (defined relative to B(sub o)), local pitch angle (defined relative to B(x)), and magnetic moment become important, some of them exhibiting periodic sloshing unrelated to the nonadiabatic processes of interest. Comparing and contrasting the various runs illustrates the phenomena that must be considered when the premise underlying quasi-linear theory are relaxed.

Pontius, D. H., Jr.↗

Heating and scattering of ring-beam distributions by turbulence

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.

Gray, P. C.↗

MHD simulations of boundary layer formation along the dayside Venus ionopause due to mass loading

A two-dimensional magnetohydrodynamic (MHD) simulation of mass-loaded solar wind flow around the dayside of Venus is presented. For conditions appropriate to a low-altitude ionopause the simulations show that mass loading from the pickup of oxygen ions produces a boundary layer of finite thickness along the ionopause. Within this layer the temperatures exhibit strong gradients normal to and away from the ionopause. Furthermore, there is a shear in the bulk flow velocity across the boundary layer, such that the (predominantly tangential) flow decreases in speed as the ionopause is approached and remains small along the ionopause, consistent with Pioneer Venus observations. The total mass density increases significantly as the flow approaches the ionopause, where the contribution of O(+) to the total number density is a few percent. Numerical simulations are carried out for various mass addition rates and demonstrate that the boundary layer develops when oxygen ion production exceeds approximately 2 x 10(exp 5)/cu m/s. For the upstream solar wind parameters and mass loading rates chosen for these simulations, the results are consistent with observations made on the dayside of Venus for average ionopause conditions near 300 km.

Mcgary, J. E.↗

Inertial currents in isotropic plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasmas, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MHD plasma. Solutions are developed by taking the MHD limit ot two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

Inertial Currents in Isotropic Plasma

The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Heinemann, M.↗

Particle drift in the Earth's plasma sheet

We generalize the derivation of the average gradient/curvature-drift for a flux tube filled with an isotropic distribution of particles at specified kinetic energy. The present treatment is restricted to a two-dimensional magnetic field with zero electric field, but it includes all chaotic and Speiser orbits, which do not correspond to the simple picture of gradient/curvature drift. We assume that particles are evenly distributed throughout the regions of phase space allowed by their energy and canonical momentum. This assumption is closely related but not exactly equivalent to the assumption of isotropic pitch-angle distribution. Our derivation assumes that the maximum Larmor radius is small compared to the scale length for equatorial variations in the flux tube volume, but it does not involve any restrictions on the curvature of the field line. The resulting expression for the drift rate is valid for situations where the particle drift velocity is comparable to the thermal speed in some regions. The apparent implication of this generalized treatment is that the existence of very complex non-adiabatic particle trajectories in the plasma sheet may not invalidate previous estimates of the average rate of particle drift out the sides of the tail, estimates that were made under the assumption of simple guiding-center drifts.

Wolf, R. A.↗

On Vasyliunas's equivalent conductivity formalism

The Vasyliunas's (1972) equivalent conductivity formalism (ECF) for representing the coupling of the ionosphere and the magnetosphere is discussed, and a new, simpler, derivation is presented of the ECF, in which certain of the underlying assumptions and their implications are made transparent. The derivation presented indicates that the only role of the ions in the ECF is to insure quasi-neutrality. It is shown that the ECF is not as robust as usually assumed and that caution must be used to insure that reasonable results are obtained.

Pontius, D. H., Jr.↗

Transient flux tubes in the terrestrial magnetosphere

Small, isolated density depletions are postulated to be introduced into the magnetotail by an uneven plasma loading process in the far tail. Such bubbles would be displaced earthward by an interchange process, and a significant depletion would move faster than either the gradient-curvature drift speed or the average convection speed. The dissolution of a bubble into the background results in a local reduction in plasma density, thereby violating ideal MHD and the subsequent conservation of plasma content. If a sufficient amount of magnetic flux in the far tail is contained in bubbles, then the net effect may be a reduction of plasma pressure in the near tail sufficient to resolve the pressure balance inconsistency.

Pontius, D. H., Jr.↗

Departure from corotation of the Io plasma torus - Local plasma production

The departure of the Jovian magnetosphere from rigid corotation is adequately explained by outward plasma transport at distances where L is greater than approximately 10. The departure of 5% observed in the Io plasma torus, however, is too large to be accounted for simply by plasma transport. Local plasma production is proposed as the main factor determining the corotation lag in the torus. The outward pick-up current provided by ionization of neutral atoms is calculated and related to the current produced in the ionosphere by the corotation lag. This leads to an expression giving the corotation lag of the torus as a function of radial distance. Charge transfer is found to be an important process, allowing the majority of the torus mass to be ejected from the magnetosphere in a neutral state. Thus, the mass loading rate is found to be several times that inferred from examination of the corotation lag associated with outward plasma transport.

Pontius, D. H., Jr.↗