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.

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.

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.