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.

23 pages, 16 figures, submitted to ApJ

References in corpus (18)

• Structure and Evolution of Giant Cells in Global Models of Solar Convection
• Meridional Circulation in Solar and Stellar Convection Zones
• A simulation of convective dynamo in the solar convective envelope: maintenance of the solar-like differential rotation and emerging flux
• Rapidly Rotating Suns and Active Nests of Convection
• Magnetically controlled stellar differential rotation near the transition from solar to anti-solar profiles
• High-resolution calculation of the solar global convection with the reduced speed of sound technique: II. Near surface shear layer with the rotation
• Reconciling solar and stellar magnetic cycles with nonlinear dynamo simulations
• Reynolds stress and heat flux in spherical shell convection
• Helioseismic Imaging of Fast Convective Flows Throughout the Near-Surface Shear Layer
• The Emergence of Solar Supergranulation as a Natural Consequence of Rotationally-Constrained Interior Convection
• Solar differential rotation reproduced with high-resolution simulation
• A Comprehensive Method to Measure Solar Meridional Circulation and Center-to-Limb Effect Using Time-Distance Helioseismology
• The Solar Meridional Circulation and Sunspot Cycle Variability
• Understanding solar torsional oscillations from global dynamo models
• Helioseismology challenges models of solar convection
• Supergranulation and multiscale flows in the solar photosphere: Global observations vs. a theory of anisotropic turbulent convection
• Spatial Scales and Time Variation of Solar Subsurface Convection
• Numerical convergence of 2D solar convection in implicit large-eddy simulations