Numerical calculations relevant to the initial expansion of the polar wind
The results of a numerical simulation of the expansion of an H(+)-0(+) electron plasma into a vacuum are reported. Variations were introduced in the ionic density ratios, the scale lengths for the density gradient at the plasma-vacuum interface, and the initial electron-ion temperature ratio. Conditions of a gravitationally bound O(+) plasma were also considered. It was found that if H(+) is a major ion, the H(+) density profile in the expansion region is always concave, with a phase-space divided into three regions, i.e., an undisturbed plasma, a rarefaction region, and an expansion region in which ions are accelerated. When the H(+) is a minor ion, the O(+) ions accelerate the H(+) ions, the H(+) phase-space is divided into five regions, and the increased electron temperatures produce enhanced electrostatic potentials and H(+) drifts in the expansion region.