BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daniel Lecoanet

Daniel Lecoanet appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRphysics.flu-dynastro-ph.IMphysics.ao-phastro-ph.EPastro-ph.HEphysics.geo-phastro-phastro-ph.GABiological Physicsmath.NAphysics.comp-phphysics.plasm-phphysics.soc-ph

Catalog footprint

What is connected

27works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

27 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
preprint2016arXiv

Tensor calculus in polar coordinates using Jacobi polynomials

Spectral methods are an efficient way to solve partial differential equations on domains possessing certain symmetries. The utility of a method depends strongly on the choice of spectral basis. In this paper we describe a set of bases built out of Jacobi polynomials, and associated operators for solving scalar, vector, and tensor partial differential equations in polar coordinates on a unit disk. By construction, the bases satisfy regularity conditions at r=0 for any tensorial field. The coordinate singularity in a disk is a prototypical case for many coordinate singularities. The work presented here extends to other geometries. The operators represent covariant derivatives, multiplication by azimuthally symmetric functions, and the tensorial relationship between fields. These arise naturally from relations between classical orthogonal polynomials, and form a Heisenberg algebra. Other past work uses more specific polynomial bases for solving equations in polar coordinates. The main innovation in this paper is to use a larger set of possible bases to achieve maximum bandedness of linear operations. We provide a series of applications of the methods, illustrating their ease-of-use and accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulent Chemical Diffusion in Convectively Bounded Carbon Flames

It has been proposed that mixing induced by convective overshoot can disrupt the inward propagation of carbon deflagrations in super-asymptotic giant branch stars. To test this theory, we study an idealized model of convectively bounded carbon flames with 3D hydrodynamic simulations of the Boussinesq equations using the pseudospectral code Dedalus. Because the flame propagation timescale is much longer than the convection timescale, we approximate the flame as fixed in space, and only consider its effects on the buoyancy of the fluid. By evolving a passive scalar field, we derive a {\it turbulent} chemical diffusivity produced by the convection as a function of height, $D_{\rm t}(z)$. Convection can stall a flame if the chemical mixing timescale, set by the turbulent chemical diffusivity, $D_{\rm t}$, is shorter than the flame propagation timescale, set by the thermal diffusivity, $κ$, i.e., when $D_{\rm t}>κ$. However, we find $D_{\rm t}<κ$ for most of the flame because convective plumes are not dense enough to penetrate into the flame. Extrapolating to realistic stellar conditions, this implies that convective mixing cannot stall a carbon flame and that "hybrid carbon-oxygen-neon" white dwarfs are not a typical product of stellar evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Experimental study of internal wave generation by convection in water

We experimentally investigate the dynamics of water cooled from below at 0^oC and heated from above. Taking advantage of the unusual property that water's density maximum is at about 4^oC, this set-up allows us to simulate in the laboratory a turbulent convective layer adjacent to a stably stratified layer, which is representative of atmospheric and stellar conditions. High precision temperature and velocity measurements are described, with a special focus on the convectively excited internal waves propagating in the stratified zone. Most of the convective energy is at low frequency, and corresponding waves are localized to the vicinity of the interface. However, we show that some energy radiates far from the interface, carried by shorter horizontal wavelength, higher frequency waves. Our data suggest that the internal wave field is passively excited by the convective fluctuations, and the wave propagation is correctly described by the dissipative linear wave theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical Simulations of Internal Wave Generation by Convection in Water

Water's density maximum at 4C makes it well suited to study internal gravity wave excitation by convection: an increasing temperature profile is unstable to convection below 4C, but stably stratified above 4C. We present numerical simulations of a water-like fluid near its density maximum in a two dimensional domain. We successfully model the damping of waves in the simulations using linear theory, provided we do not take the weak damping limit typically used in the literature. In order to isolate the physical mechanism exciting internal waves, we use the novel spectral code Dedalus to run several simplified model simulations of our more detailed simulation. We use data from the full simulation as source terms in two simplified models of internal wave excitation by convection: bulk excitation by convective Reynolds stresses, and interface forcing via the mechanical oscillator effect. We find excellent agreement between the waves generated in the full simulation and the simplified simulation implementing the bulk excitation mechanism. The interface forcing simulations over excite high frequency waves because they assume the excitation is by the "impulsive" penetration of plumes, which spreads energy to high frequencies. However, we find the real excitation is instead by the "sweeping" motion of plumes parallel to the interface. Our results imply that the bulk excitation mechanism is a very accurate heuristic for internal wave generation by convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The spin rate of pre-collapse stellar cores: wave-driven angular momentum transport in massive stars

The core rotation rates of massive stars have a substantial impact on the nature of core-collapse supernovae and their compact remnants. We demonstrate that internal gravity waves (IGW), excited via envelope convection during a red supergiant phase or during vigorous late time burning phases, can have a significant impact on the rotation rate of the pre-SN core. In typical ($10 \, M_\odot \lesssim M \lesssim 20 \, M_\odot$) supernova progenitors, IGW may substantially spin down the core, leading to iron core rotation periods $P_{\rm min,Fe} \gtrsim 30 \, {\rm s}$. Angular momentum (AM) conservation during the supernova would entail minimum NS rotation periods of $P_{\rm min,NS} \gtrsim 3 \, {\rm ms}$. In most cases, the combined effects of magnetic torques and IGW AM transport likely lead to substantially longer rotation periods. However, the stochastic influx of AM delivered by IGW during shell burning phases inevitably spin up a slowly rotating stellar core, leading to a maximum possible core rotation period. We estimate maximum iron core rotation periods of $P_{\rm max,Fe} \lesssim 5 \times 10^3 \, {\rm s}$ in typical core-collapse supernova progenitors, and a corresponding spin period of $P_{\rm max, NS} \lesssim 500 \, {\rm ms}$ for newborn neutron stars. This is comparable to the typical birth spin periods of most radio pulsars. Stochastic spin-up via IGW during shell O/Si burning may thus determine the initial rotation rate of most neutron stars. For a given progenitor, this theory predicts a Maxwellian distribution in pre-collapse core rotation frequency that is uncorrelated with the spin of the overlying envelope.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Zombie Vortex Instability I: A Purely Hydrodynamic Instability to Resurrect the Dead Zones of Protoplanetary Disks

There is considerable interest in hydrodynamic instabilities in dead zones of protoplanetary disks as a mechanism for driving angular momentum transport and as a source of particle-trapping vortices to mix chondrules and incubate planetesimal formation. We present simulations with a pseudo-spectral anelastic code and with the compressible code Athena, showing that stably stratified flows in a shearing, rotating box are violently unstable and produce space-filling, sustained turbulence dominated by large vortices with Rossby numbers of order 0.2-0.3. This Zombie Vortex Instability (ZVI) is observed in both codes and is triggered by Kolmogorov turbulence with Mach numbers less than 0.01. It is a common view that if a given constant density flow is stable, then stable vertical stratification should make the flow even more stable. Yet, we show that sufficient vertical stratification can be unstable to ZVI. ZVI is robust and requires no special tuning of boundary conditions, or initial radial entropy or vortensity gradients (though we have studied ZVI only in the limit of infinite cooling time). The resolution of this paradox is that stable stratification allows for a new avenue to instability: baroclinic critical layers. ZVI has not been seen in previous studies of flows in rotating, shearing boxes because those calculations frequently lacked vertical density stratification and/or sufficient numerical resolution. Although we do not expect appreciable angular momentum transport from ZVI in the small domains in this study, we hypothesize that ZVI in larger domains with compressible equations may lead to angular transport via spiral density waves.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular Momentum Transport via Internal Gravity Waves in Evolving Stars

Recent asteroseismic advances have allowed for direct measurements of the internal rotation rates of many sub-giant and red giant stars. Unlike the nearly rigidly rotating Sun, these evolved stars contain radiative cores that spin faster than their overlying convective envelopes, but slower than they would in the absence of internal angular momentum transport. We investigate the role of internal gravity waves in angular momentum transport in evolving low mass stars. In agreement with previous results, we find that convectively excited gravity waves can prevent the development of strong differential rotation in the radiative cores of Sun-like stars. As stars evolve into sub-giants, however, low frequency gravity waves become strongly attenuated and cannot propagate below the hydrogen burning shell, allowing the spin of the core to decouple from the convective envelope. This decoupling occurs at the base of the sub-giant branch when stars have surface temperatures of roughly 5500 K. However, gravity waves can still spin down the upper radiative region, implying that the observed differential rotation is likely confined to the deep core near the hydrogen burning shell. The torque on the upper radiative region may also prevent the core from accreting high-angular momentum material and slow the rate of core spin-up. The observed spin-down of cores on the red giant branch cannot be totally attributed to gravity waves, but the waves may enhance shear within the radiative region and thus increase the efficacy of viscous/magnetic torques.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conduction in low Mach number flows: Part I Linear & weakly nonlinear regimes

Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law - the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion - is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entropy. This paper investigates the differences between temperature and entropy diffusion for both linear internal gravity waves and weakly nonlinear convection. In addition to simulating the two thermal conduction models with the fully compressible Navier-Stokes equations, we also study their effects in the reduced, "sound-proof" anelastic and pseudo-incompressible equations. We find that in the linear and weakly nonlinear regimes, temperature and entropy diffusion give quantitatively similar results, although there are some larger errors in the pseudo-incompressible equations with temperature diffusion due to inaccuracies in the equation of state. Extrapolating our weakly nonlinear results, we speculate that differences between temperature and entropy diffusion might become more important for strongly turbulent convection.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy conservation and gravity waves in sound-proof treatments of stellar interiors: Part II Lagrangian constrained analysis

The speed of sound greatly exceeds typical flow velocities in many stellar and planetary interiors. To follow the slow evolution of subsonic motions, various sound-proof models attempt to remove fast acoustic waves whilst retaining stratified convection and buoyancy dynamics. In astrophysics, anelastic models typically receive the most attention in the class of sound-filtered stratified models. Generally, anelastic models remain valid in nearly adiabatically stratified regions like stellar convection zones, but may break down in strongly sub-adiabatic, stably stratified layers common in stellar radiative zones. However, studying stellar rotation, circulation, and dynamos requires understanding the complex coupling between convection and radiative zones, and this requires robust equations valid in both regimes. Here we extend the analysis of equation sets begun in Brown Vasil & Zweibel 2012, which studied anelastic models, to two types of pseudo-incompressible models. This class of models has received attention in atmospheric applications, and more recently in studies of white-dwarf supernovae progenitors. We demonstrate that one model conserves energy but the other does not. We use Lagrangian variational methods to extend the energy conserving model to a general equation of state, and dub the resulting equation set the Generalized Pseudo-Incompressible (GPI) model. We show that the GPI equations suitably capture low frequency phenomena in both convection and radiative zones in stars and other stratified systems, and we provide recommendations for converting low-Mach number codes to this equation set.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Internal Gravity Wave Excitation by Turbulent Convection

We calculate the flux of internal gravity waves (IGWs) generated by turbulent convection in stars. We solve for the IGW eigenfunctions analytically near the radiative-convective interface in a local, Boussinesq, and cartesian domain. We consider both discontinuous and smooth transitions between the radiative and convective regions and derive Green's functions to solve for the IGWs in the radiative region. We find that if the radiative-convective transition is smooth, the IGW flux depends on the exact form of the buoyancy frequency near the interface. IGW excitation is most efficient for very smooth interfaces, which gives an upper bound on the IGW flux of ~ F_conv (d/H), where F_conv is the flux carried by the convective motions, d is the width of the transition region, and H is the pressure scale height. This can be much larger than the standard result in the literature for a discontinuous radiative-convective transition, which gives a wave flux ~ F_conv M$, where M is the convective Mach number. However, in the smooth transition case, the most efficiently excited perturbations will break in the radiative zone. The flux of IGWs which do not break and are able to propagate in the radiative region is at most F_conv M^(5/8) (d/H)^(3/8), larger than the discontinuous transition result by (MH/d)^(-3/8). The transition region in the Sun is smooth for the energy-bearing waves; as a result, we predict that the IGW flux is a few to five times larger than previous estimates. We discuss the implications of our results for several astrophysical applications, including IGW driven mass loss and the detectability of convectively excited IGWs in main sequence stars.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Dynamics of Rayleigh-Taylor Stable and Unstable Contact Discontinuities with Anisotropic Thermal Conduction

We study the effects of anisotropic thermal conduction along magnetic field lines on an accelerated contact discontinuity in a weakly collisional plasma. We first perform a linear stability analysis similar to that used to derive the Rayleigh-Taylor instability (RTI) dispersion relation. We find that anisotropic conduction is only important for compressible modes, as incompressible modes are isothermal. Modes grow faster in the presence of anisotropic conduction, but growth rates do not change by more than a factor of order unity. We next run fully non-linear numerical simulations of a contact discontinuity with anisotropic conduction. The non-linear evolution can be thought of as a superposition of three physical effects: temperature diffusion due to vertical conduction, the RTI, and the heat flux driven buoyancy instability (HBI). In simulations with RTI-stable contact discontinuities, the temperature discontinuity spreads due to vertical heat conduction. This occurs even for initially horizontal magnetic fields due to the initial vertical velocity perturbation and numerical mixing across the interface. The HBI slows this temperature diffusion by reorienting initially vertical magnetic field lines to a more horizontal geometry. In simulations with RTI-unstable contact discontinuities, the dynamics are initially governed by temperature diffusion, but the RTI becomes increasingly important at late times. We discuss the possible application of these results to supernova remnants, solar prominences, and cold fronts in galaxy clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Violation of Richardson's Criterion via Introduction of a Magnetic Field

Shear flow instabilities can profoundly affect the diffusion of momentum in jets, stars, and disks. The Richardson criterion gives a sufficient condition for instability of a shear flow in a stratified medium. The velocity gradient $V'$ can only destabilize a stably stratified medium with squared Brunt-Vaisala frequency $N^2$ if $V'^2/4>N^2$. We find this is no longer true when the medium is a magnetized plasma. We investigate the effect of stable stratification on magnetic field and velocity profiles unstable to magneto-shear instabilities, i.e., instabilities which require the presence of both magnetic field and shear flow. We show that a family of profiles originally studied by Tatsuno & Dorland (2006) remain unstable even when $V'^2/4<N^2$, violating the Richardson criterion. However, not all magnetic fields can result in a violation of the Richardson criterion. We consider a class of flows originally considered by Kent (1968), which are destabilized by a constant magnetic field, and show that they become stable when $V'^2/4<N^2$, as predicted by the Richardson criterion. This suggests that magnetic free energy is required to violate the Richardson criterion. This work implies that the Richardson criterion cannot be used when evaluating the ideal stability of a sheared, stably stratified, and magnetized plasma. We briefly discuss the implications for astrophysical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Mean Motion Resonances in Extrasolar Planetary Systems with Turbulence, Interactions, and Damping

This paper continues previous work on the effects of turbulence on mean motion resonances in extrasolar planetary systems. Turbulence is expected to arise in the disks that form planets, and these fluctuations act to compromise resonant configurations. This paper extends previous work by considering how interactions between the planets and possible damping effects imposed by the disk affect the outcomes. These physical processes are studied using three approaches: numerical integrations of the 3-body problem with additional forcing due to turbulence, model equations that reduce the problem to stochastically driven oscillators, and Fokker-Planck equations that describe the time evolution of an ensemble of systems. With this combined approach, we elucidate the physics of how turbulence can remove extrasolar planetary systems from mean motion resonance. As expected, systems with sufficiently large damping (dissipation) can maintain resonance, in spite of turbulent forcing. In the absence of strong damping, ensembles of these systems exhibit two regimes of behavior, where the fraction of the bound states decreases as a power-law or as an exponential. Both types of behavior can be understood through the model developed herein. For systems with weak interactions between planets, the model reduces to a stochastic pendulum, and the fraction of bound states decreases as a power-law. For highly interactive systems, the dynamics are more complicated and the fraction of bound states decreases exponentially. We show how planetary interactions lead to drift terms in the Fokker-Planck equation and account for this exponential behavior. In addition to clarifying the physical processes involved, this paper strengthens the finding that turbulence implies that mean motions resonances should be rare.

0 citations0 reviewsNo code linkedOpen access