skip to main content

Title: Stokes waves with constant vorticity: I. Numerical computation

Periodic traveling waves are numerically computed in a constant vorticity flow subject to the force of gravity. The Stokes wave problem is formulated via a conformal mapping as a nonlinear pseudodifferential equation, involving a periodic Hilbert transform for a strip, and solved by the Newton‐GMRES method. For strong positive vorticity, in the finite or infinite depth, overhanging profiles are found as the amplitude increases and tend to a touching wave, whose surface contacts itself at the trough line, enclosing an air bubble; numerical solutions become unphysical as the amplitude increases further and make a gap in the wave speed versus amplitude plane; another touching wave takes over and physical solutions follow along the fold in the wave speed versus amplitude plane until they ultimately tend to an extreme wave, which exhibits a corner at the crest. Touching waves connected to zero amplitude are found to approach the limiting Crapper wave as the strength of positive vorticity increases unboundedly, while touching waves connected to the extreme waves approach the rigid body rotation of a fluid disk.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Studies in Applied Mathematics
Page Range / eLocation ID:
p. 162-189
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Numerically computed with high accuracy are periodic travelling waves at the free surface of a two-dimensional, infinitely deep, and constant vorticity flow of an incompressible inviscid fluid, under gravity, without the effects of surface tension. Of particular interest is the angle the fluid surface of an almost extreme wave makes with the horizontal. Numerically found are the following. (i) There is a boundary layer where the angle rises sharply from $0^\circ$ at the crest to a local maximum, which converges to $30.3787\ldots ^\circ$ , independently of the vorticity, as the amplitude increases towards that of the extreme wave, which displays a corner at the crest with a $30^\circ$ angle. (ii) There is an outer region where the angle descends to $0^\circ$ at the trough for negative vorticity, while it rises to a maximum, greater than $30^\circ$ , and then falls sharply to $0^\circ$ at the trough for large positive vorticity. (iii) There is a transition region where the angle oscillates about $30^\circ$ , resembling the Gibbs phenomenon. Numerical evidence suggests that the amplitude and frequency of the oscillations become independent of the vorticity as the wave profile approaches the extreme form. 
    more » « less
  2. Abstract

    This paper considers two-dimensional steady solitary waves with constant vorticity propagating under the influence of gravity over an impermeable flat bed. Unlike in previous works on solitary waves, we allow for both internal stagnation points and overhanging wave profiles. Using analytic global bifurcation theory, we construct continuous curves of large-amplitude solutions. Along these curves, either the wave amplitude approaches the maximum possible value, the dimensionless wave speed becomes unbounded, or a singularity develops in a conformal map describing the fluid domain. This is stronger than what one would expect from a straightforward generalization of existing results for periodic waves. We also show that an arbitrary solitary wave of elevation with constant vorticity must be supercritical. The existence proof relies on a novel reformulation of the problem as an elliptic system for two scalar functions in a fixed domain, one describing the conformal map of the fluid region and the other the flow beneath the wave.

    more » « less
  3. Abstract

    Recent work found evidence using aquaplanet experiments that tropical cyclone (TC) size on Earth is limited by the Rhines scale, which depends on the planetary vorticity gradientβ. This study aims to examine how the Rhines scale limits the size of an individual TC. The traditional Rhines scale is first reexpressed as a Rhines speed to characterize how the effect ofβvaries with radius in a vortex whose wind profile is known. The framework is used to define the vortex Rhines scale, which is the transition radius that divides the vortex into a vortex-dominant region at smaller radii, where the axisymmetric circulation is steady, and a wave-dominant region at larger radii, where the circulation stimulates planetary Rossby waves and dissipates. Experiments are performed using a simple barotropic model on aβplane initialized with a TC-like axisymmetric vortex defined using a recently developed theoretical TC wind profile model. The gradientβand initial vortex size are each systematically varied to investigate the detailed responses of the TC-like vortex toβ. Results show that the vortex shrinks toward an equilibrium size that closely follows the vortex Rhines scale. A larger initial vortex relative to its vortex Rhines scale will shrink faster. The shrinking time scale is well described by the vortex Rhines time scale, which is defined as the overturning time scale of the circulation at the vortex Rhines scale and is shown to be directly related to the Rossby wave group velocity. The relationship between our idealized results and the real Earth is discussed. Results may generalize to other eddy circulations, such as the extratropical cyclone.

    Significance Statement

    Tropical cyclones vary in size significantly on Earth, but how large a tropical cyclone could potentially be is still not understood. The variation of the Coriolis parameter with latitude is known to limit the size of turbulent circulations, but its effect on tropical cyclones has not been studied. This study derives a new parameter related to this concept called the “vortex Rhines scale” and shows in a simple model how and why storms will tend to shrink toward this size. These results help explain why tropical cyclone size tends to increase slowly with latitude on Earth and can help us understand what sets the size of tropical cyclones on Earth in general.

    more » « less
  4. Abstract

    Lidar and radar observations of persistent atmospheric wave activity in the Antarctic atmosphere motivate investigation of generation of acoustic‐gravity waves (AGWs) by vibrations of ice shelves and exploiting their possible ionospheric manifestations as a source of information about the ice shelves' conditions and stability. A mathematical model of the waves radiated by vibrations of a finite area of the lower boundary of the atmosphere is developed in this paper by extending to AGWs an efficient, numerically exact approach that was originally developed in seismology and underwater acoustics. The model represents three‐dimensional wave fields as Fourier integrals of numerical or analytical solutions of a one‐dimensional wave equation and accounts for the source directionality, AGW refraction and diffraction, and the wind‐induced anisotropy of wave dissipation. Application of the model to the generation of atmospheric waves in Antarctica by free vibrations of the Ross Ice Shelf reveals a complex three‐dimensional structure of the AGW field and elucidates the impact of various environmental factors on the wave field. The intricate variation of the wave amplitude with altitude and in the horizontal plane is shaped by the spatial spectrum of the ice surface vibrations and the temperature and wind velocity stratification from the troposphere to the mesosphere. It is found that the waves due to the low‐order modes of the free oscillations of the Ross Ice Shelf, which have periods of the order of several hours, can transport energy to the middle and upper atmosphere in a wide range of directions from near‐horizontal to near‐vertical.

    more » « less
  5. 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