No abstract

The first results of particle-in-cell simulations of the electrostatic sheath and magnetic pre-sheath of thermionically emitting planar tungsten surfaces in fusion plasmas are presented. Plasma conditions during edge localized modes (ELMs) and during inter-ELM periods have been considered for various inclinations of the magnetic field and for selected surface temperatures. All runs have been performed under two assumptions for the sheath potential drop; fixed or floating. The primary focus lies on the evaluation of the escaping thermionic current and the quantification of the suppression due to the combined effects of space-charge and Larmor gyration. When applicable, the results are compared with the predictions of analytical models. The heat balance in the presence of thermionic emission as well as the contribution of the escaping thermionic current to surface cooling are also investigated. Regimes are identified where emission needs to be considered in the energy budget.

Plasma-facing components receive enormous heat fluxes under steady state and especially during transient conditions that can even lead to tungsten (W) melting. Under these conditions, the unimpeded thermionic current density emitted from the W surfaces can exceed the incident plasma current densities by several orders of magnitude triggering a replacement current which drives melt layer motion via the J × B force. However, in tokamaks, the thermionic current is suppressed by space-charge effects and prompt re-deposition due to gyrorotation. We present comprehensive results of particle-in-cell modelling using the 2D3V code SPICE2 for the thermionic emissive sheath of tungsten. Simulations have been performed for various surface temperatures and selected inclinations of the magnetic field corresponding to the leading edge and sloped exposures. The surface temperature dependence of the escaping thermionic current and its limiting value are determined for various plasma parameters; For the leading edge geometry, the results agree remarkably well with the Takamura analytical model. For the sloped geometry, the limiting value is observed to be proportional to the thermal electron current and a simple analytical expression is proposed that accurately reproduces the numerical results.