Anisotropic turbulent thermal diffusion and thermal convection in a rapidly rotating fluid sphere

Estimates of the molecular values of magnetic, viscous and thermal diffusion suggest that the state of the Earth’s core is turbulent and that complete numerical simulation of the geodynamo is not realizable at present. Large eddy simulation of the geodynamo with modelling of the sub-grid scale turbulence must be used. Current geodynamo models effectively model the sub-grid scale turbulence with isotropic diffusivities larger than the molecular values appropriate for the core. In the Braginsky and Meytlis (1990) picture of core turbulence the thermal and viscous diffusivities are enhanced up to the molecular magnetic diffusivity in the directions of the rotation axis and mean magnetic field. We neglect the mean magnetic field herein to isolate the effects of anisotropic thermal diffusion, enhanced or diminished along the rotation axis, and explore the instability of a steady conductive basic state with zero mean flow in the Boussinesq approximation. This state is found to be more stable (less stable) as the thermal diffusion parallel to the rotation axis is increased (decreased), if the transverse thermal diffusion is fixed. To examine the effect of simultaneously varying the diffusion along and transverse to the rotation axis, the Frobenius norm is used to control for the total thermal diffusion. When the Frobenius norm of the thermal diffusion tensor is fixed, it is found that increasing the thermal diffusion parallel to the rotation axis is destabilising. This result suggests that, for a fixed total thermal diffusion, geodynamo codes with anisotropic thermal diffusion may operate at lower modified Rayleigh numbers.