Search

We investigate the linear and non-linear properties of the Goldreich–Schubert–Fricke (GSF) instability in stellar radiative zones with arbitrary local (radial and latitudinal) differential rotation. This instability may lead to turbulence that contributes to the redistribution of angular momentum and chemical composition in stars. In our local Boussinesq model, we investigate varying the orientation of the shear with respect to the 'effective gravity', which we describe using the angle ϕ. We first perform an axisymmetric linear analysis to explore the effects of varying ϕ on the local stability of arbitrary differential rotations. We then explore the non-linear hydrodynamical evolution in three dimensions using a modified shearing box. The model exhibits both diffusive GSF instability and a non-diffusive instability that occurs when the Solberg-Høiland criteria are violated. We observe the non-linear development of strong zonal jets ('layering' in the angular momentum) with a preferred orientation in both cases, which can considerably enhance turbulent transport. By varying ϕ, we find instability with mixed radial and latitudinal shears transports angular momentum more efficiently (particularly if adiabatically unstable) than cases with purely radial shear (ϕ = 0). By exploring the dependence on box size, we find the transport properties of the GSF instability to be largely insensitive to this, implying we can meaningfully extrapolate our results to stars. However, there is no preferred length-scale for adiabatic instability, which therefore exhibits strong box-size dependence. These instabilities may contribute to the missing angular momentum transport required in red giant and subgiant stars and drive turbulence in the solar tachocline. [ABSTRACT FROM AUTHOR]

We continue our investigation into the non-linear evolution of the Goldreich–Schubert–Fricke (GSF) instability in differentially rotating radiation zones. This instability may be a key player in transporting angular momentum in stars and giant planets, but its non-linear evolution remains mostly unexplored. In a previous paper we considered the equatorial instability, whereas here we simulate the instability at a general latitude for the first time. We adopt a local Cartesian Boussinesq model in a modified shearing box for most of our simulations, but we also perform some simulations with stress-free, impenetrable, radial boundaries. We first revisit the linear instability and derive some new results, before studying its non-linear evolution. The instability is found to behave very differently compared with its behaviour at the equator. In particular, here we observe the development of strong zonal jets ('layering' in the angular momentum), which can considerably enhance angular momentum transport, particularly in axisymmetric simulations. The jets are, in general, tilted with respect to the local gravity by an angle that corresponds initially with that of the linear modes, but which evolves with time and depends on the strength of the flow. The instability transports angular momentum much more efficiently (by several orders of magnitude) than it does at the equator, and we estimate that the GSF instability could contribute to the missing angular momentum transport required in both red giant and subgiant stars. It could also play a role in the long-term evolution of the solar tachocline and the atmospheric dynamics of hot Jupiters. [ABSTRACT FROM AUTHOR]

We present an investigation into the non-linear evolution of the Goldreich–Schubert–Fricke (GSF) instability using axisymmetric and 3D simulations near the equator of a differentially rotating radiation zone. This instability may provide an important contribution to angular momentum transport in stars and planets. We adopt a local Boussinesq Cartesian shearing box model, which represents a small patch of a differentially rotating stellar radiation zone. Complementary simulations are also performed with stress-free, impenetrable boundaries in the local radial direction. The linear and non-linear evolution of the equatorial axisymmetric instability is formally equivalent to the salt fingering instability. This is no longer the case in 3D, but we find that the instability behaves non-linearly in a similar way to salt fingering. Axisymmetric simulations – and those in 3D with short dimensions along the local azimuthal direction – quickly develop strong jets along the rotation axis, which inhibit the instability and lead to predator-prey-like temporal dynamics. In 3D, the instability initially produces homogeneous turbulence and enhanced momentum transport, though in some cases jets form on a much longer time-scale. We propose and validate numerically a simple theory for non-linear saturation of the GSF instability and its resulting angular momentum transport. This theory is straightforward to implement in stellar evolution codes incorporating rotation. We estimate that the GSF instability could contribute towards explaining the missing angular momentum transport required in red giant stars, and play a role in the long-term evolution of the solar tachocline. [ABSTRACT FROM AUTHOR]

Stellar magnetic activity – which has been observed in a diverse set of stars including the Sun – originates via a magnetohydrodynamic dynamo mechanism working in stellar interiors. The full set of magnetohydrodynamic equations governing stellar dynamos is highly complex, and so direct numerical simulation is currently out of reach computationally. An understanding of the bifurcation structure, likely to be found in the partial differential equations governing such dynamos, is vital if we are to understand the activity of solar-like stars and its evolution with varying stellar parameters such as rotation rate. Low-order models are an important aid to this understanding, and can be derived either as approximations of the governing equations themselves or by using bifurcation theory to obtain systems with the desired structure. We use normal-form theory to derive a third-order model with robust behaviour. The model is able to reproduce many of the basic types of behaviour found in observations of solar-type stars. In the appropriate parameter regime, a chaotic modulation of the basic cycle is present, together with varying periods of low activity such as that observed during the solar Maunder minima. [ABSTRACT FROM AUTHOR]