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.Comment: 13 pages, 9 figures, to appear in Astronomy and Astrophysic

The Sun's Supergranulation

Supergranulation is a fluid-dynamical phenomenon taking place in the solar photosphere, primarily detected in the form of a vigorous cellular flow pattern with a typical horizontal scale of approximately 30--35~megameters, a dynamical evolution time of 24--48~h, a strong 300--400~m/s (rms) horizontal flow component and a much weaker 20--30~m/s vertical component. Supergranulation was discovered more than sixty years ago, however, explaining its physical origin and most important observational characteristics has proven extremely challenging ever since, as a result of the intrinsic multiscale, nonlinear dynamical complexity of the problem concurring with strong observational and computational limitations. Key progress on this problem is now taking place with the advent of 21st-century supercomputing resources and the availability of global observations of the dynamics of the solar surface with high spatial and temporal resolutions. This article provides an exhaustive review of observational, numerical and theoretical research on supergranulation, and discusses the current status of our understanding of its origin and dynamics, most importantly in terms of large-scale nonlinear thermal convection, in the light of a selection of recent findings.Comment: Major update of 2010 Liv. Rev. Sol. Phys. review. Addresses many new theoretical, numerical and observational developments. All sections, including discussion, revised extensively. Also includes previously unpublished results on nonlinear dynamics of convection in large domains, and lagrangian transport at the solar surfac

Completeness of Inertial Modes of an Incompressible Non-Viscous Fluid in a Corotating Ellipsoid

Inertial modes are the eigenmodes of contained rotating fluids restored by the Coriolis force. When the fluid is incompressible, inviscid and contained in a rigid container, these modes satisfy Poincar\'e's equation that has the peculiarity of being hyperbolic with boundary conditions. Inertial modes are therefore solutions of an ill-posed boundary-value problem. In this paper we investigate the mathematical side of this problem. We first show that the Poincar\'e problem can be formulated in the Hilbert space of square-integrable functions, with no hypothesis on the continuity or the differentiability of velocity fields. We observe that with this formulation, the Poincar\'e operator is bounded and self-adjoint and as such, its spectrum is the union of the point spectrum (the set of eigenvalues) and the continuous spectrum only. When the fluid volume is an ellipsoid, we show that the inertial modes form a complete base of polynomial velocity fields for the square-integrable velocity fields defined over the ellipsoid and meeting the boundary conditions. If the ellipsoid is axisymmetric then the base can be identified with the set of Poincar\'e modes, first obtained by Bryan (1889), and completed with the geostrophic modes.Comment: 19 pages, 1 figure, to appear in Physical Review

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.Comment: 10 pages, to appear in the proceedings of the conference "New advances in stellar physics: from microscopic to macroscopic processes", Roscoff, May 201

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 {\alpha} g_eff^{\beta}. 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 {\beta}=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 {\beta} 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 {\beta} 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

2D dynamics of the radiative core of low mass stars

Understanding the internal rotation of low mass stars all along their evolution is of primary interest when studying their rotational dynamics, internal mixing and magnetic field generation. In this context, helio- and asteroseismology probe angular velocity gradients deep within solar type stars at different evolutionary stages. Still the rotation close to the center of such stars on the main sequence is hardly detectable and the dynamical interaction of the radiative core with the surface convective envelope is not well understood. For instance, the influence of the differential rotation profile sustained by convection and applied as a boundary condition to the radiation zone is very important in the formation of tachoclines. In this work, we study a 2D hydrodynamical model of a radiative core when an imposed, solar or anti-solar, differential rotation is applied at the upper boundary. This model uses the Boussinesq approximation and we find that the shear induces a cylindrical differential rotation associated with a unique cell of meridional circulation in each hemisphere (counterclockwise when the shear is solar-like and clockwise when it is anti-solar). The results are discussed in the framework of seismic observables (internal rotation rate, core-to-surface rotation ratio) while perspectives to improve our modeling by including magnetic field or transport by internal gravity waves will be discussed.Comment: 5 pages, 4 figures. To appear in the proceedings of "Seismology of the Sun and the Distant Stars 2016, Joint TASC2 & KASC9 Workshop - SPACEINN & HELAS8 Conference" (ed. M\'ario J. P. F. G. Monteiro, Margarida S. Cunha, Jo\~ao Miguel T. Ferreira ), Azores, Portugal, 11-15 July 201