Machine-learning Closure for Vlasov-Poisson Dynamics in Fourier-Hermite Space

Abstract

Accurate reduced models of turbulence are desirable to facilitate the optimization of magnetic-confinement fusion reactor designs. As a first step toward higher-dimensional turbulence applications, we use reservoir computing, a machine-learning (ML) architecture, to develop a closure model for a limiting case of electrostatic gyrokinetics. We implement a pseudo-spectral Eulerian code to solve the one-dimensional Vlasov-Poisson system on a basis of Fourier modes in configuration space and Hermite polynomials in velocity space. When cast onto the Hermite basis, the Vlasov equation becomes an infinitely coupled hierarchy of fluid moments, presenting a closure problem. We exploit the locality of interactions in the Hermite representation to introduce an ML closure model of the small-scale dynamics in velocity space. In the linear limit, when the kinetic Fourier-Hermite solver is augmented with the reservoir closure, the closure permits a reduction of the velocity resolution, with a relative error within two percent for the Hermite moment where the reservoir closes the hierarchy. In the strongly-nonlinear regime, the ML closure model more accurately resolves the low-order Fourier and Hermite spectra when compared to a na\"{i}ve closure by truncation and reduces the required velocity resolution by a factor of sixteen.