Ambipolarity condition and possibility of multivalued steady-state solutions to transport equations in stellarators

The possibility of steady-state multivalued solutions to transport equations in stellarators is considered. It is shown that the ambipolarity condition is necessary but not sufficient to find the ambipolar electric field, because the functions entering into it (the plasma density and temperature, as well as their spatial derivatives) depend on the ambipolar field. To do this correctly, it is necessary to solve the full set of time-independent transport equations (including diffusion and heat conduction equations). The possible existence of multivalued solutions to this set of equations is analyzed numerically. It is shown that, under certain conditions that depend on the form and magnitude of particle and heat sources, such solutions can exist. Their form is determined by the initial value of the ambipolar field, the source magnitudes, and the boundary conditions. Discontinuous solutions in which the radial profile of the ambipolar field undergoes jumps are found. In this case, however, the particle and energy fluxes remain continuous, because the discontinuities of the electric field are balanced by the discontinuities of the density and temperature gradients.