Axisymmetric investigation of differential rotation in contracting stellar radiative zones

Institut de Recherche en Astrophysique et Planétologie (IRAP), Université de Toulouse, 14 Avenue Edouard Belin, 31400 Toulouse, France
e-mail: bgouhier@irap.omp.eu; flignieres@irap.omp.eu; ljouve@irap.omp.eu

Received: 24 August 2020
Accepted: 11 November 2020

Abstract

Context. Stars experience rapid contraction or expansion at different phases of their evolution. Modelling the transport of angular momentum and the transport of chemical elements occurring during these phases remains an unsolved problem.

Aims. We study a stellar radiative zone undergoing radial contraction and investigate the induced differential rotation and meridional circulation.

Methods. We consider a rotating spherical layer crossed by an imposed radial velocity field that mimics the contraction, and numerically solve the axisymmetric hydrodynamical equations in both the Boussinesq and anelastic approximations. An extensive parametric study is conducted to cover regimes of contraction, rotation, stable stratification, and density stratification that are relevant for stars.

Results. The differential rotation and the meridional circulation result from a competition between the contraction-driven inward transport of angular momentum and an outward transport dominated by either viscosity or an Eddington–Sweet-type circulation, depending on the value of the P_r(N₀/Ω₀)² parameter, where P_r is the Prandtl number, N₀ the Brunt–Väisäilä frequency, and Ω₀ the rotation rate. Taking the density stratification into account is important to study more realistic radial contraction fields, and also because the resulting flow is less affected by unwanted effects of the boundary conditions. In these different regimes and for a weak differential rotation we derive scaling laws that relate the amplitude of the differential rotation to the contraction timescale.