Numerical simulation of torus-driven plasma transport in the Jovian magnetosphere

The Rice convection model has been modified for application to the transport of Io-generated plasma through the Jovian magnetosphere. The new code, called the RCM-J, has been used for several ideal-magnetohydrodynamic (MHD) numerical simulations to study how interchange instability causes an initially assumed torus configuration to break up. In simulations that start from a realistic torus configuration but include no energetic particles, the torus disintegrates too quickly (approximately 50 hours). By adding an impounding distribution of energetic particles to suppress the interchange instability, resonable lifetimes were obtained. For cases in which impoundment is insufficient to produce ideal-MHD stability, the torus breaks up predominantly into long fingers, unless the initial condition strongly favors some other geometrical form. If the initial torus has more mass on one side of the planet than the other, fingers form predominatly on the heavy side (which we associate with the active sector). Coriolis force bends the fingers to lag corotation. The simulation results are consistent with the idea that the fingers are formed with a longitudinal thickness that is roughly equal to the latitudinal distance over which the invariant density declines at the outer edges of the initial torus. Our calculations give an average longitudinal distance between plasma fingers of about 15 deg which corresponds to 20 to 30 minutes of rotation of the torus. We point to some Voyager and Ulysses data that are consistent with this scale of torus longitudinal irregularity.