Numerical calculation of the electron distribution function for EBT electron annuli

In the ELMO Bumpy Torus (EBT), high beta relativistic electron annuli in each mirror sector provide the average minimum-B configuration necessary for stability of the toroidal plasma. A nonrelativistic bounce-averaged kinetic equation is solved for the electron distribution function, which is assumed to be isotropic in velocity space. The equation is solved in time and two dimensions and includes both velocity and radial derivatives. The Fokker-Planck collision operator includes Coulomb collisions and first- and second-harmonic microwave heating terms. Ionization of cold neutrals (assuming a cylindrical neutral transport model) provides a source of cold electrons