Gyrokinetic linear instabilities and quasilinear fluxes for variations of ITER tokamak baseline parameters
Creators
- 1. DIFFER, Eindhoven, The Netherlands
- 2. DeepMind, London, UK
- 3. Max-Planck Institut für Plasmaphysik, Garching, Germany
Description
Linear instability and quasilinear fluxes calculated with the GENE plasma microturbulence code. The input parameters correspond to variations of ITER baseline scenario parameters calculated by integrated modelling using the QuaLiKiz transport model, as described in P. Mantica et al (2019) Plasma Physics and Controlled Fusion 62 014021.
The quasilinear fluxes were calculated with a bespoke saturation rule calibrated to dedicated GENE nonlinear simulations carried out in the same ITER regime. The quasilinear flux dataset was fit with a neural network (NN) regression model, which was then used for ITER baseline integrated modelling and performance projections. Alongside the linear stability dataset, the two separate quasilinear datasets correspond to an unfiltered dataset, and a filtered and data-augmented dataset used for the NN regression. For full details, please see reference [J. Citrin et al (2023) submitted to Physics of Plasmas].
The dimensionless input and output variables correspond to the IMAS gyrokinetic IDS standards. The major radius was taken as the reference length. A key and further details are found below.
| rhoN | Normalized toroidal flux coordinate |
| ky | Binormal wavenumber, normalized to the reference (ion scale) gyroradius |
| omt_DT | Normalized logarithmic main ion temperature gradient (\(R/L_{Ti}\)) |
| omt_el | Normalized logarithmic electron temperature gradient (\(R/L_{Te}\)) |
| omn_el | Normalized logarithmic electron density gradient (\(R/L_{ne}\)) |
| s | Magnetic shear |
| q | Safety factor (q-profile) |
| gamma | Instability growth rate (gyroBohm normalisation with IMAS standard) |
| omega | Instability frequency (gyroBohm normalisation). Positive frequencies correspond to the ion diamagnetic direction |
| kperp2 | Square of perpendicular wavenumber weighted over poloidal mode structure \(\langle{k_\perp^2}\rangle\) |
| Q_DT | ky-dependent ion heat flux, normalized by the square of the electrostatic potential |
| Q_el | ky-dependent electron heat flux, normalized by the square of the electrostatic potential |
| G_el | ky-dependent electron particle flux, normalized by the square of the electrostatic potential |
| QDT | Quasilinear ion heat flux, following summation of modes and a saturation rule |
| Qe | Quasilinear electron heat flux, following summation of modes and a saturation rule |
| Ge | Quasilinear electron particle flux, following summation of modes and a saturation rule |
| QDT_ITG | Quasilinear ion heat flux, when considering ITG modes only. GyroBohm normalized with IMAS convention |
| Qe_ITG | Quasilinear electron heat flux, when considering ITG modes only. GyroBohm normalized with IMAS convention |
| Ge_ITG | Quasilinear electron particle flux, when considering ITG modes only. GyroBohm normalized with IMAS convention |
| QDT_TEM | Quasilinear ion heat flux, when considering TEM modes only. GyroBohm normalized with IMAS convention |
| Qe_TEM | Quasilinear electron heat flux, when considering TEM modes only. GyroBohm normalized with IMAS convention |
| Ge_TEM | Quasilinear electron particle flux, when considering TEM modes only. GyroBohm normalized with IMAS convention |