Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCompressible Hydrodynamic Mean-Field Equations in Spherical Geometry and their Application to Turbulent Stellar Convection Data
135
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes.
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.
Small-scale dynamo action is often held responsible for the generation of quiet-Sun magnetic fields. We aim to determine the excitation conditions and saturation level of small-scale dynamos in non-rotating turbulent convection at low magnetic Prandtl numbers. We use high resolution direct numerical simulations of weakly stratified turbulent convection. We find that the critical magnetic Reynolds number for dynamo excitation increases as the magnetic Prandtl number is decreased, which might suggest that small-scale dynamo action is not automatically evident in bodies with small magnetic Prandtl numbers as the Sun. As a function of the magnetic Reynolds number (${rm Rm}$), the growth rate of the dynamo is consistent with an ${rm Rm}^{1/2}$ scaling. No evidence for a logarithmic increase of the growth rate with ${rm Rm}$ is found.
(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.
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.
In this study we investigate the effects of turbulent convection on formation of large-scale inhomogeneous magnetic structures by means of Large-Eddy Simulation (LES) for convection in solar-type stars. The main idea of this study is the implementation of a new subgrid-scale model for the effective Lorentz force in a three-dimensional nonlinear radiative magnetohydrodynamics (MHD) code developed for simulating the upper solar convection zone and lower atmosphere. To this end we derived the energy budget equations, which include the effects of the subgrid-scale turbulence on the Lorentz-force, and implemented the new subgrid-scale turbulence model (TELF-Model) in a three-dimensional nonlinear MHD LES code. Using imposed initial vertical and horizontal uniform magnetic fields in LES with the TELF-Model, we have shown that the magnetic flux tubes formation is started when the initial mean magnetic field is larger than a threshold value (about 100 G). This is in agreement with the theoretical studies by Rogachevskii and Kleeorin (2007). We have determined the vertical profiles of the velocity and magnetic fluctuations, total MHD energy and anisotropy of turbulent magneto-convection, kinetic and current and cross helicities.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
