Implicit large eddy simulations of global solar convection: effects of numerical resolution in non-rotating and rotating cases

Abstract: Simulating deep solar convection and its coupled mean-field motions is a formidable challenge where few observational results constrain models that suffer from the non-physical influence of the grid resolution. We present hydrodynamic global Implicit Large-Eddy simulations (ILES) of deep solar convection performed with the EULAG-MHD code, and explore the effects of grid resolution on the properties of rotating and non-rotating convection. The results, based on low-order moments and turbulent spectra reveal that convergence could be achieved in non-rotating simulations provided sufficient resolution in the radial direction. The flow is highly anisotropic, with the energy contained in horizontal divergent motions exceeding by more than three orders of magnitude their radial counterpart. By contrast, in rotating simulations the largest energy is in the toroidal part of the horizontal motions. As the grid resolution increases, the turbulent correlations change in such a way that a solar-like differential rotation, obtained in the simulation with the coarser grid, transitions to the anti-solar differential rotation. The reason for this change is the contribution of the effective viscosity to the balance of the forces driving large-scale flows. As the effective viscosity decreases, the angular momentum balance improves, yet the force balance in the meridional direction lessens, favoring a strong meridional flow that advects angular momentum towards the poles. The results suggest that obtaining the correct distribution of angular momentum may not be a mere issue of numerical resolution. Accounting for additional physics, such as magnetism or the near-surface shear layer, may be necessary in simulating the solar interior.