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]

The predominant force balance in rapidly rotating planetary cores is between Coriolis, pressure, buoyancy and Lorentz forces. This magnetostrophic balance leads to a Taylor state where the spatially averaged azimuthal Lorentz force is compelled to vanish on cylinders aligned with the rotation axis. Any deviation from this state leads to a torsional oscillation, signatures of which have been observed in the Earth's secular variation and are thought to influence length of day variations via angular momentum conservation. In order to investigate the dynamics of torsional oscillations (TOs), we perform several 3-D dynamo simulations in a spherical shell. We find TOs, identified by their propagation at the correct Alfvén speed, in many of our simulations. We find that the frequency, location and direction of propagation of the waves are influenced by the choice of parameters. Torsional waves are observed within the tangent cylinder and also have the ability to pass through it. Several of our simulations display waves with core traveltimes of 4–6 yr. We calculate the driving terms for these waves and find that both the Reynolds force and ageostrophic convection acting through the Lorentz force are important in driving TOs. [ABSTRACT FROM PUBLISHER]