Multiple equilibria of divertor plasmas

Abstract

A one-dimensional, two-fluid transport model with a temperature-dependent neutral recycling coefficient is shown to give rise to multiple equilibria of divertor plasmas (bifurcation). Numerical techniques for obtaining these multiple equilibria and for examining their stability are presented. Although these numerical techniques have been well known to the scientific community, this is the first time they have been applied to divertor plasma modeling to show the existence of multiple equilibria as well as the stability of these solutions. Numerical and approximate analytical solutions of the present one-dimensional transport model both indicate that there exists three steady-state solutions corresponding to (1) a high-temperature, low-density equilibrium, (2) a low-temperature, high-density equilibrium, and (3) an intermediate-temperature equilibrium. While both the low-temperature and the high-temperature equilibria are stable, with respect to small perturbations in the plasma conditions, the intermediate-temperature equilibrium is physically unstable, i.e., any small perturbation about this equilibrium will cause a transition toward either the high-temperature or low-temperature equilibrium.