L.L. LoDestro

Bringing global gyrokinetic turbulence simulations to the transport timescale using a multiscale approach

The vast separation dividing the characteristic times of energy confinement and turbulence in the core of toroidal plasmas makes first-principles prediction on long timescales extremely challenging. Here we report the demonstration of a multiple-timescale method that enables coupling global gyrokinetic simulations with a transport solver to calculate the evolution of the self-consistent temperature profile. This method, which exhibits resiliency to the intrinsic fluctuations arising in turbulence simulations, holds potential for integrating nonlocal gyrokinetic turbulence simulations into predictive, whole-device models.

Investigation of a Multiple-Timescale Turbulence-Transport Coupling Method in the Presence of Random Fluctuations

One route to improved predictive modeling of magnetically confined fusion reactors is to couple transport solvers with direct numerical simulations (DNS) of turbulence, rather than with surrogate models. An additional challenge presented by coupling directly with DNS is the inherent fluctuations in the turbulence, which limit the convergence achievable in the transport solver. In this article, we investigate the performance of one numerical coupling method in the presence of turbulent fluctuations. To test a particular numerical coupling method for the transport solver, we use an autoregressive-moving-average model to generate stochastic fluctuations efficiently with statistical properties resembling those of a gyrokinetic simulation. These fluctuations are then added to a simple, solvable problem, and we examine the behavior of the coupling method. We find that monitoring the residual as a proxy for the error can be misleading. From a pragmatic point of view, this study aids us in the full problem of transport coupled to DNS by predicting the amount of averaging required to reduce the fluctuation error and obtain a specific level of accuracy.