Particle in cell simulation of instabilities in space and astrophysical plasmas.

摘要

Several plasma instabilities relevant to space physics are investigated using the parallel PIC plasma simulation code P3arsec. This thesis addresses electrostatic micro-instabilities relevant to ion ring distributions, proceeds to electromagnetic micro-instabilities pertinent to streaming plasmas, and then to the stability of a plasma held in the field of a current rod. The physical relevance of each of these instabilities is discussed, a phenomenological description is given, and analytic and simulation results are presented and compared.; Instability of a magnetized plasma with a portion of the ions in a velocity ring distribution around the magnetic field is investigated using simulation and analytic theory. The physics of this distribution is relevant to solar flares, x-ray emission by comets, and pulsars. Physical parameters, including the mass ratio, are near those of a solar flare in the simulation. The simulation and analytic results show agreement in the linear regime. In the nonlinear stage the simulation shows highly accelerated electrons in agreement with the observed spectrum of x-rays emitted by solar flares.; A mildly relativistic streaming electron positron plasma with no ambient magnetic field is known to be unstable to electrostatic (two-stream/beam instability) and purely electromagnetic (Weibel) modes. This instability is relevant to highly energetic interstellar phenomena, including pulsars, supernova remnants, and the early universe. It is also important for experiments in which relativistic beams penetrate a background plasma, as in fast ignitor scenarios. Cold analytic theory is presented and compared to simulations. There is good agreement in the regime where cold theory applies. The simulation and theory shows that to properly characterize the instability, directions parallel and perpendicular to propagation of the beams must be considered. A residual magnetic field is observed which may be of astro-physical significance.; The stability of a plasma in the magnetic field of a current rod is investigated for various temperature and density profiles. Such a plasma obeys similar physics to a plasma in a dipole magnetic field, while the current rod is much easier to analyze theoretically and realize in simulations. The stability properties of a plasma confined in a dipole field are important for understanding a variety of space phenomena and the Levitated Dipole eXperiment (LDX). Simple energy principle calculations and simulations with a variety of temperature and density profiles show that the plasma is stable to interchange for pressure profiles ∝ r−10/3. The simulations also show that the density profile will be stationary as long as density ∝ r −2 even though the temperature profile may not be stable.