Effect of Coulomb collisions on contrastreaming plasmas
Effect of Coulomb collisions on instability of contrastreaming plasmas - numerical analysis
SEARCH · Engineering Papers
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.
Effect of Coulomb collisions on instability of contrastreaming plasmas - numerical analysis
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.
Data obtained by OGO 5 are used to confirm IMP 6 observations of an inverse dependence of the helium-to-hydrogen temperature ratio in the solar wind on the ratio of solar-wind expansion time to the Coulomb-collision equipartition time. The analysis is then extended to determine the relation of the difference between the hydrogen and helium bulk velocities (the differential flow vector) with the ratio between the solar-wind expansion time and the time required for Coulomb collisions to slow down a beam of ions passing through a plasma. It is found that the magnitude of the differential flow vector varies inversely with the time ratio when the latter is small and approaches zero when it is large. These results are shown to suggest a model of continuous preferential heating and acceleration of helium (or cooling and deceleration of hydrogen), which is cancelled or limited by Coulomb collisions by the time the plasma has reached 1 AU. Since the average dependence of the differential flow vector on the time ratio cannot explain all the systematic variations of the vector observed in corotating high-velocity streams, it is concluded that additional helium acceleration probably occurs on the leading edge of such streams.
Instability of contrastreaming plasmas investigated by taking into account Coulomb collisions via Fokker-Planck coefficients in Boltzmann equation
Expressions are given for the drag due to Coulomb collisions on a charged dust grain moving through a plasma. The commonly used Chandrasekhar (1943) result does not include large scattering angles or plasma collective effects. An equation given by Morfill et al. (1980) which does include large scattering angles, is limited to one-dimensional plasma particle motion and can give at least an order of magnitude too much drag. This paper also makes use of an analogy between the drag problem and problems in electrostatics. This analogy permits generalization to any isotropic plasma of an observation made by Chandrasekhar for a Maxwellian, namely, that the drag is independent of the presence or absence of plasma particles moving faster than the grain. Finally, the contribution of plasma collective effects to the drag is studied with the inclusion of large scattering angles.
NASA’s Endurance rocket mission resolved a subtle energy shift in photoelectrons escaping from Earth’s polar caps, from which the first successful measurement of Earth’s ambipolar electrostatic field was made. Preflight models predicted that coulomb collisions with thermal electrons could add an extra energy shift of up to 1 eV. However, no such additional shift was reported in the initial analysis. Here, we identify a previously unreported effect related to the spacecraft wake, which introduced a systematic +0.25eV energy shift that countered and masked the shift from coulomb collisions. After applying an in-flight calibration for this wake effect, we find a total photopeak energy shift consistent with the combined influence of the ambipolar field and the subtler -0.37eV coulomb shift, predicted by the fully kinetic STET model. Further analysis of these recalibrated measurements yields an ambipolar field estimate in close agreement with the initial analysis.
Here, we introduce a new approach for coupling cumulative and single-scattering distributions, extending Monte Carlo methods for Coulomb collisions in weakly coupled plasmas into the moderately coupled regime. Using time steps typical of cumulative-scattering methods, we show that this approach not only recovers leading-order corrections proportional to the inverse Coulomb logarithm but also solves the full Boltzmann collision integral. A key feature of the method is its simplicity: it can be integrated into cumulative-scattering methods for weakly coupled plasmas with only ~ 10 additional lines of code.
A major improvement of the present investigation over previous studies of the subject is related to the use of helium temperatures obtained from helium ion measurements uncontaminated by the high-velocity tail of the proton distribution. More observations, covering a large parameter range, were employed, and the effects of interspecies drift were taken into account. It is shown in a more definite way than has been done previously, that Coulomb collisions provide the most important mechanism bringing about equilibrium between helium and protons in the solar wind. Other mechanisms may play some part in restricted regions, but Coulomb collisions are dominant on the macroscale.
Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. In conclusion, the efficacy of the model is illustrated with various test problems.
The evolution of nonthermal electrons accelerated in magnetic loops is determined by solving the kinetic equation, including magnetic field convergence and Coulomb collisions in order to determine the effects of these interactions on the induced cyclotron maser and plasma wave growth. It is found that the growth rates are larger and the possibility of cyclotron maser action is stronger for smaller loop column density, for larger magnetic field convergence, for a more isotropic injected electron pitch angle distribution, and for more impulsive acceleration. For modest values of the column density in the coronal portion of a flaring loop, the growth rates of instabilities are significantly reduced, and the reduction is much larger for the cyclotron modes than for the plasma wave modes. The rapid decrease in the growth rates with increasing loop column density suggests that, in flare loops when such phenomena occur, the densities are lower than commonly accepted.
The additional energy requirements of the topside ionosphere during a magnetic storm are less than one quarter of the ring current energy. This energy is supplied largely by Coulomb collisions of ring current protons of energy less than about 20 keV with background thermal electrons which conduct the heat to the ionosphere. Past criticisms are discussed of this mechanism for the supply of energy to the SAR-arc and neighboring regions of the ionosphere.
Standard numerical integrators work well for many-body Coulomb repulsion problems when the timestep is much shorter than the timescale of relative position changes. However, for ‘hard’ collisions in which two particles have a near miss and exchange a lot of momentum within one timestep, they understandably struggle. This note proposes using the exact solution of Keplerian two-body orbits (usually hyperbolic) to calculate the momentum exchange with other particles: either a selection of the ‘closest’ ones or all of them.
A collisionless, time-dependent, kinetic plasma model is applied to the problem of baseline plasmasphere refilling of an initially depleted flux tube, without regard for the effects of wave-particle interactions. Refilling calculations for various flux tubes and for different ionospheric plasma fluxes and temperatures are performed. In each case considered, the same set of events occurs. Initially, two polar wind outflows develop from each hemisphere and set up counterstreaming beams. With time the vacant phase space region between these beams fills, primarily because of collision-induced particle diffusion but also because of lowering ambipolar potential drops from the increasing density in the plasmasphere. In contrast to all previous hydrodynamic approaches, no formation of shocks was found. The plasma first evolves an isotropic, nearly Maxwellian velocity distribution in a region that starts near the ionosphere and moves outward toward the equator. For reasonable topside ionospheric temperatures and fluxes, the thermal plasma all along an L shell is found to become nearly isotropic in 6 to 30 hr, consistent with the observations of Horwitz et al. (1984).
Finite range effects in Coulomb collisions for high atomic number targets at low bombardment particle energies - theoretical physics
A kinetic theory for the velocity distribution of solar wind electrons which illustrates the global and local properties of the solar wind expansion is proposed. By means of the Boltzmann equation with the Krook collision operator accounting for Coulomb collisions, it is found that Coulomb collisions determine the population and shape of the electron distribution function in both the thermal and suprathermal energy regimes. For suprathermal electrons, the cumulative effects of Coulomb interactions are shown to take place on the scale of the heliosphere itself, whereas the Coulomb interactions of thermal electrons occur on a local scale near the point of observation (1 AU). The bifurcation of the electron distribution between thermal and suprathermal electrons is localized to the deep solar corona (1 to 10 solar radii).
Coulomb Born approximation derived for ionization of hydrogen atom by electron collision when atom is in any given initial state
The decay of the major ion species which constitute the ring current is studied by solving the time evolution of their distribution functions during the recovery phase of a moderate geomagnetic storm. In this work, only equatorially mirroring particles are considered. Particles are assumed to move subject to E x B and gradient drifts. They also experience loses along their drift paths. Two loss mechanisms are considered: charge exchange with neutral hydrogen atoms and Coulomb collisions with thermal plasma in the plasmasphere. Thermal plasma densities are calculated with a plasmaspheric model employing a time-dependent convection electric field model. The drift-loss model successfully reproduces a number of important and observable features in the distribution function. Charge exchange is found to be the major loss mechanism for the ring current ions; however the important effects of Coulomb collisions on both the ring current and thermal populations are also presented. The model predicts the formation of a low-energy (less than 500 eV) ion population as a result of energy degradation caused by Coulomb collision of the ring current ions with the plasmaspheric electrons; this population may be one source of the low-energy ions observed during active and quiet periods in the inner magnetosphere. The energy transferred to plasmaspheric electrons through Coulomb collisions with ring current ions is believed to be the energy source for the electron temperature enhancement and the associated 6300 A (stable auroral red (SAR) arc) emission in the subauroral region. The calculated energy deposition rate is sufficient to produce a subauroral electron temperature enhancement and SAR arc emissions that are consistent with observations of these quantities during moderate magnetic activity levels.
We investigate the thermodynamic structure of interplanetary coronal mass ejections (ICMEs) using combined surveys of the ejecta between 0.3 and 20 AU. ICMEs are shown to have a moderate expansion in the solar wind compared with theoretical predictions. The expansion seems to be governed by a polytrope with gamma approx. 1.3 in this distance range. We find that Coulomb collisions are important contributors to the ion-ion equilibration process in the ICME plasma. The alpha-proton differential speed quickly drops to below 10 km/s due to strong Coulomb collisions. However, the two species of particles are far from thermal equilibrium with a temperature ratio T(sub alpha/T(sub p) = 4-6, suggestive of a preferential heating of alpha particles. The plasma heating rate as a function of heliocentric &stance required for the temperature profile is deduced by taking into account the expansion and energy transfer between protons and alphas via Coulomb collisions. The turbulence dissipation rate is also inferred from the inertial range power spectrum of magnetic fluctuations within ICMEs. Comparison of the turbulence dissipation rate with the required heating rate shows that turbulence dissipation seems sufficient to explain the ICME heating. Sources powering the turbulence are also investigated by examining the instabilities induced by temperature anisotropies and energy deposition by pickup ions.