Data-driven multi-species heat flux closures for two-stream-unstable plasmas with nonlinear sparse regression

TL;DR

This paper develops data-driven multi-species heat flux closures for plasma models using nonlinear sparse regression, improving accuracy in capturing kinetic physics in collisionless plasmas.

Contribution

It extends previous single-species closure models to multi-species systems and introduces a new nonlinear sparse regression framework for estimating closure coefficients.

Findings

Models predict heat flux with 80-90% accuracy.

Closure models account for 85-95% of pressure change rate.

Framework compares favorably with linear collisionless theory.

Abstract

The dual aims of accuracy and computational efficiency in computational plasma physics lend themselves well to the use of fluid models. The first of these goals, however, is only satisfied for such models insofar as the utilized closure can capture the neglected kinetic physics -- something which has proven challenging for multi-scale collisionless processes. In a recent article [E. R. Ingelsten et al. (2025) J. Plasma Phys. 91 E64], we used the data-driven method of sparse regression to discover a novel heat flux closure for electrostatic phenomena. Here, we generalize the six-term closure model found in that work from single- to multi-species modeling. Using data from OSIRIS particle-in-cell simulations over a range of initial conditions, we then demonstrate how the unknown coefficients in front of the three most important terms in the closure can be estimated from box-averaged fluid quantities. Both neural networks and a newly developed framework for nonlinear sparse regression are showcased. The resulting models predict the heat flux for each species with a typical accuracy of 80-90 % and regularly account for 85-95 % of the rate of change in the pressure. The models are also compared with results from multi-species linear collisionless theory.