Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Abstract Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">e</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> ) and ion temperature ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">i</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">e</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> by 46%, with much larger variance than MMM, and underpredicts <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">i</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> by 25%. As the ratio of kinetic to magnetic field pressure ( β ) is increased across NSTX discharges, TGLF predicts lower <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">e</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> and significant flattening of the <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">i</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">e</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:msub> <mml:mi>T</mml:mi> <mml:mrow> <mml:mi mathvariant="normal">i</mml:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Altogether, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.