Frequent Links
Trochoidal wave
In fluid dynamics, a trochoidal wave or Gerstner wave is an exact solution of the Euler equations for periodic surface gravity waves. It describes a progressive wave of permanent form on the surface of an incompressible fluid of infinite depth. The free surface of this wave solution is an inverted (upsidedown) trochoid – with sharper crests and flat troughs. This wave solution was discovered by Gerstner in 1802, and rediscovered independently by Rankine in 1863.
The flow field associated with the trochoidal wave is not irrotational: it has vorticity. The vorticity is of such a specific strength and vertical distribution that the trajectories of the fluid parcels are closed circles. This is in contrast with the usual experimental observation of Stokes drift associated with the wave motion. Also the phase speed is independent of the trochoidal wave's amplitude, in contrast with other nonlinear wave theories (like those of the Stokes wave and cnoidal wave) and observations. For these reasons – and for the fact that solutions for finite fluid depth are lacking – trochoidal waves are of limited use for engineering applications.
In computer graphics, the rendering of realisticlooking ocean waves can be done by use of socalled Gerstner waves. This is a multicomponent and multidirectional extension of the traditional Gerstner wave, often using fast Fourier transforms to make (realtime) animation feasible.^{[1]}
Description of classical trochoidal wave
Using a Lagrangian specification of the flow field, the motion of fluid parcels is:^{[2]}
 <math>
\begin{align} X(a,b,t) &= a + \frac{\mbox{e}^{kb}}{k} \sin k(a+ct), \\ Y(a,b,t) &= b  \frac{\mbox{e}^{kb}}{k} \cos k(a+ct), \end{align}
</math> where <math>x=X(a,b,t)</math> and <math>y=Y(a,b,t)</math> are the positions of the fluid parcels in the <math>(x,y)</math> plane at time <math>t</math>, with <math>x</math> the horizontal coordinate and <math>y</math> the vertical coordinate (positive upward, in the direction opposing gravity). The Lagrangian coordinates <math>(a,b)</math> label the fluid parcels, with <math>(x,y)=(a,b)</math> the centres of the circular orbits – around which the corresponding fluid parcel moves with constant speed <math>c\,\exp(kb).</math> Further <math>k=2\pi/\lambda</math> is the wavenumber (and <math>\lambda</math> the wavelength), while <math>c</math> is the phase speed with which the wave propagates in the <math>x</math>direction. The phase speed satisfies the dispersion relation:
 <math>c^2=\frac{g}{k},</math>
which is independent of the wave nonlinearity (i.e. does not depend on the wave height <math>H</math>), and this phase speed <math>c</math> the same as for Airy's linear waves in deep water.
The free surface is a line of constant pressure, and is found to correspond with a line <math>b=b_s</math>, where <math>b_s</math> is a (nonpositive) constant. For <math>b_s=0</math> the highest waves occur, with a cuspshaped crest. Note that the highest (irrotational) Stokes wave has a crest angle of 120°, instead of the 0° for the rotational trochoidal wave.^{[3]}
The wave height of the trochoidal wave is <math>H=(2/k)\exp(kb_s).</math> The wave is periodic in the <math>x</math>direction, with wavelength <math>\lambda;</math> and also periodic in time with period <math>T=\lambda/c=\sqrt{2\pi\lambda/g}.</math>
In computer graphics
A multicomponent and multidirectional extension of the Lagrangian description of the freesurface motion – as used in Gerstner's trochoidal wave – is used in computer graphics for the simulation of ocean waves.^{[1]} For the classical Gerstner wave the fluid motion exactly satisfies the nonlinear, incompressible and inviscid flow equations below the free surface. However, the extended Gerstner waves do in general not satisfy these flow equations exactly (although they satisfy them approximately, i.e. for the linearised Lagrangian description by potential flow). This description of the ocean can be programmed very efficiently by use of the fast Fourier transform (FFT). Moreover, the resulting ocean waves look realistically, as a result of the nonlinear deformation of the free surface (resulting from the Lagrangian specification of the motion).
The mathematical description of the freesurface in these Gerstner waves can be as follows:^{[1]} the horizontal coordinates are denoted as <math>x</math> and <math>z</math>, and the vertical coordinate is <math>y</math>. The mean level of the free surface is at <math>y=0</math> and the positive <math>y</math>direction is upward, opposing the Earth's gravity of strength <math>g.</math> The free surface is described parametrically as a function of the parameters <math>\alpha</math> and <math>\beta,</math> as well as of time <math>t.</math> The parameters are connected to the meansurface points <math>(x,y,z)=(\alpha,0,\beta)</math> around which the fluid parcels at the wavy surface orbit. The free surface is specified through <math>x=\xi(\alpha,\beta,t),</math> <math>y=\zeta(\alpha,\beta,t)</math> and <math>z=\eta(\alpha,\beta,t)</math> with:
 <math>
\begin{align}
\xi &= \alpha  \sum_{m=1}^M \frac{k_{x,m}}{k_m}\, \frac{a_m}{\tanh\left( k_m\, h \right)}\, \sin\left( \theta_m \right), \\ \eta &= \beta  \sum_{m=1}^M \frac{k_{z,m}}{k_m}\, \frac{a_m}{\tanh\left( k_m\, h \right)}\, \sin\left( \theta_m \right), \\ \zeta &= \sum_{m=1}^M a_m\, \cos\left( \theta_m \right) \quad\text{and} \\ \theta_m &= k_{x,m}\, \alpha + k_{z,m}\, \beta  \omega_m\, t  \phi_m,
\end{align} </math>
where <math>\tanh</math> is the hyperbolic tangent function, <math>M</math> is the number of wave components considered, <math>a_m</math> is the amplitude of component <math>{m=1\dots M}</math> and <math>\phi_m</math> its phase. Further <math>{k_m=\scriptstyle\sqrt{(k_{x,m}^2+k_{z,m}^2)}}</math> is its wavenumber and <math>\omega_m</math> its angular frequency. The latter two, <math>k_m</math> and <math>\omega_m,</math> can not be chosen independently but are related through the dispersion relation:
 <math>\omega_m^2 = g\, k_m\, \tanh \left( k_m\, h \right),</math>
with <math>h</math> the mean water depth. In deep water (<math>h\to\infty</math>) the hyperbolic tangent goes to one: <math>{\tanh(k_m\,h)\to 1.}</math> The components <math>k_{x,m}</math> and <math>k_{z,m}</math> of the horizontal wavenumber vector <math>\boldsymbol{k}_m</math> determine the wave propagation direction of component <math>m.</math>
The choice of the various parameters <math>a_m, k_{x,m}, k_{z,m}</math> and <math>\phi_m</math> for <math>{m=1\dots{M},}</math> and a certain mean depth <math>h</math> determines the form of the ocean surface. A clever choice is needed in order to exploit the possibility of fast computation by means of the FFT. See e.g. Tessendorf (2001) for a description how to do this. Most often, the wavenumbers are chosen on a regular grid in <math>(k_x,k_z)</math>space. Thereafter, the amplitudes <math>a_m</math> and phases <math>\phi_m</math> are chosen randomly in accord with the variancedensity spectrum of a certain desired sea state. Finally, by FFT, the ocean surface can be constructed in such a way that it is periodic both in space and time, enabling tiled rendering.
In rendering, also the normal vector <math>\boldsymbol{n}</math> to the surface is often needed. These can be computed using the cross product (<math>\times</math>) as:
 <math>
\boldsymbol{n} = \frac{\partial\boldsymbol{s}}{\partial\alpha} \times \frac{\partial\boldsymbol{s}}{\partial\beta} \quad \text{with} \quad \boldsymbol{s}(\alpha,\beta,t) = \begin{pmatrix} \xi(\alpha,\beta,t) \\ \zeta(\alpha,\beta,t) \\ \eta(\alpha,\beta,t) \end{pmatrix}.
</math>
The unit normal vector then is <math>\boldsymbol{e}_n=\boldsymbol{n}/\\boldsymbol{n}\,</math> with <math>\\boldsymbol{n}\</math> the norm of <math>\boldsymbol{n}.</math>
Notes
 ^ ^{a} ^{b} ^{c} Tessendorf (2001)
 ^ Lamb (1994, §251)
 ^ Stokes, G.G. (1880), "Supplement to a paper on the theory of oscillatory waves", Mathematical and Physical Papers, Volume I, Cambridge University Press, pp. 314–326, OCLC 314316422
References
 Gerstner, F.J. (1802), "Theorie der Wellen", Abhandlunger der Königlichen Böhmischen Geselschaft der Wissenschaften (Prague). Reprinted in: Annalen der Physik 32(8), pp. 412–445, 1809.
 Craik, A.D.D. (2004), "The origins of water wave theory", Annual Review of Fluid Mechanics 36: 1–28, Bibcode:2004AnRFM..36....1C, doi:10.1146/annurev.fluid.36.050802.122118
 Lamb, H. (1994), Hydrodynamics (6th ed.), Cambridge University Press, §251, ISBN 9780521458689, OCLC 30070401 Originally published in 1879, the 6th extended edition appeared first in 1932.
 Rankine, W.J.M. (1863), "On the exact form of waves near the surface of deep water", Philosophical Transactions of the Royal Society of London 153: 127–138, doi:10.1098/rstl.1863.0006
 Tessendorf, J. (2001), "Simulating ocean water" (PDF), SIGGRAPH 2001 (Simulating Nature: Realistic and Interactive Techniques), 26 pp.
