Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCalibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
104
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
As part of a larger program aimed at better quantifying the uncertainties in stellar computations, we attempt to calibrate the extent of convective overshooting in low to intermediate mass stars by means of eclipsing binary systems. We model 12 such systems, with component masses between 1.3 and 6.2 solar masses, using the detailed binary stellar evolution code STARS, producing grids of models in both metallicity and overshooting parameter. From these, we determine the best fit parameters for each of our systems. For three systems, none of our models produce a satisfactory fit. For the remaining systems, no single value for the convective overshooting parameter fits all the systems, but most of our systems can be well described with an overshooting parameter between 0.09 and 0.15, corresponding to an extension of the mixed region above the core of about 0.1-0.3 pressure scale heights. Of the nine systems where we are able to obtain a good fit, seven can be reasonably well fit with a single parameter of 0.15. We find no evidence for a trend of the extent of overshooting with either mass or metallicity, though the data set is of limited size. We repeat our calculations with a second evolution code, MESA, and we find general agreement between the two codes. For the extension of the mixed region above the convective core required by the MESA models is about 0.15-0.4 pressure scale heights. For the system EI Cep, we find that MESA gives an overshooting region that is larger than the STARS one by about 0.1 pressure scale heights for the primary, while for the secondary the difference is only 0.05 pressure scale heights.
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.
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 novel parallel technique for Fourier-Galerkin pseudo-spectral methods with applications to two-dimensional Navier-Stokes equations and inviscid Boussinesq approximation equations is presented. It takes the advantage of the programming structure of the phase-shift de-aliased scheme for pseudo-spectral codes, and combines the task-distribution strategy [Yin, Clercx and Montgomery, Comput. Fluids, 33, 509 (2004)] and parallelized Fast Fourier Transform scheme. The performances of the resulting MPI Fortran90 codes with the new procedure on SGI 3800 are reported. For fixed resolution of the same problem, the peak speed of the new scheme can be twice as fast as the old parallel methods. The parallelized codes are used to solve some challenging numerical problems governed by the Navier-Stokes equations and the Boussinesq equations. Two interesting physical problems, namely, the double-valued $omega$-$psi $ structure in two-dimensional decaying turbulence and the collapse of the bubble cap in the Boussinesq simulation, are solved by using the proposed parallel algorithms.
Origin of hydrodynamical instability and turbulence in the Keplerian accretion disc as well as similar laboratory shear flows, e.g. plane Couette flow, is a long standing puzzle. These flows are linearly stable. Here we explore the evolution of perturbation in such flows in the presence of an additional force. Such a force, which is expected to be stochastic in nature hence behaving as noise, could be result of thermal fluctuations (however small be), Brownian ratchet, grain-fluid interactions and feedback from outflows in astrophysical discs etc. We essentially establish the evolution of nonlinear perturbation in the presence of Coriolis and external forces, which is modified Landau equation. We show that even in the linear regime, under suitable forcing and Reynolds number, the otherwise least stable perturbation evolves to a very large saturated amplitude, leading to nonlinearity and plausible turbulence. Hence, forcing essentially leads a linear stable mode to unstable. We further show that nonlinear perturbation diverges at a shorter timescale in the presence of force, leading to a fast transition to turbulence. Interestingly, emergence of nonlinearity depends only on the force but not on the initial amplitude of perturbation, unlike original Landau equation based solution.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
