Prediction of divertor heat flux width for ITER Pre-Fusion Power Operation using BOUT++ transport code

Prediction of divertor heat flux width is performed for the first and the second Pre-Fusion Power Operation (PFPO) phases specified in the new ITER Research Plan using BOUT++ transport code [Li N.M. et al 2018 Comput. Phys. Commun. 228 69–82]. The initial plasma profiles inside the separatrix are taken from CORSICA scenario studies. Transport coefficients in transport code are calculated by inverting the plasma profiles inside the separatrix and are assumed to be constants in the scrape-off-layer (SOL). An anomalous thermal diffusivity scan is performed with E×B and magnetic drifts. The results in two scenarios identify two distinct regimes: a drift dominant regime when diffusivity is smaller than the respective critical diffusivity χc and a turbulence dominant regime when diffusivity is larger than it. The Goldston heuristic drift model and the ITPA multi-machine experimental scaling yield a lower limit of the width λq. From transport simulations, we obtain the critical diffusivity χc = 0.5 m2⁄ s in 5MA/1.77T PFPO-1 scenario and χc = 0.3 m2⁄ s in 7.5MA/2.65T PFPO-2 scenario. Separatrix temperature and collisionality also have a significant impact on the heat flux width in the drift dominant regime. The investigation clearly yields a scaling for critical thermal diffusivity χc ∝ A½ ⁄ ((Z(1+Z)½ Bp 2)) using ITER scenarios with fixed safety factor q95, major radius R, aspect ratio R/a, and the separatrix temperature T, as well as established the connection with CFETR and C-Mod discharges. This scaling implies that for a given tokamak device with q95, R, R/a, and T fixed, a reduction of poloidal magnetic field by a factor of 3 leads to a 9 times higher critical value of thermal diffusivity χc, possibly yielding a transition from turbulence to drift dominant regime.