A key aspect in the determination of stellar properties is the comparison of observational constraints with predictions from stellar models. Asteroseismic Inference on a Massive Scale (AIMS) is an open source code that uses Bayesian statistics and a Markov Chain Monte Carlo approach to find a representative set of models that reproduce a given set of classical and asteroseismic constraints. These models are obtained by interpolation on a pre-calculated grid, thereby increasing computational efficiency. We test the accuracy of the different operational modes within AIMS for grids of stellar models computed with the Liège stellar evolution code (main sequence and red giants) and compare the results to those from another asteroseismic analysis pipeline, PARAM. Moreover, using artificial inputs generated from models within the grid (assuming the models to be correct), we focus on the impact on the precision of the code when considering different combinations of observational constraints (individual mode frequencies, period spacings, parallaxes, photospheric constraints,...). Our tests show the absolute limitations of precision on parameter inferences using synthetic data with AIMS, and the consistency of the code with expected parameter uncertainty distributions. Interpolation testing highlights the significance of the underlying physics to the analysis performance of AIMS and provides caution as to the upper limits in parameter step size. All tests demonstrate the flexibility and capability of AIMS as an analysis tool and its potential to perform accurate ensemble analysis with current and future asteroseismic data yields.
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⊙ 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, RiB, varies with mass and to a smaller extent with time. The variation of RiB 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 RiB with time. An entrainment law method can therefore explain the apparent mass dependence of convective boundary mixing through RiB. New models including entrainment can better reproduce the mass dependence of the main sequence width using entrainment law parameters A ∼ 2 × 10−4 and n = 1. We compare these empirically constrained values to the results of 3D hydrodynamics simulations and discuss implications.
Renzini (1987) wrote an influential critique of mixing-length theory (Böhm-Vitense 1958, MLT) as used in stellar evolution codes, and concluded that three-dimensional (3D) fluid dynamical simulations were needed to clarify several important issues. We have critically explored the limitations of the numerical methods and conclude that they are approaching the required accuracy Woodward et al. 2015). Implicit large eddy simulations (ILES) automatically connect large scale turbulence to a Kolmogorov cascade below the grid scale, allowing turbulent boundary layers to remove singularities that appear in the theory. Interactions between coherent structures give multi-modal behavior, driving intermittency and fluctuations. Reynolds averaging (RA) allows us to abstract the essential features of this dynamical behavior of boundaries which are appropriate to stellar evolution, and consider how they relate static boundary conditions (Richardson, Schwarzschild or Ledoux). We clarify several questions concerning when and why MLT works, and does not work, using both analytical theory and 3D high resolution numerical simulations. The composition gradients and boundary layer structure which are produced by our simulations suggest a self-consistent approach to boundary layers, removing the need for ad hoc procedures for 'convective overshooting' and 'semi-convection'. In a companion paper we quantify the adequacy of our numerical resolution, determine of the length scale of dissipation (the 'mixing length') without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and extend MLT to deal with strong stratification.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.