We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

This article is focused on the dynamics of a rotating electrically conducting fluid in a turbulent state. As inside the Earth's core or in various industrial processes, a flow is altered by the presence of both background rotation and a large scale magnetic field. In this context, we present a set of 3D direct numerical simulations of incompressible decaying turbulence. We focus on parameters similar to the ones encountered in geophysical and astrophysical flows, so that the Rossby number is small, the interaction parameter is large, but the Elsasser number, defining the ratio between Coriolis and Lorentz forces, is about unity. These simulations allow to quantify the effect of rotation and thus inertial waves on the growth of magnetic fluctuations due to Alfv\'en waves. Rotation prevents the occurrence of equipartition between kinetic and magnetic energies, with a reduction of magnetic energy at decreasing Elsasser number {\Lambda}. It also causes a decrease of energy transfer mediated by cubic correlations. In terms of flow structure, a decrease of {\Lambda} corresponds to an increase in the misalignment of velocity and magnetic field.