Source author record

Michel Rieutord

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2021arXiv

Predictions for gravity-mode periods and surface abundances in intermediate-mass dwarfs from shear mixing and radiative levitation

The treatment of chemical mixing in the radiative envelopes of intermediate-mass stars has hardly been calibrated so far. Recent asteroseismic studies demonstrated that a constant diffusion coefficient in the radiative envelope is not able to explain the periods of trapped gravity modes in the oscillation spectra of $γ$ Doradus pulsators. We present a new generation of MESA stellar models with two major improvements. First, we present a new implementation for computing radiative accelerations and Rosseland mean opacities that requires significantly less CPU time. Second, the inclusion of shear mixing based on rotation profiles computed with the 2D stellar structure code ESTER is considered. We show predictions for the mode periods of these models covering stellar masses from 1.4 to 3.0${\rm M_\odot}$ across the main sequence (MS), computed for different metallicities. The morphology of the chemical mixing profile resulting from shear mixing in combination with atomic diffusion and radiative levitation does allow for mode trapping, while the diffusion coefficient in the outer envelope is large ($>10^{6}\,{\rm cm^2\,s^{-1}}$). Furthermore, we make predictions for the evolution of surface abundances for which radiative accelerations can be computed. We find that the N/C and C/O abundance ratios correlate with stellar age. We predict that these correlations are observable with precisions $\lesssim 0.1$ dex on these ratios, given that a precise age estimate can be made.

preprint2017arXiv

Gravity darkening in late-type stars. The Coriolis effect

Recent interferometric data have been used to constrain the brightness distribution at the surface of nearby stars, in particular the so-called gravity darkening that makes fast rotating stars brighter at their poles than at their equator. However, good models of gravity darkening are missing for stars that posses a convective envelope. In order to better understand how rotation affects the heat transfer in stellar convective envelopes, we focus on the heat flux distribution in latitude at the outer surface of numerical models. We carry out a systematic parameter study of three-dimensional, direct numerical simulations of anelastic convection in rotating spherical shells. Restricting our investigations to hydrodynamical models with a thermal Prandtl number fixed to unity, we consider both thick and thin (solar-like) shells, and vary the stratification over three orders of magnitude. We measure the heat transfer efficiency in terms of the Nusselt number, defined as the output luminosity normalised by the conductive state luminosity. We report diverse Nusselt number profiles in latitude, ranging from brighter (usually at the onset of convection) to darker equator and uniform profiles. We find that the variations of the surface brightness are mainly controlled by the surface value of the local Rossby number: when the Coriolis force dominates the dynamics, the heat flux is weakened in the equatorial region by the zonal wind and enhanced at the poles by convective motions inside the tangent cylinder. In the presence of a strong background density stratification however, as expected in real stars, the increase of the local Rossby number in the outer layers leads to uniformisation of the surface heat flux distribution.

preprint2016arXiv

2D dynamics of the radiation zone of low mass stars

The internal rotation of low mass stars all along their evolution is of primary interest when studying their rotational dynamics, internal mixing and magnetic fields generation. In this context, helio- and asteroseismology probe angular velocity gradients deep within solar type stars. Still the rotation of the close center of such stars on the main sequence is hardly detectable and the dynamical interactions of the radiative core with the surface convective envelope is not well understood. Among them, the influence of the differential rotation profile sustained by convection and applied as a boundary condition to the radiation zone may be very important leading to the formation of tachoclines. In the solar convective region, the equator is rotating faster than the pole while numerical simulations predict either a solar or an anti-solar rotation in other low mass stars envelopes depending on their convective Rossby number. In this work, we therefore build for the first time 2D steady hydrodynamical models of low mass stars radiation zone providing a full 2D description of their dynamics and studying the influence of a general shear boundary condition accounting for a solar or an anti-solar differential rotation in the convective envelope. We compute coherently differential rotation and the associated meridional circulation using the Boussinesq approximation.

preprint2016arXiv

Gravito-inertial waves in a differentially rotating spherical shell

The gravito-inertial waves propagating over a shellular baroclinic flow inside a rotating spherical shell are analysed using the Boussinesq approximation. The wave properties are examined by computing paths of characteristics in the non-dissipative limit, and by solving the full dissipative eigenvalue problem using a high-resolution spectral method. Gravito-inertial waves are found to obey a mixed-type second-order operator and to be often focused around short-period attractors of characteristics or trapped in a wedge formed by turning surfaces and boundaries. We also find eigenmodes that show a weak dependence with respect to viscosity and heat diffusion just like truly regular modes. Some axisymmetric modes are found unstable and likely destabilized by baroclinic instabilities. Similarly, some non-axisymmetric modes that meet a critical layer (or corotation resonance) can turn unstable at sufficiently low diffusivities. In all cases, the instability is driven by the differential rotation. For many modes of the spectrum, neat power laws are found for the dependence of the damping rates with diffusion coefficients, but the theoretical explanation for the exponent values remains elusive in general. The eigenvalue spectrum turns out to be very rich and complex, which lets us suppose an even richer and more complex spectrum for rotating stars or planets that own a differential rotation driven by baroclinicity.

preprint2016arXiv

Tidal dissipation by inertial waves in differentially rotating convective envelopes of low-mass stars

Tidal interactions in close star-planet or binary star systems may excite inertial waves (their restoring force is the Coriolis force) in the convective region of the stars. The dissipation of these waves plays a prominent role in the long-term orbital and rotational evolution of the bodies involved. If the primary star rotates as a solid body, inertial waves have a Doppler-shifted frequency restricted to the range $[-2Ω, 2Ω]$ ($Ω$ being the angular velocity of the star), and they can propagate in the entire convective region. However, turbulent convection can sustain differential rotation with an equatorial acceleration (as in the Sun) or deceleration that modifies the frequency range and propagation domain of inertial waves and allows corotation resonances for non-axisymmetric oscillations. In this work, we perform numerical simulations of tidally excited inertial waves in a differentially rotating convective envelope with a conical (or latitudinal) rotation profile. The tidal forcing that we adopt contains spherical harmonics that correspond to the case of a circular and coplanar orbit. We study the viscous dissipation of the waves as a function of tidal frequency for various stellar masses and differential rotation parameters, as well as its dependence on the turbulent viscosity coefficient. We compare our results with previous studies assuming solid-body rotation and point out the potential key role of corotation resonances in the dynamical evolution of close-in star-planet or binary systems.

preprint2015arXiv

Physical processes leading to surface inhomogeneities: the case of rotation

In this lecture I discuss the bulk surface heterogeneity of rotating stars, namely gravity darkening. I especially detail the derivation of the omega-model of Espinosa Lara & Rieutord (2011), which gives the gravity darkening in early-type stars. I also discuss the problem of deriving gravity darkening in stars owning a convective envelope and in those that are members of a binary system.

preprint2015arXiv

Two-dimensional models of early-type fast rotating stars: the ESTER project

In this talk I present the latest results of the ESTER project that has taken up the challenge of building two dimensional (axisymmetric) models of stars rotating at any rotation rate. In particular, I focus on main sequence massive and intermediate mass stars. I show what should be expected in such stars as far as the differential rotation and the associated meridional circulation are concerned, notably the emergence of a Stewartson layer along the tangent cylinder of the core. I also indicate what may be inferred about the evolution of an intermediate-mass star at constant angular momentum and how Be stars may form. I finally give some comparisons between models and observations of the gravity darkening on some nearby fast rotators as it has been derived from interferometric observations. In passing, I also discuss how 2D models can help to recover the fundamental parameters of a star.

preprint2014arXiv

Dynamics of the envelope of a rapidly rotating star or giant planet in gravitational contraction

We wish to understand the processes that control the fluid flows of a gravitationally contracting and rotating star or giant planet. We consider a spherical shell containing an incompressible fluid that is slowly absorbed by the core so as to mimick gravitational contraction. We also consider the effects of a stable stratification that may also modify the dynamics of a pre-main sequence star of intermediate mass. This simple model reveals the importance of both the Stewartson layer attached to the core and the boundary conditions met by the fluid at the surface of the object. In the case of a pre-main sequence star of intermediate mass where the envelope is stably stratified, shortly after the birth line, the spin-up flow driven by contraction overwhelms the baroclinic flow that would take place otherwise.This model also shows that for a contracting envelope, a self-similar flow of growing amplitude controls the dynamics. It suggests that initial conditions on the birth line are most probably forgotten. Finally, the model shows that the near (Stewartson) layer that lies on the tangent cylinder of the core is likely a key feature of the dynamics that is missing in 1D models.This layer can explain the core and envelope rotational coupling that is required to explain the slow rotation of cores in giant and subgiants stars.

preprint2014arXiv

Dynamics of the radiative envelope of rapidly rotating stars: Effects of spin-down driven by mass loss

(abridged) This paper aims at deciphering the dynamics of the envelope of a rotating star when some angular momentum loss due to mass loss is present. We especially wish to know when the spin-down flow forced by the mass loss supersedes the baroclinic flows that pervade the radiative envelope of rotating stars. We consider a Boussinesq fluid enclosed in a rigid sphere whose flows are forced both by the baroclinic torque, the spin-down of an outer layer, and an outward mass flux. The spin-down forcing is idealized in two ways: either by a rigid layer that imposes its spinning down velocity at some interface or by a turbulent layer that imposes a stress at this same interface to the interior of the star. In the case where the layer is rigid and imposes its velocity, we find that, as the mass-loss rate increases, the flow inside the star shows two transitions: the meridional circulation associated with baroclinic flows is first replaced by its spin-down counterpart, while at much stronger mass-loss rates the baroclinic differential rotation is superseded by the spin-down differential rotation. In fact, we find three wind regimes: weak (or no wind), moderate, and strong. In the first case, the flow in the radiative envelope is of baroclinic origin. In the moderate case, the circulation results from the spin-down while the differential rotation may either be of baroclinic or of spin-down origin, depending on the coupling between mass and angular momentum losses. For fast rotating stars, our model says that the moderate wind regime starts when mass loss is higher than ~1e-11 Msun/yr. In the strong wind case, the flow in the radiative envelope is mainly driven by angular momentum advection. This latter transition depends on the mass and the rotation rate of the star, being around 1e-8 Msun/yr for a 3 Msun ZAMS star rotating at 200 km/s according to our model.

preprint2014arXiv

Gravito-inertial modes in a differentially rotating spherical shell

While many intermediate- and high-mass main sequence stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study of axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. We consider a simplified model where the radiative zone of the star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral computations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We describe various geometries for the propagation domains, conditioning the surface visibility of the corresponding modes.

preprint2014arXiv

Gravito-inertial modes in a differentially rotating spherical shell

Oscillations have been detected in a variety of stars, including intermediate- and high-mass main sequence stars. While many of these stars are rapidly and differentially rotating, the effects of rotation on oscillation modes are poorly known. In this communication we present a first study on axisymmetric gravito-inertial modes in the radiative zone of a differentially rotating star. These modes probe the deep layers of the star around its convective core. We consider a simplified model where the radiative zone of a star is a linearly stratified rotating fluid within a spherical shell, with differential rotation due to baroclinic effects. We solve the eigenvalue problem with high-resolution spectral simulations and determine the propagation domain of the waves through the theory of characteristics. We explore the propagation properties of two kinds of modes: those that can propagate in the entire shell and those that are restricted to a subdomain. Some of the modes that we find concentrate kinetic energy around short-period shear layers known as attractors. We characterise these attractors by the dependence of their Lyapunov exponent with the \BV frequency of the background and the oscillation frequency of the mode. Finally, we note that, as modes associated with short-period attractors form dissipative structures, they could play an important role for tidal interactions but should be dismissed in the interpretation of observed oscillation frequencies.

preprint2014arXiv

Magnetohydrodynamics and Solar Physics

In this short review, I present some of the recent progresses on the pending questions of solar physics. These questions let us revisit the solar wind, the solar dynamo problem, the dynamics of the photosphere and finally have a glimpse at other solar type stars. Discussing the use of direct numerical simulations in solar physics, I show that the full numerical calculation of the flow in a single supergranule would require more electric power than the luminosity of the sun itself with present computer technology.

preprint2013arXiv

Asteroseismology of fast-rotating stars: the example of alpha Ophiuchi

Many early-type stars have been measured with high angular velocities. In such stars, mode identification is difficult as the effects of fast and differential rotation are not well known. Using fundamental parameters measured by interferometry, the ESTER structure code and the TOP oscillation code, we investigate the oscillation spectrum of Rasalhague (alpha Ophiuchi), for which observations by the MOST satellite found 57 oscillations frequencies. Results do not show a clear identification of the modes and highlight the difficulties of asteroseismology for such stars with a very complex oscillation spectrum.

preprint2013arXiv

Pulsations of rapidly rotating stars with compositional discontinuities

Recent observations of rapidly rotating stars have revealed the presence of regular patterns in their pulsation spectra. This has raised the question as to their physical origin, and in particular, whether they can be explained by an asymptotic frequency formula for low-degree acoustic modes, as recently discovered through numerical calculations and theoretical considerations. In this context, a key question is whether compositional/density gradients can adversely affect such patterns to the point of hindering their identification. To answer this question, we calculate frequency spectra using two-dimensional ESTER stellar models. These models use a multi-domain spectral approach, allowing us to easily insert a compositional discontinuity while retaining a high numerical accuracy. We analyse the effects of such discontinuities on both the frequencies and eigenfunctions of pulsation modes in the asymptotic regime. We find that although there is more scatter around the asymptotic frequency formula, the semi-large frequency separation can still be clearly identified in a spectrum of low-degree acoustic modes.

preprint2013arXiv

Two-dimensional models of early-type fast rotating stars: new challenges in stellar physics

Two-dimensional models of rapidly rotating stars are already unavoidable for the interpretation of interferometric or asteroseismic data of this kind of stars. When combined with time evolution, they will allow the including of a more accurate physics for the computation of element transport and the determination of surface abundances. In addition, modeling the evolution of rotation will improve gyrochronology. Presently, two-dimensional ESTER models predict the structure and the large-scale flows (differential rotation and meridional circulation) of stars with mass larger than 1.7Msun at any rotation rate. Main sequence evolution can be mimicked by varying the hydrogen content of the convective core. Models have been successfully tested on half a dozen of nearby fast rotating stars observed with optical or infra-red interferometers. They are now the right tool to investigate the oscillation spectrum of early-type fast rotators.

preprint2012arXiv

Ab initio modelling of steady rotating stars

Modelling isolated rotating stars at any rotation rate is a challenge for the next generation of stellar models. These models will couple dynamical aspects of rotating stars, like angular momentum and chemicals transport, with classical chemical evolution, gravitational contraction or mass-loss. Such modelling needs to be achieved in two dimensions, combining the calculation of the structure of the star, its mean flows and the time-evolution of the whole. We present here a first step in this challenging programme. It leads to the first self-consistent two-dimensional models of rotating stars in a steady state generated by the ESTER code. In these models the structure (pressure, density and temperature) and the flow fields are computed in a self-consistent way allowing the prediction of the differential rotation and the associated meridian circulation of the stars. After a presentation of the physical properties of such models and the numerical methods at work, we give the first grid of such models describing massive and intermediate-mass stars for a selection of rotation rates up to 90% of the breakup angular velocity.

preprint2012arXiv

Gravity Darkening in Binary Stars

Context.Interpretation of light curves of many types of binary stars requires the inclusion of the (cor)relation between surface brightness and local effective gravity. Until recently, this correlation has always been modeled by a power law relating the flux or the effective temperature and the effective gravity, namely T_eff α g_eff^β. Aims. We look for a simple model that can describe the variations of the flux at the surface of stars belonging to a binary system. Methods. This model assumes that the energy flux is a divergence-free vector anti-parallel to the effective gravity. The effective gravity is computed from the Roche model. Results. After explaining in a simple manner the old result of Lucy (1967), which says that β=0.08 for solar type stars, we first argue that one-dimensional models should no longer be used to evaluate gravity darkening laws. We compute the correlation between log T_eff and log g_eff using a new approach that is valid for synchronous, weakly magnetized, weakly irradiated binaries. We show that this correlation is approximately linear, validating the use of a power law relation between effective temperature and effective gravity as a first approximation. We further show that the exponent β of this power law is a slowly varying function, which we tabulate, of the mass ratio of the binary star and the Roche lobe filling factor of the stars of the system. The exponent β remains mostly in the interval (0.20, 0.25) if extreme mass ratios are eliminated. Conclusions. For binary stars that are synchronous, weakly magnetized and weakly irradiated, the gravity darkening exponent is well constrained and may be removed from the free parameters of the models.

preprint2012arXiv

On the excitation of inertial modes in an experimental spherical Couette flow

Spherical Couette flow (flow between concentric rotating spheres) is one of flows under consideration for the laboratory magnetic dynamos. Recent experiments have shown that such flows may excite Coriolis restored inertial modes. The present work aims to better understand the properties of the observed modes and the nature of their excitation. Using numerical solutions describing forced inertial modes of a uniformly rotating fluid inside a spherical shell, we first identify the observed oscillations of the Couette flow with non-axisymmetric, retrograde, equatorially anti-symmetric inertial modes, confirming first attempts using a full sphere model. Although the model has no differential rotation, identification is possible because a large fraction of the fluid in a spherical Couette flow rotates rigidly. From the observed sequence of the excited modes appearing when the inner sphere is slowed down by step, we identify a critical Rossby number associated with a given mode and below which it is excited. The matching between this critical number and the one derived from the phase velocity of the numerically computed modes shows that these modes are excited by an instability likely driven by the critical layer that develops in the shear layer staying along the tangent cylinder of the inner sphere.

preprint2011arXiv

Gravity darkening in rotating stars

Interpretation of interferometric observations of rapidly rotating stars requires a good model of their surface effective temperature. Until now, laws of the form T_eff \propto g_eff^β have been used, but they are only valid for slowly rotating stars. We propose a simple model that can describe the latitudinal variations in the flux of rotating stars at any rotation rate. This model assumes that the energy flux is a divergence-free vector that is antiparallel to the effective gravity. When mass distribution can be described by a Roche model, the latitudinal variations in the effective temperature only depend on a single parameter, namely the ratio of the equatorial velocity to the Keplerian velocity. We validate this model by comparing its predictions to those of the most realistic two-dimensional models of rotating stars issued from the ESTER code. The agreement is very good, as it is with the observations of two rapidly rotating stars, α Aql and α Leo. We suggest that as long as a gray atmosphere can be accepted, the inversion of data on flux distribution coming from interferometric observations of rotating stars uses such a model, which has just one free parameter.