An official website of the United States government
We study stimulated generation – the transfer of energy from balanced flows to existing internal waves – using an asymptotic model that couples barotropic quasi-geostrophic flow and near-inertial waves with $$\text{e}^{\text{i}mz}$$ vertical structure, where $$m$$ is the vertical wavenumber and $$z$$ is the vertical coordinate. A detailed description of the conservation laws of this vertical-plane-wave model illuminates the mechanism of stimulated generation associated with vertical vorticity and lateral strain. There are two sources of wave potential energy, and corresponding sinks of balanced kinetic energy: the refractive convergence of wave action density into anti-cyclones (and divergence from cyclones); and the enhancement of wave-field gradients by geostrophic straining. We quantify these energy transfers and describe the phenomenology of stimulated generation using numerical solutions of an initially uniform inertial oscillation interacting with mature freely evolving two-dimensional turbulence. In all solutions, stimulated generation co-exists with a transfer of balanced kinetic energy to large scales via vortex merging. Also, geostrophic straining accounts for most of the generation of wave potential energy, representing a sink of 10 %–20 % of the initial balanced kinetic energy. However, refraction is fundamental because it creates the initial eddy-scale lateral gradients in the near-inertial field that are then enhanced by advection. In these quasi-inviscid solutions, wave dispersion is the only mechanism that upsets stimulated generation: with a barotropic balanced flow, lateral straining enhances the wave group velocity, so that waves accelerate and rapidly escape from straining regions. This wave escape prevents wave energy from cascading to dissipative scales.
more »
« less
An improved model of near-inertial wave dynamics
The YBJ equation (Young & Ben Jelloul, J. Marine Res. , vol. 55, 1997, pp. 735–766) provides a phase-averaged description of the propagation of near-inertial waves (NIWs) through a geostrophic flow. YBJ is obtained via an asymptotic expansion based on the limit $$\mathit{Bu}\rightarrow 0$$ , where $$\mathit{Bu}$$ is the Burger number of the NIWs. Here we develop an improved version, the YBJ + equation. In common with an earlier improvement proposed by Thomas, Smith & Bühler ( J. Fluid Mech. , vol. 817, 2017, pp. 406–438), YBJ + has a dispersion relation that is second-order accurate in $$\mathit{Bu}$$ . (YBJ is first-order accurate.) Thus both improvements have the same formal justification. But the dispersion relation of YBJ + is a Padé approximant to the exact dispersion relation and with $$\mathit{Bu}$$ of order unity this is significantly more accurate than the power-series approximation of Thomas et al. (2017). Moreover, in the limit of high horizontal wavenumber $$k\rightarrow \infty$$ , the wave frequency of YBJ + asymptotes to twice the inertial frequency $2f$ . This enables solution of YBJ + with explicit time-stepping schemes using a time step determined by stable integration of oscillations with frequency $2f$ . Other phase-averaged equations have dispersion relations with frequency increasing as $$k^{2}$$ (YBJ) or $$k^{4}$$ (Thomas et al. 2017): in these cases stable integration with an explicit scheme becomes impractical with increasing horizontal resolution. The YBJ + equation is tested by comparing its numerical solutions with those of the Boussinesq and YBJ equations. In virtually all cases, YBJ + is more accurate than YBJ. The error, however, does not go rapidly to zero as the Burger number characterizing the initial condition is reduced: advection and refraction by geostrophic eddies reduces in the initial length scale of NIWs so that $$\mathit{Bu}$$ increases with time. This increase, if unchecked, would destroy the approximation. We show, however, that dispersion limits the damage by confining most of the wave energy to low $$\mathit{Bu}$$ . In other words, advection and refraction by geostrophic flows does not result in a strong transfer of initially near-inertial energy out of the near-inertial frequency band.
more »
« less
- Award ID(s):
- 1657041
- PAR ID:
- 10169189
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 876
- ISSN:
- 0022-1120
- Page Range / eLocation ID:
- 428 to 448
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Studies of internal wave-driven mixing in the coastal ocean have been mainly focused on internal tides, while wind-driven near-inertial waves (NIWs) have received less attention in this regard. This study demonstrates a scenario of NIW-driven mixing over the Texas-Louisiana shelf. Supported by a high-resolution simulation over the shelf, the NIWs driven by land-sea breeze radiate downward at a sharp front and enhance the mixing in the bottom boundary layer where the NIWs are focused due to slantwise critical reflection. The criterion for slantwise critical reflection of NIWs is (where ω is the wave frequency, S bot is the bottom slope, and S p is the isopycnal slope) under the assumption that the mean flow is in a thermal wind balance and only varies in the slope-normal direction. The mechanism driving the enhanced mixing is explored in an idealized simulation. During slantwise critical reflection, NIWs are amplified with enhanced shear and periodically destratify a bottom boundary layer via differential buoyancy advection, leading to periodically enhanced mixing. Turbulent transport of tracers is also enhanced during slantwise critical reflection of NIWs, which has implications for bottom hypoxia over the Texas-Louisiana shelf.more » « less
-
Wave turbulence provides a conceptual framework for weakly nonlinear interactions in dispersive media. Dating from five decades ago, applications of wave turbulence theory to oceanic internal waves assigned a leading-order role to interactions characterized by a near equivalence between the group velocity of high-frequency internal waves with the phase velocity of near-inertial waves. This scale-separated interaction leads to a Fokker–Planck (generalized diffusion) equation. More recently, starting four decades ago, this scale-separated paradigm has been investigated using ray tracing methods. These ray methods characterize spectral transport of energy by counting the amplitude and net velocity of wave packets in phase space past a high-wavenumber gate prior to ‘breaking’. This explicitly advective characterization is based on an intuitive assignment and lacks theoretical underpinning. When one takes an estimate of the net spectral drift from the wave turbulence derivation and makes the corresponding assessment, one obtains a prediction of spectral transport that is an order of magnitude larger than either observations or reported ray tracing estimates. Motivated by this contradiction, we report two parallel derivations for transport equations describing the refraction of high-frequency internal waves in a sea of random inertial waves. The first uses standard wave turbulence techniques and the second is an ensemble-averaged packet transport equation characterized by the dispersion of wave packets about a mean drift in the spectral domain. The ensemble-averaged transport equation for ray tracing differs in that it contains the intuitively motivated advective term. We conclude that the aforementioned contradiction between theory, numerics and observations needs to be taken at face value and present a pathway for resolving this contradiction.more » « less
-
Inertia-gravity waves in the atmosphere and ocean are transported and refracted by geostrophic turbulent currents. Provided that the wave group velocity is much greater than the speed of geostrophic turbulent currents, kinetic theory can be used to obtain a comprehensive statistical description of the resulting interaction (Savva et al. , J. Fluid Mech. , vol. 916, 2021, A6). The leading-order process is scattering of wave energy along a surface of constant frequency, $$\omega$$ , in wavenumber space. The constant- $$\omega$$ surface corresponding to the linear dispersion relation of inertia-gravity waves is a cone extending to arbitrarily high wavenumbers. Thus, wave scattering by geostrophic turbulence results in a cascade of wave energy to high wavenumbers on the surface of the constant- $$\omega$$ cone. Solution of the kinetic equations shows establishment of a wave kinetic energy spectrum $$\sim k_h^{-2}$$ , where $$k_h$$ is the horizontal wavenumber.more » « less
-
null (Ed.)Prior mathematical work of Constantin & Iyer ( Commun. Pure Appl. Maths , vol. 61, 2008, pp. 330–345; Ann. Appl. Probab. , vol. 21, 2011, pp. 1466–1492) has shown that incompressible Navier–Stokes solutions possess infinitely many stochastic Lagrangian conservation laws for vorticity, backward in time, which generalize the invariants of Cauchy ( Sciences mathématiques et physique , vol. I, 1815, pp. 33–73) for smooth Euler solutions. We reformulate this theory for the case of wall-bounded flows by appealing to the Kuz'min ( Phys. Lett. A , vol. 96, 1983, pp. 88–90)–Oseledets ( Russ. Math. Surv. , vol. 44, 1989, p. 210) representation of Navier–Stokes dynamics, in terms of the vortex-momentum density associated to a continuous distribution of infinitesimal vortex rings. The Constantin–Iyer theory provides an exact representation for vorticity at any interior point as an average over stochastic vorticity contributions transported from the wall. We point out relations of this Lagrangian formulation with the Eulerian theory of Lighthill (Boundary layer theory. In Laminar Boundary Layers (ed. L. Rosenhead), 1963, pp. 46–113)–Morton ( Geophys. Astrophys. Fluid Dyn. , vol. 28, 1984, pp. 277–308) for vorticity generation at solid walls, and also with a statistical result of Taylor ( Proc. R. Soc. Lond. A , vol. 135, 1932, pp. 685–702)–Huggins ( J. Low Temp. Phys. , vol. 96, 1994, pp. 317–346), which connects dissipative drag with organized cross-stream motion of vorticity and which is closely analogous to the ‘Josephson–Anderson relation’ for quantum superfluids. We elaborate a Monte Carlo numerical Lagrangian scheme to calculate the stochastic Cauchy invariants and their statistics, given the Eulerian space–time velocity field. The method is validated using an online database of a turbulent channel-flow simulation (Graham et al. , J. Turbul. , vol. 17, 2016, pp. 181–215), where conservation of the mean Cauchy invariant is verified for two selected buffer-layer events corresponding to an ‘ejection’ and a ‘sweep’. The variances of the stochastic Cauchy invariants grow exponentially backward in time, however, revealing Lagrangian chaos of the stochastic trajectories undergoing both fluid advection and viscous diffusion.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2019.557
