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
- NSF-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
-
-
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
-
more » « less
Abstract We evaluate energetic electron scattering in pitch angle and energy using realistic magnetic field and density models due to whistler mode chorus waves in Jupiter's magnetosphere and study their dependences on various wave and background parameters. We calculate the bounce‐averaged diffusion coefficients by considering the latitudinal variation of total electron density and ambient magnetic field intensity, using the VIP4 internal magnetic field and CAN current sheet model. The electron phase space density evolution due to chorus waves is simulated at
M shell of 10, using the central wave frequency at0.1 and wave amplitude of 30 pT. Under the typical values of the ratio between the plasma frequency and electron cyclotron frequency, chorus waves could cause fast pitch angle scattering loss of energetic electrons from tens to several hundred keV in several hours, and gradual acceleration of relativistic electrons at several MeV in several days. The electron pitch angle scattering at ~500 keV and the acceleration at several MeV are both enhanced using the latitudinally varying density and VIP4 + CAN magnetic field model compared to the electron evolution using the constant density and dipole magnetic field model. Our sensitivity study indicates that the electron scattering at higher energy is caused by waves at lower frequencies or in a lower‐density background plasma, and the scattering is faster for waves at smaller wave normal angles. The electron diffusion is mainly caused by waves at lower latitudes, but the waves at higher latitudes (>30°) contribute to the electron loss at higher energies (>2 MeV).f ce -
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
-
A theoretical study is made of the stability of propagating internal gravity wave modes along a horizontal stratified fluid layer bounded by rigid walls. The analysis is based on the Floquet eigenvalue problem for infinitesimal perturbations to a wave mode of small amplitude. The appropriate instability mechanism hinges on how the perturbation spatial scale relative to the basic-state wavelength, controlled by a parameter $\mu$ , compares to the basic-state amplitude parameter, $\epsilon \ll 1$ . For $\mu ={O}(1)$ , the onset of instability arises due to perturbations that form resonant triads with the underlying wave mode. For short-scale perturbations such that $\mu \ll 1$ but $\alpha =\mu /\epsilon \gg 1$ , this triad resonance instability reduces to the familiar parametric subharmonic instability (PSI), where triads comprise fine-scale perturbations with half the basic-wave frequency. However, as $\mu$ is further decreased holding $\epsilon$ fixed, higher-frequency perturbations than these two subharmonics come into play, and when $\alpha ={O}(1)$ Floquet modes feature broadband spectrum. This broadening phenomenon is a manifestation of the advection of small-scale perturbations by the basic-wave velocity field. By working with a set of ‘streamline coordinates’ in the frame of the basic wave, this advection can be ‘factored out’. Importantly, when $\alpha ={O}(1)$ PSI is replaced by a novel, multi-mode resonance mechanism which has a stabilising effect that provides an inviscid short-scale cut-off to PSI. The theoretical predictions are supported by numerical results from solving the Floquet eigenvalue problem for a mode-1 basic state.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2020.837
