A Consistent One-Dimensional Model for the Turbulent Tachocline

K. Petrovay

DOI: https://doi.org/10.1023/A%3A1024880723754

2003-07-31

Abstract:The first consistent model for the turbulent tachocline is presented, with the turbulent diffusivity computed within the model instead of being specified arbitrarily. For the origin of the 3D turbulence a new mechanism is proposed. Owing to the strongly stable stratification, the mean radial shear is stable, while the horizontal shear is expected to drive predominantly horizontal, quasi-2D motions in thin slabs. Here I suggest that a major source of 3D overturning turbulent motions in the tachocline is the secondary shear instability due to the strong, random vertical shear arising between the uncorrelated horizontal flows in neighbouring slabs. A formula for the vertical diffusivity due to this turbulence, equation (9), is derived and applied in a simplified 1D model of the tachocline. It is found that Maxwell stresses due to an oscillatory poloidal magnetic field of a few hundred gauss are able to confine the tachocline to a thickness below 5 Mm. The integral scale of the 3D overturning turbulence is the buoyancy scale, on the order of 10 km and its velocity amplitude is a few m/s, yielding a vertical turbulent diffusivity on the order of 10^8 cm^2/s.

Astrophysics

What problem does this paper attempt to address?

This paper attempts to solve some key problems in the solar tachocline, especially how to calculate the turbulent diffusivity self - consistently in the theoretical model rather than specifying it arbitrarily. Specifically:

1. **Self - consistent calculation of turbulent diffusivity**: Previous studies usually ignored turbulence or assumed a fixed turbulent diffusivity. However, this paper proposes a new mechanism to explain the origin of 3D turbulence and obtains a formula (formula (9)) to calculate the turbulent diffusivity in the vertical direction through theoretical derivation. This enables the model to describe the physical processes inside the sun more accurately.

2. **New mechanism of 3D turbulence**: Due to strong stable stratification, the radial shear is stable, while the horizontal shear drives a thin layer with mainly horizontal, quasi - two - dimensional motion. The author proposes that the strong vertical shear caused by the random horizontal flow between adjacent thin layers is the main source of 3D overturning turbulence. This secondary shear instability can generate enough turbulence to affect the structure and evolution of the solar tachocline.

3. **Effect of Maxwell stress on the thickness of the tachocline**: By introducing the Maxwell stress generated by the oscillating poloidal magnetic field, the study found that a poloidal magnetic field of several hundred gauss is sufficient to confine the tachocline within 5 megameters. This helps to understand the relationship between the magnetic field inside the sun and the thickness of the tachocline.

4. **Turbulence characteristics**: According to theoretical analysis, the integral scale of 3D overturning turbulence is about 10 kilometers, the velocity amplitude is several meters per second, resulting in a vertical turbulent diffusivity of the order of \(10^8 \, \text{cm}^2/\text{s}\). These results are very important for understanding the dynamic processes inside the sun.

In summary, this paper aims to establish a self - consistent one - dimensional model to describe the solar turbulent tachocline and proposes a new mechanism to explain the origin of 3D turbulence. By introducing Maxwell stress and the oscillating poloidal magnetic field, the model can explain why the tachocline can be maintained within a relatively small thickness range. This work provides an important theoretical basis for understanding the complex physical processes inside the sun.

### Formula display

- **Formula for magnetic diffusion depth**:
\[
H_{\text{skin}}=\left(\frac{2\eta}{\omega_{\text{cyc}}}\right)^{1/2}
\]
where \(\eta\) is the magnetic diffusivity, \(\omega_{\text{cyc}} = 2\pi/P\), and \(P = 22\) years is the solar cycle.

- **Square of Brunt - Väisälä frequency**:
\[
N_{\text{BV}}^2 [s^{-2}]=-\frac{\Delta \nabla g}{H_P}\sim 1.5\times 10^{-7}z [\text{Mm}]
\]
where \(\Delta \nabla\) is the thermal stratification parameter, \(g\) is the acceleration due to gravity, and \(H_P\) is the pressure scale height.

- **Critical value of vertical turbulent diffusivity**:
\[
\nu_{\text{cr}}=\frac{v_h^3}{4l_hN_{\text{BV}}^2}\sim 10^{-4}v_h^3 [\text{CGS}]z [\text{Mm}]
\]

These formulas help us understand the physical processes and characteristics in the solar turbulent tachocline.