  The inclusion of kinetic effects into fluid models has been a long standing
problem in magnetic reconnection and plasma physics. Generally the pressure
tensor is reduced to a scalar which is an approximation used to aid in the
modeling of large scale global systems such as the Earth's magnetosphere. This
unfortunately omits important kinetic physics which have been shown to play a
crucial role in collisionless regimes. The multi-fluid 10-moment model on the
other-hand retains the full symmetric pressure tensor. The 10-moment model is
constructed by taking moments of the Vlasov equation up to second order, and
includes the scalar density, the vector bulk-flow, and the symmetric pressure
tensor for a total of 10 separate components. Use of the multi-fluid 10-moment
model requires a closure which truncates the cascading system of equations.
Here we look to leverage data-driven methodologies to seek a closure which may
improve physical fidelity of the 10-moment multi-fluid model in collisionless
regimes. Specifically we use the Sparse Identification of Nonlinear Dynamics
(SINDy) method for symbolic equation discovery to seek the truncating closure
from fully kinetic particle-in-cell simulation data, which inherently retains
the relevant kinetic physics. We verify our method by reproducing the 10-moment
model from the PIC particle data and use the method to generate a closure
truncating the 10-moment model which is analyzed through the nonlinear phase of
reconnection.
