Stellar convective penetration: parameterized theory and dynamical simulations

Abstract

Most stars host convection zones in which heat is transported directly by fluid motion, but the behavior of convective boundaries is not well understood. Here we present 3D numerical simulations which exhibit penetration zones: regions where the entire luminosity \emph{could} be carried by radiation, but where the temperature gradient is approximately adiabatic and convection is present. To parameterize this effect, we define the "penetration parameter" $\mathcal{P}$ which compares how far the radiative gradient deviates from the adiabatic gradient on either side of the Schwarzschild convective boundary. Following Roxburgh (1989) and Zahn (1991), we construct an energy-based theoretical model in which $\mathcal{P}$ controls the extent of penetration. We test this theory using 3D numerical simulations which employ a simplified Boussinesq model of stellar convection. The convection is driven by internal heating and we use a height-dependent radiative conductivity; this allows us to separately specify $\mathcal{P}$ and the stiffness $\mathcal{S}$ of the radiative-convective boundary. We find significant convective penetration in all simulations. Our simple theory describes the simulations well. Penetration zones can take thousands of overturn times to develop, so long simulations or accelerated evolutionary techniques are required. In stars, we expect $\mathcal{P}\approx 1$ and in this regime our results suggest that convection zones may extend beyond the Schwarzschild boundary by up to $\sim$20-30% of a mixing length. We present a MESA stellar model of the Sun which employs our parameterization of convective penetration as a proof of concept. We discuss prospects for extending these results to more realistic stellar contexts.

Appears as accepted for publication in ApJ. Supplemental materials available at https://doi.org/10.5281/zenodo.5131118

References in corpus (20)

• Modules for Experiments in Stellar Astrophysics (MESA)
• Modules for Experiments in Stellar Astrophysics (MESA): Pulsating Variable Stars, Rotation, Convective Boundaries, and Energy Conservation
• The chemical make-up of the Sun: A 2020 vision
• Updated Electron-Conduction Opacities: The Impact on Low-Mass Stellar Models
• Mind the gap: The location of the lower edge of the pair instability supernovae black hole mass gap
• A new Generation of Standard Solar Models
• Internal mixing of rotating stars inferred from dipole gravity modes
• On differential rotation and overshooting in solar-like stars
• Global Seismology of the Sun
• Numerical Simulations of Penetration and Overshoot in the Sun
• The Li-age correlation: the Sun is unusually Li deficient for its age
• Solar overshoot region and small-scale dynamo with realistic energy flux
• Lithium depletion and angular momentum transport in solar-type stars
• Extreme value statistics for two-dimensional convective penetration in a pre-Main Sequence star
• Convective core sizes in rotating massive stars: I. Constraints from solar metallicity OB field stars
• Calibrating Core Overshooting Parameters With Two-dimensional Hydrodynamical Simulations
• Dynamics of Mixed Convective--Stably-Stratified Fluids
• Multidimensional low-Mach number time-implicit hydrodynamic simulations of convective helium shell burning in a massive star
• The magnitude of viscous dissipation in strongly stratified two-dimensional convection
• Prandtl number dependence of compressible convection: Flow statistics and convective energy transport

Cited by in corpus (24)

• Rejuvenated accretors have less bound envelopes: Impact of Roche lobe overflow on subsequent common envelope events
• Convective boundary mixing in main-sequence stars: theory and empirical constraints
• A study of convective core overshooting as a function of stellar mass based on two-dimensional hydrodynamical simulations
• 3D hydrodynamic simulations of massive main-sequence stars. I. Dynamics and mixing of convection and internal gravity waves
• An Atlas of Convection in Main-Sequence Stars
• Schwarzschild and Ledoux are equivalent on evolutionary timescales
• Layer formation in a stably-stratified fluid cooled from above. Towards an analog for Jupiter and other gas giants
• Magnetic field breakout from white dwarf crystallization dynamos
• Is Thermohaline Mixing the Full Story? Evidence for Separate Mixing Events near the Red Giant Branch Bump
• 3D hydrodynamics simulations of internal gravity waves in red giant branch stars
• Reverse Algols and hydrogen-rich Wolf-Rayet stars from very massive binaries
• Mixed Mode Asteroseismology of Red Giant Stars Through the Luminosity Bump
• Dissecting the microphysics behind the metallicity-dependence of massive stars radii
• Euclid detectability of pair instability supernovae in binary population synthesis models consistent with merging binary black holes
• Multidimensional Simulations of Core Convection
• The landscape of binary core-collapse supernova progenitors and the late emergence of Wolf-Rayet winds
• Glitches in solar-like oscillating F-type stars: Theoretical signature of the base of the convective envelope on the ratios $r_{010}$
• Stellar age determination in the Mass-Luminosity Plane
• Implications of a turbulent convection model for classical Cepheids
• Characterising the AGB bump and its potential to constrain mixing processes in stellar interiors
• Penetrative magneto-convection of a rotating Boussinesq flow in $f$-planes
• The shape of convection in 2D and 3D global simulations of stellar interiors
• Convective Penetration in Early-Type Stars
• Discerning internal conditions of pulsating hot subdwarf B stars. Modelling multiple trapped modes in KIC 10001893