Showing 1–27 of 27 results for author: Anders, E H

Evolution of Semi-convective Staircases in Rotating Flows: Consequences for Fuzzy Cores in Giant Planets

Authors: J. R. Fuentes, Bradley W. Hindman, Adrian E. Fraser, Evan H. Anders

Abstract: Recent observational constraints on the internal structure of Jupiter and Saturn suggest that these planets have ``fuzzy" cores, i.e., gradients of the concentration of heavy elements that might span a large fraction of the planet's radius. These cores could be composed of a semi-convective staircase, i.e., multiple convective layers separated by diffusive interfaces arising from double-diffusive… ▽ More Recent observational constraints on the internal structure of Jupiter and Saturn suggest that these planets have ``fuzzy" cores, i.e., gradients of the concentration of heavy elements that might span a large fraction of the planet's radius. These cores could be composed of a semi-convective staircase, i.e., multiple convective layers separated by diffusive interfaces arising from double-diffusive instabilities. However, to date, no study has demonstrated how such staircases can avoid layer mergers and persist over evolutionary time scales. In fact, previous work has found that these mergers occur rapidly, leading to only a single convective layer. Using 3D simulations, we demonstrate that rotation prolongs the lifetime of a convective staircase by increasing the timescale for both layer merger and erosion of the interface between the final two layers. We present an analytic model for the erosion phase, predicting that rotation increases the erosion time by a factor of approximately $\mathrm{Ro}^{-1/2}$, where $\mathrm{Ro}$ is the Rossby number of the convective flows (the ratio of the rotation period to the convective turnover time). For Jovian conditions at early times after formation (when convection is vigorous enough to mix a large fraction of the planet), we find the erosion time to be roughly $10^{9}~\mathrm{yrs}$ in the non-rotating case and $10^{11}~\mathrm{yrs}$ in the rotating case. If these timescales are confirmed with a larger suite of numerical simulations, the existence of convective staircases within the deep interiors of giant planets is a strong possibility, and rotation could be an important factor in the preservation of their fuzzy cores. △ Less

Submitted 2 October, 2024; v1 submitted 20 August, 2024; originally announced August 2024.

Comments: Accepted for publication in ApJ Letters

Journal ref: ApJL 975 L1 (2024)

Iterative Methods for Navier--Stokes Inverse Problems

Authors: Liam O'Connor, Daniel Lecoanet, Evan H. Anders, Kyle C. Augustson, Keaton J. Burns, Geoffrey M. Vasil, Jeffrey S. Oishi, Benjamin P. Brown

Abstract: Even when the partial differential equation underlying a physical process can be evolved forward in time, the retrospective (backward in time) inverse problem often has its own challenges and applications. Direct Adjoint Looping (DAL) is the defacto approach for solving retrospective inverse problems, but it has not been applied to deterministic retrospective Navier--Stokes inverse problems in 2D… ▽ More Even when the partial differential equation underlying a physical process can be evolved forward in time, the retrospective (backward in time) inverse problem often has its own challenges and applications. Direct Adjoint Looping (DAL) is the defacto approach for solving retrospective inverse problems, but it has not been applied to deterministic retrospective Navier--Stokes inverse problems in 2D or 3D. In this paper, we demonstrate that DAL is ill-suited for solving retrospective 2D Navier--Stokes inverse problems. Alongside DAL, we study two other iterative methods: Simple Backward Integration (SBI) and the Quasi-Reversible Method (QRM). Our iterative SBI approach is novel while iterative QRM has previously been used. Using these three iterative methods, we solve two retrospective inverse problems: 1D Korteweg--de Vries--Burgers (decaying nonlinear wave) and 2D Navier--Stokes (unstratified Kelvin--Helmholtz vortex). In both cases, SBI and QRM reproduce the target final states more accurately and in fewer iterations than DAL. We attribute this performance gap to additional terms present in SBI and QRM's respective backward integrations which are absent in DAL. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Modelling time-dependent convective penetration in 1D stellar evolution

Authors: Cole Johnston, Mathias Michielsen, Evan H. Anders, Mathieu Renzo, Matteo Cantiello, P. Marchant, Jared A. Goldberg, Richard H. D. Townsend, Gautham Sabhahit, Adam S. Jermyn

Abstract: 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundar… ▽ More 1D stellar evolution calculations produce uncertain predictions for quantities like the age, core mass, core compactness, and nucleo-synthetic yields; a key source of uncertainty is the modeling of interfaces between regions that are convectively stable and those that are not. Theoretical and numerical work has demonstrated that there should be numerous processes adjacent to the convective boundary that induce chemical and angular momentum transport, as well as modify the thermal structure of the star. One such process is called convective penetration, wherein vigorous convection extends beyond the nominal convective boundary and alters both the composition and thermal structure. In this work, we incorporate the process of convective penetration in stellar evolution calculations using the stellar evolution software instrument mesa. We implement convective penetration according to the description presented by Anders et al. (2022a) to calculate a grid of models from the pre main sequence to He core depletion. The extent of the convective penetration zone is self-consistently calculated at each time step without introducing new free parameters. We find both a substantial penetration zone in all models with a convective core and observable differences to global stellar properties such as the luminosity and radius. We preset how the predicted radial extent of the penetration zone scales with the total stellar mass, age and the metallicity of the star. We discuss our results in the context of existing numerical and observational studies. △ Less

Submitted 13 December, 2023; originally announced December 2023.

Comments: Accepted for publication in The Astrophysical Journal, 16 pages, 8 figures, accompanying repository (https://github.com/mathren/dbconvpen)

Internally heated and fully compressible convection: flow morphology and scaling laws

Authors: Whitney T. Powers, Evan H. Anders, Benjamin P. Brown

Abstract: In stars and planets natural processes heat convective flows in the bulk of a convective region rather than at hard boundaries. By characterizing how convective dynamics are determined by the strength of an internal heating source we can gain insight into the processes driving astrophysical convection. Internally heated convection has been studied extensively in incompressible fluids, but the effe… ▽ More In stars and planets natural processes heat convective flows in the bulk of a convective region rather than at hard boundaries. By characterizing how convective dynamics are determined by the strength of an internal heating source we can gain insight into the processes driving astrophysical convection. Internally heated convection has been studied extensively in incompressible fluids, but the effects of stratification and compressibility have not been examined in detail. In this work, we study fully compressible convection driven by a spatially uniform heating source in 2D and 3D Cartesian, hydrodynamic simulations. We use a fixed temperature upper boundary condition which results in a system that is internally heated in the bulk and cooled at the top. We find that the flow speed, as measured by the Mach number, and turbulence, as measured by the Reynolds number, can be independently controlled by separately varying the characteristic temperature gradient from internal heating and the diffusivities. 2D simulations at a fixed Mach number (flow speed) demonstrate consistent power at low wavenumber as diffusivities are decreased. We observe convection where the velocity distribution is skewed towards cold, fast downflows, and that the flow speed is related to the length scale and entropy gradient of the upper boundary where the downflows are driven. We additionally find a heat transport scaling law which is consistent with prior incompressible work. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Comments: 22 pages, 12 figures, submitted to Phys. Rev. Fluids

Journal ref: Phys. Rev. Fluids 9 (2024), 043501

Nonlinear mode coupling and energetics of driven magnetized shear-flow turbulence

Authors: B. Tripathi, A. E. Fraser, P. W. Terry, E. G. Zweibel, M. J. Pueschel, E. H. Anders

Abstract: To comprehensively understand saturation of two-dimensional ($2$D) magnetized Kelvin-Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a signif… ▽ More To comprehensively understand saturation of two-dimensional ($2$D) magnetized Kelvin-Helmholtz-instability-driven turbulence, energy transfer analysis is extended from the traditional interaction between scales to include eigenmode interactions, by using the nonlinear couplings of linear eigenmodes of the ideal instability. While both kinetic and magnetic energies cascade to small scales, a significant fraction of turbulent energy deposited by unstable modes in the fluctuation spectrum is shown to be re-routed to the conjugate-stable modes at the instability scale. They remove energy from the forward cascade at its inception. The remaining cascading energy flux is shown to attenuate exponentially at a small scale, dictated by the large-scale stable modes. Guided by a widely used instability-saturation assumption, a general quasilinear model of instability is tested by retaining all nonlinear interactions except those that couple to the large-scale stable modes. These complex interactions are analytically removed from the magnetohydrodynamic equations using a novel technique. Observations are: an explosive large-scale vortex separation instead of the well-known merger of $2$D, a dramatic enhancement in turbulence level and spectral energy fluxes, and a reduced small-scale dissipation length-scale. These show critical role of the stable modes in instability saturation. Possible reduced-order turbulence models are proposed for fusion and astrophysical plasmas, based on eigenmode-expanded energy transfer analyses. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: Selected by the editors of Physics of Plasmas as a Featured article. https://doi.org/10.1063/5.0156560

Journal ref: Physics of Plasmas 30, 072107 (2023)

The photometric variability of massive stars due to gravity waves excited by core convection

Authors: Evan H. Anders, Daniel Lecoanet, Matteo Cantiello, Keaton J. Burns, Benjamin A. Hyatt, Emma Kaufman, Richard H. D. Townsend, Benjamin P. Brown, Geoffrey M. Vasil, Jeffrey S. Oishi, Adam S. Jermyn

Abstract: Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometri… ▽ More Massive stars die in catastrophic explosions, which seed the interstellar medium with heavy elements and produce neutron stars and black holes. Predictions of the explosion's character and the remnant mass depend on models of the star's evolutionary history. Models of massive star interiors can be empirically constrained by asteroseismic observations of gravity wave oscillations. Recent photometric observations reveal a ubiquitous red noise signal on massive main sequence stars; a hypothesized source of this noise is gravity waves driven by core convection. We present the first 3D simulations of massive star convection extending from the star's center to near its surface, with realistic stellar luminosities. Using these simulations, we make the first prediction of photometric variability due to convectively-driven gravity waves at the surfaces of massive stars, and find that gravity waves produce photometric variability of a lower amplitude and lower characteristic frequency than the observed red noise. We infer that the photometric signal of gravity waves excited by core convection is below the noise limit of current observations, so the red noise must be generated by an alternative process. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: As accepted for publication in Nature Astronomy except for final editorial revisions. Supplemental materials available online at https://doi.org/10.5281/zenodo.7764997 . We have also sonified our results to make them more accessible, see https://github.com/evanhanders/gmode_variability_paper/blob/main/sound/gmode_sonification.pdf

Journal ref: Nature Astronomy, 2023, 10.1038/s41550-023-02040-7

Rotation reduces convective mixing in Jupiter and other gas giants

Authors: J. R. Fuentes, Evan H. Anders, Andrew Cumming, Bradley W. Hindman

Abstract: Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers.… ▽ More Recent measurements of Jupiter's gravitational moments by the Juno spacecraft and seismology of Saturn's rings suggest that the primordial composition gradients in the deep interior of these planets have persisted since their formation. One possible explanation is the presence of a double-diffusive staircase below the planet's outer convection zone, which inhibits mixing across the deeper layers. However, hydrodynamic simulations have shown that these staircases are not long-lasting and can be disrupted by overshooting convection. In this paper we suggests that planetary rotation could be another factor for the longevity of primordial composition gradients. Using rotational mixing-length theory and 3D hydrodynamic simulations, we demonstrate that rotation significantly reduces both the convective velocity and the mixing of primordial composition gradients. In particular, for Jovian conditions at $t\sim10^{8}~\mathrm{yrs}$ after formation, rotation reduces the convective velocity by a factor of 6, and in turn, the kinetic energy flux available for mixing gets reduced by a factor of $6^3\sim 200$. This leads to an entrainment timescale that is more than two orders of magnitude longer than without rotation. We encourage future hydrodynamic models of Jupiter and other gas giants to include rapid rotation, because the decrease in the mixing efficiency could explain why Jupiter and Saturn are not fully mixed. △ Less

Submitted 22 May, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: Accepted for publication in the Astrophysical Journal Letters

Journal ref: ApJL 950 L4 (2023)

Convective boundary mixing in main-sequence stars: theory and empirical constraints

Authors: Evan H. Anders, May G. Pedersen

Abstract: The convective envelopes of solar-type stars and the convective cores of intermediate- and high-mass stars share boundaries with stable radiative zones. Through a host of processes we collectively refer to as "convective boundary mixing" (CBM), convection can drive efficient mixing in these nominally stable regions. In this review, we discuss the current state of CBM research in the context of mai… ▽ More The convective envelopes of solar-type stars and the convective cores of intermediate- and high-mass stars share boundaries with stable radiative zones. Through a host of processes we collectively refer to as "convective boundary mixing" (CBM), convection can drive efficient mixing in these nominally stable regions. In this review, we discuss the current state of CBM research in the context of main-sequence stars through three lenses. (1) We examine the most frequently implemented 1D prescriptions of CBM -- exponential overshoot, step overshoot, and convective penetration -- and we include a discussion of implementation degeneracies and how to convert between various prescriptions. (2) Next, we examine the literature of CBM from a fluid dynamical perspective, with a focus on three distinct processes: convective overshoot, entrainment, and convective penetration. (3) Finally, we discuss observational inferences regarding how much mixing should occur in the cores of intermediate- and high-mass stars, and the implied constraints that these observations place on 1D CBM implementations. We conclude with a discussion of pathways forward for future studies to place better constraints on this difficult challenge in stellar evolution modeling. △ Less

Submitted 15 April, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: Appears roughly as published; invited review in the Galaxies special issue, "The Structure and Evolution of Stars". Supplementary materials can be found online in a Zenodo repository at https://zenodo.org/record/7735067

Journal ref: Galaxies 2023, 11(2), 56

Heat transport and convective velocities in compositionally-driven convection in neutron star and white dwarf interiors

Authors: J. R. Fuentes, A. Cumming, M. Castro-Tapia, E. H. Anders

Abstract: We investigate heat transport associated with compositionally-driven convection driven by crystallization at the ocean-crust interface in accreting neutron stars, or growth of the solid core in cooling white dwarfs. We study the effect of thermal diffusion and rapid rotation on the convective heat transport, using both mixing length theory and numerical simulations of Boussinesq convection. We det… ▽ More We investigate heat transport associated with compositionally-driven convection driven by crystallization at the ocean-crust interface in accreting neutron stars, or growth of the solid core in cooling white dwarfs. We study the effect of thermal diffusion and rapid rotation on the convective heat transport, using both mixing length theory and numerical simulations of Boussinesq convection. We determine the heat flux, composition gradient and Péclet number, $\mathrm{Pe}$ (the ratio of thermal diffusion time to convective turnover time) as a function of the composition flux. We find two regimes of convection with a rapid transition between them as the composition flux increases. At small Pe, the ratio between the heat flux and composition flux is independent of Pe,, because the loss of heat from convecting fluid elements due to thermal diffusion is offset by the smaller composition gradient needed to overcome the reduced thermal buoyancy. At large Pe, the temperature gradient approaches the adiabatic gradient, saturating the heat flux. We discuss the implications for neutron star and white dwarf cooling. Convection in neutron stars spans both regimes. We find rapid mixing of neutron star oceans, with a convective turnover time of order weeks to minutes depending on rotation. Except during the early stages of core crystallization, white dwarf convection is in the thermal-diffusion-dominated fingering regime. We find convective velocities much smaller than recent estimates for crystallization-driven dynamos. The small fraction of energy carried as kinetic energy calls into question the effectiveness of crystallization-driven dynamos as an explanation for observed white dwarf magnetic fields. △ Less

Submitted 6 April, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Comments: Accepted for Publication in ApJ (Astrophysical Journal)

Journal ref: ApJ 950 73 (2023)

The Stability of Prendergast Magnetic Fields

Authors: Emma Kaufman, Daniel Lecoanet, Evan H. Anders, Benjamin P. Brown, Geoffrey M. Vasil, Jeffrey S. Oishi, Keaton J. Burns

Abstract: Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any eviden… ▽ More Convection in massive main sequence stars generates large scale magnetic fields in their cores which persists as they evolve up the red giant branch. The remnants of these fields may take the form of the Prendergast magnetic field, a combination of poloidal and toroidal field components which are expected to stabilize each other. Previous analytic and numerical calculations did not find any evidence for instability of the Prendergast field over short timescales. In this paper, we present numerical simulations which show a long timescale, linear instability of this magnetic field. We find the instability to be robust to changes in boundary conditions and it is not stabilized by strong stable stratification. The instability is a resistive instability, and the growth rate has a power-law dependence on the resistivity, in which the growth rate decreases as the resistivity decreases. We estimate the growth rate of the instability in stars by extrapolating this power-law to stellar values of the resistivity. The instability is sufficiently rapid to destabilize the magnetic field on timescales shorter than the stellar evolution timescale, indicating that the Prendergast field is not a good model to use in studies of magnetic fields in stars. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 9 pages, 8 figures

An Atlas of Convection in Main-Sequence Stars

Authors: Adam S. Jermyn, Evan H. Anders, Daniel Lecoanet, Matteo Cantiello

Abstract: Convection is ubiquitous in stars and occurs under many different conditions. Here we explore convection in main-sequence stars through two lenses: dimensionless parameters arising from stellar structure and parameters which emerge from the application of mixing length theory. We first define each quantity in terms familiar both to the 1D stellar evolution community and the hydrodynamics community… ▽ More Convection is ubiquitous in stars and occurs under many different conditions. Here we explore convection in main-sequence stars through two lenses: dimensionless parameters arising from stellar structure and parameters which emerge from the application of mixing length theory. We first define each quantity in terms familiar both to the 1D stellar evolution community and the hydrodynamics community. We then explore the variation of these quantities across different convection zones, different masses, and different stages of main-sequence evolution. We find immense diversity across stellar convection zones. Convection occurs in thin shells, deep envelopes, and nearly-spherical cores; it can be efficient of inefficient, rotationally constrained or not, transsonic or deeply subsonic. This atlas serves as a guide for future theoretical and observational investigations by indicating which regimes of convection are active in a given star, and by describing appropriate model assumptions for numerical simulations. △ Less

Submitted 31 May, 2022; originally announced June 2022.

Comments: 62 pages, 76 figures. All your favorite dimensionless numbers and then some. Submitted to ApJS

Layer formation in a stably-stratified fluid cooled from above. Towards an analog for Jupiter and other gas giants

Authors: J. R. Fuentes, A. Cumming, E. H. Anders

Abstract: In 1D evolution models of gas giant planets, an outer convection zone advances into the interior as the surface cools, and multiple convective layers form beneath that convective front. To study layer formation below an outer convection zone in a similar scenario, we investigate the evolution of a stably-stratified fluid with a linear composition gradient that is constantly being cooled from above… ▽ More In 1D evolution models of gas giant planets, an outer convection zone advances into the interior as the surface cools, and multiple convective layers form beneath that convective front. To study layer formation below an outer convection zone in a similar scenario, we investigate the evolution of a stably-stratified fluid with a linear composition gradient that is constantly being cooled from above. We use the Boussinesq approximation in a series of 2D simulations at low and high Prandtl numbers ($\mathrm{Pr} = 0.5$ and 7), initialized with constant temperature everywhere, and cooled at different rates. We find that multiple convective layers form at $\mathrm{Pr} = 7$, {as the result of an instability in the} diffusive thermal boundary layer below the outer convection zone. At low Pr, layers do not form because the temperature gradient within the boundary layer is much smaller than at large Pr and, consequently, is not large enough to overcome the stabilizing effect of the composition gradient. For the stratification used in this study, on the long-term the composition gradient is an ineffective barrier against the propagation of the outer convection zone and the entire fluid becomes fully-mixed, whether layers form or not. Our results challenge 1D evolutionary models of gas giant planets, which predict that layers are long-lived and that the outer convective envelope stops advancing inwards. We discuss what is needed for future work to build more realistic models. △ Less

Submitted 28 November, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Accepted for publication in Physical review fluids

Journal ref: Phys. Rev. Fluids 7, 124501, 2022 (Editor's suggestion)

Characterizing Observed Extra Mixing Trends in Red Giants using the Reduced Density Ratio from Thermohaline Models

Authors: Adrian E. Fraser, Meridith Joyce, Evan H. Anders, Jamie Tayar, Matteo Cantiello

Abstract: Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we p… ▽ More Observations show an almost ubiquitous presence of extra mixing in low-mass upper giant branch stars. The most commonly invoked explanation for this is thermohaline mixing. One-dimensional stellar evolution models include various prescriptions for thermohaline mixing, but the use of observational data directly to discriminate between thermohaline prescriptions has thus far been limited. Here, we propose a new framework to facilitate direct comparison: Using carbon-to-nitrogen measurements from the SDSS-IV APOGEE survey as a probe of mixing and a fluid parameter known as the reduced density ratio from one-dimensional stellar evolution programs, we compare the observed amount of extra mixing on the upper giant branch to predicted trends from three-dimensional fluid dynamics simulations. Using this method, we are able to empirically constrain how mixing efficiency should vary with the reduced density ratio. We find the observed amount of extra mixing is strongly correlated with the reduced density ratio and that trends between reduced density ratio and fundamental stellar parameters are robust across choices for modeling prescription. We show that stars with available mixing data tend to have relatively low density ratios, which should inform the regimes selected for future simulation efforts. Finally, we show that there is increased mixing at low reduced density ratios, which is consistent with current hydrodynamical models of thermohaline mixing. The introduction of this framework sets a new standard for theoretical modeling efforts, as validation for not only the amount of extra mixing, but trends between the degree of extra mixing and fundamental stellar parameters is now possible. △ Less

Submitted 11 October, 2022; v1 submitted 18 April, 2022; originally announced April 2022.

Comments: 22 pages, 7 figures, accepted for publication in ApJ. MESA histories available on Zenodo at https://doi.org/10.5281/zenodo.7152746 . MESA inlists available on GitHub at https://github.com/afraser3/Empirical-Magnetic-Thermohaline

Moosinesq Convection in the Cores of Moosive Stars

Authors: Evan H. Anders, Evan B. Bauer, Adam S. Jermyn, Samuel J. Van Kooten, Benjamin P. Brown, Eric W. Hester, Mindy Wilkinson, Jared A. Goldberg, Tania Varesano, Daniel Lecoanet

Abstract: Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure… ▽ More Stars with masses $\gtrsim 4 \times 10^{27}M_{\rm{moose}} \approx 1.1 M_\odot$ have core convection zones during their time on the main sequence. In these moosive stars, convection introduces many uncertainties in stellar modeling. In this Letter, we build upon the Boussinesq approximation to present the first-ever simulations of Moosinesq convection, which captures the complex geometric structure of the convection zones of these stars. These flows are bounded in a manner informed by the majestic terrestrial Alces alces (moose) and could have important consequences for the evolution of these stars. We find that Moosinesq convection results in very interesting flow morphologies and rapid heat transfer, and posit this as a mechanism of biomechanical thermoregulation. △ Less

Submitted 30 March, 2022; originally announced April 2022.

Comments: 7 pages, 3 figures, 1 moose, no logic in this place

Convective Penetration in Early-Type Stars

Authors: Adam S. Jermyn, Evan H. Anders, Daniel Lecoanet, Matteo Cantiello

Abstract: Observations indicate that the convective cores of stars must ingest a substantial amount of material from the overlying radiative zone, but the extent of this mixing and the mechanism which causes it remain uncertain. Recently, Anders et al. (2021) developed a theory of convective penetration and calibrated it with 3D numerical hydrodynamics simulations. Here we employ that theory to predict the… ▽ More Observations indicate that the convective cores of stars must ingest a substantial amount of material from the overlying radiative zone, but the extent of this mixing and the mechanism which causes it remain uncertain. Recently, Anders et al. (2021) developed a theory of convective penetration and calibrated it with 3D numerical hydrodynamics simulations. Here we employ that theory to predict the extent of convective boundary mixing in early-type main-sequence stars. We find that convective penetration produces enough mixing to explain core masses inferred from asteroseismology and eclipsing binary studies, and matches observed trends in mass and age. While there are remaining uncertainties in the theory, this agreement suggests that most convective boundary mixing in early-type main-sequence stars arises from convective penetration. Finally, we provide a fitting formula for the extent of core convective penetration for main-sequence stars in the mass range from $1.1-60 M_\odot$. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: 10 pages, 8 figures. Accepted in ApJ

Schwarzschild and Ledoux are equivalent on evolutionary timescales

Authors: Evan H. Anders, Adam S. Jermyn, Daniel Lecoanet, Adrian E. Fraser, Imogen G. Cresswell, Meridith Joyce, J. R. Fuentes

Abstract: Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are co… ▽ More Stellar evolution models calculate convective boundaries using either the Schwarzschild or Ledoux criterion, but confusion remains regarding which criterion to use. Here we present a 3D hydrodynamical simulation of a convection zone and adjacent radiative zone, including both thermal and compositional buoyancy forces. As expected, regions which are unstable according to the Ledoux criterion are convective. Initially, the radiative zone adjacent to the convection zone is Schwarzschild-unstable but Ledoux-stable due to a composition gradient. Over many convective overturn timescales the convection zone grows via entrainment. The convection zone saturates at the size originally predicted by the Schwarzschild criterion, although in this final state the Schwarzschild and Ledoux criteria agree. Therefore, the Schwarzschild criterion should be used to determine the size of stellar convection zones, except possibly during short-lived evolutionary stages in which entrainment persists. △ Less

Submitted 11 March, 2022; originally announced March 2022.

Comments: Manuscript accepted in ApJ Letters. Supplemental materials are in a Zenodo repository (https://zenodo.org/record/6344868). Code is in a Github repository (https://github.com/evanhanders/schwarzschild_or_ledoux)

A Transparent Window into Early-Type Stellar Variability

Authors: Adam S. Jermyn, Evan H. Anders, Matteo Cantiello

Abstract: Subsurface convection zones are ubiquitous in early-type stars. Driven by narrow opacity peaks, these thin convective regions transport little heat but play an important role in setting the magnetic properties and surface variability of stars. Here we demonstrate that these convection zones are \emph{not} present in as wide a range of stars as previously believed. In particular, there are regions… ▽ More Subsurface convection zones are ubiquitous in early-type stars. Driven by narrow opacity peaks, these thin convective regions transport little heat but play an important role in setting the magnetic properties and surface variability of stars. Here we demonstrate that these convection zones are \emph{not} present in as wide a range of stars as previously believed. In particular, there are regions which 1D stellar evolution models report to be convectively unstable but which fall below the critical Rayleigh number for onset of convection. For sub-solar metallicity this opens up a \emph{stability window} in which there are no subsurface convection zones. For LMC metallicity this surface stability region extends roughly between $8M_\odot$ and $16M_\odot$, increasing to $8M_\odot$ -- $35M_\odot$ for SMC metallicity. Such windows are then an excellent target for probing the relative influence of subsurface convection and other sources of photometric variability in massive stars. △ Less

Submitted 25 January, 2022; originally announced January 2022.

Comments: 9 pages, 4 figures. Accepted in ApJ

Stellar convective penetration: parameterized theory and dynamical simulations

Authors: Evan H. Anders, Adam S. Jermyn, Daniel Lecoanet, Benjamin P. Brown

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 parame… ▽ More 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. △ Less

Submitted 11 March, 2022; v1 submitted 21 October, 2021; originally announced October 2021.

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

Journal ref: ApJ 926 169 (2022)

Marginally-Stable Thermal Equilibria of Rayleigh-Bénard Convection

Authors: Liam O'Connor, Daniel Lecoanet, Evan H. Anders

Abstract: Natural convection is ubiquitous throughout the physical sciences and engineering, yet many of its important properties remain elusive. To study convection in a novel context, we derive and solve a quasilinear form of the Rayleigh-Bénard problem by representing the perturbations in terms of marginally-stable eigenmodes. The amplitude of each eigenmode is determined by requiring that the background… ▽ More Natural convection is ubiquitous throughout the physical sciences and engineering, yet many of its important properties remain elusive. To study convection in a novel context, we derive and solve a quasilinear form of the Rayleigh-Bénard problem by representing the perturbations in terms of marginally-stable eigenmodes. The amplitude of each eigenmode is determined by requiring that the background state maintains marginal stability. The background temperature profile evolves due to the advective flux of every marginally-stable eigenmode, as well as diffusion. To ensure marginal stability and to obtain the eigenfunctions at every timestep, we perform a one-dimensional eigenvalue solve on each of the allowable wavenumbers. The background temperature field evolves to an equilibrium state, where the advective flux from the marginally-stable eigenmodes and the diffusive flux sum to a constant. These marginally-stable thermal equilibria (MSTE) are exact solutions of the quasilinear equations. The mean temperature profile has thinner boundary layers and larger Nusselt numbers than thermally-equilibrated 2D and 3D simulations of the full nonlinear equations. We find the Nusselt number scales like $\rm{Nu} \sim\rm{Ra}^{1/3}$. When an MSTE is used as initial conditions for a 2D simulation, we find that Nu quickly equilibrates without the burst of turbulence often induced by purely conductive initial conditions, but we also find that the kinetic energy is too large and viscously attenuates on a long viscous time scale. This is due to the thin temperature boundary layers which diffuse heat very effectively, thereby requiring high-velocity advective flows to reach an equilibrium. △ Less

Submitted 27 February, 2024; v1 submitted 28 May, 2021; originally announced May 2021.

Journal ref: Phys. Rev. Fluids 6, 093501 -- Published 21 September 2021

Surface Manifestation of Stochastically Excited Internal Gravity Waves

Authors: Daniel Lecoanet, Matteo Cantiello, Evan H. Anders, Eliot Quataert, Louis-Alexandre Couston, Mathieu Bouffard, Benjamin Favier, Michael Le Bars

Abstract: Recent photometric observations of massive stars show ubiquitous low-frequency "red-noise" variability, which has been interpreted as internal gravity waves (IGWs). Simulations of IGWs generated by convection show smooth surface wave spectra, qualitatively matching the observed red-noise. On the other hand, theoretical calculations by Shiode et al (2013) and Lecoanet et al (2019) predict IGWs shou… ▽ More Recent photometric observations of massive stars show ubiquitous low-frequency "red-noise" variability, which has been interpreted as internal gravity waves (IGWs). Simulations of IGWs generated by convection show smooth surface wave spectra, qualitatively matching the observed red-noise. On the other hand, theoretical calculations by Shiode et al (2013) and Lecoanet et al (2019) predict IGWs should manifest at the surface as regularly-spaced peaks associated with standing g-modes. In this work, we compare these theoretical approaches to simplified 2D numerical simulations. The simulations show g-mode peaks at their surface, and are in good agreement with Lecoanet et al (2019). The amplitude estimates of Shiode et al (2013) did not take into account the finite width of the g-mode peaks; after correcting for this finite width, we find good agreement with simulations. However, simulations need to be run for hundreds of convection turnover times for the peaks to become visible; this is a long time to run a simulation, but a short time in the life of a star. The final spectrum can be predicted by calculating the wave energy flux spectrum in much shorter simulations, and then either applying the theory of Shiode et al (2013) or Lecoanet et al (2019). △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: submitted to MNRAS

A Refined Model of Convectively-Driven Flicker in Kepler Light Curves

Authors: Samuel J. Van Kooten, Evan H. Anders, Steven R. Cranmer

Abstract: Light curves produced by the Kepler mission demonstrate stochastic brightness fluctuations (or "flicker") of stellar origin which contribute to the noise floor, limiting the sensitivity of exoplanet detection and characterization methods. In stars with surface convection, the primary driver of these variations on short (sub-eight-hour) timescales is believed to be convective granulation. In this w… ▽ More Light curves produced by the Kepler mission demonstrate stochastic brightness fluctuations (or "flicker") of stellar origin which contribute to the noise floor, limiting the sensitivity of exoplanet detection and characterization methods. In stars with surface convection, the primary driver of these variations on short (sub-eight-hour) timescales is believed to be convective granulation. In this work, we improve existing models of this granular flicker amplitude, or $F_8$, by including the effect of the Kepler bandpass on measured flicker, by incorporating metallicity in determining convective Mach numbers, and by using scaling relations from a wider set of numerical simulations. To motivate and validate these changes, we use a recent database of convective flicker measurements in Kepler stars, which allows us to more fully detail the remaining model--prediction error. Our model improvements reduce the typical misprediction of flicker amplitude from a factor of 2.5 to 2. We rule out rotation period and strong magnetic activity as possible explanations for the remaining model error, and we show that binary companions may affect convective flicker. We also introduce an "envelope" model which predicts a range of flicker amplitudes for any one star to account for some of the spread in numerical simulations, and we find that this range covers 78% of observed stars. We note that the solar granular flicker amplitude is lower than most Sun-like stars. This improved model of convective flicker amplitude can better characterize this source of noise in exoplanet studies as well as better inform models and simulations of stellar granulation. △ Less

Submitted 13 April, 2021; originally announced April 2021.

Comments: Accepted for publication in the Astrophysical Journal. 23 pages, 12 figures. Code and data available at https://doi.org/10.5281/zenodo.4444282

Convective dynamics with mixed temperature boundary conditions: why thermal relaxation matters and how to accelerate it

Authors: Evan H. Anders, Geoffrey M. Vasil, Benjamin P. Brown, Lydia Korre

Abstract: Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary a… ▽ More Astrophysical simulations of convection frequently impose different thermal boundary conditions at the top and the bottom of the domain in an effort to more accurately model natural systems. In this work, we study Rayleigh-Benard convection (RBC) under the Boussinesq approximation. We examine simulations with mixed temperature boundary conditions in which the flux is fixed at the bottom boundary and the temperature is fixed at the top ("FT"). We aim to understand how FT boundaries change the nature of the convective solution compared to the traditional choice of thermal boundaries, in which the temperature is fixed at the top and bottom of the domain ("TT"). We demonstrate that the timescale of thermal relaxation for FT simulations is dependent upon the initial conditions. "Classic" initial conditions which employ a hydrostatically -- and thermally -- balanced linear temperature profile exhibit a long thermal relaxation. This long relaxation is not seen in FT simulations which use a TT simulation's nonlinear state as initial conditions ("TT-to-FT"). In the thermally relaxed, statistically stationary state, the mean behavior of an FT simulation corresponds to an equivalent simulation with TT boundaries, and time- and volume-averaged flow statistics like the Nusselt number and the Peclet number are indistinguishable between FT and TT simulations. FT boundaries are fundamentally asymmetric, and we examine the asymmetries that these boundaries produce in the flow. We find that the fixed-flux boundary produces more extreme temperature events than the fixed-temperature boundary. However, these near-boundary asymmetries do not measurably break the symmetry in the convective interior. We briefly explore rotating RBC to demonstrate that our findings with respect to thermal relaxation carry over to this more complex case, and to show the power of TT-to-FT initial conditions. △ Less

Submitted 12 September, 2020; v1 submitted 28 February, 2020; originally announced March 2020.

Comments: Published open access in Physical Review Fluids. The arXiv v3 appears roughly as published, but the published version is better proofread

Journal ref: Phys. Rev. Fluids 5, 083501 (2020)

Entropy Rain: Dilution and Compression of Thermals in Stratified Domains

Authors: Evan H. Anders, Daniel Lecoanet, Benjamin P. Brown

Abstract: Large-scale convective flows called giant cells were once thought to transport the Sun's luminosity in the solar convection zone, but recent observations have called their existence into question. In place of large-scale flows, some authors have suggested the solar luminosity may instead be transported by small droplets of rapidly falling, low entropy fluid. This "entropy rain" could propagate as… ▽ More Large-scale convective flows called giant cells were once thought to transport the Sun's luminosity in the solar convection zone, but recent observations have called their existence into question. In place of large-scale flows, some authors have suggested the solar luminosity may instead be transported by small droplets of rapidly falling, low entropy fluid. This "entropy rain" could propagate as dense vortex rings, analogous to rising buoyant thermals in the Earth's atmosphere. In this work, we develop an analytical theory describing the evolution of dense, negatively buoyant thermals. We verify the theory with 2D cylindrical and 3D cartesian simulations of laminar, axisymmetric thermals in highly stratified atmospheres. Our results show that dense thermals fall in two categories: a stalling regime in which the droplets slow down and expand, and a falling regime in which the droplets accelerate and shrink as they propagate downwards. We estimate that solar downflows are in the falling regime and maintain their entropy perturbation against diffusion until they reach the base of the convection zone. This suggests that entropy rain may be an effective nonlocal mechanism for transporting the solar luminosity. △ Less

Submitted 16 October, 2019; v1 submitted 5 June, 2019; originally announced June 2019.

Comments: As published in ApJ. 12 pages, including 6 figures, 2 tables, and 3 appendices

Predicting the Rossby number in convective experiments

Authors: Evan H. Anders, Cathryn M. Manduca, Benjamin P. Brown, Jeffrey S. Oishi, Geoffrey M. Vasil

Abstract: The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of… ▽ More The Rossby number is a crucial parameter describing the degree of rotational constraint on the convective dynamics in stars and planets. However, it is not an input to computational models of convection but must be measured ex post facto. Here, we report the discovery of a new quantity, the Predictive Rossby number, which is both tightly correlated with the Rossby number and specified in terms of common inputs to numerical models. The Predictive Rossby number can be specified independent of Rayleigh number, allowing suites of numerical solutions to separate the degree of rotational constraint from the strength of the driving of convection. We examine the scaling of convective transport in terms of the Nusselt number and the degree of turbulence in terms of the Reynolds number of the flow, and we find scaling laws nearly identical to those in nonrotational convection at low Rossby number when the Predictive Rossby number is held constant. Finally, we describe the boundary layers as a function of increasing turbulence at constant Rossby number. △ Less

Submitted 15 January, 2019; v1 submitted 11 December, 2018; originally announced December 2018.

Comments: As accepted in ApJ. 9 pages (including references). 4 figures, 2 tables

Journal ref: The Astrophysical Journal. Volume 872, Number 2 (2019)

Accelerated evolution of convective simulations

Authors: Evan H. Anders, Benjamin P. Brown, Jeffrey S. Oishi

Abstract: High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in… ▽ More High Peclet number, turbulent convection is a classic system with a large timescale separation between flow speeds and the thermal relaxation time. In this paper, we present a method of fast-forwarding through the long thermal relaxation of convective simulations, and we test the validity of this method. This accelerated evolution (AE) method involves measuring the dynamics of convection early in a simulation and using its characteristics to adjust the mean thermodynamic profile within the domain towards its evolved state. We study Rayleigh-Bénard convection as a test case for AE. Evolved flow properties of AE solutions are measured to be within a few percent of solutions which are reached through standard evolution (SE) over a full thermal timescale. At the highest values of the Rayleigh number at which we compare SE and AE, we find that AE solutions require roughly an order of magnitude fewer computing hours to evolve than SE solutions. △ Less

Submitted 21 August, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: As published in PRFluids. 14 pages total with 6 figures, 2 tables, and 2 appendices

Journal ref: Phys. Rev. Fluids 3, 083502 (2018)

Convective heat transport in stratified atmospheres at low and high Mach number

Authors: Evan H. Anders, Benjamin P. Brown

Abstract: We study fully compressible convection in the context of plane-parallel, polytropically stratified atmospheres. We perform a suite of 2D and 3D simulations in which we vary the initial superadiabaticity ($ε$) and the Rayleigh number (Ra) while fixing the initial density stratification, aspect ratio, and Prandtl number. The evolved heat transport, quantified by the Nusselt number (Nu), follows scal… ▽ More We study fully compressible convection in the context of plane-parallel, polytropically stratified atmospheres. We perform a suite of 2D and 3D simulations in which we vary the initial superadiabaticity ($ε$) and the Rayleigh number (Ra) while fixing the initial density stratification, aspect ratio, and Prandtl number. The evolved heat transport, quantified by the Nusselt number (Nu), follows scaling relationships similar to those found in the well-studied, incompressible Rayleigh-Bénard problem. This scaling holds up in both 2D and 3D and is not appreciably affected by the magnitude of $ε$. △ Less

Submitted 29 August, 2017; v1 submitted 20 November, 2016; originally announced November 2016.

Comments: 8 pages, 4 figures, as published in PRFluids

Journal ref: Phys. Rev. Fluids 2, 083501 (2017)

The Advanced LIGO Photon Calibrators

Authors: S. Karki, D. Tuyenbayev, S. Kandhasamy, B. P. Abbott, T. D. Abbott, E. H. Anders, J. Berliner, J. Betzwieser, H. P. Daveloza, C. Cahillane, L. Canete, C. Conley, J. R. Gleason, E. Goetz, J. S. Kissel, K. Izumi, G. Mendell, V. Quetschke, M. Rodruck, S. Sachdev, T. Sadecki, P. B. Schwinberg, A. Sottile, M. Wade, A. J. Weinstein , et al. (2 additional authors not shown)

Abstract: The two interferometers of the Laser Interferometry Gravitaional-wave Observatory (LIGO) recently detected gravitational waves from the mergers of binary black hole systems. Accurate calibration of the output of these detectors was crucial for the observation of these events, and the extraction of parameters of the sources. The principal tools used to calibrate the responses of the second-generati… ▽ More The two interferometers of the Laser Interferometry Gravitaional-wave Observatory (LIGO) recently detected gravitational waves from the mergers of binary black hole systems. Accurate calibration of the output of these detectors was crucial for the observation of these events, and the extraction of parameters of the sources. The principal tools used to calibrate the responses of the second-generation (Advanced) LIGO detectors to gravitational waves are systems based on radiation pressure and referred to as Photon Calibrators. These systems, which were completely redesigned for Advanced LIGO, include several significant upgrades that enable them to meet the calibration requirements of second-generation gravitational wave detectors in the new era of gravitational-wave astronomy. We report on the design, implementation, and operation of these Advanced LIGO Photon Calibrators that are currently providing fiducial displacements on the order of $10^{-18}$ m/$\sqrt{\textrm{Hz}}$ with accuracy and precision of better than 1 %. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: 14 pages, 19 figures