An official website of the United States government
We use a multiple-scale expansion to average the wave action balance equation over an ensemble of sea-surface velocity fields characteristic of the ocean mesoscale and submesoscale. Assuming that the statistical properties of the flow are stationary and homogeneous, we derive an expression for a diffusivity tensor of surface-wave action density. The small parameter in this expansion is the ratio of surface current speed to gravity wave group speed. For isotropic currents, the action diffusivity is expressed in terms of the kinetic energy spectrum of the flow. A Helmholtz decomposition of the sea-surface currents into solenoidal (vortical) and potential (divergent) components shows that, to leading order, the potential component of the surface velocity field has no effect on the diffusivity of wave action: only the vortical component of the sea-surface velocity results in diffusion of surface-wave action. We validate our analytic results for the action diffusivity by Monte Carlo ray-tracing simulations through an ensemble of stochastic velocity fields.
more »
« less
Frequency diffusion of waves by unsteady flows
The production of broadband frequency spectra from narrowband wave forcing in geophysical flows remains an open problem. Here we consider a related theoretical problem that points to the role of time-dependent vortical flow in producing this effect. Specifically, we apply multi-scale analysis to the transport equation of wave action density in a homogeneous stationary random background flow under the Wentzel–Kramers–Brillouin approximation. We find that, when some time dependence in the mean flow is retained, wave action density diffuses both along and across surfaces of constant frequency in wavenumber–frequency space; this stands in contrast to previous results showing that diffusion occurs only along constant-frequency surfaces when the mean flow is steady. A self-similar random background velocity field is used to show that the magnitude of this frequency diffusion depends non-monotonically on the time scale of variation of the velocity field. Numerical solutions of the ray-tracing equations for rotating shallow water illustrate and confirm our theoretical predictions. Notably, the mean intrinsic wave frequency increases in time, which by wave action conservation implies a concomitant increase of wave energy at the expense of the energy of the background flow.
more »
« less
- Award ID(s):
- 1813891
- PAR ID:
- 10309191
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 905
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Turbulence and mixing in a near-bottom convectively driven flow are examined by numerical simulations of a model problem: a statically unstable disturbance at a slope with inclination $$\unicode[STIX]{x1D6FD}$$ in a stable background with buoyancy frequency $$N$$ . The influence of slope angle and initial disturbance amplitude are quantified in a parametric study. The flow evolution involves energy exchange between four energy reservoirs, namely the mean and turbulent components of kinetic energy (KE) and available potential energy (APE). In contrast to the zero-slope case where the mean flow is negligible, the presence of a slope leads to a current that oscillates with $$\unicode[STIX]{x1D714}=N\sin \unicode[STIX]{x1D6FD}$$ and qualitatively changes the subsequent evolution of the initial density disturbance. The frequency, $$N\sin \unicode[STIX]{x1D6FD}$$ , and the initial speed of the current are predicted using linear theory. The energy transfer in the sloping cases is dominated by an oscillatory exchange between mean APE and mean KE with a transfer to turbulence at specific phases. In all simulated cases, the positive buoyancy flux during episodes of convective instability at the zero-velocity phase is the dominant contributor to turbulent kinetic energy (TKE) although the shear production becomes increasingly important with increasing $$\unicode[STIX]{x1D6FD}$$ . Energy that initially resides wholly in mean available potential energy is lost through conversion to turbulence and the subsequent dissipation of TKE and turbulent available potential energy. A key result is that, in contrast to the explosive loss of energy during the initial convective instability in the non-sloping case, the sloping cases exhibit a more gradual energy loss that is sustained over a long time interval. The slope-parallel oscillation introduces a new flow time scale $$T=2\unicode[STIX]{x03C0}/(N\sin \unicode[STIX]{x1D6FD})$$ and, consequently, the fraction of initial APE that is converted to turbulence during convective instability progressively decreases with increasing $$\unicode[STIX]{x1D6FD}$$ . For moderate slopes with $$\unicode[STIX]{x1D6FD}<10^{\circ }$$ , most of the net energy loss takes place during an initial, short ( $$Nt\approx 20$$ ) interval with periodic convective overturns. For steeper slopes, most of the energy loss takes place during a later, long ( $Nt>100$ ) interval when both shear and convective instability occur, and the energy loss rate is approximately constant. The mixing efficiency during the initial period dominated by convectively driven turbulence is found to be substantially higher (exceeds 0.5) than the widely used value of 0.2. The mixing efficiency at long time in the present problem of a convective overturn at a boundary varies between 0.24 and 0.3.more » « less
-
Abstract The generation of broadband wave energy frequency spectra from narrowband wave forcing in geophysical flows remains a conundrum. In contrast to the long-standing view that nonlinear wave–wave interactions drive the spreading of wave energy in frequency space, recent work suggests that Doppler-shifting by geostrophic flows may be the primary agent. We investigate this possibility by ray tracing a large number of inertia–gravity wave packets through three-dimensional, geostrophically turbulent flows generated either by a quasigeostrophic (QG) simulation or by synthetic random processes. We find that, in all cases investigated, a broadband quasi-stationary inertia–gravity wave frequency spectrum forms, irrespective of the initial frequencies and wave vectors of the packets. The frequency spectrum is well represented by a power law. A possible theoretical explanation relies on the analogy between the kinematic stretching of passive tracer gradients and the refraction of wave vectors. Consistent with this hypothesis, the spectrum of eigenvalues of the background flow velocity gradients predicts a frequency spectrum that is nearly identical to that found by integration of the ray tracing equations.more » « less
-
Wave-Powered Reverse Osmosis (WPRO) represents a promising convergence of ocean energy harvesting and advanced Reverse Osmosis (RO) desalination techniques. The significant fluctuations in pressures and flow rates within the integrated WPRO system present a critical challenge, necessitating an accurate transient model for effective performance estimation. This study presents a two-dimensional transient model based on pressure-correction algorithm to simulate channel flow with membrane boundary conditions under varying inlet conditions. The coupled dynamics of pressure, velocity, and salt concentration are addressed iteratively by decoupling and updating each term separately. The model investigates the performance of RO systems under different input conditions, including constant, sinusoidal, and irregular flow. The results indicate that constant input with higher pressure and lower flow rate achieves a better Recovery Ratio (RR). It is emphasized that for WPRO systems, a fair comparison requires choosing the same average power or pressure when evaluating different inputs. Under equivalent input power, sinusoidal waves result in a lower RR compared to constant inputs due to reduced average pressure. Conversely, under equivalent inlet pressure and flow rate, sinusoidal waves achieve a higher RR than constant inputs due to the phase difference between pressuredriven permeate velocity and diffusion-driven Concentration Polarization (CP). Specifically, sinusoidal inputs with higher frequency and higher amplitude display a higher RR. Additionally, irregular input yields a higher RR than constant inputs, as mean pressure and power can be maintained at levels comparable to those of constant input. The model’s adaptability to diverse flow regimes — from steady to sinusoidal and irregular fluctuations — highlights its potential as a critical tool for optimizing RO desalination processes powered by renewable ocean energy.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2020.837
