A three-dimensional numerical model for current collection in a magnetized plasma is reported. The model is based on an electrostatic particle-in-cell code. The model yields self-consistent sheath structure including distributions of plasma and the electric potential around the body and the collection of electrons. The analytical theory of current collection by a body in a magnetized plasma yields an upper bound for the collected current determined by the conservation of energy and canonical angular momentum. The theory shows that the collected charged particles come from a cylindrical volume aligned with the magnetic shadow of the body; the maximum radius r0 of this volume is determined by the body size, body potential, and the ambient magnetic field. This theory does not deal with the sheath structure around the body. The condition for the actual current to approach the upper-bound value has been a matter of debate. Our simulations reveal when and why the collected current becomes equal to its upper-bound value.

Sheath size in the radial direction perpendicular to the axial ambient magnetic field is determined by the current-limiting radius r0. Our simulation yields time-average current in good agreement with its upper bound. This feature of the current collection is explained as follows: Once electrons enter the sheath, some of them are freely accelerated perpendicular to the magnetic field because they are demagnetized by the large gradients in the perpendicular electric fields. Simulations show a large perpendicular acceleration, producing perpendicular energy as large as that determined by the potential on the body, especially in the region where perpendicular electric fields (E⊥) are the strongest. An analysis shows that the demagnetization of electrons occurs above a threshold potential on the body. This threshold condition follows from the breakdown of the adiabaticity of the electron dynamics inside the shealth. ¿ American Geophysical Union 1994