Английская Википедия:Airy wave theory

Содержание

Текст

Текст

Текст

Текст

Текст

Текст

Текст

Текст

Текст

Шаблон:Short description In fluid dynamics, Airy wave theory (often referred to as linear wave theory) gives a linearised description of the propagation of gravity waves on the surface of a homogeneous fluid layer. The theory assumes that the fluid layer has a uniform mean depth, and that the fluid flow is inviscid, incompressible and irrotational. This theory was first published, in correct form, by George Biddell Airy in the 19th century.

Airy wave theory is often applied in ocean engineering and coastal engineering for the modelling of random sea states – giving a description of the wave kinematics and dynamics of high-enough accuracy for many purposes. Further, several second-order nonlinear properties of surface gravity waves, and their propagation, can be estimated from its results. Airy wave theory is also a good approximation for tsunami waves in the ocean, before they steepen near the coast.

This linear theory is often used to get a quick and rough estimate of wave characteristics and their effects. This approximation is accurate for small ratios of the wave height to water depth (for waves in shallow water), and wave height to wavelength (for waves in deep water).

Description

Файл:Sine wave amplitude.svg

Wave characteristics.

Файл:Dispersion gravity 1.svg

Dispersion of gravity waves on a fluid surface. Phase and group velocity divided by Шаблон:Math as a function of Шаблон:Math. A: phase velocity, B: group velocity, C: phase and group velocity Шаблон:Math valid in shallow water. Drawn lines: based on dispersion relation valid in arbitrary depth. Dashed lines: based on dispersion relation valid in deep water.

Airy wave theory uses a potential flow (or velocity potential) approach to describe the motion of gravity waves on a fluid surface. The use of (inviscid and irrotational) potential flow in water waves is remarkably successful, given its failure to describe many other fluid flows where it is often essential to take viscosity, vorticity, turbulence or flow separation into account. This is due to the fact that for the oscillatory part of the fluid motion, wave-induced vorticity is restricted to some thin oscillatory Stokes boundary layers at the boundaries of the fluid domain.

Airy wave theory is often used in ocean engineering and coastal engineering. Especially for random waves, sometimes called wave turbulence, the evolution of the wave statistics – including the wave spectrum – is predicted well over not too long distances (in terms of wavelengths) and in not too shallow water. Diffraction is one of the wave effects which can be described with Airy wave theory. Further, by using the WKBJ approximation, wave shoaling and refraction can be predicted.

Earlier attempts to describe surface gravity waves using potential flow were made by, among others, Laplace, Poisson, Cauchy and Kelland. But Airy was the first to publish the correct derivation and formulation in 1841. Soon after, in 1847, the linear theory of Airy was extended by Stokes for non-linear wave motion – known as Stokes' wave theory – correct up to third order in the wave steepness. Even before Airy's linear theory, Gerstner derived a nonlinear trochoidal wave theory in 1802, which however is not irrotational.

Airy wave theory is a linear theory for the propagation of waves on the surface of a potential flow and above a horizontal bottom. The free surface elevation Шаблон:Math of one wave component is sinusoidal, as a function of horizontal position Шаблон:Mvar and time Шаблон:Mvar:

<math>\eta(x,t) = a \cos \left( kx - \omega t\right)</math>

where

• Шаблон:Mvar is the wave amplitude in metres,
• Шаблон:Math is the cosine function,
• Шаблон:Mvar is the angular wavenumber in radians per metre, related to the wavelength Шаблон:Mvar by Шаблон:Math,
• Шаблон:Mvar is the angular frequency in radians per second, related to the period Шаблон:Mvar and frequency Шаблон:Mvar by Шаблон:Math.

The waves propagate along the water surface with the phase speed Шаблон:Mvar:

<math>c_p = \frac{\omega}{k} = \frac{\lambda}{T}.</math>

The angular wavenumber Шаблон:Mvar and frequency Шаблон:Mvar are not independent parameters (and thus also wavelength Шаблон:Mvar and period Шаблон:Mvar are not independent), but are coupled. Surface gravity waves on a fluid are dispersive waves – exhibiting frequency dispersion – meaning that each wavenumber has its own frequency and phase speed.

Note that in engineering the wave height Шаблон:Mvar – the difference in elevation between crest and trough – is often used:

<math>H = 2 a \quad \text{and} \quad a = \tfrac12 H,</math>

valid in the present case of linear periodic waves.

Файл:Orbital wave motion.svg

Orbital motion under linear waves. The yellow dots indicate the momentary position of fluid particles on their (orange) orbits. The black dots are the centres of the orbits.

Файл:Wave motion-i18n-mod-with-ground.svg

Particle motion in an ocean wave at deep (A) and shallow (B) depths. 1) Propagation direction. 2) Wave crest. 3) Wave trough.

Underneath the surface, there is a fluid motion associated with the free surface motion. While the surface elevation shows a propagating wave, the fluid particles are in an orbital motion. Within the framework of Airy wave theory, the orbits are closed curves: circles in deep water and ellipses in finite depth—with the circles dying out before reaching the bottom of the fluid layer, and the ellipses becoming flatter near the bottom of the fluid layer. So while the wave propagates, the fluid particles just orbit (oscillate) around their average position. With the propagating wave motion, the fluid particles transfer energy in the wave propagation direction, without having a mean velocity. The diameter of the orbits reduces with depth below the free surface. In deep water, the orbit's diameter is reduced to 4% of its free-surface value at a depth of half a wavelength.

In a similar fashion, there is also a pressure oscillation underneath the free surface, with wave-induced pressure oscillations reducing with depth below the free surface – in the same way as for the orbital motion of fluid parcels.

Mathematical formulation of the wave motion

Flow problem formulation

The waves propagate in the horizontal direction, with coordinate Шаблон:Mvar, and a fluid domain bound above by a free surface at Шаблон:Math, with Шаблон:Mvar the vertical coordinate (positive in the upward direction) and Шаблон:Mvar being time. The level Шаблон:Math corresponds with the mean surface elevation. The impermeable bed underneath the fluid layer is at Шаблон:Math. Further, the flow is assumed to be incompressible and irrotational – a good approximation of the flow in the fluid interior for waves on a liquid surface – and potential theory can be used to describe the flow. The velocity potential Шаблон:Math is related to the flow velocity components Шаблон:Mvar and Шаблон:Mvar in the horizontal (Шаблон:Mvar) and vertical (Шаблон:Mvar) directions by:

<math> u_x = \frac{\partial\Phi}{\partial x} \quad \text{and} \quad u_z = \frac{\partial\Phi}{\partial z}.</math>

Then, due to the continuity equation for an incompressible flow, the potential Шаблон:Math has to satisfy the Laplace equation: Шаблон:NumBlk

Boundary conditions are needed at the bed and the free surface in order to close the system of equations. For their formulation within the framework of linear theory, it is necessary to specify what the base state (or zeroth-order solution) of the flow is. Here, we assume the base state is rest, implying the mean flow velocities are zero.

The bed being impermeable, leads to the kinematic bed boundary-condition: Шаблон:NumBlk

In case of deep water – by which is meant infinite water depth, from a mathematical point of view – the flow velocities have to go to zero in the limit as the vertical coordinate goes to minus infinity: Шаблон:Math.

At the free surface, for infinitesimal waves, the vertical motion of the flow has to be equal to the vertical velocity of the free surface. This leads to the kinematic free-surface boundary-condition: Шаблон:NumBlk

If the free surface elevation Шаблон:Math was a known function, this would be enough to solve the flow problem. However, the surface elevation is an extra unknown, for which an additional boundary condition is needed. This is provided by Bernoulli's equation for an unsteady potential flow. The pressure above the free surface is assumed to be constant. This constant pressure is taken equal to zero, without loss of generality, since the level of such a constant pressure does not alter the flow. After linearisation, this gives the dynamic free-surface boundary condition: Шаблон:NumBlk

Because this is a linear theory, in both free-surface boundary conditions – the kinematic and the dynamic one, equations (Шаблон:EquationNote) and (Шаблон:EquationNote) – the value of Шаблон:Math and Шаблон:Math at the fixed mean level Шаблон:Math is used.

Solution for a progressive monochromatic wave

For a propagating wave of a single frequency – a monochromatic wave – the surface elevation is of the form:

<math>\eta = a \cos ( k x - \omega t ).</math>

The associated velocity potential, satisfying the Laplace equation (1) in the fluid interior, as well as the kinematic boundary conditions at the free surface (2), and bed (3), is:

<math>\Phi = \frac{\omega}{k} a \frac{\cosh k (z+h) }{\sinh k h} \sin ( k x - \omega t),</math>

with Шаблон:Math and Шаблон:Math the hyperbolic sine and hyperbolic cosine function, respectively. But Шаблон:Mvar and Шаблон:Math also have to satisfy the dynamic boundary condition, which results in non-trivial (non-zero) values for the wave amplitude Шаблон:Mvar only if the linear dispersion relation is satisfied:

<math>\omega^2 = g k \tanh k h ,</math>

with Шаблон:Math the hyperbolic tangent. So angular frequency Шаблон:Mvar and wavenumber Шаблон:Mvar – or equivalently period Шаблон:Mvar and wavelength Шаблон:Mvar – cannot be chosen independently, but are related. This means that wave propagation at a fluid surface is an eigenproblem. When Шаблон:Mvar and Шаблон:Mvar satisfy the dispersion relation, the wave amplitude Шаблон:Mvar can be chosen freely (but small enough for Airy wave theory to be a valid approximation).

Table of wave quantities

In the table below, several flow quantities and parameters according to Airy wave theory are given. The given quantities are for a bit more general situation as for the solution given above. Firstly, the waves may propagate in an arbitrary horizontal direction in the Шаблон:Math plane. The wavenumber vector is Шаблон:Math, and is perpendicular to the cams of the wave crests. Secondly, allowance is made for a mean flow velocity Шаблон:Math, in the horizontal direction and uniform over (independent of) depth Шаблон:Mvar. This introduces a Doppler shift in the dispersion relations. At an Earth-fixed location, the observed angular frequency (or absolute angular frequency) is Шаблон:Mvar. On the other hand, in a frame of reference moving with the mean velocity Шаблон:Math (so the mean velocity as observed from this reference frame is zero), the angular frequency is different. It is called the intrinsic angular frequency (or relative angular frequency), denoted Шаблон:Mvar. So in pure wave motion, with Шаблон:Math, both frequencies Шаблон:Mvar and Шаблон:Mvar are equal. The wave number Шаблон:Mvar (and wavelength Шаблон:Mvar) are independent of the frame of reference, and have no Doppler shift (for monochromatic waves).

The table only gives the oscillatory parts of flow quantities – velocities, particle excursions and pressure – and not their mean value or drift. The oscillatory particle excursions Шаблон:Math and Шаблон:Math are the time integrals of the oscillatory flow velocities Шаблон:Math and Шаблон:Math respectively.

Water depth is classified into three regimes:

Файл:Deep and shallow water waves.pdf

Visualization of deep and shallow water waves by relating wavelength to depth to bed.

• deep water – for a water depth larger than half the wavelength, Шаблон:Math, the phase speed of the waves is hardly influenced by depth (this is the case for most wind waves on the sea and ocean surface),
• shallow water – for a water depth smaller than 5% of the wavelength, Шаблон:Math, the phase speed of the waves is only dependent on water depth, and no longer a function of period or wavelength; and
• intermediate depth – all other cases, Шаблон:Math, where both water depth and period (or wavelength) have a significant influence on the solution of Airy wave theory.

In the limiting cases of deep and shallow water, simplifying approximations to the solution can be made. While for intermediate depth, the full formulations have to be used.

Properties of gravity waves on the surface of deep water, shallow water and at intermediate depth, according to Airy wave theory
quantity	symbol	units	deep water (Шаблон:Math)	shallow water (Шаблон:Math)	intermediate depth (all Шаблон:Mvar and Шаблон:Mvar)
surface elevation	<math>\eta(\mathbf{x},t)</math>	m	<math>a \cos \theta(\mathbf{x},t)</math>
wave phase	<math>\theta(\mathbf{x},t)</math>	rad	<math>\mathbf{k}\cdot\mathbf{x} - \omega t</math>
observed angular frequency	<math>\omega</math>	rad·s−1	<math> \left( \omega - \mathbf{k}\cdot\mathbf{U} \right)^2 = \bigl( \Omega(k) \bigr)^2 \quad \text{ with } \quad k= |\mathbf{k}| </math>
intrinsic angular frequency	<math>\sigma</math>	rad·s−1	<math> \quad \sigma^2 = \bigl( \Omega(k) \bigr)^2 \quad \text{ with } \quad \sigma = \omega - \mathbf{k}\cdot\mathbf{U} </math>
unit vector in the wave propagation direction	<math> \mathbf{e}_k</math>	–	<math>\frac{\mathbf{k}}{k}</math>
Шаблон:Anchordispersion relation	<math>\Omega(k)</math>	rad·s−1	<math>\Omega(k) = \sqrt{g k}</math>	<math>\Omega(k) = k \sqrt{g h}</math>	<math>\Omega(k) = \sqrt{g k \tanh k h}</math>
phase speed	<math>c_p=\frac{\Omega(k)}{k}</math>	m·s−1	<math>\sqrt{\frac{g}{k}} = \frac{g}{\sigma}</math>	<math>\sqrt{g h}</math>	<math>\sqrt{\frac{g}{k} \tanh k h}</math>
group speed	<math>c_g = \frac{\partial\Omega}{\partial k}</math>	m·s−1	<math>\tfrac{1}{2} \sqrt{\frac{g}{k}} = \tfrac{1}{2} \frac{g}{\sigma}</math>	<math>\sqrt{g h}</math>	<math>\tfrac{1}{2} c_p \left( 1 + k h \frac{1 - \tanh^2 k h}{\tanh k h} \right)</math>
ratio	<math> \frac{c_g}{c_p}</math>	–	<math>\tfrac{1}{2}</math>	<math>1</math>	<math>\tfrac{1}{2} \left( 1 + k h \frac{1 - \tanh^2 k h}{\tanh k h} \right)</math>
horizontal velocity	<math>\mathbf{u}_x(\mathbf{x},z,t)</math>	m·s−1	<math>\mathbf{e}_k \sigma a\, e^{k z} \cos \theta</math>	<math>\mathbf{e}_k \sqrt{\frac{g}{h}} a \cos \theta</math>	<math>\mathbf{e}_k \sigma a \frac{\cosh k (z+h) }{\sinh k h} \cos \theta</math>
vertical velocity	<math>u_z(\mathbf{x},z,t)</math>	m·s−1	<math>\sigma a\, e^{k z} \sin \theta</math>	<math>\sigma a \frac{z + h}{h} \sin \theta</math>	<math>\sigma a \frac{\sinh k (z+h) }{\sinh k h} \sin \theta</math>
horizontal particle excursion	<math>\boldsymbol{\xi}_x(\mathbf{x},z,t)</math>	m	<math>-\mathbf{e}_k a\, e^{k z} \sin \theta</math>	<math>-\mathbf{e}_k \frac{1}{k h} a \sin \theta</math>	<math>-\mathbf{e}_k a \frac{\cosh k (z+h) }{\sinh k h} \sin \theta</math>
vertical particle excursion	<math>\xi_z(\mathbf{x},z,t)</math>	m	<math>a\, e^{k z} \cos \theta</math>	<math>a \frac{z + h}{h} \cos \theta</math>	<math>a \frac{\sinh k (z+h) }{\sinh k h} \cos \theta</math>
pressure oscillation	<math>p(\mathbf{x},z,t)</math>	N·m−2	<math>\rho g a\, e^{k z} \cos \theta</math>	<math>\rho g a \cos \theta</math>	<math>\rho g a \frac{\cosh k (z+h) }{\cosh k h} \cos \theta</math>

Шаблон:Anchor Surface tension effects

Файл:Dispersion capillary.svg

Dispersion of gravity–capillary waves on the surface of deep water. Phase and group velocity divided by Шаблон:Math as a function of inverse relative wavelength Шаблон:Math.
Blue lines (A): phase velocity Шаблон:Mvar, Red lines (B): group velocity Шаблон:Mvar.
Drawn lines: gravity–capillary waves.
Dashed lines: gravity waves.
Dash-dot lines: pure capillary waves.
Note: Шаблон:Math is surface tension in this graph.

Due to surface tension, the dispersion relation changes to:

<math> \Omega^2(k) = \left( g + \frac{\gamma}{\rho} k^2 \right) k\, \tanh k h ,</math>

with Шаблон:Mvar the surface tension in newtons per metre. All above equations for linear waves remain the same, if the gravitational acceleration Шаблон:Mvar is replaced by

<math>\tilde{g} = g + \frac{\gamma}{\rho} k^2.</math>

As a result of surface tension, the waves propagate faster. Surface tension only has influence for short waves, with wavelengths less than a few decimeters in case of a water–air interface. For very short wavelengths – 2 mm or less, in case of the interface between air and water – gravity effects are negligible. Note that surface tension can be altered by surfactants.

The group velocity Шаблон:Math of capillary waves – dominated by surface tension effects – is greater than the phase velocity Шаблон:Math. This is opposite to the situation of surface gravity waves (with surface tension negligible compared to the effects of gravity) where the phase velocity exceeds the group velocity.

Interfacial waves

Surface waves are a special case of interfacial waves, on the interface between two fluids of different density.

Two layers of infinite depth

Consider two fluids separated by an interface, and without further boundaries. Then their dispersion relation Шаблон:Math is given through

<math> \Omega^2(k) = |k| \left( \frac{\rho-\rho'}{\rho+\rho'} g + \frac{\gamma}{\rho+\rho'} k^2 \right),</math>

where Шаблон:Mvar and Шаблон:Math are the densities of the two fluids, below (Шаблон:Mvar) and above (Шаблон:Math) the interface, respectively. Further γ is the surface tension on the interface.

For interfacial waves to exist, the lower layer has to be heavier than the upper one, Шаблон:Math. Otherwise, the interface is unstable and a Rayleigh–Taylor instability develops.

Two layers between horizontal rigid planes

Файл:Interfacial wave rigid boundaries.svg

Wave motion on the interface between two layers of inviscid homogeneous fluids of different density, confined between horizontal rigid boundaries (at the top and bottom). The motion is forced by gravity. The upper layer has mean depth Шаблон:Math and density Шаблон:Math, while the lower layer has mean depth Шаблон:Mvar and density Шаблон:Mvar. The wave amplitude is Шаблон:Mvar, the wavelength is denoted by Шаблон:Mvar (related to the wavenumber Шаблон:Mvar by Шаблон:Math), the gravitational acceleration by Шаблон:Mvar and the phase speed by Шаблон:Mvar (with Шаблон:Math).

For two homogeneous layers of fluids, of mean thickness Шаблон:Mvar below the interface and Шаблон:Math above – under the action of gravity and bounded above and below by horizontal rigid walls – the dispersion relationship Шаблон:Math for gravity waves is provided by:

<math> \Omega^2(k) = \frac{g k (\rho - \rho')}{\rho \coth k h + \rho' \coth k h'},</math>

where again Шаблон:Mvar and Шаблон:Math are the densities below and above the interface, while Шаблон:Math is the hyperbolic cotangent function. For the case Шаблон:Math is zero this reduces to the dispersion relation of surface gravity waves on water of finite depth Шаблон:Mvar.

Two layers bounded above by a free surface

In this case the dispersion relation allows for two modes: a barotropic mode where the free surface amplitude is large compared with the amplitude of the interfacial wave, and a baroclinic mode where the opposite is the case – the interfacial wave is higher than and in antiphase with the free surface wave. The dispersion relation for this case is of a more complicated form.

Second-order wave properties

Several second-order wave properties, ones that are quadratic in the wave amplitude Шаблон:Mvar, can be derived directly from Airy wave theory. They are of importance in many practical applications, such as forecasts of wave conditions. Using a WKBJ approximation, second-order wave properties also find their applications in describing waves in case of slowly varying bathymetry, and mean-flow variations of currents and surface elevation. As well as in the description of the wave and mean-flow interactions due to time and space-variations in amplitude, frequency, wavelength and direction of the wave field itself.

Table of second-order wave properties Шаблон:Anchor

In the table below, several second-order wave properties – as well as the dynamical equations they satisfy in case of slowly varying conditions in space and time – are given. More details on these can be found below. The table gives results for wave propagation in one horizontal spatial dimension. Further on in this section, more detailed descriptions and results are given for the general case of propagation in two-dimensional horizontal space.

Second-order quantities and their dynamics, using results of Airy wave theory
quantity	symbol	units	formula
mean wave-energy density per unit horizontal area	<math>E</math>	J·m−2	<math>E = \tfrac12 \rho g a^2</math>
radiation stress or excess horizontal momentum flux due to the wave motion	<math>S_{xx}</math>	N·m−1	<math>S_{xx} = \left( 2 \frac{c_g}{c_p} - \tfrac12 \right) E</math>
wave action	<math>\mathcal{A}</math>	J·s·m−2	<math>\mathcal{A} = \frac{E}{\sigma} = \frac{E}{\omega - k U}</math>
mean mass-flux due to the wave motion or the wave pseudo-momentum	<math>M</math>	kg·m−1·s−1	<math>M = \frac{E}{c_p} = k \frac{E}{\sigma}</math>
mean horizontal mass-transport velocity	<math>\tilde{U}</math>	m·s−1	<math>\tilde{U} = U + \frac{M}{\rho h} = U + \frac{E}{\rho h c_p}</math>
Stokes drift	<math>\bar{u}_S</math>	m·s−1	<math>\bar{u}_S = \tfrac12 \sigma k a^2 \frac{\cosh 2k(z+h)}{\sinh^2 kh}</math>
wave-energy propagation		J·m−2·s−1	<math>\frac{\partial E}{\partial t} + \frac{\partial}{\partial x} \bigl( (U + c_g) E \bigr) + S_{xx} \frac{\partial U}{\partial x} =0</math>
wave action conservation		J·m−2	<math>\frac{\partial \mathcal{A}}{\partial t} + \frac{\partial}{\partial x} \bigl( (U + c_g) \mathcal{A} \bigr) = 0</math>
wave-crest conservation		rad·m−1·s−1	<math>\frac{\partial k}{\partial t} + \frac{\partial \omega}{\partial x} = 0 \quad \text{with} \quad \omega = \Omega(k) + k U</math>
mean mass conservation		kg·m−2·s−1	<math>\frac{\partial}{\partial t}( \rho h ) + \frac{\partial}{\partial x} \left( \rho h \tilde{U} \right) = 0</math>
mean horizontal-momentum evolution		N·m−2	<math>\frac{\partial}{\partial t}\left( \rho h \tilde{U} \right) + \frac{\partial}{\partial x} \left( \rho h {\tilde U}^2 + \tfrac12 \rho g h^2 + S_{xx} \right) = \rho g h \frac{\partial d}{\partial x}</math>

The last four equations describe the evolution of slowly varying wave trains over bathymetry in interaction with the mean flow, and can be derived from a variational principle: Whitham's averaged Lagrangian method. In the mean horizontal-momentum equation, Шаблон:Math is the still water depth, that is, the bed underneath the fluid layer is located at Шаблон:Math. Note that the mean-flow velocity in the mass and momentum equations is the mass transport velocity Шаблон:Mvar, including the splash-zone effects of the waves on horizontal mass transport, and not the mean Eulerian velocity (for example, as measured with a fixed flow meter).

Wave energy density

Wave energy is a quantity of primary interest, since it is a primary quantity that is transported with the wave trains. As can be seen above, many wave quantities like surface elevation and orbital velocity are oscillatory in nature with zero mean (within the framework of linear theory). In water waves, the most used energy measure is the mean wave energy density per unit horizontal area. It is the sum of the kinetic and potential energy density, integrated over the depth of the fluid layer and averaged over the wave phase. Simplest to derive is the mean potential energy density per unit horizontal area Шаблон:Math of the surface gravity waves, which is the deviation of the potential energy due to the presence of the waves:

<math>\begin{align}

E_\text{pot} &= \overline{\int_{-h}^{\eta} \rho gz\,\mathrm{d}z} - \int_{-h}^0 \rho gz\, \mathrm{d}z \\[6px]

            &= \overline{\tfrac12\rho g\eta^2} = \tfrac14 \rho ga^2.

\end{align}</math>

The overbar denotes the mean value (which in the present case of periodic waves can be taken either as a time average or an average over one wavelength in space).

The mean kinetic energy density per unit horizontal area Шаблон:Math of the wave motion is similarly found to be:

<math>\begin{align}

E_\text{kin} &= \overline{\int_{-h}^0 \tfrac12 \rho \left[ \left| \mathbf{U} + \mathbf{u}_x \right|^2 + u_z^2 \right]\, \mathrm{d}z} - \int_{-h}^0 \tfrac12 \rho \left| \mathbf{U} \right|^2\, \mathrm{d}z \\[6px]

            &= \tfrac14 \rho \frac{\sigma^2}{k \tanh k h}a^2, 

\end{align}</math>

with Шаблон:Mvar the intrinsic frequency, see the table of wave quantities. Using the dispersion relation, the result for surface gravity waves is:

<math>E_\text{kin} = \tfrac14 \rho g a^2.</math>

As can be seen, the mean kinetic and potential energy densities are equal. This is a general property of energy densities of progressive linear waves in a conservative system. Adding potential and kinetic contributions, Шаблон:Math and Шаблон:Math, the mean energy density per unit horizontal area Шаблон:Mvar of the wave motion is:

<math>E = E_\text{pot} + E_\text{kin} = \tfrac12 \rho g a^2.</math>

In case of surface tension effects not being negligible, their contribution also adds to the potential and kinetic energy densities, giving

<math> E_\text{pot} = E_\text{kin} = \tfrac14 \left( \rho g + \gamma k^2 \right) a^2, </math>

so

<math> E = E_\text{pot} + E_\text{kin} = \tfrac12 \left( \rho g + \gamma k^2 \right) a^2, </math>

with Шаблон:Mvar the surface tension.

Wave action, wave energy flux and radiation stress

In general, there can be an energy transfer between the wave motion and the mean fluid motion. This means, that the wave energy density is not in all cases a conserved quantity (neglecting dissipative effects), but the total energy density – the sum of the energy density per unit area of the wave motion and the mean flow motion – is. However, there is for slowly varying wave trains, propagating in slowly varying bathymetry and mean-flow fields, a similar and conserved wave quantity, the wave action Шаблон:Math:

<math>\frac{\partial \mathcal{A}}{\partial t} + \nabla\cdot\left[ \left(\mathbf{U}+\mathbf{c}_g\right) \mathcal{A}\right] = 0,</math>

with Шаблон:Math the action flux and Шаблон:Math the group velocity vector. Action conservation forms the basis for many wind wave models and wave turbulence models. It is also the basis of coastal engineering models for the computation of wave shoaling. Expanding the above wave action conservation equation leads to the following evolution equation for the wave energy density:

<math>\frac{\partial E}{\partial t} + \nabla\cdot\left[\left( \mathbf{U}+\mathbf{c}_g\right) E \right] + \boldsymbol{S}:\left(\nabla\mathbf{U}\right) = 0,</math>

with:

• Шаблон:Math is the mean wave energy density flux,
• Шаблон:Mvar is the radiation stress tensor and
• Шаблон:Math is the mean-velocity shear rate tensor.

In this equation in non-conservation form, the Frobenius inner product Шаблон:Math is the source term describing the energy exchange of the wave motion with the mean flow. Only in the case that the mean shear-rate is zero, Шаблон:Math, the mean wave energy density Шаблон:Mvar is conserved. The two tensors Шаблон:Mvar and Шаблон:Math are in a Cartesian coordinate system of the form:

<math>

 \begin{align}
 \boldsymbol{S} &= \begin{pmatrix} S_{xx} & S_{xy} \\ S_{yx} & S_{yy} \end{pmatrix}
   = \boldsymbol{I} \left( \frac{c_g}{c_p} - \frac12 \right) E
   + \frac{1}{k^2} \begin{pmatrix} k_x k_x & k_x k_y \\[2ex] k_y k_x & k_y k_y \end{pmatrix} \frac{c_g}{c_p} E,
 \\[6px]
 \boldsymbol{I} &= \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}, 
 \\[6px]
 \nabla \mathbf{U} &= 
   \begin{pmatrix}
     \displaystyle \frac{\partial U_x}{\partial x} & \displaystyle \frac{\partial U_y}{\partial x}
     \\[2ex]
     \displaystyle \frac{\partial U_x}{\partial y} & \displaystyle \frac{\partial U_y}{\partial y}
   \end{pmatrix},
 \end{align}

</math>

with Шаблон:Mvar and Шаблон:Mvar the components of the wavenumber vector Шаблон:Math and similarly Шаблон:Mvar and Шаблон:Mvar the components in of the mean velocity vector Шаблон:Math.

Wave mass flux and wave momentum

The mean horizontal momentum per unit area Шаблон:Math induced by the wave motion – and also the wave-induced mass flux or mass transport – is:

<math>\begin{align}

\mathbf{M} &= \overline{\int_{-h}^\eta \rho \left( \mathbf{U}+\mathbf{u}_x\right)\, \mathrm{d}z} - \int_{-h}^0 \rho \mathbf{U}\, \mathrm{d}z \\[6px]

          &= \frac{E}{c_p} \mathbf{e}_k,

\end{align}</math>

which is an exact result for periodic progressive water waves, also valid for nonlinear waves. However, its validity strongly depends on the way how wave momentum and mass flux are defined. Stokes already identified two possible definitions of phase velocity for periodic nonlinear waves:

• Stokes first definition of wave celerity (S1) – with the mean Eulerian flow velocity equal to zero for all elevations Шаблон:Mvar' below the wave troughs, and
• Stokes second definition of wave celerity (S2) – with the mean mass transport equal to zero.

The above relation between wave momentum Шаблон:Math and wave energy density Шаблон:Mvar is valid within the framework of Stokes' first definition.

However, for waves perpendicular to a coast line or in closed laboratory wave channel, the second definition (S2) is more appropriate. These wave systems have zero mass flux and momentum when using the second definition. In contrast, according to Stokes' first definition (S1), there is a wave-induced mass flux in the wave propagation direction, which has to be balanced by a mean flow Шаблон:Math in the opposite direction – called the undertow.

So in general, there are quite some subtleties involved. Therefore also the term pseudo-momentum of the waves is used instead of wave momentum.

Mass and momentum evolution equations

For slowly varying bathymetry, wave and mean-flow fields, the evolution of the mean flow can de described in terms of the mean mass-transport velocity Шаблон:Math defined as:

<math>\tilde{\mathbf U} = \mathbf{U} + \frac{\mathbf{M}}{\rho h}.</math>

Note that for deep water, when the mean depth Шаблон:Mvar goes to infinity, the mean Eulerian velocity Шаблон:Math and mean transport velocity Шаблон:Math become equal.

The equation for mass conservation is:

<math> \frac{\partial}{\partial t}\left( \rho h \right) + \nabla \cdot \left( \rho h\tilde{\mathbf{U}} \right) = 0,</math>

where Шаблон:Math is the mean water depth, slowly varying in space and time.

Similarly, the mean horizontal momentum evolves as:

<math> \frac{\partial}{\partial t}\left( \rho h \tilde{\mathbf{U}}\right) + \nabla \cdot \left( \rho h \tilde{\mathbf{U}} \otimes \tilde{\mathbf{U}} + \tfrac12\rho gh^2\boldsymbol{I} + \boldsymbol{S} \right) = \rho g h \nabla d,</math>

with Шаблон:Mvar the still-water depth (the sea bed is at Шаблон:Math), Шаблон:Mvar is the wave radiation-stress tensor, Шаблон:Mvar is the identity matrix and Шаблон:Math is the dyadic product:

<math> \tilde{\mathbf{U}} \otimes \tilde{\mathbf{U}} =

 \begin{pmatrix} 
   \tilde{U}_x \tilde{U}_x & \tilde{U}_x \tilde{U}_y \\
   \tilde{U}_y \tilde{U}_x & \tilde{U}_y \tilde{U}_y 
 \end{pmatrix}.</math>

Note that mean horizontal momentum is only conserved if the sea bed is horizontal (the still-water depth Шаблон:Mvar is a constant), in agreement with Noether's theorem.

The system of equations is closed through the description of the waves. Wave energy propagation is described through the wave-action conservation equation (without dissipation and nonlinear wave interactions):

<math> \frac{\partial}{\partial t} \left( \frac{E}{\sigma} \right) + \nabla \cdot \left[ \left( \mathbf{U} +\mathbf{c}_g \right) \frac{E}{\sigma} \right] = 0.</math>

The wave kinematics are described through the wave-crest conservation equation:

<math>\frac{\partial \mathbf{k}}{\partial t} + \nabla \omega = \mathbf{0},</math>

with the angular frequency Шаблон:Mvar a function of the (angular) wavenumber Шаблон:Math, related through the dispersion relation. For this to be possible, the wave field must be coherent. By taking the curl of the wave-crest conservation, it can be seen that an initially irrotational wavenumber field stays irrotational.

Stokes drift

When following a single particle in pure wave motion (Шаблон:Math), according to linear Airy wave theory, a first approximation gives closed elliptical orbits for water particles. However, for nonlinear waves, particles exhibit a Stokes drift for which a second-order expression can be derived from the results of Airy wave theory (see the table above on second-order wave properties). The Stokes drift velocity Шаблон:Math, which is the particle drift after one wave cycle divided by the period, can be estimated using the results of linear theory:

<math>\bar{\mathbf{u}}_S = \tfrac12 \sigma k a^2 \frac{\cosh 2k(z+h)}{\sinh^2 kh} \mathbf{e}_k,</math>

so it varies as a function of elevation. The given formula is for Stokes first definition of wave celerity. When Шаблон:Math is integrated over depth, the expression for the mean wave momentum Шаблон:Math is recovered.

See also

• Boussinesq approximation (water waves) – nonlinear theory for waves in shallow water.
• Capillary wave – surface waves under the action of surface tension
• Cnoidal wave – nonlinear periodic waves in shallow water, solutions of the Korteweg–de Vries equation
• Mild-slope equation – refraction and diffraction of surface waves over varying depth
• Ocean surface wave – real water waves as seen in the ocean and sea
• Stokes wave – nonlinear periodic waves in non-shallow water
• Wave power – using ocean and sea waves for power generation.

References

Historical

• Шаблон:Cite book Also: "Trigonometry, On the Figure of the Earth, Tides and Waves", 396 pp.
• Шаблон:Cite journal
  Reprinted in: Шаблон:Cite book

Further reading

• Шаблон:Cite journal
• Шаблон:Cite book
• Шаблон:Cite book Two parts, 967 pages.
• Шаблон:Cite book Originally published in 1879, the 6th extended edition appeared first in 1932.
• Шаблон:Cite book
• Шаблон:Cite book 504 pp.
• Шаблон:Cite book
• Шаблон:Citation

External links

• Linear theory of ocean surface waves on WikiWaves.
• Water waves at MIT.