Prandtl number dependence of stellar convection: Flow statistics and convective energy transport
Abstract
Context. The ratio of kinematic viscosity to thermal diffusivity, the Prandtl number, is much smaller than unity in stellar convection zones.
Aims: The main goal of this work is to study the statistics of convective flows and energy transport as functions of the Prandtl number.
Methods: Threedimensional numerical simulations of compressible nonrotating hydrodynamic convection in Cartesian geometry are used. The convection zone (CZ) is embedded between two stably stratified layers. The dominant contribution to the diffusion of entropy fluctuations comes in most cases from a subgridscale diffusivity whereas the mean radiative energy flux is mediated by a diffusive flux employing Kramers opacity law. Here, we study the statistics and transport properties of up and downflows separately.
Results: The volumeaveraged rms velocity increases with decreasing Prandtl number. At the same time, the filling factor of downflows decreases and leads to, on average, stronger downflows at lower Prandtl numbers. This results in a strong dependence of convective overshooting on the Prandtl number. Velocity power spectra do not show marked changes as a function of Prandtl number except near the base of the convective layer where the dominance of vertical flows is more pronounced. At the highest Reynolds numbers, the velocity power spectra are more compatible with the BolgianoObukhov k^{−11/5} than the KolmogorovObukhov k^{−5/3} scaling. The horizontally averaged convected energy flux (F̅_{conv}), which is the sum of the enthalpy (F̅_{enth}) and kinetic energy fluxes (F̅_{kin}), is independent of the Prandtl number within the CZ. However, the absolute values of F̅_{enth} and F̅_{kin} increase monotonically with decreasing Prandtl number. Furthermore, F̅_{enth} and F̅_{kin} have opposite signs for downflows and their sum F̅^{↓}_{conv} diminishes with Prandtl number. Thus, the upflows (downflows) are the dominant contribution to the convected flux at low (high) Prandtl numbers. These results are similar to those from RayleighBenárd convection in the low Prandtl number regime where convection is vigorously turbulent but inefficient at transporting energy.
Conclusions: The current results indicate a strong dependence of convective overshooting and energy flux on the Prandtl number. Numerical simulations of astrophysical convection often use a Prandtl number of unity because it is numerically convenient. The current results suggest that this can lead to misleading results and that the astrophysically relevant low Prandtl number regime is qualitatively different from the parameter regimes explored in typical contemporary simulations.
 Publication:

Astronomy and Astrophysics
 Pub Date:
 November 2021
 DOI:
 10.1051/00046361/202141337
 arXiv:
 arXiv:2105.08453
 Bibcode:
 2021A&A...655A..78K
 Keywords:

 turbulence;
 convection;
 Astrophysics  Solar and Stellar Astrophysics;
 Physics  Fluid Dynamics
 EPrint:
 15 pages, 20 figures, submitted to Astron. Astrophys