Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userImproving 1D Stellar Models with 3D Atmospheres
86
0
0.0
(
0
)
Added by
Jakob R{\\o}rsted Mosumgaard
Publication date
2016
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Stellar evolution codes play a major role in present-day astrophysics, yet they share common issues. In this work we seek to remedy some of those by the use of results from realistic and highly detailed 3D hydrodynamical simulations of stellar atmospheres. We have implemented a new temperature stratification extracted directly from the 3D simulations into the Garching Stellar Evolution Code to replace the simplified atmosphere normally used. Secondly, we have implemented the use of a variable mixing-length parameter, which changes as a function of the stellar surface gravity and temperature -- also derived from the 3D simulations. Furthermore, to make our models consistent, we have calculated new opacity tables to match the atmospheric simulations. Here, we present the modified code and initial results on stellar evolution using it.
rate research
Read More
In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.
CONTEXT: In recent years, the solar chemical abundances have been studied in considerable detail because of discrepant values of solar metallicity inferred from different indicators, i.e., on the one hand, the sub-solar photospheric abundances resulting from spectroscopic chemical composition analyses with the aid of 3D hydrodynamical models of the solar atmosphere, and, on the other hand, the high metallicity inferred by helioseismology. AIMS: After investigating the solar oxygen abundance using a CO5BOLD 3D hydrodynamical solar model in previous work, we undertake a similar approach studying the solar abundance of nitrogen, since this element accounts for a significant fraction of the overall solar metallicity, Z. METHOD: We used a selection of atomic spectral lines to determine the solar nitrogen abundance, relying mainly on equivalent width measurements in the literature. We investigate the influence on the abundance analysis, of both deviations from local thermodynamic equilibrium (NLTE effects) and photospheric inhomogeneities (granulation effects). RESULTS: We recommend use of a solar nitrogen abundance of A(N)=7.86+-0.12 whose error bar reflects the line-to-line scatter. CONCLUSION: The solar metallicity implied by the CO5BOLD-based nitrogen and oxygen abundances is in the range 0.0145<= Z <= 0.0167. This result is a step towards reconciling photospheric abundances with helioseismic constraints on Z. Our most suitable estimates are Z=0.0156 and Z/X=0.0213.
We present an overview of the current status of our efforts to derive the microturbulence and macroturbulence parameters (ximic and ximac) from the CIFIST grid of CO5BOLD 3D model atmospheres as a function of the basic stellar parameters Teff, log g, and [M/H]. The latest results for the Sun and Procyon show that the derived microturbulence parameter depends significantly on the numerical resolution of the underlying 3D simulation, confirming that `low-resolution models tend to underestimate the true value of ximic. Extending the investigation to twelve further simulations with different Teff, log g, and [M/H], we obtain a first impression of the predicted trend of ximic over the Hertzsprung-Russell diagram: in agreement with empirical evidence, microturbulence increases towards higher effective temperature and lower gravity. The metallicity dependence of ximic must be interpreted with care, since it also reflects the deviation between the 1D and 3D photospheric temperature stratifications that increases systematically towards lower metallicity.
We present the implementation of a radiative transfer solver with coherent scattering in the new BIFROST code for radiative magneto-hydrodynamical (MHD) simulations of stellar surface convection. The code is fully parallelized using MPI domain decomposition, which allows for large grid sizes and improved resolution of hydrodynamical structures. We apply the code to simulate the surface granulation in a solar-type star, ignoring magnetic fields, and investigate the importance of coherent scattering for the atmospheric structure. A scattering term is added to the radiative transfer equation, requiring an iterative computation of the radiation field. We use a short-characteristics-based Gauss-Seidel acceleration scheme to compute radiative flux divergences for the energy equation. The effects of coherent scattering are tested by comparing the temperature stratification of three 3D time-dependent hydrodynamical atmosphere models of a solar-type star: without scattering, with continuum scattering only, and with both continuum and line scattering. We show that continuum scattering does not have a significant impact on the photospheric temperature structure for a star like the Sun. Including scattering in line-blanketing, however, leads to a decrease of temperatures by about 350,K below log tau < -4. The effect is opposite to that of 1D hydrostatic models in radiative equilibrium, where scattering reduces the cooling effect of strong LTE lines in the higher layers of the photosphere. Coherent line scattering also changes the temperature distribution in the high atmosphere, where we observe stronger fluctuations compared to a treatment of lines as true absorbers.
3D hydrodynamics models of deep stellar convection exhibit turbulent entrainment at the convective-radiative boundary which follows the entrainment law, varying with boundary penetrability. We implement the entrainment law in the 1D Geneva stellar evolution code. We then calculate models between 1.5 and 60 M$_{odot}$ at solar metallicity ($Z=0.014$) and compare them to previous generations of models and observations on the main sequence. The boundary penetrability, quantified by the bulk Richardson number, $Ri_{mathrm{B}}$, varies with mass and to a smaller extent with time. The variation of $Ri_{mathrm{B}}$ with mass is due to the mass dependence of typical convective velocities in the core and hence the luminosity of the star. The chemical gradient above the convective core dominates the variation of $Ri_{mathrm{B}}$ with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through $Ri_{mathrm{B}}$. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters $A sim 2 times 10^{-4}$ and $n=1$. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
