Search for: All records

Abstract Starting from the fully compressible fluid equations in a plane-parallel atmosphere, we demonstrate that linear internal gravity waves are naturally pseudo-incompressible in the limit that the wave frequencyωis much less than that of surface gravity waves, i.e., $\omega \ll \sqrt{{\mathit{gk}}_h}$, wheregis the gravitational acceleration andk_his the horizontal wavenumber. We accomplish this by performing a formal expansion of the wave functions and the local dispersion relation in terms of a dimensionless frequency $\epsilon =\omega /\sqrt{{\mathit{gk}}_h}$. Further, we show that, in this same low-frequency limit, several forms of the anelastic approximation, including the Lantz–Braginsky–Roberts formulation, poorly reproduce the correct behavior of internal gravity waves. The pseudo-incompressible approximation is achieved by assuming that Eulerian fluctuations of the pressure are small in the continuity equation—whereas, in the anelastic approximation, Eulerian density fluctuations are ignored. In an adiabatic stratification, such as occurs in a convection zone, the two approximations become identical. However, in a stable stratification, the differences between the two approximations are stark and only the pseudo-incompressible approximation remains valid. 

Geophysical and astrophysical fluid flows are typically driven by buoyancy and strongly constrained at large scales by planetary rotation. Rapidly rotating Rayleigh–Bénard convection (RRRBC) provides a paradigm for experiments and direct numerical simulations (DNS) of such flows, but the accessible parameter space remains restricted to moderately fast rotation rates (Ekman numbers$${ {Ek}} \gtrsim 10^{-8}$$), while realistic$${Ek}$$for geo- and astrophysical applications are orders of magnitude smaller. On the other hand, previously derived reduced equations of motion describing the leading-order behaviour in the limit of very rapid rotation ($$ {Ek}\to 0$$) cannot capture finite rotation effects, and the physically most relevant part of parameter space with small but finite$${Ek}$$has remained elusive. Here, we employ the rescaled rapidly rotating incompressible Navier–Stokes equations (RRRiNSE) – a reformulation of the Navier–Stokes–Boussinesq equations informed by the scalings valid for$${Ek}\to 0$$, recently introduced by Julienet al.(2024) – to provide full DNS of RRRBC at unprecedented rotation strengths down to$$ {Ek}=10^{-15}$$and below, revealing the disappearance of cyclone–anticyclone asymmetry at previously unattainable Ekman numbers ($${Ek}\approx 10^{-9}$$). We also identify an overshoot in the heat transport as$${Ek}$$is varied at fixed$$\widetilde { {Ra}} \equiv {Ra}{Ek}^{4/3}$$, where$$Ra$$is the Rayleigh number, associated with dissipation due to ageostrophic motions in the boundary layers. The simulations validate theoretical predictions based on thermal boundary layer theory for RRRBC and show that the solutions of RRRiNSE agree with the reduced equations at very small$${Ek}$$. These results represent a first foray into the vast, largely unexplored parameter space of very rapidly rotating convection rendered accessible by RRRiNSE. 

Turbulence is a widely observed state of fluid flows, characterized by complex, nonlinear interactions between motions across a broad spectrum of length and time scales. While turbulence is ubiquitous, from teacups to planetary atmospheres, oceans, and stars, its manifestations can vary considerably between different physical systems. For instance, three-dimensional turbulent flows display a forward energy cascade from large to small scales, while in two-dimensional turbulence, energy cascades from small to large scales. In a given physical system, a transition between such disparate regimes of turbulence can occur when a control parameter reaches a critical value. The behavior of flows close to such transition points, which separate qualitatively distinct phases of turbulence, has been found to be unexpectedly rich. Here, we survey recent findings on such transitions in highly anisotropic turbulent fluid flows, including turbulence in thin layers and under the influence of rapid rotation. We also review recent work on transitions induced by turbulent fluctuations, such as random reversals and transitions between large-scale vortices and jets, among others. The relevance of these results and their ramifications for future investigations are discussed. 

Convection is a ubiquitous process driving geophysical/astrophysical fluid flows, which are typically strongly constrained by planetary rotation on large scales. A celebrated model of such flows, rapidly rotating Rayleigh-Bénard convection, has been extensively studied in direct numerical simulations (DNS) and laboratory experiments, but the parameter values attainable by state-of-the-art methods are limited to moderately rapid rotation (Ekman numbers Ek≳10−8), while realistic geophysical/astrophysical Ek are significantly smaller. Asymptotically reduced equations of motion, the nonhydrostatic quasi-geostrophic equations (NHQGE), describing the flow evolution in the limit Ek→0, do not apply at finite rotation rates. The geophysical/astrophysical regime of small but finite Ek therefore remains currently inaccessible. Here, we introduce a new, numerically advantageous formulation of the Navier-Stokes-Boussinesq equations informed by the scalings valid for Ek→0, the \textit{Rescaled Rapidly Rotating incompressible Navier-Stokes Equations} (RRRiNSE). We solve the RRRiNSE using a spectral quasi-inverse method resulting in a sparse, fast algorithm to perform efficient DNS in this previously unattainable parameter regime. We validate our results against the literature across a range of Ek and demonstrate that the algorithmic approaches taken remain accurate and numerically stable at Ek as low as 10−15. Like the NHQGE, the RRRiNSE derive their efficiency from adequate conditioning, eliminating spurious growing modes that otherwise induce numerical instabilities at small Ek. We show that the time derivative of the mean temperature is inconsequential for accurately determining the Nusselt number in the stationary state, significantly reducing the required simulation time, and demonstrate that full DNS using RRRiNSE agree with the NHQGE at very small Ek. 

We study structure formation in two-dimensional turbulence driven by an external force, interpolating between linear instability forcing and random stirring, subject to nonlinear damping. Using extensive direct numerical simulations, we uncover a rich parameter space featuring four distinct branches of stationary solutions: large-scale vortices, hybrid states with embedded shielded vortices (SVs) of either sign, and two states composed of many similar SVs. Of the latter, the first is a dense vortex gas where all SVs have the same sign and diffuse across the domain. The second is a hexagonal vortex crystal forming from this gas when the linear instability is sufficiently weak. These solutions coexist stably over a wide parameter range. The late-time evolution of the system from small-amplitude initial conditions is nearly self-similar, involving three phases: initial inverse cascade, random nucleation of SVs from turbulence and, once a critical number of vortices is reached, a phase of explosive nucleation of SVs, leading to a statistically stationary state. The vortex gas is continued in the forcing parameter, revealing a sharp transition towards the crystal state as the forcing strength decreases. This transition is analysed in terms of the diffusivity of individual vortices using ideas from statistical physics. The crystal can also decay via an inverse cascade resulting from the breakdown of shielding or insufficient nonlinear damping acting on SVs. Our study highlights the importance of the forcing details in two-dimensional turbulence and reveals the presence of non-trivial SV states in this system, specifically the emergence and melting of a vortex crystal.