A physically comprehensive and theoretically based transport model tuned to three-dimensional (3-D) ballooning mode gyrokinetic instabilities and gyrofluid nonlinear turbulence simulations is formulated with global and local magnetic shear stabilization and E×B rotational shear stabilization. Taking no fit coefficients from experiment, the model is tested against a large transport profile database with good agreement. This model is capable of describing enhanced core confinement transport barriers in negative central shear discharges based on rotational shear stabilization. The model is used to make ignition projections from relative gyroradius scaling discharges.

A new theory-based transport model with comprehensive physics (trapping, general toroidal geometry, fully electromagnetic, electron-ion collisions, impurity ions) has been developed. The core of the model is the new trapped-gyro-Landau-fluid (TGLF) equations, which provide a fast and accurate approximation to the linear eigenmodes for gyrokinetic drift-wave instabilities (trapped ion and electron modes, ion and electron temperature gradient modes, and kinetic ballooning modes). The new TGLF transport model is more accurate, and has an extended range of validity, compared to its predecessor GLF23. The TGLF model unifies trapped and passing particles in a single set of gyro-Landau-fluid equations. A model for the averaging of the Landau resonance by the trapped particles makes the equations work seamlessly over the whole drift-wave wave-number range from trapped ion modes to electron temperature gradient modes. A fast eigenmode solution method enables unrestricted magnetic geometry. The transport model uses the TGLF eigenmodes to compute quasilinear fluxes of energy and particles. A model for the saturated intensity of the turbulence completes the flux calculation. The intensity model is constructed to fit a large set of nonlinear gyrokinetic turbulence simulations with kinetic electrons. The TGLF model is valid in new physical regimes that GLF23 was not. These include the low aspect ratio spherical torus, which has both a high trapped fraction and strong shaping of magnetic flux surfaces. The TGLF model is also valid close to the magnetic separatrix so the transport physics of the H-mode pedestal region can be explored.

A new system of gyro-Landau fluid (GLF) equations for tokamak plasmas is presented. The new equations include both trapped particles, which can average the Landau resonance, and passing particles which do have a Landau resonance. The trap GLF (TGLF) model is unrestricted in trapped fraction or perpendicular wave number of the electrostatic perturbation. The linearly unstable eigenmodes of the TGLF equations include low-frequency trapped ion modes all the way up to high-frequency electron temperature gradient driftwaves. Extensive benchmarking of the linear TGLF eigenmodes with a large database of gyrokinetic linear stability calculations verifies that the TGLF model is accurate over the full range of plasma parameters tested.

One modeling framework for integrated tasks (OMFIT) is a comprehensive integrated modeling framework which has been developed to enable physics codes to interact in complicated workflows, and support scientists at all stages of the modeling cycle. The OMFIT development follows a unique bottom-up approach, where the framework design and capabilities organically evolve to support progressive integration of the components that are required to accomplish physics goals of increasing complexity. OMFIT provides a workflow for easily generating full kinetic equilibrium reconstructions that are constrained by magnetic and motional Stark effect measurements, and kinetic profile information that includes fast-ion pressure modeled by a transport code. It was found that magnetic measurements can be used to quantify the amount of anomalous fast-ion diffusion that is present in DIII-D discharges, and provide an estimate that is consistent with what would be needed for transport simulations to match the measured neutron rates. OMFIT was used to streamline edge-stability analyses, and evaluate the effect of resonant magnetic perturbation (RMP) on the pedestal stability, which have been found to be consistent with the experimental observations. The development of a five-dimensional numerical fluid model for estimating the effects of the interaction between magnetohydrodynamic (MHD) and microturbulence, and its systematic verification against analytic models was also supported by the framework. OMFIT was used for optimizing an innovative high-harmonic fast wave system proposed for DIII-D. For a parallel refractive index > ∥ n 3, the conditions for strong electron-Landau damping were found to be independent of launched ∥ n and poloidal angle. OMFIT has been the platform of choice for developing a neural-network based approach to efficiently perform a non-linear multivariate regression of local transport fluxes as a function of local dimensionless parameters. Transport predictions for thousands of DIII-D discharges showed excellent agreement with the power balance calculations across the whole plasma radius and over a broad range of operating Nuclear Fusion

The 2D spectrum of the saturated electric potential from gyrokinetic turbulence simulations that include both ion and electron scales (multi-scale) in axisymmetric tokamak geometry is analyzed. The paradigm that the turbulence is saturated when the zonal (axisymmetic) ExB flow shearing rate competes with linear growth is shown to not apply to the electron scale turbulence. Instead, it is the mixing rate by the zonal ExB velocity spectrum with the turbulent distribution function that competes with linear growth. A model of this mechanism is shown to be able to capture the suppression of electron-scale turbulence by ion-scale turbulence and the threshold for the increase in electron scale turbulence when the ion-scale turbulence is reduced. The model computes the strength of the zonal flow velocity and the saturated potential spectrum from the linear growth rate spectrum. The model for the saturated electric potential spectrum is applied to a quasilinear transport model and shown to accurately reproduce the electron and ion energy fluxes of the non-linear gyrokinetic multi-scale simulations. The zonal flow mixing saturation model is also shown to reproduce the non-linear upshift in the critical temperature gradient caused by zonal flows in ion-scale gyrokinetic simulations.