Stellar convection. III - Convection at large Rayleigh numbers

The results of a numerical, three-dimensional study of convection in a self-gravitating sphere of Boussinesq fluid with a Rayleigh number of 10-billion and a Prandtl number of 1 are presented. The velocity and temperature are computed by using spectral methods in the horizontal and finite-differencing in the radial directions. An eddy viscosity and diffusivity are needed to model the subresolution flow. For Rayleigh numbers much less than 10-billion the flows do not have well-defined inertial ranges, and an eddy viscosity and diffusivity cannot be assigned in this self-consistent manner. By computing the energy spectra as well as the detailed energy budgets as a function of wavenumber, it is shown that for Rs = 10-billion there is an inertial range for the modes corresponding to spherical harmonics with l greater than 6.