-
Effects of stratification on overshooting and waves atop the convective core of $5M_{\odot}$ main-sequence stars
Authors:
Adrien Morison,
Arthur Le Saux,
Isabelle Baraffe,
Jack Morton,
Thomas Guillet,
Dimitar Vlaykov,
Tom Goffrey,
Jane Pratt
Abstract:
As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate 2D convection in the core at three different stages of evolution of a $5M_{\odot}$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably-stratified envelope, along with the overshooting le…
▽ More
As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate 2D convection in the core at three different stages of evolution of a $5M_{\odot}$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably-stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.
△ Less
Submitted 8 July, 2024;
originally announced July 2024.
-
Two-dimensional simulations of internal gravity waves in a 5 $M_{\odot}$ Zero-Age-Main-Sequence model
Authors:
A. Le Saux,
I. Baraffe,
T. Guillet,
D. G. Vlaykov,
A. Morison,
J. Pratt,
T. Constantino,
T. Goffrey
Abstract:
Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-…
▽ More
Main-sequence intermediate-mass stars present a radiative envelope that supports internal gravity waves (IGWs). Excited at the boundary with the convective core, IGWs propagate towards the stellar surface and are suspected to impact physical processes such as rotation and chemical mixing. Using the fully compressible time-implicit code MUSIC, we study IGWs in two-dimensional simulations of a zero-age-main-sequence 5 solar mass star model up to 91\% of the stellar radius with different luminosity and radiative diffusivity enhancements. Our results show that low frequency waves excited by core convection are strongly impacted by radiative effects as they propagate. This impact depends on the radial profile of radiative diffusivity which increases by almost 5 orders of magnitude between the centre of the star and the top of the simulation domain. In the upper layers of the simulation domain, we observe an increase of the temperature. Our study suggests that this is due to heat added in these layers by IGWs damped by radiative diffusion. We show that non-linear effects linked to large amplitude IGWs may be relevant just above the convective core. Both these effects are intensified by the artificial enhancement of the luminosity and radiative diffusivity, with enhancement factors up to $10^4$ times the realistic values. Our results also highlight that direct comparison between numerical simulations with enhanced luminosity and observations must be made with caution. Finally, our work suggests that thermal effects linked to the damping of IGWs could have a non-negligible impact on stellar structure.
△ Less
Submitted 5 April, 2023;
originally announced April 2023.
-
A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations
Authors:
I. Baraffe,
J. Clarke,
A. Morison,
D. G. Vlaykov,
T. Constantino,
T. Goffrey,
T. Guillet,
A. Le Saux,
J. Pratt
Abstract:
We perform two-dimensional numerical simulations of core convection for zero-age-main-sequence stars covering a mass range from 3 $M_\odot$ to 20 $M_\odot$. The simulations are performed with the fully compressible time-implicit code MUSIC. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mas…
▽ More
We perform two-dimensional numerical simulations of core convection for zero-age-main-sequence stars covering a mass range from 3 $M_\odot$ to 20 $M_\odot$. The simulations are performed with the fully compressible time-implicit code MUSIC. We study the efficiency of overshooting, which describes the ballistic process of convective flows crossing a convective boundary, as a function of stellar mass and luminosity. We also study the impact of artificially increasing the stellar luminosity for 3 $M_\odot$ models. The simulations cover hundreds to thousands of convective turnover timescales. Applying the framework of extreme plume events previously developed for convective envelopes, we derive overshooting lengths as a function of stellar masses. We find that the overshooting distance ($d_{\rm ov}$) scales with the stellar luminosity ($L$) and the convective core radius ($r_{\rm conv}$). We derive a scaling law $d_{\rm ov} \propto L^{1/3} r_{\rm conv}^{1/2}$ which is implemented in a 1D stellar evolution code and the resulting stellar models are compared to observations. The scaling predicts values for the overshooting distance that significantly increase with stellar mass, in qualitative agreement with observations. Quantitatively, however, the predicted values are underestimated for masses $\gtrsim 10 M_\odot$. Our 2D simulations show the formation of a nearly-adiabatic layer just above the Schwarzschild boundary of the convective core, as exhibited in recent 3D simulations of convection. The most luminous models show a growth in size with time of the nearly-adiabatic layer. This growth seems to slow down as the upper edge of the nearly-adiabatic layer gets closer to the maximum overshooting length and as the simulation time exceeds the typical thermal diffusive timescale in the overshooting layer.
△ Less
Submitted 6 January, 2023;
originally announced January 2023.
-
Impact of radial truncation on global 2D hydrodynamic simulations for a Sun-like model
Authors:
D. G. Vlaykov,
I. Baraffe,
T. Constantino,
T. Goffrey,
T. Guillet,
A. Le Saux,
A. Morison,
J. Pratt
Abstract:
Stellar convection is a non-local process responsible for the transport of heat and chemical species. It can lead to enhanced mixing through convective overshooting and excitation of internal gravity waves (IGWs) at convective boundaries. The relationship between these processes is still not well understood and requires global hydrodynamic simulations to capture the important large-scale dynamics.…
▽ More
Stellar convection is a non-local process responsible for the transport of heat and chemical species. It can lead to enhanced mixing through convective overshooting and excitation of internal gravity waves (IGWs) at convective boundaries. The relationship between these processes is still not well understood and requires global hydrodynamic simulations to capture the important large-scale dynamics. The steep stratification in stellar interiors suggests that the radial extent of such simulations can affect the convection dynamics, the IGWs in the stably stratified radiative zone, and the depth of the overshooting layer. We investigate these effects using two-dimensional global simulations performed with the fully compressible stellar hydrodynamics code MUSIC. We compare eight different radial truncations of the same solar-like stellar model evolved over approximately 400 convective turnover times. We find that the location of the inner boundary has an insignificant effect on the convection dynamics, the convective overshooting and the travelling IGWs. We relate this to the background conditions at the lower convective boundary which are unaffected by the truncation, as long as a significantly deep radiative layer is included in the simulation domain. However, we find that extending the outer boundary by only a few percent of the stellar radius significantly increases the velocity and temperature perturbations in the convection zone, the overshooting depth, the power and the spectral slope of the IGWs. The effect is related to the background conditions at the outer boundary, which are determined in essence by the hydrostatic stratification and the given luminosity.
△ Less
Submitted 9 May, 2022;
originally announced May 2022.
-
Two-dimensional simulations of solar-like models with artificially enhanced luminosity. II. Impact on internal gravity waves
Authors:
A. Le Saux,
T. Guillet,
I. Baraffe,
D. G. Vlaykov,
T. Constantino,
J. Pratt,
T. Goffrey,
M. Sylvain,
V. Réville,
A. S. Brun
Abstract:
Artificially increasing the luminosity and the thermal diffusivity of a model is a common tactic adopted in hydrodynamical simulations of stellar convection. In this work, we analyse the impact of these artificial modifications on the physical properties of stellar interiors and specifically on internal gravity waves. We perform two-dimensional simulations of solar-like stars with the MUSIC code.…
▽ More
Artificially increasing the luminosity and the thermal diffusivity of a model is a common tactic adopted in hydrodynamical simulations of stellar convection. In this work, we analyse the impact of these artificial modifications on the physical properties of stellar interiors and specifically on internal gravity waves. We perform two-dimensional simulations of solar-like stars with the MUSIC code. We compare three models with different luminosity enhancement factors to a reference model. The results confirm that properties of the waves are impacted by the artificial enhancement of the luminosity and thermal diffusivity. We find that an increase in the stellar luminosity yields a decrease in the bulk convective turnover timescale and an increase in the characteristic frequency of excitation of the internal waves. We also show that a higher energy input in a model, corresponding to a larger luminosity, results in higher energy in high frequency waves. Across our tests with the luminosity and thermal diffusivity enhanced together by up to a factor of 104, our results are consistent with theoretical predictions of radiative damping. Increasing the luminosity also has an impact on the amplitude of oscillatory motions across the convective boundary. One must use caution when interpreting studies of internal gravity waves based on hydrodynamical simulations with artificially enhanced luminosity.
△ Less
Submitted 1 February, 2022;
originally announced February 2022.
-
Dynamics in a stellar convective layer and at its boundary: Comparison of five 3D hydrodynamics codes
Authors:
R. Andrassy,
J. Higl,
H. Mao,
M. Mocák,
D. G. Vlaykov,
W. D. Arnett,
I. Baraffe,
S. W. Campbell,
T. Constantino,
P. V. F. Edelmann,
T. Goffrey,
T. Guillet,
F. Herwig,
R. Hirschi,
L. Horst,
G. Leidi,
C. Meakin,
J. Pratt,
F. Rizzuti,
F. K. Roepke,
P. Woodward
Abstract:
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility o…
▽ More
Our ability to predict the structure and evolution of stars is in part limited by complex, 3D hydrodynamic processes such as convective boundary mixing. Hydrodynamic simulations help us understand the dynamics of stellar convection and convective boundaries. However, the codes used to compute such simulations are usually tested on extremely simple problems and the reliability and reproducibility of their predictions for turbulent flows is unclear. We define a test problem involving turbulent convection in a plane-parallel box, which leads to mass entrainment from, and internal-wave generation in, a stably stratified layer. We compare the outputs from the codes FLASH, MUSIC, PPMSTAR, PROMPI, and SLH, which have been widely employed to study hydrodynamic problems in stellar interiors. The convection is dominated by the largest scales that fit into the simulation box. All time-averaged profiles of velocity components, fluctuation amplitudes, and fluxes of enthalpy and kinetic energy are within $\lesssim 3σ$ of the mean of all simulations on a given grid ($128^3$ and $256^3$ grid cells), where $σ$ describes the statistical variation due to the flow's time dependence. They also agree well with a $512^3$ reference run. The $128^3$ and $256^3$ simulations agree within $9\%$ and $4\%$, respectively, on the total mass entrained into the convective layer. The entrainment rate appears to be set by the amount of energy that can be converted to work in our setup and details of the small-scale flows in the boundary layer seem to be largely irrelevant. Our results lend credence to hydrodynamic simulations of flows in stellar interiors. We provide in electronic form all outputs of our simulations as well as all information needed to reproduce or extend our study.
△ Less
Submitted 26 January, 2022; v1 submitted 1 November, 2021;
originally announced November 2021.
-
Two-dimensional simulations of solar-like models with artificially enhanced luminosity -- I. Impact on convective penetration
Authors:
I. Baraffe,
J. Pratt,
D. G. Vlaykov,
T. Guillet,
T. Goffrey,
A. Le Saux,
T. Constantino
Abstract:
We perform 2D, fully compressible, time-implicit simulations of convection in a solar-like model with the MUSIC code. Our main motivation is to explore the impact of a common tactic adopted in numerical simulations of convection that use realistic stellar conditions. This tactic is to artificially increase the luminosity and to modify the thermal diffusivity of the reference stellar model. This wo…
▽ More
We perform 2D, fully compressible, time-implicit simulations of convection in a solar-like model with the MUSIC code. Our main motivation is to explore the impact of a common tactic adopted in numerical simulations of convection that use realistic stellar conditions. This tactic is to artificially increase the luminosity and to modify the thermal diffusivity of the reference stellar model. This work focuses on the impact of these modifications on convective penetration (or overshooting) at the base of the convective envelope of a solar-like model. We explore a range of enhancement factors for the energy input and confirm the increase in the characteristic overshooting depth with the increase in the energy input. Our results highlight the importance of the impact of penetrative downflows on the thermal background below the convective boundary. This is a result of compression and shear which induce local heating and thermal mixing. The artificial increase in the energy flux intensifies the heating process by increasing the velocities in the convective zone and at the convective boundary, revealing a subtle connection between the local heating of the thermal background and the plume dynamics. This heating also increases the efficiency of heat transport by radiation which may counterbalance further heating and helps to establish a steady state. The modification of the thermal background by penetrative plumes impacts the width of the overshooting layer. Our results suggest that an artificial modification of the radiative diffusivity in the overshooting layer, rather than only accelerating the thermal relaxation, could also alter the dynamics of the penetrating plumes and thus the width of the overshooting layer. Results from simulations with an artificial modification of the energy flux and of the thermal diffusivity should be regarded with caution if used to determine an overshooting width.
△ Less
Submitted 2 September, 2021;
originally announced September 2021.
-
Suppression of lithium depletion in young low-mass stars from fast rotation
Authors:
Thomas Constantino,
Isabelle Baraffe,
Thomas Goffrey,
Jane Pratt,
Thomas Guillet,
Dimitar Vlaykov,
Louis Amard
Abstract:
We compute rotating 1D stellar evolution models that include a modified temperature gradient in convection zones and criterion for convective instability inspired by rotating 3D hydrodynamical simulations performed with the MUSIC code. In those 3D simulations we found that convective properties strongly depend on the Solberg-Høiland criterion for stability. We therefore incorporated this into 1D s…
▽ More
We compute rotating 1D stellar evolution models that include a modified temperature gradient in convection zones and criterion for convective instability inspired by rotating 3D hydrodynamical simulations performed with the MUSIC code. In those 3D simulations we found that convective properties strongly depend on the Solberg-Høiland criterion for stability. We therefore incorporated this into 1D stellar evolution models by replacing the usual Schwarzschild criterion for stability and also modifying the temperature gradient in convection zones. We computed a grid of 1D models between 0.55 and 1.2 stellar masses from the pre-main sequence to the end of main sequence in order to study the problem of lithium depletion in low-mass main sequence stars. This is an ideal test case because many of those stars are born as fast rotators and the rate of lithium depletion is very sensitive to the changes in the stellar structure. Additionally, observations show a correlation between slow rotation and lithium depletion, contrary to expectations from standard models of rotationally driven mixing. By suppressing convection, and therefore decreasing the temperature at the base of the convective envelope, lithium burning is strongly quenched in our rapidly rotating models to an extent sufficient to account for the lithium spread observed in young open clusters.
△ Less
Submitted 23 August, 2021;
originally announced August 2021.
-
Unveiling the magnetic nature of chromospheric vortices
Authors:
Mariarita Murabito,
Juie Shetye,
Marco Stangalini,
Erwin Verwichte,
Tony Arber,
Ilaria Ermolli,
Fabrizio Giorgi,
Tom Goffrey
Abstract:
Vortex structures in the Sun's chromosphere are believed to channel energy between different layers of the solar atmosphere. We investigate the nature and dynamics of two small-scale quiet-Sun rotating structures in the chromosphere. We analyse two chromospheric structures that show clear rotational patterns in spectropolarimetric observations taken with the Interferometric Bidimensional Spectrome…
▽ More
Vortex structures in the Sun's chromosphere are believed to channel energy between different layers of the solar atmosphere. We investigate the nature and dynamics of two small-scale quiet-Sun rotating structures in the chromosphere. We analyse two chromospheric structures that show clear rotational patterns in spectropolarimetric observations taken with the Interferometric Bidimensional Spectrometer (IBIS) at the Ca II 8542 Å~ line. We present the detection of spectropolarimetric signals that manifest the magnetic nature of rotating structures in the chromosphere. Our observations show two long-lived structures of plasma that each rotate clockwise inside a 10 arcsec$^{2}$~ quiet-Sun region. Their circular polarization signals are 5-10 times above the noise level. Line-of-sight Doppler velocity and horizontal velocity maps from the observations reveal clear plasma flows at and around the two structures. An MHD simulation shows these two structures are plausibly magnetically connected. Wave analysis suggests that the observed rotational vortex pattern could be due to a combination of slow actual rotation and a faster azimuthal phase speed pattern of a magneto-acoustic mode. Our results imply that the vortex structures observed in the Sun's chromosphere are magnetic in nature and that they can be connected locally through the chromosphere.
△ Less
Submitted 25 June, 2020; v1 submitted 24 June, 2020;
originally announced June 2020.
-
Comparison of two- and three-dimensional compressible convection in a pre-main sequence star
Authors:
J. Pratt,
I. Baraffe,
T. Goffrey,
C. Geroux,
T. Constantino,
D. Folini,
R. Walder
Abstract:
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 statist…
▽ More
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 star's 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.
△ Less
Submitted 10 March, 2020;
originally announced March 2020.
-
The Acceleration of Superrotation in Simulated Hot Jupiter Atmospheres
Authors:
Florian Debras,
Nathan Mayne,
Isabelle Baraffe,
Etienne Jaupart,
Pierre Mourier,
Guillaume Laibe,
Tom Goffrey,
John Thuburn
Abstract:
Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady st…
▽ More
Context. Atmospheric superrotating flows at the equator are an almost ubiquitous result of simulations of hot Jupiters, and a theory explaining how this zonally coherent flow reaches an equilibrium has been developed in the literature. However, this understanding relies on the existence of either an initial superrotating or a sheared flow, coupled with a slow evolution such that a linear steady state can be reached. Aims. A consistent physical understanding of superrotation is needed for arbitrary drag and radiative timescales, and the relevance of considering linear steady states needs to be assessed. Methods. We obtain an analytical expression for the structure, frequency and decay rate of propagating waves in hot Jupiter atmospheres around a state at rest in the 2D shallow water beta plane limit. We solve this expression numerically and confirm the robustness of our results with a 3D linear wave algorithm. We then compare with 3D simulations of hot Jupiter atmospheres and study the non linear momentum fluxes. Results. We show that under strong day night heating the dynamics does not transit through a linear steady state when starting from an initial atmosphere in solid body rotation. We further show that non linear effects favour the initial spin up of superrotation and that the acceleration due to the vertical component of the eddy momentum flux is critical to the initial development of superrotation. Conclusions. Overall, we describe the initial phases of the acceleration of superrotation, including consideration of differing radiative and drag timescales, and conclude that eddy-momentum driven superrotating equatorial jets are robust, physical phenomena in simulations of hot Jupiter atmospheres.
△ Less
Submitted 8 November, 2019;
originally announced November 2019.
-
Eigenvectors, Circulation and Linear Instabilities for Planetary Science in 3 Dimensions (ECLIPS3D)
Authors:
Florian Debras,
Nathan Mayne,
Isabelle Baraffe,
Tom Goffrey,
John Thuburn
Abstract:
Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions,…
▽ More
Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter.
△ Less
Submitted 9 September, 2019;
originally announced September 2019.
-
A well-balanced scheme for the simulation tool-kit A-MaZe: implementation, tests, and first applications to stellar structure
Authors:
M. V. Popov,
R. Walder,
D. Folini,
T. Goffrey,
I. Baraffe,
T. Constantino,
C. Geroux,
J. Pratt,
M. Viallet
Abstract:
Characterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe…
▽ More
Characterizing stellar convection in multiple dimensions is a topic at the forefront of stellar astrophysics. Numerical simulations are an essential tool for this task. We present an extension of the existing numerical tool-kit A-MaZe that enables such simulations of stratified flows in a gravitational field. The finite-volume based, cell-centered, and time-explicit hydrodynamics solver of A-MaZe was extended such that the scheme is now well-balanced in both momentum and energy. The algorithm maintains an initially static balance between gravity and pressure to machine precision. Quasi-stationary convection in slab-geometry preserves gas energy (internal plus kinetic) on average despite strong local up- and down-drafts. By contrast, a more standard numerical scheme is demonstrated to result in substantial gains of energy within a short time on purely numerical grounds. The test is further used to point out the role of dimensionality, viscosity, and Rayleigh number for compressible convection. Applications to a young sun in 2D and 3D, covering a part of the inner radiative zone as well as the outer convective zone, demonstrate that the scheme meets its initial design goal. Comparison with results obtained for a physically identical setup with a time-implicit code show qualitative agreement.
△ Less
Submitted 5 September, 2019;
originally announced September 2019.
-
Lithium depletion in solar-like stars: effect of overshooting based on realistic multi-dimensional simulations
Authors:
I. Baraffe,
J. Pratt,
T. Goffrey,
T. Constantino,
D. Folini,
M. V. Popov,
R. Walder,
M. Viallet
Abstract:
We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star…
▽ More
We study lithium depletion in low-mass and solar-like stars as a function of time, using a new diffusion coefficient describing extra-mixing taking place at the bottom of a convective envelope. This new form is motivated by multi-dimensional fully compressible, time implicit hydrodynamic simulations performed with the MUSIC code. Intermittent convective mixing at the convective boundary in a star can be modeled using extreme value theory, a statistical analysis frequently used for finance, meteorology, and environmental science. In this letter, we implement this statistical diffusion coefficient in a one-dimensional stellar evolution code, using parameters calibrated from multi-dimensional hydrodynamic simulations of a young low-mass star. We propose a new scenario that can explain observations of the surface abundance of lithium in the Sun and in clusters covering a wide range of ages, from $\sim$ 50 Myr to $\sim$ 4 Gyr. Because it relies on our physical model of convective penetration, this scenario has a limited number of assumptions. It can explain the observed trend between rotation and depletion, based on a single additional assumption, namely that rotation affects the mixing efficiency at the convective boundary. We suggest the existence of a threshold in stellar rotation rate above which rotation strongly prevents the vertical penetration of plumes and below which rotation has small effects. In addition to providing a possible explanation for the long standing problem of lithium depletion in pre-main sequence and main sequence stars, the strength of our scenario is that its basic assumptions can be tested by future hydrodynamic simulations.
△ Less
Submitted 29 July, 2017;
originally announced July 2017.
-
Extreme value statistics for two-dimensional convective penetration in a pre-Main Sequence star
Authors:
J. Pratt,
I. Baraffe,
T. Goffrey,
T. Constantino,
M. Viallet,
M. V. Popov,
R. Walder,
D. Folini
Abstract:
We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this…
▽ More
We examine a penetration layer formed between a central radiative zone and a large convection zone in the deep interior of a young low-mass star. Using the Multidimensional Stellar Implicit Code (MUSIC) to simulate two-dimensional compressible stellar convection in a spherical geometry over long times, we produce statistics that characterize the extent and impact of convective penetration in this layer. We apply extreme value theory to the maximal extent of convective penetration at any time. We compare statistical results from simulations which treat non-local convection, throughout a large portion of the stellar radius, with simulations designed to treat local convection in a small region surrounding the penetration layer. For each of these situations, we compare simulations of different resolution, which have different velocity magnitudes. We also compare statistical results between simulations that radiate energy at a constant rate to those that allow energy to radiate from the stellar surface according to the local surface temperature. Based on the frequency and depth of penetrating convective structures, we observe two distinct layers that form between the convection zone and the stable radiative zone. We show that the probability density function of the maximal depth of convective penetration at any time corresponds closely in space with the radial position where internal waves are excited. We find that the maximal penetration depth can be modeled by a Weibull distribution with a small shape parameter. Using these results, and building on established scalings for diffusion enhanced by large-scale convective motions, we propose a new form for the diffusion coefficient that may be used for one-dimensional stellar evolution calculations in the large Péclet number regime. These results should contribute to the 321D link.
△ Less
Submitted 15 June, 2017;
originally announced June 2017.
-
Benchmarking the Multi-dimensional Stellar Implicit Code MUSIC
Authors:
T. Goffrey,
J. Pratt,
M. Viallet,
I. Baraffe,
M. V. Popov,
R. Walder,
D. Folini,
C. Geroux,
T. Constantino
Abstract:
We present the results of a numerical benchmark study for the MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two- and three-dimensional compressible hydrodynamics problems relevant to stellar interiors. MUSIC is an implicit large eddy simulation code that uses implicit time integration, implemented as a Jacobian-free Newton Krylov method. A physics based preconditioning…
▽ More
We present the results of a numerical benchmark study for the MUlti-dimensional Stellar Implicit Code (MUSIC) based on widely applicable two- and three-dimensional compressible hydrodynamics problems relevant to stellar interiors. MUSIC is an implicit large eddy simulation code that uses implicit time integration, implemented as a Jacobian-free Newton Krylov method. A physics based preconditioning technique which can be adjusted to target varying physics is used to improve the performance of the solver. The problems used for this benchmark study include the Rayleigh-Taylor and Kelvin-Helmholtz instabilities, and the decay of the Taylor-Green vortex. Additionally we show a test of hydrostatic equilibrium, in a stellar environment which is dominated by radiative effects. In this setting the flexibility of the preconditioning technique is demonstrated. This work aims to bridge the gap between the hydrodynamic test problems typically used during development of numerical methods and the complex flows of stellar interiors. A series of multi-dimensional tests are performed and analysed. Each of these test cases is analysed with a simple, scalar diagnostic, with the aim of enabling direct code comparisons. As the tests performed do not have analytic solutions we verify MUSIC by comparing to established codes including ATHENA and the PENCIL code. MUSIC is able to both reproduce behaviour from established and widely-used codes as well as results expected from theoretical predictions. This benchmarking study concludes a series of papers describing the development of the MUSIC code and provides confidence in the future applications.
△ Less
Submitted 31 October, 2016;
originally announced October 2016.
-
Spherical-shell boundaries for two-dimensional compressible convection in a star
Authors:
J. Pratt,
I. Baraffe,
T. Goffrey,
C. Geroux,
M. Viallet,
D. Folini,
T. Constantino,
M. Popov,
R. Walder
Abstract:
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 MUltid…
▽ More
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.
△ Less
Submitted 23 June, 2016;
originally announced June 2016.
-
Multi-dimensional structure of accreting young stars
Authors:
C. Geroux,
I. Baraffe,
M. Viallet,
T. Goffrey,
J. Pratt,
T. Constantino,
D. Folini,
M. V. Popov,
R. Walder
Abstract:
This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young S…
▽ More
This work is the first attempt to describe the multi-dimensional structure of accreting young stars based on fully compressible time implicit multi-dimensional hydrodynamics simulations. One major motivation is to analyse the validity of accretion treatment used in previous 1D stellar evolution studies. We analyse the effect of accretion on the structure of a realistic stellar model of the young Sun. Our work is inspired by the numerical work of Kley \& Lin (1996, ApJ, 461, 933) devoted to the structure of the boundary layer in accretion disks. We analyse the redistribution of accreted material with a range of values of specific entropy relative to the bulk specific entropy of the material in the accreting object's convective envelope. A primary goal is to understand whether and how accreted energy deposited onto a stellar surface is redistributed in the interior. This study focusses on the high accretion rates characteristic of FU Ori systems. We find that the highest entropy cases produce a distinctive behaviour in the mass redistribution, rms velocities, and enthalpy flux in the convective envelope. This change in behaviour is characterised by the formation of a hot layer on the surface of the accreting object, which tends to suppress convection in the envelope. We analyse the long-term effect of such a hot buffer zone on the structure and evolution of the accreting object with 1D stellar evolution calculations. We study the relevance of the assumption of redistribution of accreted energy into the stellar interior used in the literature. One conclusion is that, for a given amount of accreted energy transferred to the accreting object, a treatment assuming accretion energy redistribution throughout the stellar interior could significantly overestimate the effects on the stellar structure, in particular, on the resulting expansion.
△ Less
Submitted 10 February, 2016;
originally announced February 2016.
-
A Jacobian-free Newton-Krylov method for time-implicit multidimensional hydrodynamics
Authors:
Maxime Viallet,
Tom Goffrey,
Isabelle Baraffe,
Doris Folini,
Chris Geroux,
Mikhail Popov,
Jane Pratt,
Rolf Walder
Abstract:
This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the invi…
▽ More
This work is a continuation of our efforts to develop an efficient implicit solver for multidimensional hydrodynamics for the purpose of studying important physical processes in stellar interiors, such as turbulent convection and overshooting. We present an implicit solver that results from the combination of a Jacobian-Free Newton-Krylov method and a preconditioning technique tailored to the inviscid, compressible equations of stellar hydrodynamics. We assess the accuracy and performance of the solver for both 2D and 3D problems for Mach numbers down to $10^{-6}$. Although our applications concern flows in stellar interiors, the method can be applied to general advection and/or diffusion-dominated flows. The method presented in this paper opens up new avenues in 3D modeling of realistic stellar interiors allowing the study of important problems in stellar structure and evolution.
△ Less
Submitted 14 December, 2015; v1 submitted 11 December, 2015;
originally announced December 2015.
