Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPrandtl number dependence of compressible convection: Flow statistics and convective energy transport
79
0
0.0
(
0
)
Added by
Petri K\\\"apyl\\\"a
Publication date
2021
fields
Physics
and research's language is
English
Authors
Petri J. Kapyla
Ask ChatGPT about the research
No Arabic abstract
(abridged) Context: The ratio of kinematic viscosity to thermal diffusivity, the Prandtl number, is much smaller than unity in stellar convection zones. Aims: To study the statistics of convective flows and energy transport as functions of the Prandtl number. Methods: Three-dimensional numerical simulations convection in Cartesian geometry are used. The convection zone (CZ) is embedded between two stably stratified layers. Statistics and transport properties of up- and downflows are studied separately. Results: The rms velocity increases with decreasing Prandtl number. At the same time the filling factor of downflows decreases and leads to stronger downflows at lower Prandtl numbers, and to a strong dependence of overshooting on the Prandtl number. Velocity power spectra do not show marked changes as a function of Prandtl number. At the highest Reynolds numbers the velocity power spectra are compatible with the Bolgiano-Obukhov $k^{-11/5}$ scaling. The horizontally averaged convected energy flux ($overline{F}_{rm conv}$) is independent of the Prandtl number within the CZ. However, the upflows (downflows) are the dominant contribution to the convected flux at low (high) Prandtl number. These results are similar to those from Rayleigh-Benard 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 Prandtl number of unity. 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 parameters regimes explored in typical simulations.
rate research
Read More
In this paper, we investigate the upward overshooting by three-dimensional numerical simulations. We find that the above convectively stable zone can be partitioned into three layers: the thermal adjustment layer (mixing both entropy and material), the turbulent dissipation layer (mixing material but not entropy), and the thermal dissipation layer (mixing neither entropy nor material). The turbulent dissipation layer is separated from the thermal adjustment layer and the thermal dissipation layer by the first and second zero points of the vertical velocity correlation. The simulation results are in good agreement with the prediction of the one-dimensional turbulent Reynolds stress model. First, the layer structure is similar. Second, the upper boundary of the thermal adjustment layer is close to the peak of the magnitude of the temperature perturbation. Third, the Peclet number at the upper boundary of the turbulent dissipation layer is close to 1. In addition, we have studied the scalings of the overshooting distance on the relative stability parameter $S$, the Prandtl number $rm Pr$, and the Peclet number $rm Pe$. The scaling on $S$ is not unique. The trend is that the overshooting distance decreases with $S$. Fitting on $rm Pr$ shows that the overshooting distance increases with $rm Pr$. Fitting on $rm Pe$ shows that the overshooting distance decreases with $rm Pe$. Finally, we calculate the ratio of the thickness of the turbulent dissipation layer to that of the thermal adjustment layer. The ratio remains almost constant, with an approximate value of 2.4.
Results from direct numerical simulation for three-dimensional Rayleigh-Benard convection in samples of aspect ratio $Gamma=0.23$ and $Gamma=0.5$ up to Rayleigh number $Ra=2times10^{12}$ are presented. The broad range of Prandtl numbers $0.5<Pr<10$ is considered. In contrast to some experiments, we do not see any increase in $Nu/Ra^{1/3}$, neither due to $Pr$ number effects, nor due to a constant heat flux boundary condition at the bottom plate instead of constant temperature boundary conditions. Even at these very high $Ra$, both the thermal and kinetic boundary layer thicknesses obey Prandtl-Blasius scaling.
(abridged) Context: Turbulent diffusion of large-scale flows and magnetic fields play major roles in many astrophysical systems. Aims: Our goal is to compute turbulent viscosity and magnetic diffusivity, relevant for diffusing large-scale flows and magnetic fields, respectively, and their ratio, the turbulent magnetic Prandtl number, ${rm Pm}_{rm t}$, for isotropically forced homogeneous turbulence. Methods: We use simulations of forced turbulence in fully periodic cubes composed of isothermal gas with an imposed large-scale sinusoidal shear flow. Turbulent viscosity is computed either from the resulting Reynolds stress or from the decay rate of the large-scale flow. Turbulent magnetic diffusivity is computed using the test-field method. The scale dependence of the coefficients is studied by varying the wavenumber of the imposed sinusoidal shear and test fields. Results: We find that turbulent viscosity and magnetic diffusivity are in general of the same order of magnitude. Furthermore, the turbulent viscosity depends on the fluid Reynolds number (${rm Re}$) and scale separation ratio of turbulence. The scale dependence of the turbulent viscosity is found to be well approximated by a Lorentzian. The results for the turbulent transport coefficients appear to converge at sufficiently high values of ${rm Re}$ and the scale separation ratio. However, a weak decreasing trend is found even at the largest values of ${rm Re}$. The turbulent magnetic Prandtl number converges to a value that is slightly below unity for large ${rm Re}$ whereas for small ${rm Re}$, we find values between 0.5 and 0.6. Conclusions: The turbulent magnetic diffusivity is in general consistently higher than the turbulent viscosity. The actual value of ${rm Pm}_{rm t}$ found from the simulations ($approx0.9ldots0.95$) at large ${rm Re}$ and scale separation ratio is higher than any of the analytic predictions.
Extending our recent studies of two-dimensional stellar convection to 3D, we compare three-dimensional hydrodynamic simulations to identically set-up two-dimensional simulations, for a realistic pre-main sequence star. We compare statistical quantities related to convective flows including: average velocity, vorticity, local enstrophy, and penetration depth beneath a convection zone. These statistics are produced during stationary, steady-state compressible convection in the stars convection zone. Our simulations with the MUSIC code confirm the common result that two-dimensional simulations of stellar convection have a higher magnitude of velocity on average than three-dimensional simulations. Boundary conditions and the extent of the spherical shell can affect the magnitude and variability of convective velocities. The difference between 2D and 3D velocities is dependent on these background points; in our simulations this can have an effect as large as the difference resulting from the dimensionality of the simulation. Nevertheless, radial velocities near the convective boundary are comparable in our 2D and 3D simulations. The average local enstrophy of the flow is lower for two-dimensional simulations than for three-dimensional simulations, indicating a different shape and structuring of 3D stellar convection. We perform a statistical analysis of the depth of convective penetration below the convection zone, using the model proposed in our recent study (Pratt et al. 2017). Here we analyze the convective penetration in three dimensional simulations, and compare the results to identically set-up 2D simulations. In 3D the penetration depth is as large as the penetration depth calculated from 2D simulations.
A turbulent transport of radiation in the solar convective zone is investigated. The mean-field equation for the irradiation intensity is derived. It is shown that due to the turbulent effects, the effective penetration length of radiation can be increased in several times in comparison with the mean penetration length of radiation (defined as an inverse mean absorption coefficient). Using the model of the solar convective zone based on the mixing length theory, where the mean penetration length of radiation is usually much smaller than the turbulent correlation length, it is demonstrated that the ratio of the effective penetration length to the mean penetration length of radiation increases in 2.5 times in the vicinity of the solar surface. The main reason are the compressibility effects that become important in the vicinity of the solar surface where temperature and density fluctuations increase towards the solar surface, enhancing fluctuations of the radiation absorption coefficient and increasing the effective penetration length of radiation.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
