Numerical Simulations of Convective 3-Dimensional Red Supergiant Envelopes
Abstract
We explore the three-dimensional properties of convective, luminous (L≈104.5-105L), Hydrogen-rich envelopes of Red Supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using Athena++. These computations comprise ≈30\% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈ 3M of the simulated M≈15M stars. This work reveals a radius, Rcorr, around which the nature of the convection changes. For r>Rcorr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well-studied in less luminous stars, but in RSGs, the near- (or above-) Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outwards rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This "halo" of material will impact predictions for both shock breakout and early lightcurves of Type II-P supernovae. Inside of Rcorr, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length-theory (MLT). Radiation pressure provides ≈1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length of α=3-4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes (≈300 days), with amplitudes consistent with those observed photometrically.
Turn this paper into a full lesson
ArcXiv compiles a staged curriculum from this paper: 8-12 lessons across beginner → advanced, synthesised section guides, visuals, flashcards, a quiz, exercises, and on-demand deep dives per section. Grounded in the abstract, never invented.