On the discontinuity in kinetic solutions of the collisionless polar wind
The use of a Maxwellian ion velocity distribution in kinetic solutions of the collisionless polar wind regime results in a velocity distribution and associated moments that display sharp discontinuities at the boundary, which quantitatively affects the solution at high altitudes. This paper investigates the causes of this discontinuity, using both analytical and numerical approaches; it is shown that the choice of appropriate boundary conditions will eliminate the discontinuity.