Abstract

Context. Some contracting or expanding stars are thought to host a large-scale magnetic field in their radiative interior. By interacting with the contraction-induced flows, such fields may significantly alter the rotational history of the star. They thus constitute a promising way to address the problem of angular momentum transport during the rapid phases of stellar evolution. Aims. In this work, we aim to study the interplay between flows and magnetic fields in a contracting radiative zone. Methods. We performed axisymmetric Boussinesq and anelastic numerical simulations in which a portion of the radiative zone was modelled by a rotating spherical layer, stably stratified and embedded in a large-scale (either dipolar or quadrupolar) magnetic field. This layer is subject to a mass-conserving radial velocity field mimicking contraction. The quasi-steady flows were studied in strongly or weakly stably stratified regimes relevant for pre-main sequence stars and for the cores of subgiant and red giant stars. The parametric study consists in varying the amplitude of the contraction velocity and of the initial magnetic field. The other parameters were fixed with the guidance of a previous study. Results. After an unsteady phase during which the toroidal field grew linearly and then back-reacted on the flow, a quasi-steady configuration was reached, characterised by the presence of two magnetically decoupled regions. In one of them, magnetic tension imposes solid-body rotation. In the other, called the dead zone, the main force balance in the angular momentum equation does not involve the Lorentz force and a differential rotation exists. In the strongly stably stratified regime, when the initial magnetic field is quadrupolar, a magnetorotational instability is found to develop in the dead zones. The large-scale structure is eventually destroyed and the differential rotation is able to build up in the whole radiative zone. In the weakly stably stratified regime, the instability is not observed in our simulations, but we argue that it may be present in stars. Conclusions. We propose a scenario that may account for the post-main sequence evolution of solar-like stars, in which quasi-solid rotation can be maintained by a large-scale magnetic field during a contraction timescale. Then, an axisymmetric instability would destroy this large-scale structure and this enables the differential rotation to set in. Such a contraction-driven instability could also be at the origin of the observed dichotomy between strongly and weakly magnetic intermediate-mass stars.

Highlights

  • Rotation is ubiquitous at every stage of stellar evolution

  • In this work we investigated the dynamics of a contracting radiative spherical layer embedded in a large-scale magnetic field

  • The contr#a»ction is modelled through an imposed radial velocity field V f and the gas dynamics is modelled using either the Boussinesq or the anelastic approximations

Read more

Summary

Introduction

Rotation is ubiquitous at every stage of stellar evolution. Yet, it is still often considered as a second order effect in stellar evolution models. In a previous work (Gouhier et al 2021), we investigated the differential rotation and meridional circulation produced in a modelled contracting stellar radiative zone This axisymmetric study included the effects of stable stratification and was conducted both under the Boussinesq and the anelastic approximations but the effects of a magnetic field were ignored. A distinct population of MS intermediate-mass stars, including the A-type star Vega and the Am-type stars Sirius, β Ursae Majoris and θ Leonis, exhibits much weaker (∼ 1 G) multi-polar magnetic fields (Lignières et al 2009; Petit et al 2010, 2011; Blazère et al 2016) This magnetic dichotomy could be explained if, during the PMS, contraction forces a differential rotation that destroys pre-existing large scale weak magnetic fields through magnetohydrodynamic (MHD) instabilities (Aurière et al 2007; Lignières et al 2014; Jouve et al 2015, 2020).

Mathematical formulation
Timescales of physical processes
Initial and boundary conditions
Numerical method
The stellar context
The numerical study
Numerical results
Unsteady evolution
Steady state in the viscous regime with a dipolar field
Region outside the dead zone
Dead zone
Effect of the density stratification
Description of the instability
Quadrupolar field in the viscous regime
Non-linear evolution
Post-instability description
Meridional circulation and differential rotation confined to the dead zone
Meridional circulation and differential rotation not confined to a dead zone
Summary and conclusions
D39 Dipole
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call