The validity of the anelastic approximation has recently been questioned in the regime of rapidly-rotating compressible convection in low Prandtl number fluids (Calkins et al. 2015). Given the broad usage and the high computational efficiency of sound-proof approaches in this astrophysically relevant regime, this paper clarifies the conditions for a safe application. The potential of the alternative pseudo-incompressible approximation is investigated, which in contrast to the anelastic approximation is shown to never break down for predicting the point of marginal stability. Its accuracy, however, decreases close to the parameters corresponding to the failure of the anelastic approach, which is shown to occur when the sound-crossing time of the domain exceeds a rotation time scale, i.e. for rotational Mach numbers greater than one. Concerning the supercritical case, which is naturally characterised by smaller rotational Mach numbers, we find that the anelastic approximation does not show unphysical behaviour. Growth rates computed with the linearised anelastic equations converge toward the corresponding fully compressible values as the Rayleigh number increases. Likewise, our fully nonlinear turbulent simulations, produced with our fully compressible and anelastic models and carried out in a highly supercritical, rotating, compressible, low Prandtl number regime show good agreement. However, this nonlinear test example is for only a moderately low convective Rossby number of 0.14.

Shear-induced turbulence could play a significant role in mixing momentum and chemical species in stellar radiation zones, as discussed by Zahn (1974). In this paper we analyze the results of direct numerical simulations of stratified plane Couette flows, in the limit of rapid thermal diffusion, to measure the turbulent diffusivity and turbulent viscosity as a function of the local shear and the local stratification. We find that the stability criterion proposed by Zahn (1974), namely that the product of the gradient Richardson number and the Prandtl number must be smaller than a critical values $(J\Pr)_c$ for instability, adequately accounts for the transition to turbulence in the flow, with $(J\Pr)_c \simeq 0.007$. This result recovers and confirms the prior findings of Prat et al. (2016). Zahn's model for the turbulent diffusivity and viscosity (Zahn 1992), namely that the mixing coefficient should be proportional to the ratio of the thermal diffusivity to the gradient Richardson number, does not satisfactorily match our numerical data when applied as is. It fails (as expected) in the limit of large stratification where the Richardson number exceeds the aforementioned threshold for instability, but it also fails in the limit of low stratification where the turbulent eddy scale becomes limited by the computational domain size. We propose a revised model for turbulent mixing by diffusive stratified shear instabilities, that now properly accounts for both limits, fits our data satisfactorily, and recovers Zahn's 1992 model in the limit of large Reynolds numbers.

Numerical simulations of turbulent Rayleigh-B\'enard convection in an ideal gas, using either the anelastic approximation or the fully compressible equations, are compared. Theoretically, the anelastic approximation is expected to hold in weakly superadiabatic systems with $\epsilon = \Delta T / T_r \ll 1$, where $\Delta T$ denotes the superadiabatic temperature drop over the convective layer and $T_r$ the bottom temperature. Using direct numerical simulations, a systematic comparison of anelastic and fully compressible convection is carried out. With decreasing superadiabaticity $\epsilon$, the fully compressible results are found to converge linearly to the anelastic solution with larger density contrasts generally improving the match. We conclude that in many solar and planetary applications, where the superadiabaticity is expected to be vanishingly small, results obtained with the anelastic approximation are in fact more accurate than fully compressible computations, which typically fail to reach small $\epsilon$ for numerical reasons. On the other hand, if the astrophysical system studied contains $\epsilon\sim O(1)$ regions, such as the solar photosphere, fully compressible simulations have the advantage of capturing the full physics. Interestingly, even in weakly superadiabatic regions, like the bulk of the solar convection zone, the errors introduced by using artificially large values for $\epsilon$ for efficiency reasons remain moderate. If quantitative errors of the order of $10\%$ are acceptable in such low $\epsilon$ regions, our work suggests that fully compressible simulations can indeed be computationally more efficient than their anelastic counterparts.

Giant planets like Jupiter and Saturn feature strong zonal wind patterns on their surfaces. Although several different mechanisms that may drive these jets have been proposed over the last decades, the origin of the zonal winds is still unclear. Here, we explore the possibility that the interplay of planetary rotation with the compression and expansion of the convecting fluid can drive multiple deep zonal jets by a compressional Rhines-type mechanism, as originally proposed by Ingersoll and Pollard (1982). In a certain limit, this deep mechanism is shown to be mathematically analogous to the classical Rhines mechanism possibly operating at cloud level. Jets are predicted to occur on a compressional Rhines length $l_R = (2 \Omega \langle H_\rho^{-1} \rangle v_{jet}^{-1} )^{-1/2}$, where $\Omega$ is the angular velocity, $\langle H_\rho^{-1} \rangle$ is the mean inverse density scale height and $v_{jet}$ is the typical jet velocity. Two-dimensional numerical simulations using the anelastic approximation reveal that this mechanism robustly generates jets of the predicted width, and that it typically dominates the dynamics in systems deeper than $O(l_R)$. Potential vorticity staircases are observed to form spontaneously and are typically accompanied by unstably stratified buoyancy staircases. The mechanism only operates at large rotation rates, exceeding those typically reached in three-dimensional simulations of deep convection in spherical shells. Applied to Jupiter and Saturn, the compressional Rhines scaling reasonably fits the available observations. Interestingly, even weak vertical density variations such as those in the Earth core can give rise to a large number of jets, leading to fundamentally different flow structures than predicted by the Boussinesq models typically used in this context.