2D divertor heat flux distribution using a 3D heat conduction solver in National Spherical Torus Experiment

Abstract

The divertor heat flux footprint in tokamaks is often observed to be non-axisymmetric due to intrinsic error fields, applied 3D magnetic fields or during transients such as edge localized modes. Typically, only 1D radial heat flux profiles are analyzed; however, analysis of the full 2D divertor measurements provides opportunities to study the asymmetric nature of the deposited heat flux. To accomplish this an improved 3D Fourier analysis method has been successfully applied in a heat conduction solver (TACO) to determine the 2D heat flux distribution at the lower divertor surface in the National Spherical Torus Experiment (NSTX) tokamak. This advance enables study of helical heat deposition onto the divertor. In order to account for heat transmission through poorly adhered surface layers on the divertor plate, a heat transmission coefficient, defined as the surface layer thermal conductivity divided by the thickness of the layer, was introduced to the solution of heat conduction equation. This coefficient is denoted as α and a range of values were tested in the model to ensure a reliable heat flux calculation until a specific value of α led to the constant total deposited energy in the numerical solution after the end of discharge. A comparison between 1D heat flux profiles from TACO and from a 2D heat flux calculation code, THEODOR, shows good agreement. Advantages of 2D heat flux distribution over the conventional 1D heat flux profile are also discussed, and examples of 2D data analysis in the study of striated heat deposition pattern as well as the toroidal degree of asymmetry of peak heat flux and heat flux width are demonstrated.