The simplest formulation of the problem of thermal convection of a fluid is called the Rayleigh–Bénard problem. It describes convection of a liquid enclosed between two differentially heated horizontal plates, each held at a fixed temperature. The mathematical model is derived under the Boussinesq approximation, which amounts to assuming that the flow is highly subsonic and that density perturbations to a uniform and constant background density are negligible everywhere except in the buoyancy term , where stands for the gravity (Chandrasekhar, 1961). The equilibrium background state is a linear temperature profile with temperature decreasing from the bottom to the top of the layer. This case is in many respects different and simpler than the strongly stratified SCZ case, which treatment requires using more general compressible fluid and energy equations than those given below (Nordlund, 1982), but is sufficient to discuss many of the important physical (Section 5.2) and numerical (Section 6.1) issues pertaining to supergranulation.
Anticipating upcoming discussions on the origin of supergranulation, we extend the simplest hydrodynamic formulation of the Rayleigh–Bénard problem to the case of an electrically conducting liquid threaded by a mean vertical magnetic field denoted by and rotating around a vertical axis, with a rotation rate . This set-up is shown on Figure 9.
In nondimensional form, the equations for momentum and energy conservation, the induction equation, the equations for mass conservation and magnetic field solenoidality read
Several important numbers appear in the equations above, starting with the Rayleigh numberChandrasekhar number Taylor number, thermal and magnetic Prandtl numbers (see Section 2.2).
In the simplest non-rotating hydrodynamic case (, no induction), when the Rayleigh number is less than a critical value that depends on the particular choice of boundary conditions, diffusive processes dominate over buoyancy: the hydrostatic solution is stable, i.e., any velocity or temperature perturbations decays. For , convection sets in as a linear instability and perturbations grow exponentially in the form of convection rolls or hexagons with a horizontal spatial periodicity comparable to the convective layer depth in most cases. The effects of magnetic fields and rotation on the linear stability analysis are discussed in the next paragraphs.
It should be noted that , , and are all extremely large numbers in the Sun, if they are computed from microscopic transport coefficients (Section 2.2). So, in principle, there is no reason why solar convection should be close to the instability threshold. However, theoretical studies of large-scale convection (such as supergranulation) commonly assume that viscous, thermal, and magnetic diffusion at such scales are determined by turbulent transport, not microscopic transport. This leads to much larger transport coefficients (which can be estimated, for instance, using the typical scale and velocity of the granulation pattern) and much smaller “effective” , , and , so the “large-scale” system is generally considered not too far away from criticality. Making this (strong) mean-field assumption serves to legitimate using the standard toolkits of linear and weakly nonlinear analysis to understand the large-scale behaviour of solar convection.
This work is licensed under a Creative Commons License.