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

Pascale Garaud

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.EPphysics.flu-dynphysics.plasm-ph

Catalog footprint

What is connected

25works
4topics
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

25 published item(s)

preprint2022arXiv

Magnetized oscillatory double-diffusive convection

We study the properties of oscillatory double-diffusive convection (ODDC) in the presence of a uniform vertical background magnetic field. ODDC takes place in stellar regions that are unstable according to the Schwarzschild criterion and stable according to the Ledoux criterion (sometimes called semiconvective regions), which are often predicted to reside just outside the core of intermediate-mass main sequence stars. Previous hydrodynamic studies of ODDC have shown that the basic instability saturates into a state of weak wave-like convection, but that a secondary instability can sometimes transform it into a state of layered convection, where layers then rapidly merge and grow until the entire region is fully convective. We find that magnetized ODDC has very similar properties overall, with some important quantitative differences. A linear stability analysis reveals that the fastest-growing modes are unaffected by the field, but that other modes are. Numerically, the magnetic field is seen to influence the saturation of the basic instability, overall reducing the turbulent fluxes of temperature and composition. This in turn affects layer formation, usually delaying it, and occasionally suppressing it entirely for sufficiently strong fields. Further work will be needed, however, to determine the field strength above which layer formation is actually suppressed in stars. Potential observational implications are briefly discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

On the dynamical interaction between overshooting convection and an underlying dipole magnetic field -- I. The non-dynamo regime

Motivated by the dynamics in the deep interiors of many stars, we study the interaction between overshooting convection and the large-scale poloidal fields residing in radiative zones. We have run a suite of 3D Boussinesq numerical calculations in a spherical shell that consists of a convection zone with an underlying stable region that initially compactly contains a dipole field. By varying the strength of the convective driving, we find that, in the less turbulent regime, convection acts as turbulent diffusion that removes the field faster than solely molecular diffusion would do. However, in the more turbulent regime, turbulent pumping becomes more efficient and partially counteracts turbulent diffusion, leading to a local accumulation of the field below the overshoot region. These simulations suggest that dipole fields might be confined in underlying stable regions by highly turbulent convective motions at stellar parameters. The confinement is of large-scale field in an average sense and we show that it is reasonably modeled by mean-field ideas. Our findings are particularly interesting for certain models of the Sun, which require a large-scale, poloidal magnetic field to be confined in the solar radiative zone in order to explain simultaneously the uniform rotation of the latter and the thinness of the solar tachocline.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The tachocline revisited

The solar tachocline is a shear layer located at the base of the solar convection zone. The horizontal shear in the tachocline is likely turbulent, and it is often assumed that this turbulence would be strongly anisotropic as a result of the local stratification. What role this turbulence plays in the tachocline dynamics, however, remains to be determined. In particular, it is not clear whether it would result in a turbulent eddy diffusivity, or anti-diffusivity, or something else entirely. In this paper, we present the first direct numerical simulations of turbulence in horizontal shear flows at low Prandtl number, in an idealized model that ignores rotation and magnetic fields. We find that several regimes exist, depending on the relative importance of the stratification, viscosity and thermal diffusivity. Our results suggest that the tachocline is in the stratified turbulence regime, which has very specific properties controlled by a balance between buoyancy, inertia, and thermal diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effect of rotation on oscillatory double-diffusive convection (semiconvection)

Oscillatory double-diffusive convection (ODDC, more traditionally called semiconvection) is a form of linear double-diffusive instability that occurs in fluids that are unstably stratified in temperature (Schwarzschild unstable), but stably stratified in chemical composition (Ledoux stable). This scenario is thought to be quite common in the interiors of stars and giant planets, and understanding the transport of heat and chemical species by ODDC is of great importance to stellar and planetary evolution models. Fluids unstable to ODDC have a tendency to form convective thermo-compositional layers which significantly enhance the fluxes of temperature and chemical composition compared with microscopic diffusion. Although a number of recent studies have focused on studying properties of both layered and non-layered ODDC, few have addressed how additional physical processes such as global rotation affect its dynamics. In this work we study first how rotation affects the linear stability properties of rotating ODDC. Using direct numerical simulations we then analyze the effect of rotation on properties of layered and non-layered ODDC, and study how the angle of the rotation axis with respect to the direction of gravity affects layering. We find that rotating systems can be broadly grouped into two categories, based on the strength of rotation. Qualitative behavior in the more weakly rotating group is similar to non-rotating ODDC, but strongly rotating systems become dominated by vortices that are invariant in the direction of the rotation vector and strongly influence transport. We find that whenever layers form, rotation always acts to reduce thermal and compositional transport.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

2D or not 2D: the effect of dimensionality on the dynamics of fingering convection at low Prandtl number

Fingering convection (otherwise known as thermohaline convection) is an instability that occurs in stellar radiative interiors in the presence of unstable compositional gradients. Numerical simulations have been used in order to estimate the efficiency of mixing induced by this instability. However, fully three-dimensional (3D) computations in the parameter regime appropriate for stellar astrophysics (i.e. low Prandtl number) are prohibitively expensive. This raises the question of whether two-dimensional (2D) simulations could be used instead to achieve the same goals. In this work, we address this issue by comparing the outcome of 2D and 3D simulations of fingering convection at low Prandtl number. We find that 2D simulations are never appropriate. However, we also find that the required 3D computational domain does not have to be very wide: the third dimension need only contain a minimum of two wavelengths of the fastest-growing linearly unstable mode to capture the essentially 3D dynamics of small-scale fingering. Narrow domains, however, should still be used with caution since they could limit the subsequent development of any large-scale dynamics typically associated with fingering convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A new model for mixing by double-diffusive convection (semi-convection). II. The transport of heat and composition through layers

Regions of stellar and planetary interiors that are unstable according to the Schwarzschild criterion, but stable according to the Ledoux criterion, are subject to a form of oscillatory double-diffusive (ODD) convection often called "semi-convection". In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and composition by ODD convection, and ultimately propose a new 1D prescription that can be used in stellar and planetary structure and evolution models. The first paper in this series demonstrated that under certain conditions ODD convection spontaneously transitions from an initially homogeneously turbulent state into a staircase of convective layers, which results in a substantial increase in the transport of heat and composition. Here, we present simulations of ODD convection in this layered regime, we describe the dynamical behavior of the layers, and we derive empirical scaling laws for the transport through layered convection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A New Model for Mixing By Double-Diffusive Convection (Semi-Convection). III. Thermal and Compositional Transport Through Non-Layered ODDC

Oscillatory double-diffusive convection (ODDC) (also known as semi- convection) refers to a type of double diffusive instability that occurs in regions of planetary and stellar interiors which have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. (2011) first showed that ODDC can either spontaneously form layers, which significantly enhance the transport of heat and chemical species compared to mi- croscopic transport, or remain in a state dominated by large scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series have focused on identifying under what condi- tions layers form (Mirouh et al. 2012), and quantifying transport through layered systems (Wood et al. 2013). Here we proceed to characterize transport through systems that are unstable to the ODDC instability, but do not undergo spon- taneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both 2D and 3D numerical simulations and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport by ODDC.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Main Sequence Evolution with Layered Semiconvection

Semiconvection - mixing that occurs in regions that are stable when considering compositional gradients, but unstable when ignoring them - is shown to have the greatest potential impact on main sequence stars with masses in the range 1.2 - 1.7 solar masses. We present the first stellar evolution calculations using a prescription for semiconvection derived from extrapolation of direct numerical simulations of double-diffusive mixing down to stellar parameters. The dominant mode of semiconvection in stars is layered semiconvection, where the layer height is an adjustable parameter analogous to the mixing length in convection. The rate of mixing across the semiconvective region is sensitively dependent on the layer height. We find that there is a critical layer height that separates weak semiconvective mixing (where evolution is well-approximated by using the Ledoux criterion) from strong semiconvective mixing (where evolution is well-approximated by using the Schwarzschild criterion). This critical layer height is much smaller than the minimum layer height expected from simulations so we predict that for realistic layer heights, the evolution is nearly the same as a model ran with the Schwarzschild criterion. We also investigate the effects of compositional gradient smoothing, finding that it causes convective cores to artificially shrink in the absence of additional mixing beyond the convective boundary. Layered semiconvection with realistic layer heights provides enough such mixing that stars will still evolve as if the Schwarzschild criterion is employed. Finally, we discuss the potential of detecting such semiconvection and its implication on convective core sizes in solar-like oscillators.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stability of stratified spatially periodic shear flows at low Péclet number

This work addresses the question of the stability of stratified, spatially periodic shear flows at low Péclet number but high Reynolds number. This little-studied limit is motivated by astrophysical systems, where the Prandtl number is often very small. Furthermore, it can be studied using a reduced set of "low-Péclet-number equations" proposed by Lignieres [Astronomy & Astrophysics, 348, 933-939, 1999]. Through a linear stability analysis, we first determine the conditions for instability to infinitesimal perturbations. We formally extend Squire's theorem to the low-Péclet-number equations, which shows that the first unstable mode is always two-dimensional. We then perform an energy stability analysis of the low-Péclet-number equations and prove that for a given value of the Reynolds number, above a critical strength of the stratification, any smooth periodic shear flow is stable to perturbations of arbitrary amplitude. In that parameter regime, the flow can only be laminar and turbulent mixing does not take place. Finding that the conditions for linear and energy stability are different, we thus identify a region in parameter space where finite-amplitude instabilities could exist. Using direct numerical simulations, we indeed find that the system is subject to such finite-amplitude instabilities. We determine numerically how far into the linearly stable region of parameter space turbulence can be sustained.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Turbulent transport in a strongly stratified forced shear layer with thermal diffusion

This work presents numerical results on the transport of heat and chemical species by shear-induced turbulence in strongly stratified but thermally diffusive environments. The shear instabilities driven in this regime are sometimes called "secular" shear instabilities, and can take place even when the gradient Richardson number of the flow (the square of the ratio of the buoyancy frequency to the shearing rate) is large, provided the Péclet number (the ratio of the thermal diffusion timescale to the turnover timescale of the turbulent eddies) is small. We have identified a set of simple criteria to determine whether these instabilities can take place or not. Generally speaking, we find that they may be relevant whenever the thermal diffusivity of the fluid is very large (typically larger than $10^{14}$cm$^2$/s), which is the case in the outer layers of high-mass stars ($M\ge 10 M_\odot$) for instance. Using a simple model setup in which the shear is forced by a spatially sinusoidal, constant-amplitude body-force, we have identified several regimes ranging from effectively unstratified to very strongly stratified, each with its own set of dynamical properties. Unless the system is in one of the two extreme regimes (effectively unstratified or completely stable), however, we find that (1) only about 10% of the input power is used towards heat transport, while the remaining 90% is viscously dissipated; (2) that the effective compositional mixing coefficient is well-approximated by the model of Zahn (1992), with $D \simeq 0.02 κ_T /J$ where $κ_T$ is the thermal diffusivity and $J$ is the gradient Richardson number. These results need to be confirmed, however, with simulations in different model setups and at higher effective Reynolds number.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Double-Diffusive Convection

Much progress has recently been made in understanding and quantifying vertical mixing induced by double-diffusive instabilities such as fingering convection (usually called thermohaline convection) and oscillatory double-diffusive convection (a process closely related to semiconvection). This was prompted in parts by advances in supercomputing, which allow us to run Direct Numerical Simulations of these processes at parameter values approaching those relevant in stellar interiors, and in parts by recent theoretical developments in oceanography where such instabilities also occur. In this paper I summarize these recent findings, and propose new mixing parametrizations for both processes that can easily be implemented in stellar evolution codes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Double-diffusive mixing in stellar interiors in the presence of horizontal gradients

We have identified an important source of mixing in stellar radiation zones, that would arise whenever two conditions are satisfied: (1) the presence of an inverse vertical compositional gradient, and (2) the presence of density-compensating horizontal gradients of temperature (alternatively, entropy) and composition. The former can be caused naturally by any off-center burning process, by atomic diffusion, or by surface accretion. The latter could be caused by rotation, tides, meridional flows, etc. The linear instability and its nonlinear development have been well-studied in the oceanographic context. It is known to drive the formation of stacks of fingering layers separated by diffusive interfaces, called intrusions. Using 3D numerical simulations of the process in the astrophysically-relevant region of parameter space, we find similar results, and demonstrate that the material transport in the intrusive regime can be highly enhanced compared with pure diffusion, even in systems which would otherwise be stable to fingering (thermohaline) convection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fingering convection induced by atomic diffusion in stars: 3D numerical computations and applications to stellar models

Iron-rich layers are known to form in the stellar subsurface through a combination of gravitational settling and radiative levitation. Their presence, nature and detailed structure can affect the excitation process of various stellar pulsation modes, and must therefore be modeled carefully in order to better interpret Kepler asteroseismic data. In this paper, we study the interplay between atomic diffusion and fingering convection in A-type stars, and its role in the establishment and evolution of iron accumulation layers. To do so, we use a combination of three-dimensional idealized numerical simulations of fingering convection, and one-dimensional realistic stellar models. Using the three-dimensional simulations, we first validate the mixing prescription for fingering convection recently proposed by Brown et al. (2013), and identify what system parameters (total mass of iron, iron diffusivity, thermal diffusivity, etc.) play a role in the overall evolution of the layer. We then implement the Brown et al. (2013) prescription in the Toulouse-Geneva Evolution code to study the evolution of the iron abundance profile beneath the stellar surface. We find, as first discussed by Théado et al. (2009), that when the concurrent settling of helium is ignored, this accumulation rapidly causes an inversion in the mean molecular weight profile, which then drives fingering convection. The latter mixes iron with the surrounding material very efficiently, and the resulting iron layer is very weak. However, taking helium settling into account partially stabilizes the iron profile against fingering convection, and a large iron overabundance can accumulate. The opacity also increases significantly as a result, and in some cases ultimately triggers dynamical convection.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Spin-down dynamics of magnetized solar-type stars

It has long been known that solar-type stars undergo significant spin-down, via magnetic braking, during their Main-Sequence lifetimes. However, magnetic braking only operates on the surface layers; it is not yet completely understood how angular momentum is transported within the star, and how rapidly the spin-down information is communicated to the deep interior. In this work, we use insight from recent progress in understanding internal solar dynamics to model the interior of other solar-type stars. We assume, following Gough and McIntyre (1998), that the bulk of the radiation zone of these stars is held in uniform rotation by the presence of an embedded large-scale primordial field, confined below a stably-stratified, magnetic-free tachocline by large-scale meridional flows downwelling from the convection zone. We derive simple equations to describe the response of this model interior to spin-down of the surface layers, that are identical to the two-zone model of MacGregor and Brenner (1991), with a coupling timescale proportional to the local Eddington-Sweet timescale across the tachocline. This timescale depends both on the rotation rate of the star and on the thickness of the tachocline, and can vary from a few hundred thousand years to a few Gyr, depending on stellar properties. Qualitative predictions of the model appear to be consistent with observations, although depend sensitively on the assumed functional dependence of the tachocline thickness on the stellar rotation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chemical Transport and Spontaneous Layer Formation in Fingering Convection in Astrophysics

A region of a star that is stable to convection according to the Ledoux criterion may nevertheless undergo additional mixing if the mean molecular weight increases with radius. This process is called fingering (thermohaline) convection and may account for some of the unexplained mixing in stars such as those that have been polluted by planetary infall and those burning helium-3. We propose a new model for mixing by fingering convection in the parameter regime relevant for stellar (and planetary) interiors. Our theory is based on physical principles and supported by three-dimensional direct numerical simulations. We also discuss the possibility of formation of thermocompositional staircases in fingering regions, and their role in enhancing mixing. Finally, we provide a simple algorithm to implement this theory in one-dimensional stellar codes, such as KEPLER and MESA.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of the solar tachocline III: Numerical solutions of the Gough and McIntyre model

We present the first numerical simulations of the solar interior to exhibit a tachocline consistent with the Gough and McIntyre (1998) model. We find nonlinear, axisymmetric, steady-state numerical solutions in which: (1) a large-scale primordial field is confined within the radiation zone by downwelling meridional flows that are gyroscopically pumped in the convection zone (2) the radiation zone is in almost-uniform rotation, with a rotation rate consistent with observations (3) the bulk of the tachocline is magnetic free, in thermal-wind balance and in thermal equilibrium and (4) the interaction between the field and the flows takes place within a very thin magnetic boundary layer, the tachopause, located at the bottom of the tachocline. We show that the thickness of the tachocline scales with the amplitude of the meridional flows exactly as predicted by Gough and McIntyre. We also determine the parameter conditions under which such solutions can be obtained, and provide a simple explanation for the failure of previous numerical attempts at reproducing the Gough and McIntyre model. Finally, we discuss the implications of our findings for future numerical models of the solar interior, and for future observations of the Sun and other stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large grains can grow in circumstellar discs

We perform coagulation & fragmentation simulations to understand grain growth in T Tauri & brown dwarf discs. We present a physically-motivated approach using a probability distribution function for the collision velocities and separating the deterministic & stochastic velocities. We find growth to larger sizes compared to other models. Furthermore, if brown dwarf discs are scaled-down versions of T Tauri discs (in terms of stellar & disc mass, and disc radius), growth at the same location with respect to the outer edge occurs to similar sizes in both discs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From dust to planetesimals: an improved model for collisional growth in protoplanetary disks

Planet formation occurs within the gas and dust rich environments of protoplanetary disks. Observations of these objects show that the growth of primordial sub micron sized particles into larger aggregates occurs at the earliest stages of the disks. However, theoretical models of particle growth that use the Smoluchowski equation to describe collisional coagulation and fragmentation have so far failed to produce large particles while maintaining a significant populations of small grains. This has been generally attributed to the existence of two barriers impeding growth due to bouncing and fragmentation of colliding particles. In this paper, we demonstrate that the importance of these barriers has been artificially inflated through the use of simplified models that do not take into account the stochastic nature of the particle motions within the gas disk. We present a new approach in which the relative velocities between two particles is described by a probability distribution function that models both deterministic motion and stochastic motion. Taking both into account can give quite different results to what has been considered recently in other studies. We demonstrate the vital effect of two "ingredients" for particle growth: the proper implementation of a velocity distribution function that overcomes the bouncing barrier and, in combination with mass transfer in high-mass-ratio collisions, boosts the growth of larger particles beyond the fragmentation barrier. A robust result of our simulations is the emergence of two particle populations (small and large), potentially explaining simultaneously a number of long-standing problems in protoplanetary disks, including planetesimal formation close to the central star, the presence of mm to cm size particles far out in the disk, and the persistence of micron-size grains for millions of years.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A new model for mixing by double-diffusive convection (semi-convection): I. The conditions for layer formation

The process referred to as "semi-convection" in astrophysics and "double-diffusive convection in the diffusive regime" in Earth and planetary sciences, occurs in stellar and planetary interiors in regions which are stable according to the Ledoux criterion but unstable according to the Schwarzschild criterion. In this series of papers, we analyze the results of an extensive suite of 3D numerical simulations of the process, and ultimately propose a new 1D prescription for heat and compositional transport in this regime which can be used in stellar or planetary structure and evolution models. In a preliminary study of the phenomenon, Rosenblum et al. (2011) showed that, after saturation of the primary instability, a system can evolve in one of two possible ways: the induced turbulence either remains homogeneous, with very weak transport properties, or transitions into a thermo-compositional staircase where the transport rate is much larger (albeit still smaller than in standard convection). In this paper, we show that this dichotomous behavior is a robust property of semi-convection across a wide region of parameter space. We propose a simple semi-analytical criterion to determine whether layer formation is expected or not, and at what rate it proceeds, as a function of the background stratification and of the diffusion parameters (viscosity, thermal diffusivity and compositional diffusivity) only. The theoretical criterion matches the outcome of our numerical simulations very adequately in the numerically accessible "planetary" parameter regime, and can easily be extrapolated to the stellar parameter regime. Subsequent papers will address more specifically the question of quantifying transport in the layered case and in the non-layered case.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Sun's meridional circulation and interior magnetic field

To date, no self-consistent numerical simulation of the solar interior has succeeded in reproducing the observed thinness of the solar tachocline, and the persistence of uniform rotation beneath it. Although it is known that the uniform rotation can be explained by the presence of a global-scale confined magnetic field, numerical simulations have thus far failed to produce any solution where such a field remains confined against outward diffusion. We argue that the problem lies in the choice of parameters for which these numerical simulations have been performed. We construct a simple analytical magneto-hydrodynamic model of the solar interior and identify several distinct parameter regimes. For realistic solar parameter values, our results are in broad agreement with the tachocline model of Gough & McIntyre. In this regime, meridional flows driven at the base of the convection zone are of sufficient amplitude to hold back the interior magnetic field against diffusion. For the parameter values used in existing numerical simulations, on the other hand, we find that meridional flows are significantly weaker and, we argue, unable to confine the interior field. We propose a method for selecting parameter values in future numerical models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A model of the entropy flux and Reynolds stress in turbulent convection

We propose a closure model for the transport of entropy and momentum in astrophysical turbulence, intended for application to rotating stellar convective regions. Our closure model is first presented in the Boussinesq formalism, and compared with laboratory and numerical experimental results on Rayleigh-Benard convection and Homogeneous Rayleigh-Benard convection. The predicted angular momentum transport properties of the turbulence in the slowly rotating case recover the well-known Lambda-effect, with an amplitude uniquely related to the convective heat flux. The model is then extended to the anelastic case as well as the fully compressible case. In the special case of spherical symmetry, the predicted radial heat flux is equivalent to that of mixing-length theory. For rotating stars, our model describes the coupled transport of heat and angular momentum, and provides a unified formalism in which to study both differential rotation and thermal inhomogeneities in stellar convection zones.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gyroscopic pumping of large-scale flows in stellar interiors, and application to Lithium Dip stars

The maintenance of large-scale differential rotation in stellar convective regions by rotationally influenced convective stresses also drives large-scale meridional flows by angular--momentum conservation. This process is an example of ``gyroscopic pumping'', and has recently been studied in detail in the solar context. An important question concerns the extent to which these gyroscopically pumped meridional flows penetrate into nearby stably stratified (radiative) regions, since they could potentially be an important source of non-local mixing. Here we present an extensive study of the gyroscopic pumping mechanism, using a combination of analytical calculations and numerical simulations both in Cartesian geometry and in spherical geometry. The various methods, when compared with one another, provide physical insight into the process itself, as well as increasingly sophisticated means of estimating the gyroscopic pumping rate. As an example of application, we investigate the effects of this large-scale mixing process on the surface abundances of the light elements Li and Be for stars in the mass range 1.3-1.5 solar masses (so-called ``Li-dip stars''). We find that gyroscopic pumping is a very efficient mechanism for circulating material between the surface and the deep interior, so much in fact that it over-estimates Li and Be depletion by orders of magnitude for stars on the hot side of the dip.However, when the diffusion of chemical species back into the surface convection zone is taken into account, a good fit with observed surface abundances of Li and Be as a function of stellar mass in the Hyades cluster can be found for reasonable choices of model parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Numerically determined transport laws for fingering ("thermohaline") convection in astrophysics

We present the first three-dimensional simulations of fingering convection performed in a parameter regime close to the one relevant for astrophysics, and reveal the existence of simple asymptotic scaling laws for turbulent heat and compositional transport. These laws can straightforwardly be extrapolated to the true astrophysical regime. Our investigation also indicates that thermocompositional "staircases," a key consequence of fingering convection in the ocean, cannot form spontaneously in the fingering regime in stellar interiors. Our proposed empirically-determined transport laws thus provide simple prescriptions for mixing by fingering convection in a variety of astrophysical situations, and should, from here on, be used preferentially over older and less accurate parameterizations. They also establish that fingering convection does not provide sufficient extra mixing to explain observed chemical abundances in RGB stars.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Turbulent mixing and layer formation in double-diffusive convection: 3D numerical simulations and theory

Double-diffusive convection, often referred to as semi-convection in astrophysics, occurs in thermally and compositionally stratified systems which are stable according to the Ledoux-criterion but unstable according to the Schwarzchild criterion. This process has been given relatively little attention so far, and its properties remain poorly constrained. In this paper, we present and analyze a set of three-dimensional simulations of this phenomenon in a Cartesian domain under the Boussinesq approximation. We find that in some cases the double-diffusive convection saturates into a state of homogeneous turbulence, but with turbulent fluxes several orders of magnitude smaller than those expected from direct overturning convection. In other cases the system rapidly and spontaneously develops closely-packed thermo-compositional layers, which later successively merge until a single layer is left. We compare the output of our simulations with an existing theory of layer formation in the oceanographic context, and find very good agreement between the model and our results. The thermal and compositional mixing rates increase significantly during layer formation, and increase even further with each merger. We find that the heat flux through the staircase is a simple function of the layer height. We conclude by proposing a new approach to studying transport by double-diffusive convection in astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

What happened to the other Mohicans? Realistic models of metallicity dilution by fingering convection and observational implications

When a planet falls onto the surface of its host star, the added high-metallicity material does not remain in the surface layers, as often assumed, but is diluted into the interior through fingering (thermohaline) convection. Until now, however, the timescale over which this process happens remained very poorly constrained. Using recently-measured turbulent mixing rates for fingering convection, I provide reliable numerical and semi-analytical estimates for the rate at which the added heavy elements drain into the interior. I find that the relative metallicity enhancement post-infall drops by a factor of ten over a timescale which depends only on the structure of the host star, and decreases very rapidly with increasing stellar mass (from about 1Gyr for a 1.3M_sun star to 10Myr for a 1.5M_sun star). This result offers an elegant explanation to the lack of observed trend between metallicity and convection zone mass in planet-bearing stars. More crucially, it strongly suggests that the statistically-higher metallicity of planet-bearing stars must be of primordial origin. Finally, the fingering region is found to extend deeply into the star, a result which would provide a simple theoretical explanation to the measurements of higher lithium depletion rates in planet-bearing stars.

0 citations0 reviewsNo code linkedOpen access