Energy cascades in rapidly rotating and stratified turbulence within elongated domains

We study forced, rapidly rotating and stably stratified turbulence in an elongated domain using an asymptotic expansion at simultaneously low Rossby number $\mathit {Ro}\ll 1$ and large domain height compared with the energy injection scale, $h=H/\ell _{in}\gg 1$ . The resulting equations depend on the parameter $\lambda =(h \mathit {Ro} )^{-1}$ and the Froude number $\mathit {Fr}$ . An extensive set of direct numerical simulations (DNS) is performed to explore the parameter space $(\lambda,\mathit {Fr})$ . We show that a forward energy cascade occurs in one region of this space, and a split energy cascade outside it. At weak stratification (large $\mathit {Fr}$ ), an inverse cascade is observed for sufficiently large $\lambda$ . At strong stratification (small $\mathit {Fr}$ ) the flow becomes approximately hydrostatic and an inverse cascade is always observed. For both weak and strong stratification, we present theoretical arguments supporting the observed energy cascade phenomenology. Our results shed light on an asymptotic region in the phase diagram of rotating and stratified turbulence, which is difficult to attain by brute-force DNS.

On the Dynamics of Overshooting Convection in Spherical Shells: Effect of Density Stratification and Rotation

Overshooting of turbulent motions from convective regions into adjacent stably stratified zones plays a significant role in stellar interior dynamics, as this process may lead to mixing of chemical species and contribute to the transport of angular momentum and magnetic fields. We present a series of fully nonlinear, three-dimensional (3D) anelastic simulations of overshooting convection in a spherical shell that are focused on the dependence of the overshooting dynamics on the density stratification and the rotation, both key ingredients in stars that however have not been studied systematically together via global simulations. We demonstrate that the overshoot lengthscale is not simply a monotonic function of the density stratification in the convective region, but instead it depends on the ratio of the density stratifications in the two zones. Additionally, we find that the overshoot lengthscale decreases with decreasing Rossby number Ro and scales as Ro0.23 while it also depends on latitude with higher Rossby cases leading to a weaker latitudinal variation. We examine the mean flows arising due to rotation and find that they extend beyond the base of the convection zone into the stable region. Our findings may provide a better understanding of the dynamical interaction between stellar convective and radiative regions, and motivate future studies particularly related to the solar tachocline and the implications of its overlapping with the overshoot region.

Large-scale Vortices in Rapidly Rotating Rayleigh–Bénard Convection at Small Prandtl number

One prominent feature in the atmospheres of Jupiter and Saturn is the appearance of large-scale vortices (LSVs). However, the mechanism that sustains these LSVs remains unclear. One possible mechanism is that these LSVs are driven by rotating convection. Here, we present numerical simulation results on a rapidly rotating Rayleigh–Bénard convection at a small Prandtl number Pr = 0.1 (close to the turbulent Prandtl numbers of Jupiter and Saturn). We identified four flow regimes in our simulation: multiple small vortices, a coexisting large-scale cyclone and anticyclone, large-scale cyclone, and turbulence. The formation of LSVs requires that two conditions be satisfied: the vertical Reynolds number is large ( Re z ≥ 400 ), and the Rossby number is small (Ro ≤ 0.4). A large-scale cyclone first appears when Ro decreases to be smaller than 0.4. When Ro further decreases to be smaller than 0.1, a coexisting large-scale cyclone and anticyclone emerges. We have studied the heat transfer in rapidly rotating convection. The result reveals that the heat transfer is more efficient in the anticyclonic region than in the cyclonic region. Besides, we find that the 2D effect increases and the 3D effect decreases in transporting convective flux as the rotation rate increases. We find that aspect ratio has an effect on the critical Rossby number in the emergence of LSVs. Our results provide helpful insights into understanding the dynamics of LSVs in gas giants.

Lifecycle of a Submesoscale Front Birthed from a Nearshore Internal Bore

The ocean is home to many different submesoscale phenomena, including internal waves, fronts, and gravity currents. Each of these processes entail complex nonlinear dynamics, even in isolation. Here we present shipboard, moored, and remote observations of a submesoscale gravity current front created by a shoaling internal tidal bore in the coastal ocean. The internal bore is observed to flatten as it shoals, leaving behind a gravity current front that propagates significantly slower than the bore. We posit that the generation and separation of the front from the bore is related to particular stratification ahead of the bore, which allows the bore to reach the maximum possible internal wave speed. After the front is calved from the bore, it is observed to propagate as a gravity current for ≈4 hours, with associated elevated turbulent dissipation rates. A strong cross-shore gradient of along-shore velocity creates enhanced vertical vorticity (Rossby number ≈ 40) that remains locked with the front. Lateral shear instabilities develop along the front and may hasten its demise.

Balanced ellipsoidal vortex equilibria in a background shear flow at finite Rossby number

We consider a uniform ellipsoid of potential vorticity (PV), where we exploit analytical solutions derived for a balanced model at the second order in the Rossby number, the next order to quasi-geostrophic (QG) theory, the so-called QG+1 model. We consider this vortex in the presence of an external background shear flow, acting as a proxy for the effect of external vortices. For the QG model the system depends on four parameters, the height-to-width aspect ratio of the vortex, $h/r$ , as well as three parameters characterising the background flow, the strain rate, $\gamma$ , the ratio of the background rotation rate to the strain, $\beta$ , and the angle from which the flow is applied, $\theta$ . However, the QG+1 model also depends on the PV, as well as the Prandtl ratio, $f/N$ ( $f$ and $N$ are the Coriolis and buoyancy frequencies, respectively). For QG and QG+1 we determine equilibria for different values of the background flow parameters for increasing values of the imposed strain rate up to the critical strain rate, $\gamma _c$ , beyond which equilibria do not exist. We also compute the linear stability of this vortex to second-order modes, determining the marginal strain $\gamma _m$ at which ellipsoidal instability erupts. The results show that for QG+1 the most resilient cyclonic ellipsoids are slightly prolate, while anticyclonic ellipsoids tend to be more oblate. The highest values of $\gamma _m$ occur as $\beta \to 1$ . For large values of $f/N$ , changes in the marginal strain rates occur, stabilising anticyclonic ellipsoids while destabilising cyclonic ellipsoids.

Coastal imbalance: generation of oceanic Kelvin waves by atmospheric perturbations

The response of a semi-infinite ocean to a slowly travelling atmospheric perturbation crossing the coast provides a simple example of the breakdown of nearly geostrophic balance induced by a boundary. We examine this response in the linear shallow-water model at small Rossby number $\varepsilon \ll 1$ . Using matched asymptotics, we show that a long Kelvin wave, with $O(\varepsilon ^{-1})$ length scale and $O(\varepsilon )$ amplitude relative to quasi-geostrophic response, is generated as the perturbation crosses the coast. Accounting for this Kelvin wave restores the conservation of mass that is violated in the quasi-geostrophic approximation.

Constraints from Invariant Subtropical Vertical Velocities on the Scalings of Hadley Cell Strength and Downdraft Width with Rotation Rate

Weak-temperature-gradient influences from the tropics and quasigeostrophic influences from the extratropics plausibly constrain the subtropical-mean static stability in terrestrial atmospheres. Because mean descent acting on this static stability is a leading-order term in the thermodynamic balance, a state-invariant static stability would impose constraints on the Hadley cells, which this paper explores in simulations of varying planetary rotation rate. If downdraft-averaged effective heating (the sum of diabatic heating and eddy heat flux convergence) too is invariant, so must be vertical velocity—an “omega governor.” In that case, the Hadley circulation overturning strength and downdraft width must scale identically—the cell can strengthen only by widening or weaken only by narrowing. Semiempirical scalings demonstrate that subtropical eddy heat flux convergence weakens with rotation rate (scales positively) while diabatic heating strengthens (scales negatively), compensating one another if they are of similar magnitude. Simulations in two idealized, dry GCMs with a wide range of planetary rotation rates exhibit nearly unchanging downdraft-averaged static stability, effective heating, and vertical velocity, as well as nearly identical scalings of the Hadley cell downdraft width and strength. In one, eddy stresses set this scaling directly (the Rossby number remains small); in the other, eddy stress and bulk Rossby number changes compensate to yield the same, ~Ω−1/3 scaling. The consistency of this power law for cell width and strength variations may indicate a common driver, and we speculate that Ekman pumping could be the mechanism responsible for this behavior. Diabatic heating in an idealized aquaplanet GCM is an order of magnitude larger than in dry GCMs and reanalyses, and while the subtropical static stability is insensitive to rotation rate, the effective heating and vertical velocity are not.

LARGE EDDY SIMULATIONS OF HIGH ROSSBY NUMBER FLOW IN THE HIGH PRESSURE COMPRESSOR INTER-DISK CAVITY

The focus of the present study is to understand the effect of Rayleigh number on a high Rossby number flow in a high pres- sure compressor (HPC) inter-disk cavity. These cavities form between the compressor disks of a gas turbine engine, and they are an integral part of the internal air cooling system. We perform highly resolved large eddy simulations for two Rayleigh numbers of 0.76 × 108 and 1.54 × 108 at a fixed Rossby number of 4.5 by solving the compressible Navier–Stokes equations. The results show a flow structure dominated by a toroidal vortex in the inner region of the cavity. In the outer region, the flow is observed to move radially outwards by Ekman layers formed on the side disks and to move radially inwards through the central core region of the cavity. An enhancement in the intensity of the radial flares is observed in the outer region of the cavity for the high Rayleigh number case with no perceivable effect in the inner region. The near shroud region is mostly dominated by the centrifugal buoyancy-induced flow and the wall Nusselt number calculated at the shroud is in close agreement with centrifugal buoyancy-induced flow without an axial bore flow.