Simulations of the stationary Q=10 and the exit phase from the flat-top of an ITER 15MA baseline scenario Predictive JINTRAC simulations with consistent treatment of D and T in the whole plasma

Frida Eriksson,Emmi Tholerus,Gerard Corrigan,Yuriy Baranov,Xavier Bonnin,Daniela Farina,Lorenzo Figini,Luca Garzotti,Sunhee Kim,Florian Koechl,Alberto Loarte,Elina Militello-Asp,cameron Olde,Vassili Parail,Simon D Pinches,Alexei R Polevoi,Pär Strand

DOI: https://doi.org/10.1088/1741-4326/ad7c63

IF: 3.3

2024-09-20

Nuclear Fusion

Abstract:To design a robust termination scenario for a burning ITER plasma is a challenge that requires extensive core plasma and divertor modelling. The presented work consists of coupled core/edge/SOL/divertor simulations, performed with the JINTRAC code, studying the Q=10 flat-top phase and exit phase of the ITER 15MA/5.3T DT scenario. The modelling utilizes the recently implemented option to treat deuterium and tritium separately in the SOL/divertor enabling a consistent treatment of deuterium and tritium in the whole plasma volume which is a unique capability of JINTRAC. In addition, these are the first JINTRAC simulations of this scenario to use a first-principle transport model, self-consistently model the ECRH power deposition and to include tungsten, keeping track of tungsten sputtering and accumulation. The flat-top simulations demonstrate the possibility of sustaining a steady state fusion Q of 10 using pure deuterium gas puff together with DT mixed pellets which is an option to make a more effective use of tritium. Simulations of the exit phase are set up sequentially, with each phase providing initial conditions to the next, starting with a density decay at full current and auxiliary power and demonstrate the possibility to reduce the density robustly within a few seconds. Following the density decay, a subsequent auxiliary power ramp down in H-mode is performed with a late H-L transition at low auxiliary power which may provide an option for optimization of the plasma termination. The final ramp-down phase consists of a current ramp-down in L-mode to 3.75MA.

physics, fluids & plasmas