We investigate the linear properties of the steady stress-driven spin-down flow of a viscous fluid inside a spherical shell, both within the incompressible and anelastic approximations, and in the asymptotic limit of small viscosities. Neglecting buoyancy force, we find a steady-state to settle in a viscous time for the incompressible flow. From boundary layer analysis, we derive an analytical geostrophic solution for the 3D incompressible steady flow, inside and outside a cylinder $\mathcal{C}$ that is tangent to the inner shell. The meridional circulation associated with the spin-down flow in these two regions scales with the Ekman number $E$. A power expansion of the fields with the small parameter $E^{1/3}$ corresponding to the width of the inner Stewartson shear layer laying along the tangent cylinder, reveals it to be the site of an upwelling flow of amplitude $O(E^{1/3})$. Despite its narrowness, this shear layer concentrates most of the global meridional kinetic energy of the system. Furthermore, a stable stratification does not perturb the spin-down flow provided the Prandtl number is small enough. If this is not the case, the Stewartson layer disappears and meridional circulation is confined within the thermal layers. The scalings for the amplitude of the anelastic secondary flow have been found to be the same as for the incompressible flow in all three regions, at the lowest order. However, because the velocity no longer conforms the Taylor-Proudman theorem, its shape differs outside the tangent cylinder $\mathcal{C}$, that is, where differential rotation takes place. Finally, we find the settling of the steady-state to be reached on a viscous time for the weakly, strongly and thermally unstratified incompressible flows. Large density variations relevant to astro- and geophysical systems, tend to slightly shorten the transient.

Fast rotation is responsible for important changes in the structure and evolution of stars. Optical long baseline interferometry now permits the study of its effects on the stellar surface, mainly gravity darkening and flattening. We aim to determine the fundamental parameters of the fast-rotating star Altair, in particular its evolutionary stage, mass, and differential rotation, using state-of-the-art stellar interior and atmosphere models together with interferometric, spectroscopic, and asteroseismic observations. We use ESTER 2D stellar models to produce the relevant surface parameters needed to create intensity maps from atmosphere models. Interferometric and spectroscopic observables are computed from these intensity maps and several stellar parameters are then adjusted using the MCMC algorithm Emcee. We determined Altair's equatorial radius to be 2.008 +/- 0.006 Rsun, the position angle 301.1 +/- 0.3 degrees, the inclination 50.7 +/- 1.2 degrees, and the equatorial angular velocity 0.74 +/- 0.01 times the Keplerian angular velocity. This angular velocity leads to a flattening of 0.220 +/- 0.003. We also deduce from the spectroscopically derived vsini ~ 243 km/s, a true equatorial velocity of ~314 km/s corresponding to a rotation period of 7h46m (~3 c/d). The data also impose a strong correlation between mass, metallicity, hydrogen abundance, and core evolution. Thanks to asteroseismic data, we constrain the mass of Altair to 1.86 +/- 0.03 Msun and further deduce its metallicity Z = 0.019 and its core hydrogen mass fraction Xc = 0.71, assuming an initial solar hydrogen mass fraction X = 0.739. These values suggest that Altair is ~100 Myrs old. Finally, the 2D ESTER model also gives the internal differential rotation of Altair, showing that its core rotates approximately 50% faster than the envelope, while the surface differential rotation does not exceed 6%.

We present the first asteroseismic results for $\delta$ Scuti and $\gamma$ Doradus stars observed in Sectors 1 and 2 of the TESS mission. We utilise the 2-min cadence TESS data for a sample of 117 stars to classify their behaviour regarding variability and place them in the Hertzsprung-Russell diagram using Gaia DR2 data. Included within our sample are the eponymous members of two pulsator classes, $\gamma$ Doradus and SX Phoenicis. Our sample of pulsating intermediate-mass stars observed by TESS also allows us to confront theoretical models of pulsation driving in the classical instability strip for the first time and show that mixing processes in the outer envelope play an important role. We derive an empirical estimate of 74% for the relative amplitude suppression factor as a result of the redder TESS passband compared to the Kepler mission using a pulsating eclipsing binary system. Furthermore, our sample contains many high-frequency pulsators, allowing us to probe the frequency variability of hot young $\delta$ Scuti stars, which were lacking in the Kepler mission data set, and identify promising targets for future asteroseismic modelling. The TESS data also allow us to refine the stellar parameters of SX Phoenicis, which is believed to be a blue straggler.

We consider a binary system composed of a pulsar and a massive, fast rotating, highly distorted main sequence star as a potential scenario to dynamically put to the test certain post-Keplerian effects of both Newtonian and post-Newtonian nature. We numerically produce time series of the perturbations $\Delta\left(\delta\tau\right)$ of the R{\o}mer-like, orbital component of the pulsar's time delay $\delta\tau$ induced over 10 years by the pN gravitoelectric mass monopole, quadrupole, gravitomagnetic spin dipole and octupole accelerations along with the Newtonian quadrupolar one. We do not deal with the various propagation time delays due to the travelling electromagnetic waves. It turns out that, for a Be-type star with $M = 15\ \textrm{M}_\odot$, $R_\textrm{e} = 5.96\ \textrm{R}_\odot$, $\nu = 0.203$, $S = 3.41\times 10^{45}\ \textrm{J}\ \textrm{s}$, $J_2 = 1.92\times 10^{-3}$ orbited by a pulsar with an orbital period $P_\textrm{b}\simeq 40-70\ \textrm{d}$, the classical oblateness-driven effects are at the $\lesssim 4-150\ \textrm{s}$ level, while the pN shifts are of the order of $\lesssim 1.5-20\ \textrm{s}\ \left(GMc^{-2}\right)$, $\lesssim 10-40\ \textrm{ms}\ \left(GMR^2_\textrm{e} J_2 c^{-2}\right)$, $\lesssim 0.5 - 6\ \textrm{ms}\ \left(GSc^{-2}\right)$, $\lesssim 5 - 20\ \mu\textrm{s}\ \left(GSR^2_\textrm{e} \varepsilon^2 c^{-2}\right)$, depending on their orbital configuration. The root-mean-square (rms) timing residuals $\sigma_{\tau}$ of almost all the existing non-recycled, non-millisecond pulsars orbiting massive, fast rotating main sequence stars are $\lesssim\textrm{ms}$. Thus, such kind of binaries have the potential to become interesting laboratories to measure, or, at least, constrain, some Newtonian and post-Newtonian key features of the distorted gravitational fields of the fast rotating stars hosted by them [Abridged].

The understanding of the evolution of early-type stars is tightly related to that of the effects of rapid rotation. For massive stars, rapid rotation combines with their strong radiation-driven wind. The aim of this paper is to investigate two questions that are prerequisite to the study of the evolution of massive rapidly rotating stars: (i) What is the critical angular velocity of a star when radiative acceleration is significant in its atmosphere? (ii) How do mass and angular momentum loss depend on the rotation rate? To investigate fast rotation, which makes stars oblate, we used the 2D ESTER models and a simplified approach, the $\omega$-model, which gives the latitudinal dependence of the radiative flux in a centrifugally flattened radiative envelope. We find that radiative acceleration only mildly influences the critical angular velocity, at least for stars with masses lower than 40 Msun. We explain this mild reduction of the critical angular velocity compared to the classical Keplerian angular velocity by the combined effects of gravity darkening and a reduced equatorial opacity that is due to the centrifugal acceleration. To answer the second question, we first devised a model of the local surface mass flux, which we calibrated with previously developed 1D models. The discontinuity (the so-called bi-stability jump) included in the $\dot{M}-T_{\rm eff}$ relation of 1D models means that the mass flux of a fast-rotating star is controlled by either a single wind or a two-wind regime. Mass and angular momentum losses are strong around the equator if the star is in the two-wind regime. We also show that the difficulty of selecting massive stars that are viewed pole-on makes detecting the discontinuity in the relation between mass loss and effective temperature also quite challenging.

The understanding of the rotational evolution of early-type stars is deeply related to that of anisotropic mass and angular momentum loss. In this paper, we aim to clarify the rotational evolution of rapidly rotating early-type stars along the main sequence (MS). We have used the 2D ESTER code to compute and evolve isolated rapidly rotating early-type stellar models along the MS, with and without anisotropic mass loss. We show that stars with $Z=0.02$ and masses between $5$ and $7~M_\odot$ reach criticality during the main sequence provided their initial angular velocity is larger than 50% of the Keplerian one. More massive stars are subject to radiation-driven winds and to an associated loss of mass and angular momentum. We find that this angular momentum extraction from the outer layers can prevent massive stars from reaching critical rotation and greatly reduce the degree of criticality at the end of the MS. Our model includes the so-called bi-stability jump of the $\dot{M}-T_{\rm eff}$ relation of 1D-models. This discontinuity now shows up in the latitude variations of the mass-flux surface density, endowing rotating massive stars with either a single-wind regime (no discontinuity) or a two-wind regime (a discontinuity). In the two-winds-regime, mass loss and angular momentum loss are strongly increased at low latitudes inducing a faster slow-down of the rotation. However, predicting the rotational fate of a massive star is difficult, mainly because of the non-linearity of the phenomena involved and their strong dependence on uncertain prescriptions. Moreover, the very existence of the bi-stability jump in mass-loss rate remains to be substantiated by observations.

Gravity darkening (GD) and flattening are important consequences of stellar rotation. The precise characterization of these effects across the HRD is crucial to a deeper understanding of stellar structure and evolution. We seek to characterize such important effects on Sargas, an evolved, fast-rotating, intermediate-mass star, located in a region of the HRD where they have never been directly measured as far as we know. We use our numerical model CHARRON to analyze interferometric (VLTI/PIONIER) and spectroscopic (VLT/UVES) observations through a MCMC model-fitting procedure. The visibilities and closure phases from the PIONIER data are particularly sensitive to rotational flattening and GD. Adopting the Roche approximation, we investigate two GD models: (1) the beta-model (classical von Zeipel's law), and (2) the omega-model. Using this approach we measure several physical parameters of Sargas, namely, equatorial radius, mass, equatorial rotation velocity, mean Teff, inclination and position angle of the rotation axis, and beta. In particular, we show that the measured beta leads to a surface flux distribution equivalent to the one given by the omega-model. Thanks to our results, we also show that Sargas is most probably located in a rare and interesting region of the H-R diagram: within the Hertzsprung gap and over the hot edge of the instability strip. These results show once more the power of optical/IR long-baseline interferometry, combined with high-resolution spectroscopy, to directly measure fast-rotation effects and stellar parameters, in particular GD. As was the case for a few fast rotators previously studied by interferometry, the omega-model provides a physically more profound description of Sargas' GD, without the need of a beta exponent.

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.

Apart from recent progress in Gamma Dor stars, identifying modes in rapidly rotating stars is a formidable challenge due to the lack of simple, easily identifiable frequency patterns. As a result, it is necessary to look to observational methods for identifying modes. Two popular techniques are spectroscopic mode identification based on line profile variations (LPVs) and photometric mode identification based on amplitude ratios and phase differences between multiple photometric bands. In this respect, the BRITE constellation is particularly interesting as it provides space-based multi-colour photometry. The present contribution describes the latest developments in obtaining theoretical predictions for amplitude ratios and phase differences for pulsation modes in rapidly rotating stars. These developments are based on full 2D non-adiabatic pulsation calculations, using models from the ESTER code, the only code to treat in a self-consistent way the thermal equilibrium of rapidly rotating stars. These predictions are then specifically applied to the BRITE photometric bands to explore the prospects of identifying modes based on BRITE observations.

We investigate the asymptotic properties of axisymmetric inertial modes propagating in a spherical shell when viscosity tends to zero. We identify three kinds of eigenmodes whose eigenvalues follow very different laws as the Ekman number $E$ becomes very small. First are modes associated with attractors of characteristics that are made of thin shear layers closely following the periodic orbit traced by the characteristic attractor. Second are modes made of shear layers that connect the critical latitude singularities of the two hemispheres of the inner boundary of the spherical shell. Third are quasi-regular modes associated with the frequency of neutral periodic orbits of characteristics. We thoroughly analyse a subset of attractor modes for which numerical solutions point to an asymptotic law governing the eigenvalues. We show that three length scales proportional to $E^{1/6}$, $E^{1/4}$ and $E^{1/3}$ control the shape of the shear layers that are associated with these modes. These scales point out the key role of the small parameter $E^{1/12}$ in these oscillatory flows. With a simplified model of the viscous Poincar\'e equation, we can give an approximate analytical formula that reproduces the velocity field in such shear layers. Finally, we also present an analysis of the quasi-regular modes whose frequencies are close to $\sin(\pi/4)$ and explain why a fluid inside a spherical shell cannot respond to any periodic forcing at this frequency when viscosity vanishes.

Large scale flows in the Sun play an important role in the dynamo process linked to the solar cycle. The important large scale flows are the differential rotation and the meridional circulation with an amplitude of km/s and few m/s , respectively. These flows also have a cycle related components, namely the torsional oscillations. Our attempt is to determine large-scale plasma flows on the solar surface by deriving horizontal flow velocities using the techniques of solar granule tracking, dopplergrams, and time distance helioseismology. Coherent structure tracking (CST) and time distance helioseismology were used to investigate the solar differential rotation and meridional circulation at the solar surface on a 30 day HMI SDO sequence. The influence of a large sunspot on these large scale flows with a specific 7 day HMI SDO sequence has been also studied. The large scale flows measured by the CST on the solar surface and the same flow determined from the same data with the helioseismology in the first 1 Mm below the surface are in good agreement in amplitude and direction. The torsional waves are also located at the same latitudes with amplitude of the same order. We are able to measure the meridional circulation correctly using the CST method with only three days of data and after averaging between + and -15 degrees in longitude. We conclude that the combination of CST and Doppler velocities allows us to detect properly the differential solar rotation and also smaller amplitude flows such as the meridional circulation and torsional waves. The results of our methods are in good agreement with helioseismic measurements.

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.

In the context of secular evolution, we describe the dynamics of the radiative core of low-mass stars to understand the internal transport of angular momentum in such stars which results in a solid rotation in the Sun from 0.7R_sun to 0.2R_sun and a weak radial core-envelope differential rotation in solar-type stars. This study requires at least a 2D description to capture the latitudinal variations of the differential rotation. We build 2D numerical models of a radiative core on the top of which we impose a latitudinal shear so as to reproduce a cylindrical differential rotation in a convective envelope. We perform a systematic study over the Rossby number measuring the latitudinal differential rotation at the radiative-convective interface. The imposed shear generates a geostrophic flow implying a cylindrical differential rotation. When compared to the baroclinic flow that arises from the stable stratification, we find that the geostrophic flow is dominant when the Rossby number is high enough with a cylindrical rotation profile. For low Rossby numbers, the baroclinic solution dominates with a quasi-shellular rotation profile. Using scaling laws from 3D simulations, we show that slow rotators are expected to have a cylindrical rotation profile. Fast rotators may have a shellular profile at the beginning of the main-sequence in stellar radiative zones. This study enables us to predict different types of differential rotation and emphasizes the need of a new generation of 2D rotating stellar models developed in synergy with 3D numerical simulations. The shear induced by a surface convective zone has a strong impact on the dynamics of the underlying radiative zone in low-mass stars. But, it cannot produce a flat internal rotation profile in a solar configuration calling for additional processes for the transport of angular momentum in both radial and latitudinal directions.

Early-type stars generally tend to be fast rotators. In these stars, mode identification is very challenging as the effects of rotation are not well known. We consider here the example of $\alpha$ Ophiuchi, for which dozens of oscillation frequencies have been measured. We model the star using the two-dimensional structure code ESTER, and we compute both adiabatic and non-adiabatic oscillations using the TOP code. Both calculations yield very complex spectra, and we used various diagnostic tools to try and identify the observed pulsations. While we have not reached a satisfactory mode-to-mode identification, this paper presents promising early results.

One of the greatest challenges in interpreting the pulsations of rapidly rotating stars is mode identification, i.e. correctly matching theoretical modes to observed pulsation frequencies. Indeed, the latest observations as well as current theoretical results show the complexity of pulsation spectra in such stars, and the lack of easily recognisable patterns. In the present contribution, the latest results on non-adiabatic effects in such pulsations are described, and we show how these come into play when identifying modes. These calculations fully take into account the effects of rapid rotation, including centrifugal distortion, and are based on models from the ESTER project, currently the only rapidly rotating models in which the energy conservation equation is satisfied, a prerequisite for calculating non-adiabatic effects. Non-adiabatic effects determine which modes are excited and play a key role in the near-surface pulsation-induced temperature variations which intervene in multi-colour amplitude ratios and phase differences, as well as line profile variations.

The interpretation of stellar apparent fundamental parameters must be treated consistently with the characteristics of their surface rotation law. We develop a model to determine the distribution of the effective temperature and gravity, which depend on the surface differential rotation law and on the stellar external geometry. The basic assumptions are: a) the atmosphric layers are in radiative equilibrium; b) the bolometric flux is anti-parallel with the effective gravity; c) the angular velocity in the surface is Omega=Omega(theta,alpha,k)] where (alpha,k) are free parameters. The effective temperature varies with co-latitude `theta' and on the differential-rotation law through the effective gravity and the gravity-darkening function (GDF). Although the derived expressions can be treated numerically, for some low integer values of k, analytical forms of the `integral ofcharacteristic curves', on which the determination of the GDF relies, are obtained. The effects of the quantities (eta,alpha,k) (eta=ratio between centrifugal and gravitational accelerations at the equator) on the determination of the Vsini parameter and on the `gravity-darkening exponent' are studied. Depending on the values of (eta,alpha,k) the velocity V in the derived Vsini may strongly deviate from the equatorial rotational velocity. It is shown that the von Zeipel's-like gravity-darkening exponent $\beta_1$ depends on all parameters $(\eta,\alpha,k)$ and that its value also depends on the viewing-angle i. Hence, there no unique interpretation of this exponent determined empirically in terms of (i,alpha). The data on rotating stars should be analyzed by taking into account the rotational effects through the GDF, by assuming k=2 as a first approximation. Instead of the classical pair (eta,beta1), it would be more useful to determine the quantities (eta,alpha,i) to characterize stellar rotation.

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.

The Sun provides us with the only spatially well-resolved astrophysical example of turbulent thermal convection. While various aspects of solar photospheric turbulence, such as granulation (one-Megameter horizontal scale), are well understood, the questions of the physical origin and dynamical organization of larger-scale flows, such as the 30-Megameters supergranulation and flows deep in the solar convection zone, remain largely open in spite of their importance for solar dynamics and magnetism. Here, we present a new critical global observational characterization of multiscale photospheric flows and subsequently formulate an anisotropic extension of the Bolgiano-Obukhov theory of hydrodynamic stratified turbulence that may explain several of their distinctive dynamical properties. Our combined analysis suggests that photospheric flows in the horizontal range of scales between supergranulation and granulation have a typical vertical correlation scale of 2.5 to 4 Megameters and operate in a strongly anisotropic, self-similar, nonlinear, buoyant dynamical regime. While the theory remains speculative at this stage, it lends itself to quantitative comparisons with future high-resolution acoustic tomography of subsurface layers and advanced numerical models. Such a validation exercise may also lead to new insights into the asymptotic dynamical regimes in which other, unresolved turbulent anisotropic astrophysical fluid systems supporting waves or instabilities operate.

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\Omega, 2\Omega]$ ($\Omega$ 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.

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.

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.

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.

Stars may be understood as self-gravitating masses of a compressible fluid whose radiative cooling is compensated by nuclear reactions or gravitational contraction. The understanding of their time evolution requires the use of detailed models that account for a complex microphysics including that of opacities, equation of state and nuclear reactions. The present stellar models are essentially one-dimensional, namely spherically symmetric. However, the interpretation of recent data like the surface abundances of elements or the distribution of internal rotation have reached the limits of validity of one-dimensional models because of their very simplified representation of large-scale fluid flows. In this article, we describe the ESTER code, which is the first code able to compute in a consistent way a two-dimensional model of a fast rotating star including its large-scale flows. Compared to classical 1D stellar evolution codes, many numerical innovations have been introduced to deal with this complex problem. First, the spectral discretization based on spherical harmonics and Chebyshev polynomials is used to represent the 2D axisymmetric fields. A nonlinear mapping maps the spheroidal star and allows a smooth spectral representation of the fields. The properties of Picard and Newton iterations for solving the nonlinear partial differential equations of the problem are discussed. It turns out that the Picard scheme is efficient on the computation of the simple polytropic stars, but Newton algorithm is unsurpassed when stellar models include complex microphysics. Finally, we discuss the numerical efficiency of our solver of Newton iterations. This linear solver combines the iterative Conjugate Gradient Squared algorithm together with an LU-factorization serving as a preconditionner of the Jacobian matrix.

Context: The determination of the underlying mechanisms of the magnetic elements diffusion over the solar surface is still a challenge. Understanding the formation and evolution of the solar network (NE) is a challenge, because it provides a magnetic flux over the solar surface comparable to the flux of active regions at solar maximum. Aims: We investigate the structure and evolution of interior cells of solar supergranulation. From Hinode observations, we explore the motions on solar surface at high spatial and temporal resolution. We derive the main organization of the flows inside supergranules and their effect on the magnetic elements. Method: To probe the superganule interior cell, we used the Trees of Fragmenting Granules (TFG) evolution and their relations to horizontal Results: Evolution of TFG and their mutual interactions result in cumulative effects able to build horizontal coherent flows with longer lifetime than granulation (1 to 2 hours) over a scale up to 12\arcsec. These flows clearly act on the diffusion of the intranetwork (IN) magnetic elements and also on the location and shape of the network. Conclusions: From our analysis during 24 hours, TFG appear as one of the major elements of the supergranules which diffuse and advect the magnetic field on the Sun's surface. The strongest supergranules contribute the most to magnetic flux diffusion in the solar photosphere.

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. In low-mass stars, their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. Turbulent convection can sustain differential rotation in their envelope, with an equatorial acceleration (as in the Sun) or deceleration, which can modify the waves' propagation properties. We explore in this first paper the general propagation properties of free linear inertial waves in a differentially rotating homogeneous fluid inside a spherical shell. We assume that the angular velocity background flow depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We use i) an analytical approach in the inviscid case to get the dispersion relation, from which we compute the characteristic trajectories along which energy propagates. This allows us to study the existence of attractor cycles and infer the different families of inertial modes; ii) high-resolution numerical calculations based on a spectral method for the viscous problem. We find that modes that propagate in the whole shell (D modes) behave the same way as with solid-body rotation. However, another family of inertial modes exists (DT modes), which can propagate only in a restricted part of the convective zone. Our study shows that they are less common than D modes and that the characteristic rays and shear layers often focus towards a wedge - or point-like attractor. More importantly, we find that for non-axisymmetric oscillation modes, shear layers may cross a corotation resonance with a local accumulation of kinetic energy. Their damping rate scales very differently from what we obtain for standard D modes and we show an example where it is independent of viscosity (Ekman number) in the astrophysical regime in which it is small.

Tidally-excited inertial waves in stellar convective regions are a key mechanism for tidal dissipation in stars and therefore the evolution of close-in binary or planetary systems. As a first step, we explore here the impact of latitudinal differential rotation on the properties of free inertial modes and identify the different families of modes. We show that they differ from the case of solid-body rotation. Using an analytical approach as well as numerical calculations, we conclude that critical layers (where the Doppler-shifted frequency vanishes) could play a very important role for tidal dissipation.

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.

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.

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.

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.

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.

Star-planet tidal interactions may result in the excitation of inertial waves in the convective region of stars. Their dissipation plays a prominent role in the long-term orbital evolution of short-period planets. If the star is assumed to be rotating as a solid-body, the waves' Doppler-shifted frequency is restricted to $[-2 \Omega, 2 \Omega]$ ($\Omega$ 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 may modify waves propagation. We thus explore the properties of inertial modes of oscillation in a conically differentially rotating background flow whose angular velocity depends on the latitudinal coordinate only, close to what is expected in the external convective envelope of low-mass stars. We find that their frequency range is broadened by differential rotation, and that they may propagate only in a restricted part of the envelope. In some cases, inertial waves form shear layers around short-period attractor cycles. In others, they exhibit a remarkable behavior when a turning surface or a corotation layer exists in the star. We discuss how all these cases can impact tidal dissipation in stars.

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.

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

Context: Studying the motions on the solar surface is fundamental for understanding how turbulent convection transports energy and how magnetic fields are distributed across the solar surface. Aims: From horizontal velocity measurements all over the visible disc of the Sun and using data from the Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI), we investigate the structure and evolution of solar supergranulation. Methods: Horizontal velocity fields were measured by following the proper motions of solar granules using a newly developed version of the coherent structure tracking (CST) code. With this tool, maps of horizontal divergence were computed. We then segmented and identified supergranular cells and followed their histories by using spatio-temporal labelling. With this dataset we derived the fundamental properties of supergranulation, including their motion. Results: We find values of the fundamental parameters of supergranulation similar to previous studies: a mean lifetime of 1.5 days and a mean diameter of 25~Mm. The tracking of individual supergranular cells reveals the solar differential rotation and a poleward circulation trend of the meridional flow. The shape of the derived differential rotation and meridional flow does not depend on the cell size. If there is a background magnetic field, the diverging flows in supergranules are weaker. Conclusions: This study confirms that supergranules are suitable tracers that may be used to investigate the large-scale flows of the solar convection as long as they are detectable enough on the surface.

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.

ESTER two-dimensional models solve the steady state structure of fast rotating early-type stars including the large scale flows associated with the baroclinicity of the radiative zones. Models are compared successfully to the fundamental parameters of the two main components of the triple system $\delta$ Velorum that have been derived from interferometric and orbit measurements. Testing the models on the Be star Achernar ($\alpha$ Eri), we cannot reproduce the data and conclude that this star has left the main sequence and is likely crossing the Herzsprung gap. Computing main sequence evolution of fast rotating stars at constant angular momentum shows that their criticality increases with time suggesting that the Be phenomenon and the ensuing mass ejections is the result of evolution.

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.

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.

We compare measurements of horizontal flows on the surface of the Sun using helioseismic time--distance inversions and coherent structure tracking of solar granules. Tracking provides 2D horizontal flows on the solar surface, whereas the time--distance inversions estimate the full 3-D velocity flows in the shallow near-surface layers. Both techniques use HMI observations as an input. We find good correlations between the various measurements resulting from the two techniques. Further, we find a good agreement between these measurements and the time-averaged Doppler line-of-sight velocity, and also perform sanity checks on the vertical flow that resulted from the 3-D time--distance inversion.

In this lecture I present the way stars can be modeled in two dimensions and especially the fluid flows that are driven by rotation. I discuss some of the various ways of taking into account turbulence and conclude this contribution by a short presentation of some of the first results obtained with the ESTER code on the modeling of interferometrically observed fast rotating early-type stars.

This work aims at presenting the first two-dimensional models of an isolated rapidly rotating star that include the derivation of the differential rotation and meridional circulation in a self-consistent way.We use spectral methods in multidomains, together with a Newton algorithm to determine the steady state solutions including differential rotation and meridional circulation for an isolated non-magnetic, rapidly rotating early-type star. In particular we devise an asymptotic method for small Ekman numbers (small viscosities) that removes the Ekman boundary layer and lifts the degeneracy of the inviscid baroclinic solutions.For the first time, realistic two-dimensional models of fast-rotating stars are computed with the actual baroclinic flows that predict the differential rotation and the meridional circulation for intermediate-mass and massive stars. These models nicely compare with available data of some nearby fast-rotating early-type stars like Ras Alhague ($\alpha$ Oph), Regulus ($\alpha$ Leo), and Vega ($\alpha$ Lyr). It is shown that baroclinicity drives a differential rotation with a slow pole, a fast equator, a fast core, and a slow envelope. The differential rotation is found to increase with mass, with evolution (here measured by the hydrogen mass fraction in the core), and with metallicity. The core-envelope interface is found to be a place of strong shear where mixing will be efficient.Two-dimensional models offer a new view of fast-rotating stars, especially of their differential rotation, which turns out to be strong at the core-envelope interface. They also offer more accurate models for interpreting the interferometric and spectroscopic data of early-type stars.

The measurement of the Sun's surface motions with a high spatial and temporal resolution is still a challenge. We wish to validate horizontal velocity measurements all over the visible disk of the Sun from Solar Dynamics Observatory/ Helioseismic and Magnetic Imager (SDO/HMI) data. Horizontal velocity fields are measured by following the proper motions of solar granules using a newly developed version of the Coherent Structure Tracking (CST) code. The comparison of the surface flows measured at high spatial resolution (Hinode, 0.1 arcsec) and low resolution (SDO/HMI, 0.5 arcsec) allows us to determine corrections to be applied to the horizontal velocity measured from HMI white light data. We derive horizontal velocity maps with spatial and temporal resolutions of respectively 2.5 Mm and 30 min. From the two components of the horizontal velocity Vx and Vy measured in the sky plane and the simultaneous line of sight component from SDO/HMI dopplergrams v_D, we derive the spherical velocity components (Vr, Vtheta, Vphi). The azimuthal component Vphi gives the solar differential rotation with a high precision (+-0.037km/s) from a temporal sequence of only three hours. By following the proper motions of the solar granules, we can revisit the dynamics of the solar surface at high spatial and temporal resolutions from hours to months and years with the SDO data.

We investigate the properties of small-amplitude inertial waves propagating in a differentially rotating incompressible fluid contained in a spherical shell. For cylindrical and shellular rotation profiles and in the inviscid limit, inertial waves obey a second-order partial differential equation of mixed type. Two kinds of inertial modes therefore exist, depending on whether the hyperbolic domain where characteristics propagate covers the whole shell or not. The occurrence of these two kinds of inertial modes is examined, and we show that the range of frequencies at which inertial waves may propagate is broader than with solid-body rotation. Using high-resolution calculations based on a spectral method, we show that, as with solid-body rotation, singular modes with thin shear layers following short-period attractors still exist with differential rotation. They exist even in the case of a full sphere. In the limit of vanishing viscosities, the width of the shear layers seems to weakly depend on the global background shear, showing a scaling in E^{1/3} with the Ekman number E, as in the solid-body rotation case. There also exist modes with thin detached layers of width scaling with E^{1/2} as Ekman boundary layers. The behavior of inertial waves with a corotation resonance within the shell is also considered. For cylindrical rotation, waves get dramatically absorbed at corotation. In contrast, for shellular rotation, waves may cross a critical layer without visible absorption, and such modes can be unstable for small enough Ekman numbers.

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.

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.

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.

We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet sun using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard \textit{Hinode}. Line profiles of Fe \textsc{i} 557.6 nm were recorded by the Narrow band Filter Imager (NFI) on a $82" \times 82"$ FOV during 75 min with a time step of 28.75 s and 0.08$"$ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130 and 180 km) using the bisector technique allowing the determination of energy flux in the range 3-10 mHz using two complementary methods (Hilbert transform and Fourier power spectra). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. The net energy flux is upward. In the range 3-10 mHz, a full FOV space and time averaged flux of 2700 W m$^{-2}$ (lower layer 80-130 km) and 2000 W m$^{-2}$ (upper layer 130-180 km) is concentrated in less than 1% of the solar surface in the form of narrow (0.3$"$) AE. Their total duration (including rise and decay) is of the order of $10^{3}$ s. Inside each AE, the mean flux is $1.6 10^{5}$ W m$^{-2}$ (lower layer) and $1.2 10^{5}$ W m$^{-2}$ (upper). Each event carries an average energy (flux integrated over space and time) of $2.5 10^{19}$ J (lower layer) to $1.9 10^{19}$ J (upper). More than $10^{6}$ events could exist permanently on the Sun, with a birth and decay rate of 3500 s$^{-1}$. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.

We briefly review the observations of the solar photosphere and pinpoint some open questions related to the magnetohydrodynamics of this layer of the Sun. We then discuss the current modelling efforts, addressing among other problems, that of the origin of supergranulation.

For the first time, the motion of granules (solar plasma on the surface on scales larger than 2.5 Mm) has been followed over the entire visible surface of the Sun, using SDO/HMI white-light data. Horizontal velocity fields are derived from image correlation tracking using a new version of the coherent structure tracking algorithm.The spatial and temporal resolutions of the horizontal velocity map are 2.5 Mm and 30 min respectively . From this reconstruction, using the multi-resolution analysis, one can obtain to the velocity field at different scales with its derivatives such as the horizontal divergence or the vertical component of the vorticity. The intrinsic error on the velocity is ~0.25 km/s for a time sequence of 30 minutes and a mesh size of 2.5 Mm.This is acceptable compared to the granule velocities, which range between 0.3 km/s and 1.8 km/s. A high correlation between velocities computed from Hinode and SDO/HMI has been found (85%). From the data we derive the power spectrum of the supergranulation horizontal velocity field, the solar differential rotation, and the meridional velocity.

We present complete 2D computations of g modes in distorted polytropic models of stars performed with the Two-dimensional Oscillation Program (TOP). We computed low-degree modes (l=1 modes with radial order n=-1...-14, and l=2,3 modes with n=-1...-5 and -16...-20) of a nonrotating model and followed them by slowly increasing the rotation rate up to 70 % of the Keplerian break-up velocity. We use these computations to determine the domain of validity of perturbative methods up to the 3rd order. We study the evolution of the regularities of the spectrum and show quantitative agreement with the traditional approximation for not too large values of the ratio of the rotation rate to the pulsation frequency. We also show the appearance of new types of modes, called "rosette" modes due to their spatial structure. Thanks to the ray theory for gravito-inertial waves that we developed, we can associate these modes with stable periodic rays.

The impact of the combination of rotation and magnetic fields on oscillations of stellar fluids is still not well known theoretically. It mixes Alfv\'en and inertial waves. Neutron stars are a place where both effects may be at work. We wish to decipher the solution of this problem in the context of $r$-modes instability in neutron stars, as it appears when these modes are coupled to gravitational radiation. We consider a rotating spherical shell filled with a viscous fluid but of infinite electrical conductivity and analyze propagation of modal perturbations when a dipolar magnetic field is bathing the fluid layer. We perform an extensive numerical analysis and find that the $m=2$ $r$-mode oscillation is influenced by the magnetic field when the Lehnert number (ratio of Alfv\'en speed to rotation speed) exceeds a value proportional to the one-fourth power of the Ekman number (non-dimensional measure of viscosity). This scaling is interpreted as the coincidence of the width of internal shear layers of inertial modes and the wavelength of the Alfv\'en waves. Applied to the case of rotating magnetic neutron stars, we find that dipolar magnetic fields above $10^{14}$ G are necessary to perturb the $r$-modes instability.

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^{\beta} 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, {\alpha} Aql and {\alpha} 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.

This work aims at determining the impact of slow to moderate rotation on the regular patterns often present in solar-like oscillation spectra. We focus on the well-known asteroseismic diagnostic echelle diagrams, examining how rotation may modify the estimates of the large and small spacings, as well as the identification of modes. We illustrate the work with a real case: the solar-like star $\eta$ Bootis. The modeling takes into account rotation effects on the equilibrium models through an effective gravity and on the oscillation frequencies through both perturbative and non-perturbative calculations. We compare the results of both type of calculations in the context of the regular spacings (like the small spacings and the scaled small spacings) and echelle diagrams. We show that for echelle diagrams the perturbative approach remains valid for rotational velocities up to 40-50 km/s. We show that for the rotational velocities measured in solar-like stars, theoretical oscillation frequencies must be corrected up to the second-order in terms of rotation rate, including near degeneracy effects. For rotational velocities of about 16 km/S and higher, diagnostics on large spacings and on modal identification through echelle diagrams can be significantly altered by the presence of the $m\neq0$ components of the rotationally split modes. We found these effects to be detectable in the observed frequency range. Analysis of the effects of rotation on small spacings and scaled small spacings reveals that these can be of the order of, or even larger than surface effects, typically turbulence, microscopic diffusion, etc. Furthermore, we show that scaled spacings are significantly affected by stellar distortion even for small stellar rotational velocities (from 10-15 km/s) and therefore some care must be taken when using them as indicators for probing deep stellar interiors.

CoRoT and Kepler missions are now providing high-quality asteroseismic data for a large number of stars. Among intermediate-mass and massive stars, fast rotators are common objects. Taking the rotation effects into account is needed to correctly understand, identify, and interpret the observed oscillation frequencies of these stars. A classical approach is to consider the rotation as a perturbation. In this paper, we focus on gravity modes, such as those occurring in gamma Doradus, slowly pulsating B (SPB), or Be stars. We aim to define the suitability of perturbative methods. With the two-dimensional oscillation program (TOP), we performed complete computations of gravity modes -including the Coriolis force, the centrifugal distortion, and compressible effects- in 2-D distorted polytropic models of stars. We started with the modes l=1, n=1-14, and l=2-3, n=1-5,16-20 of a nonrotating star, and followed these modes by increasing the rotation rate up to 70% of the break-up rotation rate. We then derived perturbative coefficients and determined the domains of validity of the perturbative methods. Second-order perturbative methods are suited to computing low-order, low-degree mode frequencies up to rotation speeds ~100 km/s for typical gamma Dor stars or ~150 km/s for B stars. The domains of validity can be extended by a few tens of km/s thanks to the third-order terms. For higher order modes, the domains of validity are noticeably reduced. Moreover, perturbative methods are inefficient for modes with frequencies lower than the Coriolis frequency 2Omega. We interpret this failure as a consequence of a modification in the shape of the resonant cavity that is not taken into account in the perturbative approach.

The aim of this work is to give new observational constraints on solar surface flows by determining the horizontal scale dependence of the velocity and intensity fields, as represented by their power spectra, and to offer some theoretical guidelines to interpret these spectra. We use long time series of images taken by SOT/Hinode and reconstruct both horizontal (by granule tracking) and vertical (by Doppler effect) velocity fields in a field of view 75x75Mm^2. At small sub-granulation scales, the kinetic energy spectral density associated with vertical motions exhibits a k^{-10/3}-like spectrum, while the intensity fluctuation spectrum follows a k^{-3} or k^{-17/3}-like spectrum at the two continuum levels investigated (525 and 450 nm respectively). We discuss the physical origin of these scalings and argue that they provide a direct observational signature of buoyancy-driven turbulent dynamics in a strongly thermally diffusive regime. In the mesogranulation range and up to a scale of 25Mm, we find that the vertical velocity field amplitude decreases like L^{-3/2} with the horizontal scale L. This behaviour corresponds to a k^2 spectral power law. Still in the mesoscale range, [2.5, 10]Mm, we find that intensity fluctuations in the blue continuum follow a power law in k^2. We show that granule tracking cannot sample scales below 2.5Mm. We locate the supergranulation energy peak around 30Mm and show that the emergence of a pore erases this spectral peak. Thanks to a scale height estimate, we find that supergranular flows are shallow. (abridged abstract).

Using the Two-dimensional Oscillation Program (TOP), we have explored the effects of rapid rotation on gravity modes in polytropic stars. Coriolis force, centrifugal distortion as well as compressible effects have been taken into account. Thanks to our complete calculation, we have first studied the validity domain of perturbative methods and started to explore properties of these modes. We focus on l=1 in this analysis.

We investigate the properties of forced inertial modes of a rotating fluid inside a spherical shell. Our forcing is tidal like, but its main property is that it is on the large scales. Our solutions first confirm some analytical results obtained on a two-dimensional model by Ogilvie (2005). We also note that as the frequency of the forcing varies, the dissipation varies drastically if the Ekman number E is low (as is usually the case). We then investigate the three-dimensional case and compare the results to the foregoing model. These solutions show, like their 2D counterpart, a spiky dissipation curve when the frequency of the forcing is varied; they also display small frequency intervals where the viscous dissipation is independent of viscosity. However, we show that the response of the fluid in these frequency intervals is crucially dominated by the shear layer that is emitted at the critical latitude on the inner sphere. The asymptotic regime is reached when an attractor has been excited by this shear layer. This property is not shared by the two-dimensional model. Finally, resonances of the three-dimensional model correspond to some selected least-damped eigenmodes. Unlike their two-dimensional counter parts these modes are not associated with simple attractors; instead, they show up in frequency intervals with a weakly contracting web of characteristics. Besides, we show that the inner core is negligible when its relative radius is less than the critical value 0.4E^{1/5}. For these spherical shells, the full sphere solutions give a good approximation of the flows (abridged abstract).

Aims: The interactions of velocity scales on the Sun's surface, from granulation to supergranulation are still not understood, nor are their interaction with magnetic fields. We thus aim at giving a better description of dynamics in the mesoscale range which lies between the two scales mentioned above. Method: We analyse a 48h high-resolution time sequence of the quiet Sun photosphere at the disk center obtained with the Solar Optical Telescope onboard Hinode. The observations, which have a field of view of 100 \arcsec$\times$ 100 \arcsec, typically contain four supergranules. We monitor in detail the motion and evolution of granules as well as those of the radial magnetic field. Results: This analysis allows us to better characterize Trees of Fragmenting Granules issued from repeated fragmentation of granules, especially their lifetime statistics. Using floating corks advected by measured velocity fields, we show their crucial role in the advection of the magnetic field and in the build up of the network. Finally, thanks to the long duration of the time series, we estimate that the turbulent diffusion coefficient induced by horizontal motion is approximately $430 \mathrm{km}^2 \mathrm{s}^{-1}$. Conclusions: These results demonstrate that the long living families contribute to the formation of the magnetic network and suggest that supergranulation could be an emergent length scale building up as small magnetic elements are advected and concentrated by TFG flows. Our estimate for the magnetic diffusion associated with this horizontal motion might provide a useful input for mean-field dynamo models.

Context: Rapid rotation modifies the structure of the frequency spectrum of pulsating stars, thus making mode identification difficult. Aims: We look for new forms of organisation for the frequency spectrum that can provide a basis for mode identification at high rotation rates. Methods: Acoustic modes in uniformly rotating polytropic models of stars are computed using a numerical code that fully takes the effects of rotation (centrifugal distortion and Coriolis acceleration) into account. All low-degree modes, l=0 to 3, with radial orders n=1-10 and 21-25 for N=3 polytropic models and n=1-10 for N=1.5 polytropic models are followed from a zero rotation rate up to 59 % of the break-up velocity. Results: We find an empirical formula that gives a good description of the high-frequency range of the computed acoustic spectrum for high rotation rates. Differences between this formula and complete eigenmode calculations are shown to be substantially smaller than those obtained with a third order perturbative method valid at low rotation rates.

Context: Supergranulation is a pattern of the velocity field at the surface of the Sun, which has been known about for more than fifty years, however, no satisfactory explanation of its origin has been proposed. Aims: New observational constraints are therefore needed to guide theoretical approaches which hesitate between scenarios that either invoke a large-scale instability of the surface turbulent convection or a direct forcing by buoyancy. Method: Using the 14-Mpixel CALAS camera at the Pic-du-Midi observatory, we obtained a 7.5h-long sequence of high resolution images with unprecedented field size. Tracking granules, we have determined the velocity field at the Sun's surface in great detail from a scale of 2.5Mm up to 250Mm. Results: The kinetic energy density spectrum shows that supergranulation peaks at 36Mm and spans on scales ranging between 20Mm and 75Mm. The decrease of supergranular flows in the small scales is close to a $k^{-2}$-power law, steeper than the equipartition Kolmogorov one. The probability distribution function of the divergence field shows the signature of intermittency of the supergranulation and thus its turbulent nature.

We study the influence of large-scale photospheric motions on the estabilization of an eruptive filament, observed on October 6, 7, and 8, 2004, as part of an international observing campaign (JOP 178). Large-scale horizontal flows were invetigated from a series of MDI full-disc Dopplergrams and magnetograms. From the Dopplergrams, we tracked supergranular flow patterns using the local correlation tracking (LCT) technique. We used both LCT and manual tracking of isolated magnetic elements to obtain horizontal velocities from magnetograms. We find that the measured flow fields obtained by the different methods are well-correlated on large scales. The topology of the flow field changed significantly during the filament eruptive phase, suggesting a possible coupling between the surface flow field and the coronal magnetic field. We measured an increase in the shear below the point where the eruption starts and a decrease in shear after the eruption. We find a pattern in the large-scale horizontal flows at the solar surface that interact with differential rotation. We conclude that there is probably a link between changes in surface flow and the disappearance of the eruptive filament.

Determination of horizontal velocity fields on the solar surface is crucial for understanding the dynamics of structures like mesogranulation or supergranulation or simply the distribution of magnetic fields. We pursue here the development of a method called CST for coherent structure tracking, which determines the horizontal motion of granules in the field of view. We first devise a generalization of Strous method for the segmentation of images and show that when segmentation follows the shape of granules more closely, granule tracking is less effective for large granules because of increased sensitivity to granule fragmentation. We then introduce the multi-resolution analysis on the velocity field, based on Daubechies wavelets, which provides a view of this field on different scales. An algorithm for computing the field derivatives, like the horizontal divergence and the vertical vorticity, is also devised. The effects from the lack of data or from terrestrial atmospheric distortion of the images are also briefly discussed.

The determination of horizontal velocity fields at the solar surface is crucial to understanding the dynamics and magnetism of the convection zone of the sun. These measurements can be done by tracking granules. Tracking granules from ground-based observations, however, suffers from the Earth's atmospheric turbulence, which induces image distortion. The focus of this paper is to evaluate the influence of this noise on the maps of velocity fields. We use the coherent structure tracking algorithm developed recently and apply it to two independent series of images that contain the same solar signal. We first show that a k-\omega filtering of the times series of images is highly recommended as a pre-processing to decrease the noise, while, in contrast, using destretching should be avoided. We also demonstrate that the lifetime of granules has a strong influence on the error bars of velocities and that a threshold on the lifetime should be imposed to minimize errors. Finally, although solar flow patterns are easily recognizable and image quality is very good, it turns out that a time sampling of two images every 21 s is not frequent enough, since image distortion still pollutes velocity fields at a 30% level on the 2500 km scale, i.e. the scale on which granules start to behave like passive scalars. The coherent structure tracking algorithm is a useful tool for noise control on the measurement of surface horizontal solar velocity fields when at least two independent series are available.

Recent results from interferometry and asteroseismology require models of rapidly rotating stars that are more and more precise. We describe the basic structure and the hydrodynamics of a fully radiative star as a preliminary step towards more realistic models of rotating stars. We consider a solar mass of perfect gas enclosed in a spherical container. The gas is self-gravitating and rotating, and is the seat of nuclear heating, and heat diffusion is due to radiative diffusion with Kramers type opacities. Equations are solved numerically with spectral methods in two dimensions with a radial Gauss-Lobatto grid and spherical harmonics. We computed the centrifugally flattened structure of such a star: the von Zeipel model, which says that the energy flux is proportional to the local effective gravity is tested. We show that it overestimates the ratio of the polar to the equatorial energy flux by almost a factor 2. We also determine the Brunt-Vaisala frequency distribution and show that outer equatorial regions in a radiative zone are convectively unstable when the rotation is fast enough. We compute the differential rotation and meridional circulation stemming from the baroclinicity of the star and show that, in such radiative zones, equatorial regions rotate faster than polar ones. The surface differential rotation is also shown to reach a universal profile when rotation is slow enough (less than 36% of the breakup one), as long as viscosity and Prandlt number remain small.

We review the quest of modeling rapidly rotating stars during the past 40 years and detail the challenges to be taken up by models facing new data from interferometry, seismology, spectroscopy... We then present the progress of the ESTER project aimed at giving a physically self-consistent model for the structure and evolution of rapidly rotating stars.

We present the results of a very large aspect ratio (42.6) numerical simulation of fully compressible turbulent convection in a polytropic atmosphere, and focus on the properties of large-scale flows. Mesoscale patterns dominate the turbulent energy spectrum. We show that these structures, which had already been observed in Boussinesq simulations by Cattaneo et al. (2001), have a genuine convective origin and do not result directly from collective interactions of the smaller scales of the flow, even though their growth is strongly affected by nonlinear transfers. If this result is relevant to the solar photosphere, it suggests that the dominant convective mode below the Sun's surface may be at mesoscales.

Context: With the launch of space missions devoted to asteroseismology (like COROT), the scientific community will soon have accurate measurements of pulsation frequencies in many rapidly rotating stars. Aims: The present work focuses on the effects of rotation on pulsations of rapidly rotating stars when both the Coriolis and centrifugal accelerations require a non-perturbative treatment. Method: We develop a 2-dimensional spectral numerical approach which allows us to compute acoustic modes in centrifugally distorted polytropes including the full influence of the Coriolis force. This method is validated through comparisons with previous studies, and the results are shown to be highly accurate. Results: In the frequency range considered and with COROT's accuracy, we establish a domain of validity for perturbative methods, thus showing the need for complete calculations beyond v.sin i = 50 km/s for a R = 2.3 R_\odot, M = 1.9 M_\odot polytropic star. Furthermore, it is shown that the main differences between complete and perturbative calculations come essentially from the centrifugal distortion.

In this lecture I try to explain the basic dynamical processes at work in a radiative zone of a rotating star. In particular, the notion of baroclinicity is thoroughly discussed. Attention is specially directed to the case of circulations and the key role of angular momentum conservation is stressed. The specific part played by viscosity is also explained. The old approach of Eddington and Sweet is reviewed and criticized in the light of the seminal papers of Busse 1981 and Zahn 1992. Other examples taken in the recent literature are also presented; finally, I summarize the important points.

An efficient way of solving 2D stability problems in fluid mechanics is to use, after discretization of the equations that cast the problem in the form of a generalized eigenvalue problem, the incomplete Arnoldi-Chebyshev method. This method preserves the banded structure sparsity of matrices of the algebraic eigenvalue problem and thus decreases memory use and CPU-time consumption. The errors that affect computed eigenvalues and eigenvectors are due to the truncation in the discretization and to finite precision in the computation of the discretized problem. In this paper we analyze those two errors and the interplay between them. We use as a test case the two-dimensional eigenvalue problem yielded by the computation of inertial modes in a spherical shell. This problem contains many difficulties that make it a very good test case. It turns out that that single modes (especially most-damped modes i.e. with high spatial frequency) can be very sensitive to round-off errors, even when apparently good spectral convergence is achieved. The influence of round-off errors is analyzed using the spectral portrait technique and by comparison of double precision and extended precision computations. Through the analysis we give practical recipes to control the truncation and round-off errors on eigenvalues and eigenvectors.

A new non-perturbative method to compute accurate oscillation modes in rapidly rotating stars is presented. In this paper, the effect of the centrifugal force is fully taken into account while the Coriolis force is neglected. This assumption is valid when the time scale of the oscillation is much shorter than the inverse of the rotation rate and is expected to be suitable for high radial order p-modes of $\delta$ Scuti stars. Axisymmetric p-modes have been computed in uniformly rotating polytropic models of stars. In the frequency and rotation range considered, we found that as rotation increases (i) the asymptotic structure of the non-rotating frequency spectrum is first destroyed then replaced by a new form of organization (ii) the mode amplitude tends to concentrate near the equator (iii) differences with perturbative methods become significant as soon as the rotation rate exceeds about fifteen percent of the Keplerian limit. The implications for the seismology of rapidly rotating stars are then discussed.

Context: The observations of rapidly rotating stars are increasingly detailed and precise thanks to interferometry and asteroseismology; two-dimensional models taking into account the hydrodynamics of these stars are very much needed. Aims: A model for studying the dynamics of baroclinic stellar envelope is presented. Methods: This models treats the stellar fluid at the Boussinesq approximation and assumes that it is contained in a rigid spherical domain. The temperature field along with the rotation of the system generate the baroclinic flow. Results: We manage to give an analytical solution to the asymptotic problem at small Ekman and Prandtl numbers. We show that, provided the Brunt-Vaisala frequency profile is smooth enough, differential rotation of a stably stratified envelope takes the form a fast rotating pole and a slow equator while it is the opposite in a convective envelope. We also show that at low Prandtl numbers and without $\mu$-barriers, the jump in viscosity at the core-envelope boundary generates a shear layer staying along the tangential cylinder of the core. Its role in mixing processes is discussed. Conclusions: Such a model provides an interesting tool for investigating the fluid dynamics of rotating stars in particular for the study of the various instabilities affecting baroclinic flows or, even more, of a dynamo effect.

The solar magnetism is no more considered as a purely superficial phenomenon. The SoHO community has shown that the length of the solar cycle depends on the transition region between radiation and convection. Nevertheless, the internal solar (stellar) magnetism stays poorly known. Starting in 2008, the American instrument HMI/SDO and the European microsatellite PICARD will enrich our view of the Sun-Earth relationship. Thus obtaining a complete MHD solar picture is a clear objective for the next decades and it requires complementary observations of the dynamics of the radiative zone. For that ambitious goal, space prototypes are being developed to improve gravity mode detection. The Sun is unique to progress on the topology of deep internal magnetic fields and to understand the complex mechanisms which provoke photospheric and coronal magnetic changes and possible longer cycles important for human life. We propose the following roadmap in Europe to contribute to this "impressive" revolution in Astronomy and in our Sun-Earth relationship: SoHO (1995-2007), PICARD (2008-2010), DynaMICS (2009-2017) in parallel to SDO (2008-2017) then a world-class mission located at the L1 orbit or above the solar poles.

We carry out numerical and mathematical investigations of shear Alfven waves inside of a spherical shell filled with an incompressible conducting fluid, and bathed in a strong dipolar magnetic field. We focus on axisymmetric toroidal and non-axisymmetric modes, in continuation of a previous work by Rincon & Rieutord (2003, A&A, 398, 663). Analytical expressions are obtained for toroidal eigenmodes and their corresponding frequencies at low diffusivities. These oscillations behave like magnetic shear layers, in which the magnetic poles play a key role, and hence become singular when diffusivities vanish. It is also demonstrated that non-axisymmetric modes are split into two categories, namely poloidal or toroidal types, following similar asymptotic behaviours as their axisymmetric counterparts when the diffusivities become arbitrarily small.

We compute statistical equilibrium states of rotating self-gravitating systems enclosed within a box by maximizing the Fermi-Dirac entropy at fixed mass, energy and angular momentum. We increase the rotation up to the Keplerian limit and describe the flattening of the configuration until mass shedding occurs. At the maximum rotation, the system develops a cusp at the equator. We draw the equilibrium phase diagram of the rotating self-gravitating Fermi gas and discuss the structure of the caloric curve as a function of degeneracy parameter and angular velocity. We argue that systems described by the Fermi-Dirac distribution in phase space do not bifurcate to non-axisymmetric structures, in continuity with the case of polytropes with index n>0.808 (the Fermi gas at T=0 corresponds to n=3/2). This contrasts with the study of Votyakov et al. (2002) who consider a Fermi-Dirac distribution in configuration space and find ``double star'' structures (their model at T=0 corresponds to n=0). We also discuss the influence of rotation on the onset of the gravothermal catastrophe for globular clusters. On general grounds, we complete previous investigations concerning the nature of phase transitions in self-gravitating systems. We emphasize the inequivalence of statistical ensembles regarding the formation of binaries (or low-mass condensates) in the microcanonical ensemble and Dirac peaks (or massive condensates) in the canonical ensemble. We also describe an hysteretic cycle between the gaseous phase and the condensed phase that are connected by a ``collapse'' or an ``explosion''. This notion of hysteresis in self-gravitating systems is new.

In this review, we describe the physical processes driving the dynamical evolution of binary stars, namely the circularization of the orbit and the synchronization of their spin and orbital rotation. We also discuss the possible role of the elliptic instability which turns out to be an unavoidable ingredient of the evolution of binary stars.

We present here numerical simulations of surface solar convection which cover a box of 30$\times30\times$3.2 Mm$^3$ with a resolution of 315$\times315\times$82, which is used to investigate the dynamics of scales larger than granulation. No structure resembling supergranulation is present; possibly higher Reynolds numbers (i.e. higher numerical resolution), or magnetic fields, or greater depth are necessary. The results also show interesting aspects of granular dynamics which are briefly presented, like extensive p-mode ridges in the k-$\omega$ diagram and a ringlike distribution of horizontal vorticity around granules. At large scales, the horizontal velocity is much larger than the vertical velocity and the vertical motion is dominated by p-mode oscillations.

Two approximations, namely the subseismic approximation and the anelastic approximation, are presently used to filter out the acoustic modes when computing low frequency modes of a star (gravity modes or inertial modes). In a precedent paper (Dintrans & Rieutord 2001), we observed that the anelastic approximation gave eigenfrequencies much closer to the exact ones than the subseismic approximation. Here, we try to clarify this behaviour and show that it is due to the different physical approach taken by each approximation: On the one hand, the subseismic approximation considers the low frequency part of the spectrum of (say) gravity modes and turns out to be valid only in the central region of a star; on the other hand, the anelastic approximation considers the Brunt-Vaisala frequency as asymptotically small and makes no assumption on the order of the modes. Both approximations fail to describe the modes in the surface layers but eigenmodes issued from the anelastic approximation are closer to those including acoustic effects than their subseismic equivalent. We conclude that, as far as stellar eigenvalue problems are concerned, the anelastic approximation is better suited for simplifying the eigenvalue problem when low-frequency modes of a star are considered, while the subseismic approximation is a useful concept when analytic solutions of high order low-frequency modes are needed in the central region of a star.

We present results concerning the linear (radial and non-radial) oscillations of non-rotating superfluid neutron stars in Newtonian physics. We use a simple two-fluid model to describe the superfluid neutron star, where one fluid consists of the superfluid neutrons,while the second fluid contains all the comoving constituents (protons, electrons). The two fluids are assumed to be "free" in the sense of absence of vortex-mediated forces like mutual friction or pinning, but they can be coupled by the equation of state, in particular by entrainment.We calculate numerically the eigen-frequencies and -modes of adiabatic oscillations, neglecting beta-reactions that would lead to dissipation.We find a doubling of all acoustic-type modes (f-modes, p-modes), and confirm the absence of g-modes in these superfluid models. We show analytically and numerically that only in the case of non-stratified background models (i.e. with no composition gradient) can these doublets of acoustic modes be separated into two distinct families, which are characterised by either co- or counter-moving fluids respectively, and which are sometimes referred to as "ordinary" and "superfluid" modes. In the general, stratified case, however, this separation is not possible, and these acoustic modes can not be classified as being either purely "ordinary" or "superfluid". We show how the properties of the two-fluid modes change as functions of the coupling by entrainment. We find avoided mode-crossings for the stratified models, while the crossings are not avoided in the non-stratified, separable case. The oscillations of normal-fluid neutron stars are recovered as a special case simply by locking the two fluids together.In this effective one-fluid case we find the usual singlet f- and p-modes, and we also find the expected g-modes.

We present numerical simulations of gravito-inertial waves propagating in radiative zones of rapidly rotating stars. A first model, using the Boussinesq approximation, allows us to study the oscillations of a quasi-incompressible stratified fluid embedded in a rapidly rotating sphere or spherical shell. In a second step, we investigate the case of a Gamma-Doradus star using the anelastic approximation. Some fascinating features of rapidly rotating fluids, such as wave attractors, appear in both cases.

After showing that the anelastic approximation is better than the subseismic one to filter out acoustic waves when studying low-frequency stellar oscillations, we compute gravito-inertial modes of a typical Gamma-Doradus star using this approximation. We show that eigenmodes can be regular or singular, according to the possible focusing towards attractors of the underlying characteristics. Consequences on the oscillations spectrum are then discussed.

Oscillation modes of rapidly rotating stars have not yet been calculated with precision, rotational effects being generally approximated by perturbation methods. We are developing a numerical method able to account for the deformation of the star by the centrifugal force and, as a first step, we determined the acoustic modes of a uniform density spheroid and studied how the spheroid flatness affects these modes.

Using a numerical simulation of compressible convection with radiative transfer mimicking the solar photosphere, we compare the velocity field derived from granule motions to the actual velocity field of the plasma. We thus test the idea that granules may be used to trace large-scale velocity fields at the sun's surface. Our results show that this is indeed the case provided the scale separation is sufficient. We thus estimate that neither velocity fields at scales less than 2500 km nor time evolution at scales shorter than 0.5 hr can be faithfully described by granules. At larger scales the granular motions correlate linearly with the underlying fluid motions with a slope of ~< 2 reaching correlation coefficients up to ~0.9.

A comparative study of the properties of the anelastic and subseismic approximations is presented. The anelastic approximation is commonly used in astrophysics in compressible convection studies whereas the subseismic approximation comes from geophysics where it is used to describe long-period seismic oscillations propagating in the Earth's outer fluid core. Both approximations aim at filtering out the acoutic waves while retaining the density variations of the equilibrium configuration. However, they disagree on the form of the equation of mass conservation. We show here that the anelastic approximation is in fact the only consistent approximation as far as stellar low-frequency oscillations are concerned. We also show that this approximation implies Cowling's approximation which neglects perturbations of the gravity field. Examples are given to assert the efficiency of the anelastic approximation.

Recently, eigenmodes of rotating fluids, namely inertial modes, have received much attention in relation to their destabilization when coupled to gravitational radiation within neutron stars. However, these modes have been known for a long time by fluid dynamicists. We give a short account of their history and review our present understanding of their properties. Considering the case of a spherical container, we then give the exact solution of the boundary (Ekman) layer flow associated with inertial r-modes and show that previous estimations all underestimated the dissipation by these layers. We also show that the presence of an inner core has little influence on this dissipation. As a conclusion, we compute the window of instability in the Temperature/rotation plane for a crusted neutron star when it is modeled by an incompressible fluid.

In the limit of low viscosity, we show that the amplitude of the modes of oscillation of a rotating fluid, namely inertial modes, concentrate along an attractor formed by a periodic orbit of characteristics of the underlying hyperbolic Poincar\'e equation. The dynamics of characteristics is used to elaborate a scenario for the asymptotic behaviour of the eigenmodes and eigenspectrum in the physically relevant r\'egime of very low viscosities which are out of reach numerically. This problem offers a canonical ill-posed Cauchy problem which has applications in other fields.

We investigate the asymptotic properties of inertial modes confined in a spherical shell when viscosity tends to zero. We first consider the mapping made by the characteristics of the hyperbolic equation (Poincar\'e's equation) satisfied by inviscid solutions. Characteristics are straight lines in a meridional section of the shell, and the mapping shows that, generically, these lines converge towards a periodic orbit which acts like an attractor. We then examine the relation between this characteristic path and eigensolutions of the inviscid problem and show that in a purely two-dimensional problem, convergence towards an attractor means that the associated velocity field is not square-integrable. We give arguments which generalize this result to three dimensions. We then consider the viscous problem and show how viscosity transforms singularities into internal shear layers which in general betray an attractor expected at the eigenfrequency of the mode. We find that there are nested layers, the thinnest and most internal layer scaling with $E^{1/3}$-scale, $E$ being the Ekman number. Using an inertial wave packet traveling around an attractor, we give a lower bound on the thickness of shear layers and show how eigenfrequencies can be computed in principle. Finally, we show that as viscosity decreases, eigenfrequencies tend towards a set of values which is not dense in $[0,2\Omega]$, contrary to the case of the full sphere ($\Omega$ is the angular velocity of the system). Hence, our geometrical approach opens the possibility of describing the eigenmodes and eigenvalues for astrophysical/geophysical Ekman numbers ($10^{-10}-10^{-20}$), which are out of reach numerically, and this for a wide class of containers.

We present arguments which show that in all likelihood mesogranulation is not a true scale of solar convection but the combination of the effects of both highly energetic granules, which give birth to strong positive divergences (SPDs) among which we find exploders, and averaging effects of data processing. The important role played by SPDs in horizontal velocity fields appears in the spectra of these fields where the scale $\sim$4 Mm is most energetic; we illustrate the effect of averaging with a one-dimensional toy model which shows how two independent non-moving (but evolving) structures can be transformed into a single moving structure when time and space resolution are degraded. The role of SPDs in the formation of the photospheric network is shown by computing the advection of floating corks by the granular flow. The coincidence of the network bright points distribution and that of the corks is remarkable. We conclude with the possibility that supergranulation is not a proper scale of convection but the result of a large-scale instability of the granular flow, which manifests itself through a correlation of the flows generated by SPDs.