Engineering topics

Quest, K. B.

Publications and source records attributed to Quest, K. B..

At least 19 records

NASA NTRS1995

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗

NASA NTRS1994

On the energy principle and ion tearing in the magnetotail

We re-examine the use of the energy principle as applied to the tearing instability in the magnetotail. We demonstrate that when a magnetic field component normal to the current sheet is present, electron pitch-angle diffusion (PAD) either by micro-turbulence or by chaotic orbits cannot remove the strong stabilization of the tearing mode caused by electron compressibility. We find that our conclusions are in agreement with those of Pellat et al. (1991), who argued on the basis of canonical P(sub y) conservation that the stabilization of the ion tearing mode cannot be removed by the introduction of PAD. Our results are at variance with those of Kuznetsova and Zelenyi (1991), who argued that the application of the energy principle used by Pellat et al. (1991) is incorrect, and that tearing is in fact unstable in the limit of strong PAD. We show that the disagreement between these two studies can be traced to an incorrect orbit evaluation first introduced by Coroniti (1980) and subnsequently used by Kuznetsova and Zelenyi (1991).

Brittnacher, M.↗

NASA NTRS1994

Mode properties of low-frequency waves: Kinetic theory versus Hall-MHD

In fluid theory, the ordering of low-frequency modes in a homogeneous plasma is based on the phase velocity, since modes do not intersect each other in dispersion diagrams as a function of wavenumber or other parameters. In linear kinetic theory, modes cross each other. Thus a consistent and useful classification should be based on the physical properties of the modes instead. This paper attempts such a classification by documeting the dispersion and general mode properties of the low-frequency waves (omega much less than (OMEGA(sub ci) OMEGA(sub ce) (exp 1/2)), where OMEGA(sub ci), OMEGA(sub ce) are the cyclotron frequencies of the ions and electrons, respectively) in kinetic theory, and by comparing them to the results of two-fluid theory. Kinetic theory gives a seperate Alfven/ion-cyclotron (A/IC) wave with phase speed Omega/k approximately = v(sub A) cos theta for omega much less than OMEGA(sub ci), where v(sub A) is the Alfven velocity and theta the angle of propagation between wave vector k and background magnetic field B(sub o). For a given wavenumber, the magnetosonic mode is a double-valued solution with a singular point in theta, beta parameter space, where beta is the ratio of thermal pressure to magnetic pressure. It is shown that a branch cut starting at the singular point theta approximately 30 deg, beta approximately 3 and leading to larger beta gives a practical and consitent seperation of this double-valued magnetosonic solution. Selection of this branch cut results in a moderately damped fast/magnetos onic and a heavily damped slow/sound wave. A comprehensive review of the polarization, compressibility and other mode properties is given and shown to be consistent with the selected branch cut. At small wavenumbers, the kinetic mode properties typically start to deviate significantly from their fluid counterparts at beta approximately 0.5. At larger beta, there is no longer a consistent correspondence between the fluid and kinetic modes. Kinetic theory also dictates the use of different mode properties to distinguish between them in observational data. For example, the phase between the density and magnetic field perturbation may become useless at high beta, whereas the direction of the magnetic field perturbations with respect to k and B(sub o) remains a useful characteristic. Two quantities based on this characteristic are suggested and are shown to be useful also to distinguish between the mirror mode and A/IC waves in a plasma with temperature anisotropy.

Krauss-Varban, D.↗

NASA NTRS1991

Weak, quasiparallel profiles of earth's bow shock - A comparison between numerical simulations and ISEE 3 observations on the far flank

Over 200 crossings of the distant downwind flanks of earth's magnetosonic bow shock by ISEE 3 included many cases of weak, or low Mach number, quasi-parallel shocks. A consistent feature of the magnetic field profiles was the presence of large amplitude, near periodic to irregular transverse oscillations downstream from even the weakest Q-parallel shocks. Large downstream perturbations with whistler-like features similar to those of the observations appear in 1D simulations when the Alfven Mach number M(A) is greater than 2.5 but not when M(A) = 2.1. The observed cases with downstream waves also occurred when M(A) is greater than about 2.5, suggesting the importance of the Alfven as opposed to magnetosonic Mach number in determining the signature of weak, Q-parallel shocks.

Greenstadt, E. W.↗

NASA NTRS1991

Two-dimensional simulations of the ion/ion acoustic instability and electrostatic shocks

A newly developed 2D electrostatic code with particle ions and Boltzmann electrons is used to investigate the details of the ion/ion acoustic instability and the structure of electrostatic shocks. The simulation results show that, for the parameters relevant to the plasma sheet boundary layer, the saturation mechanism of the ion/ion acoustic instability is ion trapping. It is also shown that the 2D structure of electrostatic shocks is considerably different from that suggested by previous 1D simulations. The main reason for this difference is the presence of shock reflected ions, which through the ion/ion acoustic instability lead to the generation of large amplitude waves in the upstream region propagating obliquely to the shock normal. These waves play an important role in the shock dissipation process.

Karimabadi, H.↗

NASA NTRS1990

Low Mach number parallel and quasi-parallel shocks

The properties of low-Mach-number parallel and quasi-parallel shocks are studied using the results of one-dimensional hybrid simulations. It is shown that both the structure and ion dissipation at the shocks differ considerably. In the parallel limit, the shock remains coupled to the piston and consists of large-amplitude magnetosonic-whistler waves in the upstream, through the shock and into the downstream region, where the waves eventually damp out. These waves are generated by an ion beam instability due to the interaction between the incident and piston-reflected ions. The excited waves decelerate the plasma sufficiently that it becomes stable far into the downstream. The increase in ion temperature along the shock normal in the downstream region is due to superposition of incident and piston-rflected ions. These two populations of ions remain distinct through the downstream region. While they are both gyrophase-bunched, their counterstreaming nature results in a 180-deg phase shift in their perpendicular velocities.

Omidi, N.↗

NASA NTRS1990

Re-forming supercritical quasi-parallel shocks. II - Mechanism for wave generation and front re-formation

This paper continues the study of Thomas et al. (1990) in which hybrid simulations of quasi-parallel shocks were performed in one and two spatial dimensions. To identify the wave generation processes, the electromagnetic structure of the shock is examined by performing a number of one-dimensional hybrid simulations of quasi-parallel shocks for various upstream conditions. In addition, numerical experiments were carried out in which the backstreaming ions were removed from calculations to show their fundamental importance in reformation process. The calculations show that the waves are excited before ions can propagate far enough upstream to generate resonant modes. At some later times, the waves are regenerated at the leading edge of the interface, with properties like those of their initial interactions.

Winske, D.↗

NASA NTRS1988

On the origin of hot diamagnetic cavities near the earth's bow shock

The origin of hot diamagnetic cavities (HDCs) observed occasionally upstream from the earth's bow shock is investigated by examining the results of November 16, 1977, observation, when four of these events occurred on a single day, as well as plasma and field data from that day. The results suggest that HDCs may form as a result of an unusually strong interaction between shock-reflected ions and the incoming solar wind. It is proposed that this interaction stems from a temporary and localized reflection of a larger-than-normal fraction of the incident ions, which is stimulated by sudden changes in the upstream field orientation; the consequences of such a temporary overreflection are found to be consistent with many of the observed features of HDCs, including the strong slowing, deflection, and heating of the flow, as well as the localization, internal recoveries, and occasional formation upstream from the shock itself.

Thomsen, M. F.↗

NASA NTRS1988

Theory and simulation of collisionless parallel shocks

This paper presents a self-consistent theoretical model for collisionless parallel shock structure, based on the hypothesis that shock dissipation and heating can be provided by electromagnetic ion beam-driven instabilities. It is shown that shock formation and plasma heating can result from parallel propagating electromagnetic ion beam-driven instabilities for a wide range of Mach numbers and upstream plasma conditions. The theoretical predictions are compared with recently published observations of quasi-parallel interplanetary shocks. It was found that low Mach number interplanetary shock observations were consistent with the explanation that group-standing waves are providing the dissipation; two high Mach number observations confirmed the theoretically predicted rapid thermalization across the shock.

Quest, K. B.↗

NASA NTRS1987

On the noncoplanarity of the magnetic field within a fast collisionless shock

Within the magnetic ramp of fast collisionless plasma shocks observed with spacecraft instruments and simulated numerically, the magnetic field undergoes an excursion out of the plane of coplanarity. This rotation is consistently in the direction such that the electrostatic potential jump across the shock, as measured in the de Hoffman-Teller frame of the reference (HTF), is about 2-6 times smaller than the electrostatic potential jump measured in the normal incidence frame. The preferred direction is consistent with a basic whistler mode transition between the upstream and downstream orientations. The potential jump in the HTF is considerably smaller than the change in bulk flow energy across the shock, confirming the recent suggestion that magnetic forces contribute importantly to the slowing of the plasma in that frame. A further consequence is that suprathermal particles leaking back into the upstream region across the shock do not gain much energy from the cross-shock electric field.

Thomsen, M. F.↗

NASA NTRS1986

Electromagnetic ion beam instabilities - Comparison of oneand two-dimensional simulations

Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principal features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.

Winske, D.↗

NASA NTRS1986

Simulations of high Mach number perpendicular shocks with resistive electrons

A simulation code which models the ions as microparticles and the electrons as a resistive massless fluid is employed to study the structure of high Mach number perpendicular shocks. It is found that stable stationary shock solutions can be obtained for Alfven Mach numbers (M sub A) between 5 and 60 for upstream plasmas where the ratio of the plasma pressure to the magnetic pressure is 1, providing that the upstream resistive diffusion length is much smaller than the ion inertial length. For much larger resistive diffusion lengths, the magnetic field overshoot is damped, and the imbalance in the electron momentum equation results in a periodic fluctuation of the fraction of reflected ions. In the limit of M sub A of less than 10, the magnetic overshoot and the fraction of reflected ions increase with increasing M sub A, while at higher Mach numbers the fraction of reflected ions peaks at about 40 percent and the magnetic field overshoot increases at a much slower rate. Electron inertial effects are also considered.

Quest, K. B.↗

NASA NTRS1986

Nonadiabatic electron heating at high-Mach-number perpendicular shocks

Fully kinetic simulations of high-Mach-number (HMN) perpendicular collisionless shocks are described. It is shown that electron acceleration in the cross-shock electron field can produce downstream electron temperature significantly higher than those expected for adiabatic compression. The momentum space for test electrons at Mach 6 is illustrated.

Tokar, R. L.↗

NASA NTRS1986

Simulations of collisionless shocks

A problem of critical importance to space and astrophysics is the existence and properties of high-Mach-number shocks. In this letter the results of simulations of perpendicular shocks with Alfven Mach number 22 are presented. It is shown that the shock structure is a sensitive function of resistivity, becoming turbulent when the resistivity is too low. The problem of electron heating, and the extension of the results to higher Mach numbers are discussed.

Quest, K. B.↗

NASA NTRS1985

Simulations of high-Mach-number collisionless perpendicular shocks in astrophysical plasmas

A problem of critical importance to space physics and astrophysics is the existence and properties of high-Mach-number shocks. Preliminary results of a simulation of a perpendicular shock with Alfven Mach number 22 are reported. It is shown that for sufficiently small electron resistivity the dissipation for this shock is provided by a periodic rather than time-stationary reflection of ions. The problem of electron heating and the extension to higher Mach numbers are discussed.

Quest, K. B.↗

NASA NTRS1985

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗

NASA NTRS1985

Collisional tearing in a field-reversed sheet pinch assuming nonparallel propagation

Linear tearing in a collisional reversed-field sheet pinch is examined assuming that the wave vector k is not parallel to the equilibrium magnetic field. Equilibrium and magnetic geometry are defined, and a set of perturbed moment equations is derived assuming quasi-parallel propagation. It is shown that the usual expression for collisional growth is recovered, assuming that k sub y = 0. It is shown that the y component of momentum balance requires the generation of nonzero dJ sub x well away from the null, and an interial coupling when z not equal to 0. The effects of k sub y not equal to 0 on the growth rate are discussed.

Quest, K. B.↗

NASA NTRS1985

Plasma heating at collisionless shocks due to the kinetic cross-field streaming instability

Heating at collisionless shocks due to the kinetic cross-field streaming instability, which is the finite beta (ratio of plasma to magnetic pressure) extension of the modified two stream instability, is studied. Heating rates are derived from quasi-linear theory and compared with results from particle simulations to show that electron heating relative to ion heating and heating parallel to the magnetic field relative to perpendicular heating for both the electrons and ions increase with beta. The simulations suggest that electron dynamics determine the saturation level of the instability, which is manifested by the formation of a flattop electron distribution parallel to the magnetic field. As a result, both the saturation levels of the fluctuations and the heating rates decrease sharply with beta. Applications of these results to plasma heating in simulations of shocks and the earth's bow shock are described.

Winske, D.↗