Transport of gyration-dominated space plasmas of thermal origin. II: Numerical solution (Q1316024)

[For part I, see \textit{T. I. Gombosi} and \textit{C. E. Rasmussen}, J. Geophys. Res. 96, 7759 (1991)]. A numerical solution has been found for the hyperbolic initial-value problem, which follows from the 16-moment set of equations describing the transport of a plasma in a strong magnetic field. In addition to the time-dependent evolution of mass and momentum, these equations describe the variations in anisotropic temperatures and heat fluxes for both ion and electron species. The numerical solution employs a high-order Godunov method to achieve third-order accuracy in space and second-order accuracy in time. The method includes an approximate Riemann solver.