Encoder–decoder neural network for solving the nonlinear Fokker–Planck–Landau collision operator in XGC

• Published in: Journal of plasma physics Vol. 87; no. 2
• Main Authors: Miller, M. A., Churchill, R. M., Dener, A., Chang, C. S., Munson, T., Hager, R.
• Format: Journal Article
• Language: English
• Published: United States Cambridge University Press 24.03.2021
• Subjects: 70 PLASMA PHYSICS AND FUSION TECHNOLOGY; fusion plasma; plasma simulation
• ISSN: 0022-3778; 1469-7807; 1469-7807
• DOI: 10.1017/S0022377821000155

An encoder–decoder neural network has been used to examine the possibility for acceleration of a partial integro-differential equation, the Fokker–Planck–Landau collision operator. This is part of the governing equation in the massively parallel particle-in-cell code XGC, which is used to study turbulence in fusion energy devices. The neural network emphasizes physics-inspired learning, where it is taught to respect physical conservation constraints of the collision operator by including them in the training loss, along with the $\ell _2$ loss. In particular, network architectures used for the computer vision task of semantic segmentation have been used for training. A penalization method is used to enforce the ‘soft’ constraints of the system and integrate error in the conservation properties into the loss function. During training, quantities representing the particle density, momentum and energy for all species of the system are calculated at each configuration vertex, mirroring the procedure in XGC. This simple training has produced a median relative loss, across configuration space, of the order of $10^{-4}$ , which is low enough if the error is of random nature, but not if it is of drift nature in time steps. The run time for the current Picard iterative solver of the operator is $O(n^2)$ , where $n$ is the number of plasma species. As the XGC1 code begins to attack problems including a larger number of species, the collision operator will become expensive computationally, making the neural network solver even more important, especially since its training only scales as $O(n)$ . A wide enough range of collisionality has been considered in the training data to ensure the full domain of collision physics is captured. An advanced technique to decrease the losses further will be subject of a subsequent report. Eventual work will include expansion of the network to include multiple plasma species.