Magnetic confinement of the solar tachocline
  We study the physics of the solar tachocline and related MHD instabilities. We have performed 3-D MHD simulations of the solar radiative interior to check whether a fossil magnetic field is able to prevent the spread of the tachocline. Starting with a purely poloidal magnetic field and a latitudinal shear meant to be imposed by the convection zone at the top of the radiation zone, we have investigated the interactions between magnetic fields, rotation and shear, using the spectral code ASH on massive parallel supercomputers. In all cases we have explored, the fossil field diffuses outward and ends up connecting with the convection zone, whose differential rotation is then imprinted at latitudes above 40 deg throughout the radiative interior, according to Ferraro's law of isorotation. Rotation remains uniform in the lower latitude region which is contained within closed field lines. We find that the meridional flow cannot stop the inward progression of the differential rotation. Further, we observe the development of non-axisymmetric magnetohydrodynamic instabilities, first due to the initial poloidal configuration of the fossil field, and later to the toroidal field produced by shearing the poloidal field through the differential rotation. We do not find dynamo action as such in the radiative interior, since the mean poloidal field is not regenerated. But the instability persists during the whole evolution, while slowly decaying with the mean poloidal field. According to our numerical simulations, a fossil magnetic field cannot prevent the radiative spread of the tachocline, and thus it is unable to enforce uniform rotation in the radiation zone. Neither can the observed thinness of that layer be invoked as a proof for such an internal fossil magnetic field. 
A. S. Brun & J.-P. Zahn