skip to main content


Title: 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
NSF-PAR ID:
10169189
Author(s) / Creator(s):
;
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
  1. We investigate the effect of constant-vorticity background shear on the properties of wavetrains in deep water. Using the methodology of Fokas ( A Unified Approach to Boundary Value Problems , 2008, SIAM), we derive a higher-order nonlinear Schrödinger equation in the presence of shear and surface tension. We show that the presence of shear induces a strong coupling between the carrier wave and the mean-surface displacement. The effects of the background shear on the modulational instability of plane waves is also studied, where it is shown that shear can suppress instability, although not for all carrier wavelengths in the presence of surface tension. These results expand upon the findings of Thomas et al.  ( Phys. Fluids , vol. 24 (12), 2012, 127102). Using a modification of the generalized Lagrangian mean theory in Andrews & McIntyre ( J. Fluid Mech. , vol. 89, 1978, pp. 609–646) and approximate formulas for the velocity field in the fluid column, explicit, asymptotic approximations for the Lagrangian and Stokes drift velocities are obtained for plane-wave and Jacobi elliptic function solutions of the nonlinear Schrödinger equation. Numerical approximations to particle trajectories for these solutions are found and the Lagrangian and Stokes drift velocities corresponding to these numerical solutions corroborate the theoretical results. We show that background currents have significant effects on the mean transport properties of waves. In particular, certain combinations of background shear and carrier wave frequency lead to the disappearance of mean-surface mass transport. These results provide a possible explanation for the measurements reported in Smith ( J. Phys. Oceanogr. , vol. 36, 2006, pp. 1381–1402). Our results also provide further evidence of the viability of the modification of the Stokes drift velocity beyond the standard monochromatic approximation, such as recently proposed in Breivik et al.  ( J. Phys. Oceanogr. , vol. 44, 2014, pp. 2433–2445) in order to obtain a closer match to a range of complex ocean wave spectra. 
    more » « less
  2. 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
  3. Abstract

    Fronts and near-inertial waves (NIWs) are energetic motions in the upper ocean that have been shown to interact and provide a route for kinetic energy (KE) dissipation of balanced oceanic flows. In this paper, we study these KE exchanges using an idealized model consisting of a two-dimensional geostrophically balanced front undergoing strain-induced semigeostrophic frontogenesis and internal wave (IW) vertical modes. The front–IW KE exchanges are quantified separately during two frontogenetic stages: an exponential sharpening stage that is characterized by a low Rossby number and is driven by the imposed strain (i.e., mesoscale frontogenesis), followed by a superexponential sharpening stage that is characterized by anRossby number and is driven by the convergence of the secondary circulation (i.e., submesoscale frontogenesis). It is demonstrated that high-frequency IWs quickly escape the frontal zone and are very efficient at extracting KE from the imposed geostrophic strain field through the deformation shear production (DSP). Part of the extracted KE is then converted to wave potential energy. On the contrary, NIWs remain locked to the frontal zone and readily exchange energy with the ageostrophic frontal circulation. During the exponential stage, NIWs extract KE from the geostrophic strain through DSP and transfer it to the frontal secondary circulation via the ageostrophic shear production (AGSP) mechanism. During the superexponential stage, a newly identified mechanism, convergence production (CP), plays an important role in the NIW KE budget. The CP transfers KE from the convergent ageostrophic secondary circulation to the NIWs and largely cancels out the KE loss due to the AGSP. This CP may explain previous findings of KE transfer enhancement from balanced motions to IWs in frontal regions of realistic ocean models. We provide analytical estimates for the aforementioned energy exchange mechanisms that match well the numerical results. This highlights that the strength of the exchanges strongly depends on the frontal Rossby and Richardson numbers.

    Significance Statement

    Fronts with large horizontal density and velocity gradients are ubiquitous in the upper ocean. They are generated by a process known as frontogenesis, which is often initialized by straining motions of mesoscale balanced circulations. Here we examine the energy exchanges between fronts and internal waves in an idealized configuration, aiming to elucidate the mechanisms that can drain energy from oceanic balanced circulations. We identify a new mechanism for energy transfers from the frontal circulation to near-inertial internal waves called convergence production. This mechanism is especially effective during the later stages of frontogenesis when the convergent ageostrophic secondary circulation that develops is strong.

     
    more » « less
  4. 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
  5. We present the results of an integrated experimental and theoretical investigation of the vertical motion of millimetric droplets bouncing on a vibrating fluid bath. We characterize experimentally the dependence of the phase of impact and contact force between a drop and the bath on the drop’s size and the bath’s vibrational acceleration. This characterization guides the development of a new theoretical model for the coupling between a drop’s vertical and horizontal motion. Our model allows us to relax the assumption of constant impact phase made in models based on the time-averaged trajectory equation of Moláček and Bush ( J. Fluid Mech. , vol. 727, 2013b, pp. 612–647) and obtain a robust horizontal trajectory equation for a bouncing drop that accounts for modulations in the drop’s vertical dynamics as may arise when it interacts with boundaries or other drops. We demonstrate that such modulations have a critical influence on the stability and dynamics of interacting droplet pairs. As the bath’s vibrational acceleration is increased progressively, initially stationary pairs destabilize into a variety of dynamical states including rectilinear oscillations, circular orbits and side-by-side promenading motion. The theoretical predictions of our variable-impact-phase model rationalize our observations and underscore the critical importance of accounting for variability in the vertical motion when modelling droplet–droplet interactions. 
    more » « less