In this study, we consider plasmas evolving in a high Reynolds number regime leading to characteristic length scales associated with viscous effects that are generally much smaller than any chosen grid size. In this context, an artificial subgrid-scale model of these viscous terms \(\tau _{\mathfrak {a}\mathcal {I}},\,\mathfrak {a}\mathcal {I}\in \mathcal {M}\)  or \(\mathcal {M}_{\text{TF}}\)  has been considered in order to represent shocks and turbulence. This approach is inspired from [1]}, [2]} and has been extended for MFMS and TF models.
