Hybrid simulations of quasineutral phenomena in magnetized plasma

A new class of numerical algorithms for computer simulation of low frequency electromagnetic and electrostatic phenomena in magnetized plasma is presented. Maxwell's equations are solved in the limits of quasineutrality and negligible transverse displacement current (Darwin's model). Electrons are modeled as a fluid with polarization effects ignored. Ions are described as particles. A novel feature of these algorithms is the use of the electron fluid equation of motion to determine the electric field, which renders these numerical schemes remarkably simple and direct. The simulation plasma is either periodic or bounded by particle reflecting conducting walls. Both fully nonlinear codes with spatial grids and linearized gridless codes have been implemented