A computer simulation of the midlatitude plasmasphere and ionosphere

A computer model has been developed to simulate species density, temperature, and plasma flow in the ionosphere and plasmasphere. A new approach, the flux preserving method involving a Newton iteration, is used to solve a system of governing equations comprising four second-order partial differentials. Tests are performed to demonstrate that the simulation converges to a stable steady state solution, then steady state ion fluxes are analyzed. Finally, simulations of the collapse of the sunset ionosphere are presented. Comparisons with satellite and radar data show good agreement in a number of cases.