Application of Gaussian process regression to plasma turbulent transport model validation via integrated modelling

• Published in: arXiv.org
• Main Authors: Ho, Aaron, Citrin, Jonathan, Auriemma, Fulvio, Bourdelle, Clarisse, Casson, Francis J, Hyun-Tae, Kim, Manas, Pierre, Szepesi, Gabor, Weisen, Henri, Contributors, JET
• Format: Paper Journal Article
• Language: English
• Published: Ithaca Cornell University Library, arXiv.org 03.05.2021
• Subjects: Boundary conditions; Deuterium; Discharge; Fluxes; Gaussian process; Joint European Torus; Modelling; Momentum; Physics - Plasma Physics; Simulation; Tokamak devices; Tritium; Turbulence
• ISSN: 2331-8422
• DOI: 10.48550/arxiv.2105.01177

This paper outlines an approach towards improved rigour in tokamak turbulence transport model validation within integrated modelling. Gaussian process regression (GPR) techniques were applied for profile fitting during the preparation of integrated modelling simulations allowing for rigourous sensitivity tests of prescribed initial and boundary conditions as both fit and derivative uncertainties are provided. This was demonstrated by a JETTO integrated modelling simulation of the JET ITER-like-wall H-mode baseline discharge #92436 with the QuaLiKiz quasilinear turbulent transport model, which is the subject of extrapolation towards a deuterium-tritium plasma. The simulation simultaneously evaluates the time evolution of heat, particle, and momentum fluxes over \(\sim10\) confinement times, with a simulation boundary condition at \(\rho_{tor} = 0.85\). Routine inclusion of momentum transport prediction in multi-channel flux-driven transport modelling is not standard and is facilitated here by recent developments within the QuaLiKiz model. Excellent agreement was achieved between the fitted and simulated profiles for \(n_e\), \(T_e\), \(T_i\), and \(\Omega_{tor}\) within \(2\sigma\), but the simulation underpredicts the mid-radius \(T_i\) and overpredicts the core \(n_e\) and \(T_e\) profiles for this discharge. Despite this, it was shown that this approach is capable of deriving reasonable inputs, including derivative quantities, to tokamak models from experimental data. Furthermore, multiple figures-of-merit were defined to quantitatively assess the agreement of integrated modelling predictions to experimental data within the GPR profile fitting framework.