Scientific Machine Learning Based Reduced-Order Models for Plasma Turbulence Simulations

Abstract

This paper investigates non-intrusive Scientific Machine Learning (SciML) Reduced-Order Models (ROMs) for plasma turbulence simulations. In particular, we focus on Operator Inference (OpInf) to build low-cost physics-based ROMs from data for such simulations. As a representative example, we consider the (classical) Hasegawa-Wakatani (HW) equations used for modeling two-dimensional electrostatic drift-wave turbulence. For a comprehensive perspective of the potential of OpInf to construct predictive ROMs, we consider three setups for the HW equations by varying a key parameter, namely the adiabaticity coefficient. These setups lead to the formation of complex and nonlinear dynamics, which makes the construction of predictive ROMs of any kind challenging. We generate the training datasets by performing direct numerical simulations of the HW equations and recording the computed state data and outputs the over a time horizon of $100$ time units in the turbulent phase. We then use these datasets to construct OpInf ROMs for predictions over $400$ additional time units, that is, $400\%$ more than the training horizon. Our results show that the OpInf ROMs capture important statistical features of the turbulent dynamics and generalize beyond the training time horizon while reducing the computational effort of the high-fidelity simulation by up to five orders of magnitude. In the broader context of fusion research, this shows that non-intrusive SciML ROMs have the potential to drastically accelerate numerical studies, which can ultimately enable tasks such as the design of optimized fusion devices.

14 pages in double column format, 9 figures, 8 tables

References in corpus (14)

• Array Programming with NumPy
• On closures for reduced order models $-$ A spectrum of first-principle to machine-learned avenues
• Bifurcation in electrostatic resistive drift wave turbulence
• Operator inference for non-intrusive model reduction with quadratic manifolds
• Characterizing Magnetized Plasmas with Dynamic Mode Decomposition
• Neural Galerkin Schemes with Active Learning for High-Dimensional Evolution Equations
• Multi-scale data reconstruction of turbulent rotating flows with Gappy POD, Extended POD and Generative Adversarial Networks
• Predicting Solar Wind Streams from the Inner-Heliosphere to Earth via Shifted Operator Inference
• Turbulence model reduction by deep learning
• Domain decomposition for data-driven reduced modeling of large-scale systems
• Fourier Neural Operator for Plasma Modelling
• Synthesizing impurity clustering in the edge plasma of tokamaks using neural networks
• Physics-Preserving AI-Accelerated Simulations of Plasma Turbulence
• Non-Intrusive Reduced Models based on Operator Inference for Chaotic Systems