Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userUpward Overshooting in Turbulent Compressible Convection. III. Calibrate Parameters for One-dimensional Reynolds Stress Model
117
0
0.0
(
0
)
Authors
Tao Cai
Ask ChatGPT about the research
No Arabic abstract
In this paper, we calibrate the coefficients for the one-dimensional Reynolds stress model with the data generated from the three-dimensional numerical simulations of upward overshooting in turbulent compressible convection. It has been found that the calibrated convective and isotropic coefficients are almost the same as those calibrated in the pure convection zone. However, the calibrated diffusive coefficients differ significantly from those calibrated in the pure convection zone. We suspect that the diffusive effect induced by the boundary is stronger than by the adjacent stable zone. We have checked the validity of the downgradient approximation. We find that the prediction of the downgradient approximation on the third-order moments is unsatisfactory. However, the prediction on their derivatives is much better. It explains why the performance of the Reynolds stress model is reasonable in application to the real stars. With the calibrated coefficients, we have solved the full set of nonlocal turbulent equations on Reynolds stress model. We find that the Reynolds stress model has successfully produced the thermal adjustment layer and turbulent dissipation layer, which were identified in the three-dimensional numerical simulations. We suggest to use the inflection point of the auto-correlation of temperature perturbation and the Peclet number as the indicators on measuring the extents of the thermal adjustment layer and turbulent dissipation layer, respectively. This result may offer a practical guidance on the application of the Reynolds stress model in 1D stellar structure and evolution models.
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.
The Reynolds-averaged Navier-Stokes (RANS) equations are widely used in turbulence applications. They require accurately modeling the anisotropic Reynolds stress tensor, for which traditional Reynolds stress closure models only yield reliable results in some flow configurations. In the last few years, there has been a surge of work aiming at using data-driven approaches to tackle this problem. The majority of previous work has focused on the development of fully-connected networks for modeling the anisotropic Reynolds stress tensor. In this paper, we expand upon recent work for turbulent channel flow and develop new convolutional neural network (CNN) models that are able to accurately predict the normalized anisotropic Reynolds stress tensor. We apply the new CNN model to a number of one-dimensional turbulent flows. Additionally, we present interpretability techniques that help drive the model design and provide guidance on the model behavior in relation to the underlying physics.
Context: We study the impact of two-dimensional spherical shells on compressible convection. Realistic profiles for density and temperature from a one-dimensional stellar evolution code are used to produce a model of a large stellar convection zone representative of a young low-mass star. Methods: We perform hydrodynamic implicit large-eddy simulations of compressible convection using the MUltidimensional Stellar Implicit Code (MUSIC). Because MUSIC has been designed to use realistic stellar models produced from one-dimensional stellar evolution calculations, MUSIC simulations are capable of seamlessly modeling a whole star. Simulations in two-dimensional spherical shells that have different radial extents are performed over hundreds of convective turnover times, permitting the collection of well-converged statistics. Results: We evaluate basic statistics of the convective turnover time, the convective velocity, and the overshooting layer. These quantities are selected for their relevance to one-dimensional stellar evolution calculations, so that our results are focused toward the 321D link. The inclusion in the spherical shell of the boundary between the radiative and convection zones decreases the amplitude of convective velocities in the convection zone. The inclusion of near-surface layers in the spherical shell can increase the amplitude of convective velocities, although the radial structure of the velocity profile established by deep convection is unchanged. The impact from including the near-surface layers depends on the speed and structure of small-scale convection in the near-surface layers. Larger convective velocities in the convection zone result in a commensurate increase in the overshooting layer width and decrease in the convective turnover time. These results provide support for non-local aspects of convection.
The extent of mixed regions around convective zones is one of the biggest uncertainties in stellar evolution. 1D overshooting descriptions introduce a free parameter ($f_{ov}$) that is in general not well constrained from observations. Especially in small central convective regions the value is highly uncertain due to its tight connection to the pressure scale height. Long-term multi-dimensional hydrodynamic simulations can be used to study the size of the overshooting region and the involved mixing processes. Here we show how one can calibrate an overshooting parameter by performing 2D Maestro simulations of Zero-Age-Main-Sequence stars ranging from $1.3$ to $3.5 M_odot$. The simulations cover the convective cores of the stars and a large fraction of the surrounding radiative envelope. We follow the convective flow for at least 20 convective turnover times, while the longest simulation covers 430 turnover time scales. This allows us to study how the mixing as well as the convective boundary evolve with time, and how the resulting entrainment can be interpreted in terms of overshooting parameters. We find that increasing the overshooting parameter $f_{ov}$ beyond a certain value in the initial model of our simulations, changes the mixing behaviour completely. This result can be used to put limits on the overshooting parameter. We find $0.010 < f_{ov} < 0.017$ to be in good agreement with our simulations of a $3.5 M_odot$ mass star. We also identify a diffusive mixing component due to internal gravity waves (IGW) that is active throughout the convectively stable layer, but likely overestimated in our simulations. Furthermore, applying our calibration method to simulations of less massive stars suggests a need for a mass-dependent overshooting description where the mixing in terms of the pressure scale height is reduced for small convective cores.
We present numerical simulations of hydrodynamic overshooting convection in local Cartesian domains. We find that a substantial fraction of the lower part of the convection zone (CZ) is stably stratified according to the Schwarzschild criterion while the enthalpy flux is outward directed. This occurs when the heat conduction profile at the bottom of the CZ is smoothly varying, based either on a Kramers-like opacity prescription as a function of temperature and density or a static profile of a similar shape. We show that the subadiabatic layer arises due to nonlocal energy transport by buoyantly driven downflows in the upper parts of the CZ. Analysis of the force balance of the upflows and downflows confirms that convection is driven by cooling at the surface. We find that the commonly used prescription for the convective enthalpy flux being proportional to the negative entropy gradient does not hold in the stably stratified layers where the flux is positive. We demonstrate the existence of a non-gradient contribution to the enthalpy flux, which is estimated to be important throughout the convective layer. A quantitative analysis of downflows indicates a transition from a tree-like structure where smaller downdrafts merge into larger ones in the upper parts to a structure in the deeper parts where a height-independent number of strong downdrafts persist. This change of flow topology occurs when a substantial subadiabatic layer is present in the lower part of the CZ.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
