BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Daniel Lecoanet

Daniel Lecoanet contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.SRphysics.flu-dynastro-ph.IMphysics.ao-phastro-ph.HEastro-ph.EPastro-ph.GABiological Physicsphysics.comp-phphysics.geo-phphysics.plasm-phphysics.soc-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

An Atlas of Convection in Main-Sequence Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Asteroseismic inference of the near-core magnetic field strength in the main sequence B star HD 43317

About 10 per cent of intermediate- and high-mass dwarf stars are observed to host a strong large-scale magnetic field at their surface, which is thought to be of fossil field origin. However, there are few inferences as to the magnetic field strength and geometry within the deep interiors of stars. Considering that massive stars harbour a convective core whilst on the main sequence, asteroseismology of gravity (g) modes is able to provide constraints on their core masses, yet it has so far not been used to probe the strength of their interior magnetic fields. Here we use asteroseismology to constrain an upper limit for the magnetic field strength in the near-core region of the pulsating and magnetic B star HD 43317, based on the expected interaction of a magnetic field and its g modes. We find a magnetic field strength of order $5 \times 10^5$ G is sufficient to suppress high-radial order g modes and reproduce the observed frequency spectrum of HD 43317, which contains only high-frequency g modes. This result is the first inference of the magnetic field strength inside a main-sequence star.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Convective Penetration in Early-Type Stars

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gender Disparity in Publishing Six Months after the KITP Workshop "Probes of Transport in Stars"

Conferences and workshops shape scientific discourse. The Kavli Institute for Theoretical Physics (KITP) hosts long-term workshops to stimulate scientific collaboration that would not otherwise have taken place. One goal of KITP programs is to increase diversity in the next generation of scientists. In this analysis, we examine gender trends in authorship of papers generated as a result of the KITP program \textit{Probes of Transport in Stars}, which ran from October 11th, 2021 to December 17th, 2021. While 38\% of workshop participants were women, only 19\% of publications produced between December 1st, 2021 and June 3rd, 2022 had female first-authors. Further, of these early publications, 61\% had all-male author lists. Among publications resulting from the KITP program, the proportions of both male first-author papers and papers with all-male author lists are higher than predicted by models that take into account the gender distribution of the KITP participants. These results motivate more thorough investigations of collaboration networks at scientific conferences and workshops. Importantly, they also suggest that programs, conferences, and workshops of any kind need to take steps beyond those implemented in this KITP program to enable more diverse collaborations and address gender disparities in science.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Moosinesq Convection in the Cores of Moosive Stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Schwarzschild and Ledoux are equivalent on evolutionary timescales

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar convective penetration: parameterized theory and dynamical simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the Origin of Stochastic, Low-Frequency Photometric Variability in Massive Stars

High-precision photometric observations have revealed ubiquitous stochastic low-frequency photometric variability in early type stars. It has been suggested that this variability arises due to either subsurface convection or internal gravity waves launched by the convective core. Here we show that relevant properties of convection in subsurface convective layers correlate very well with the timescale and amplitude of stochastic low-frequency photometric variability, as well as with the amplitude of macroturbulence. We suggest that low-frequency, stochastic photometric variability and surface turbulence in massive stars are caused by the the presence of subsurface convection. We show that an explanation for the observed surface photometric variability and macroturbulence relying on convective core driven internal gravity waves encounters a number of difficulties and seems unlikely to be able to explain the observed trends.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shape and size of large-scale vortices: A generic fluid pattern in geophysical fluid dynamics

Planetary rotation organizes fluid motions into coherent, long-lived swirls, known as large-scale vortices (LSVs), which play an important role in the dynamics and long-term evolution of geophysical and astrophysical fluids. Here, using direct numerical simulations, we show that LSVs in rapidly rotating mixed convective and stably stratified fluids, which approximates the two-layer, turbulent-stratified dynamics of many geophysical and astrophysical fluids, have a generic shape and that their size can be predicted. We show that LSVs emerge in the convection zone from upscale energy transfers and can penetrate into the stratified layer. At the convective-stratified interface, the LSV cores have a positive buoyancy anomaly. Due to the thermal wind constraint, this buoyancy anomaly leads to winds in the stratified layer that decay over a characteristic vertical length scale. Thus LSVs take the shape of a depth-invariant cylinder with a finite-size radius in the turbulent layer and of a penetrating half dome in the stratified layer. Importantly, we demonstrate that when LSVs penetrate all the way through the stratified layer and reach a boundary that is no-slip, they saturate by boundary friction. We provide a prediction for the penetration depth and maximum radius of LSVs as a function of the LSV vorticity, the stratified layer depth, and the stratification. Our results, which apply for cyclonic LSVs, suggest that LSVs in slowly rotating stars and Earth's liquid core are confined to the convective layer, while in Earth's atmosphere and oceans they can penetrate far into the stratified layer.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Buoyancy-Driven Entrainment in Dry Thermals

\citet{turner1957} proposed that dry thermals entrain because of buoyancy (via a constraint which requires an increase in the radius $a$). This however, runs counter to the scaling arguments commonly used to derive the entrainment rate, which rely on either the self-similarity of \citet{scorer1957} or the turbulent entrainment hypothesis of \citet{morton1956}. The assumption of turbulence-driven entrainment was investigated by \citet{lecoanet2018}, who found that the entrainment efficiency $e$ varies by less than $20\%$ between laminar (Re = 630) and turbulent (Re = 6300) thermals. This motivated us to utilize Turner's argument of buoyancy-controlled entrainment in addition to the thermal's vertical momentum equation to build a model for thermal dynamics which does not invoke turbulence or self-similarity. We derive simple expressions for the thermals' kinematic properties and their fractional entrainment rate $ε$ and find close quantitative agreement with the values in direct numerical simulations. In particular, our expression for entrainment rate is consistent with the parameterization $ε\sim B/w^2$, for Archimedean buoyancy $B$ and vertical velocity $w$. We also directly validate the role of buoyancy-driven entrainment by running simulations where gravity is turned off midway through a thermal's rise. The entrainment efficiency $e$ is observed to drop to less than 1/3 of its original value in both the laminar and turbulent cases when $g=0$, affirming the central role of buoyancy in entrainment in dry thermals.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Dedalus: A Flexible Framework for Numerical Simulations with Spectral Methods

Numerical solutions of partial differential equations enable a broad range of scientific research. The Dedalus Project is a flexible, open-source, parallelized computational framework for solving general partial differential equations using spectral methods. Dedalus translates plain-text strings describing partial differential equations into efficient solvers. This paper details the numerical method that enables this translation, describes the design and implementation of the codebase, and illustrates its capabilities with a variety of example problems. The numerical method is a first-order generalized tau formulation that discretizes equations into banded matrices. This method is implemented with an object-oriented design. Classes for spectral bases and domains manage the discretization and automatic parallel distribution of variables. Discretized fields and mathematical operators are symbolically manipulated with a basic computer algebra system. Initial value, boundary value, and eigenvalue problems are efficiently solved using high-performance linear algebra, transform, and parallel communication libraries. Custom analysis outputs can also be specified in plain text and stored in self-describing portable formats. The performance of the code is evaluated with a parallel scaling benchmark and a comparison to a finite-volume code. The features and flexibility of the codebase are illustrated by solving several examples: the nonlinear Schrodinger equation on a graph, a supersonic magnetohydrodynamic vortex, quasigeostrophic flow, Stokes flow in a cylindrical annulus, normal modes of a radiative atmosphere, and diamagnetic levitation. The Dedalus code and the example problems are available online at http://dedalus-project.org/.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Entrainment in Resolved, Dry Thermals

Entrainment in cumulus convection remains ill-understood and difficult to quantify. For instance, entrainment is widely believed to be a fundamentally turbulent process, even though Turner (1957) pointed out that dry thermals entrain primarily because of buoyancy (via a dynamical constraint requiring an increase in radius $r$), rather than turbulence. Furthermore, entrainment has been postulated to obey a $1/r$ scaling, but this scaling has not been firmly established. Here, we study the classic case of dry, turbulent thermals in a neutrally stratified environment using fully resolved direct numerical simulation. We combine this with a thermal tracking algorithm which defines a control volume for the thermal at each time, allowing us to directly measure entrainment. We test Turner's argument by varying the Reynolds number Re of our thermals between laminar (Re~600) and turbulent (Re~6000) regimes, finding only a 20% variation in entrainment rate $ε$, supporting the claim that turbulence is not necessary for entrainment. We also directly verify the postulated $ε\sim 1/r$ scaling law.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Low-frequency variability in massive stars: Core generation or surface phenomenon?

Bowman et al. (2019) reported low-frequency photometric variability in 164 O- and B-type stars observed with K2 and TESS. They interpret these motions as internal gravity waves, which could be excited stochastically by convection in the cores of these stars. The detection of internal gravity waves in massive stars would help distinguish between massive stars with convective or radiative cores, determine core size, and would provide important constraints on massive star structure and evolution. In this work, we study the observational signature of internal gravity waves generated by core convection. We calculate the \textit{wave transfer function}, which links the internal gravity wave amplitude at the base of the radiative zone to the surface luminosity variation. This transfer function varies by many orders of magnitude for frequencies $\lesssim 1 \, {\rm d}^{-1}$, and has regularly-spaced peaks near $1 \, {\rm d}^{-1}$ due to standing modes. This is inconsistent with the observed spectra which have smooth ``red noise'' profiles, without the predicted regularly-spaced peaks. The wave transfer function is only meaningful if the waves stay predominately linear. We next show that this is the case: low frequency traveling waves do not break unless their luminosity exceeds the radiative luminosity of the star; and, the observed luminosity fluctuations at high frequencies are so small that standing modes would be stable to nonlinear instability. These simple calculations suggest that the observed low-frequency photometric variability in massive stars is not due to internal gravity waves generated in the core of these stars. We finish with a discussion of (sub)surface convection, which produces low-frequency variability in low-mass stars, very similar to that observed in Bowman et al. (2019) in higher mass stars.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Magnetorotational Instability Prefers Three Dimensions

The magnetorotational instability (MRI) occurs when a weak magnetic field destabilises a rotating, electrically conducting fluid with inwardly increasing angular velocity. The MRI is essential to astrophysical disk theory where the shear is typically Keplerian. Internal shear layers in stars may also be MRI unstable, and they take a wide range of profiles, including near-critical. We show that the fastest growing modes of an ideal magnetofluid are three-dimensional provided the shear rate, $S$, is near the two-dimensional onset value, $S_c$. For a Keplerian shear, three-dimensional modes are unstable above $S\approx0.10S_c$, and dominate the two-dimensional modes until $S\approx2.05S_{c}$. These three-dimensional modes dominate for shear profiles relevant to stars and at magnetic Prandtl numbers relevant to liquid-metal laboratory experiments. Significant numbers of rapidly growing three-dimensional modes remain well past $2.05S_{c}$. These finding are significant in three ways. First, weakly nonlinear theory suggests that the MRI saturates by pushing the shear rate to its critical value. This can happen for systems, like stars and laboratory experiments, that can rearrange their angular velocity profiles. Second, the non-normal character and large transient growth of MRI modes should be important whenever three-dimensionality exists. Finally, three-dimensional growth suggests direct dynamo action driven from the linear instability.

0 citations0 reviewsNo code linkedOpen access