An official website of the United States government
The two-dimensional stability of vertically sheared inertial oscillations at ocean fronts is explored through a linear stability analysis and nonlinear simulations. Baroclinic effects reduce the minimum frequency of inertia-gravity waves to an extent determined by the balanced Richardson number$${{Ri}}$$of the front. Below a critical value of$${{Ri}}$$, which depends on the strength of the inertial shear, the inertial oscillations become unstable to parametric subharmonic instability (PSI) resulting in growing perturbations that oscillate at half the inertial frequency$$f$$. Since the critical value is always greater than 1, PSI can occur at fronts stable to symmetric instability. Although modest in weak inertial shear, growth rates exceeding$$f/2$$can be achieved for inertial shear greater than or equal to the thermal wind shear. Our formulation allows for non-hydrostatic perturbations and can be applied to initially unstratified geostrophic adjustment problems. We find that PSI will almost totally damp the transient oscillations that arise during geostrophic adjustment. The perturbations gain energy at the expense of the inertial oscillations through ageostrophic shear production. The perturbations then themselves become unstable to secondary Kelvin–Helmholtz instabilities creating a pathway by which the inertial oscillations can be dissipated rapidly. In contrast to symmetric and baroclinic instabilities that draw on a front's kinetic or potential energy, PSI acts to increase the energy stored in the balanced front as the convergence and divergence of the eddy-momentum fluxes set up a secondary circulation in the sense to stand up the front.
more »
« less
The influence of front strength on the development and equilibration of symmetric instability. Part 1. Growth and saturation
Submesoscale fronts with large horizontal buoyancy gradients and$$O(1)$$Rossby numbers are common in the upper ocean. These fronts are associated with large vertical transport and are hotspots for biological activity. Submesoscale fronts are susceptible to symmetric instability (SI) – a form of stratified inertial instability which can occur when the potential vorticity is of the opposite sign to the Coriolis parameter. Here, we use a weakly nonlinear stability analysis to study SI in an idealised frontal zone with a uniform horizontal buoyancy gradient in thermal wind balance. We find that the structure and energetics of SI strongly depend on the front strength, defined as the ratio of the horizontal buoyancy gradient to the square of the Coriolis frequency. Vertically bounded non-hydrostatic SI modes can grow by extracting potential or kinetic energy from the balanced front and the relative importance of these energy reservoirs depends on the front strength and vertical stratification. We describe two limiting behaviours as ‘slantwise convection’ and ‘slantwise inertial instability’ where the largest energy source is the buoyancy flux and geostrophic shear production, respectively. The growing linear SI modes eventually break down through a secondary shear instability, and in the process transport considerable geostrophic momentum. The resulting breakdown of thermal wind balance generates vertically sheared inertial oscillations and we estimate the amplitude of these oscillations from the stability analysis. We finally discuss broader implications of these results in the context of current parameterisations of SI.
more »
« less
- Award ID(s):
- 1851450
- PAR ID:
- 10631544
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 926
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In Part 1 (Wienkers, Thomas & Taylor,J. Fluid Mech., vol. 926, 2021, A6), we described the theory for linear growth and weakly nonlinear saturation of symmetric instability (SI) in the Eady model representing a broad frontal zone. There, we found that both the fraction of the balanced thermal wind mixed down by SI and the primary source of energy are strongly dependent on the front strength, defined as the ratio of the horizontal buoyancy gradient to the square of the Coriolis frequency. Strong fronts with steep isopycnals develop a flavour of SI we call ‘slantwise inertial instability’ by extracting kinetic energy from the background flow and rapidly mixing down the thermal wind profile. In contrast, weak fronts extract more potential energy from the background density profile, which results in ‘slantwise convection.’ Here, we extend the theory from Part 1 using nonlinear numerical simulations to focus on the adjustment of the front following saturation of SI. We find that the details of adjustment and amplitude of the induced inertial oscillations depend on the front strength. While weak fronts develop narrow frontlets and excite small-amplitude vertically sheared inertial oscillations, stronger fronts generate large inertial oscillations and produce bore-like gravity currents that propagate along the top and bottom boundaries. The turbulent dissipation rate in these strong fronts is large, highly intermittent and intensifies during periods of weak stratification. We describe each of these mechanisms and energy pathways as the front evolves towards the final adjusted state, and in particular focus on the effect of varying the dimensionless front strength.more » « less
-
This work studies two-dimensional fixed-flux Rayleigh–Bénard convection with periodic boundary conditions in both horizontal and vertical directions and analyses its dynamics using numerical continuation, secondary instability analysis and direct numerical simulation. The fixed-flux constraint leads to time-independent elevator modes with a well-defined amplitude. Secondary instability of these modes leads to tilted elevator modes accompanied by horizontal shear flow. For$$Pr=1$$, where$$Pr$$is the Prandtl number, a subsequent subcritical Hopf bifurcation leads to hysteresis behaviour between this state and a time-dependent direction-reversing state, followed by a global bifurcation leading to modulated travelling waves without flow reversal. Single-mode equations reproduce this moderate Rayleigh number behaviour well. At high Rayleigh numbers, chaotic behaviour dominated by modulated travelling waves appears. These transitions are characteristic of high wavenumber elevator modes since the vertical wavenumber of the secondary instability is linearly proportional to the horizontal wavenumber of the elevator mode. At a low$$Pr$$, relaxation oscillations between the conduction state and the elevator mode appear, followed by quasi-periodic and chaotic behaviour as the Rayleigh number increases. In the high$$Pr$$regime, the large-scale shear weakens, and the flow shows bursting behaviour that can lead to significantly increased heat transport or even intermittent stable stratification.more » « less
-
Wind turbines operate in the atmospheric boundary layer (ABL), where Coriolis effects are present. As wind turbines with larger rotor diameters are deployed, the wake structures that they create in the ABL also increase in length. Contemporary utility-scale wind turbines operate at rotor diameter-based Rossby numbers, the non-dimensional ratio between inertial and Coriolis forces, of$$\mathcal {O}(100)$$where Coriolis effects become increasingly relevant. Coriolis forces provide a direct forcing on the wake, but also affect the ABL base flow, which indirectly influences wake evolution. These effects may constructively or destructively interfere because both the magnitude and sign of the direct and indirect Coriolis effects depend on the Rossby number, turbulence and buoyancy effects in the ABL. Using large eddy simulations, we investigate wake evolution over a wide range of Rossby numbers relevant to offshore wind turbines. Through an analysis of the streamwise and lateral momentum budgets, we show that Coriolis effects have a small impact on the wake recovery rate, but Coriolis effects induce significant wake deflections which can be parsed into two regimes. For high Rossby numbers (weak Coriolis forcing), wakes deflect clockwise in the northern hemisphere. By contrast, for low Rossby numbers (strong Coriolis forcing), wakes deflect anti-clockwise. Decreasing the Rossby number results in increasingly anti-clockwise wake deflections. The transition point between clockwise and anti-clockwise deflection depends on the direct Coriolis forcing, pressure gradients and turbulent fluxes in the wake. At a Rossby number of 125, Coriolis deflections are comparable to wake deflections induced by$${\sim} 20^{\circ }$$of yaw misalignment.more » « less
-
We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity$$\tilde {w}_b$$compared with the rise velocity in quiescent flow$$w_b$$. The decrease in mean rise velocity of bubbles$$\tilde {w}_b/w_b$$is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number$$Fr=u'/\sqrt {gd}$$, where$$u'$$is the root-mean-square velocity fluctuations,$$g$$is gravity and$$d$$is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For$$Fr > 0.5$$, we confirm the scaling$$\tilde {w}_b / w_b \propto 1 / Fr$$, as proposed previously by Ruthet al.(J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2021.680
