capillary wave

{{short description|Wave on the surface of a fluid, dominated by surface tension}}

{{redirect|Rippled waves|the insect|Rippled wave|the electrical phenomenon|Ripple waveform||Ripple (disambiguation)}}

File:2006-01-14 Surface waves.jpg

File:Multy droplets impact.JPG impacts on the interface between water and air.]]

A capillary wave is a wave traveling along the phase boundary of a fluid, whose dynamics and phase velocity are dominated by the effects of surface tension.

Capillary waves are common in nature, and are often referred to as ripples. The wavelength of capillary waves on water is typically less than a few centimeters, with a phase speed in excess of 0.2–0.3 meter/second.

A longer wavelength on a fluid interface will result in gravity–capillary waves which are influenced by both the effects of surface tension and gravity, as well as by fluid inertia. Ordinary gravity waves have a still longer wavelength.

Light breezes upon the surface of water which stir up such small ripples are also sometimes referred to as 'cat's paws'. On the open ocean, much larger ocean surface waves (seas and swells) may result from coalescence of smaller wind-caused ripple-waves.

Dispersion relation

The dispersion relation describes the relationship between wavelength and frequency in waves. Distinction can be made between pure capillary waves – fully dominated by the effects of surface tension – and gravity–capillary waves which are also affected by gravity.

=Capillary waves, proper=

The dispersion relation for capillary waves is

: $\omega^2=\frac{\sigma}{\rho+\rho'}\, |k|^3,$

where $\omega$ is the angular frequency, $\sigma$ the surface tension, $\rho$ the density of the

heavier fluid, $\rho'$ the density of the lighter fluid and $k$ the wavenumber. The wavelength is

$\lambda=\frac{2 \pi}{k}.$

For the boundary between fluid and vacuum (free surface), the dispersion relation reduces to

: $\omega^2=\frac{\sigma}{\rho}\, |k|^3.$

=Gravity–capillary waves=

File:Dispersion capillary.svg

When capillary waves are also affected substantially by gravity, they are called gravity–capillary waves. Their dispersion relation reads, for waves on the interface between two fluids of infinite depth:Lamb (1994), §267, page 458–460.Dingemans (1997), Section 2.1.1, p. 45.
Phillips (1977), Section 3.2, p. 37.

: $\omega^2=|k|\left( \frac{\rho-\rho'}{\rho+\rho'}g+\frac{\sigma}{\rho+\rho'}k^2\right),$

where $g$ is the acceleration due to gravity, $\rho$ and $\rho'$ are the densities of the two fluids $(\rho > \rho')$ . The factor $(\rho-\rho')/(\rho+\rho')$ in the first term is the Atwood number.

==Gravity wave regime==

For large wavelengths (small $k = 2\pi/\lambda$ ), only the first term is relevant and one has gravity waves.

In this limit, the waves have a group velocity half the phase velocity: following a single wave's crest in a group one can see the wave appearing at the back of the group, growing and finally disappearing at the front of the group.

==Capillary wave regime==

Shorter (large $k$ ) waves (e.g. 2 mm for the water–air interface), which are proper capillary waves, do the opposite: an individual wave appears at the front of the group, grows when moving towards the group center and finally disappears at the back of the group. Phase velocity is two thirds of group velocity in this limit.

==Phase velocity minimum==

Between these two limits is a point at which the dispersion caused by gravity cancels out the dispersion due to the capillary effect. At a certain wavelength, the group velocity equals the phase velocity, and there is no dispersion. At precisely this same wavelength, the phase velocity of gravity–capillary waves as a function of wavelength (or wave number) has a minimum. Waves with wavelengths much smaller than this critical wavelength $\lambda_{m}$ are dominated by surface tension, and much above by gravity. The value of this wavelength and the associated minimum phase speed $c_{m}$ are:

: $\lambda_m = 2 \pi \sqrt{ \frac{\sigma}{(\rho-\rho') g}}
\quad \text{and} \quad
c_m = \sqrt{ \frac{2 \sqrt{ (\rho - \rho') g \sigma }}{\rho+\rho'} }.$

For the air–water interface, $\lambda_{m}$ is found to be {{convert|1.7|cm|in|abbr=on}}, and $c_{m}$ is {{convert|0.23|m/s|ft/s|abbr=on}}.

If one drops a small stone or droplet into liquid, the waves then propagate outside an expanding circle of fluid at rest; this circle is a caustic which corresponds to the minimal group velocity.{{cite book |last=Falkovich |first=G. |title=Fluid Mechanics, a short course for physicists |publisher=Cambridge University Press |year=2011 |isbn=978-1-107-00575-4 |no-pp=yes |pages=Section 3.1 and Exercise 3.3}}

==Derivation==

As Richard Feynman put it, "[water waves] that are easily seen by everyone and which are usually used as an example of waves in elementary courses [...] are the worst possible example [...]; they have all the complications that waves can have."R.P. Feynman, R.B. Leighton, and M. Sands (1963). The Feynman Lectures on Physics. Addison-Wesley. Volume I, Chapter 51-4. The derivation of the general dispersion relation is therefore quite involved.See e.g. Safran (1994) for a more detailed description.

There are three contributions to the energy, due to gravity, to surface tension, and to hydrodynamics. The first two are potential energies, and responsible for the two terms inside the parenthesis, as is clear from the appearance of $g$ and $\sigma$ . For gravity, an assumption is made of the density of the fluids being constant (i.e., incompressibility), and likewise $g$ (waves are not high enough for gravitation to change appreciably). For surface tension, the deviations from planarity (as measured by derivatives of the surface) are supposed to be small. For common waves both approximations are good enough.

The third contribution involves the kinetic energies of the fluids. It is the most complicated and calls for a hydrodynamic framework. Incompressibility is again involved (which is satisfied if the speed of the waves is much less than the speed of sound in the media), together with the flow being irrotational – the flow is then potential. These are typically also good approximations for common situations.

The resulting equation for the potential (which is Laplace equation) can be solved with the proper boundary conditions. On one hand, the velocity must vanish well below the surface (in the "deep water" case, which is the one we consider, otherwise a more involved result is obtained, see Ocean surface waves.) On the other, its vertical component must match the motion of the surface. This contribution ends up being responsible for the extra $k$ outside the parenthesis, which causes all regimes to be dispersive, both at low values of $k$ , and high ones (except around the one value at which the two dispersions cancel out.)

class="toccolours collapsible collapsed" width="90%" style="text-align:left"

!Dispersion relation for gravity–capillary waves on an interface between two semi–infinite fluid domains

Consider two fluid domains, separated by an interface with surface tension. The mean interface position is horizontal. It separates the upper from the lower fluid, both having a different constant mass density,

\rho

and

\rho'

for the lower and upper domain respectively. The fluid is assumed to be inviscid and incompressible, and the flow is assumed to be irrotational. Then the flows are potential, and the velocity in the lower and upper layer can be obtained from

\nabla \phi

and

\nabla \phi'

, respectively. Here

\phi(x,y,z,t)

and

\phi'(x,y,z,t)

are velocity potentials.

Three contributions to the energy are involved: the potential energy $V_{g}$ due to gravity, the potential energy $V_{st}$ due to the surface tension and the kinetic energy $T$ of the flow. The part $V_{g}$ due to gravity is the simplest: integrating the potential energy density due to gravity, $\rho g z$ (or $\rho' g z$ ) from a reference height to the position of the surface, $z = \eta(x,y,t)$ :Lamb (1994), §174 and §230.

: $V_\mathrm{g}
= \iint dx\, dy\; \int_0^\eta dz\; (\rho - \rho') g z
= \frac{1}{2} (\rho-\rho') g \iint dx\, dy\; \eta^2,$

assuming the mean interface position is at $z=0$ .

An increase in area of the surface causes a proportional increase of energy due to surface tension:Lamb (1994), §266.

: $V_\mathrm{st}
= \sigma \iint dx\, dy\;
\left[
\sqrt{ 1 + \left( \frac{\partial \eta}{\partial x} \right)^2
+ \left( \frac{\partial \eta}{\partial y} \right)^2}
- 1
\right]
\approx \frac{1}{2} \sigma \iint dx\, dy\;
\left[
\left( \frac{\partial \eta}{\partial x} \right)^2
+
\left( \frac{\partial \eta}{\partial y} \right)^2
\right],$

where the first equality is the area in this (Monge's) representation, and the second

applies for small values of the derivatives (surfaces not too rough).

The last contribution involves the kinetic energy of the fluid:Lamb (1994), §61.

: $T=
\frac{1}{2} \iint dx\, dy\;
\left[
\int_{-\infty}^\eta dz\; \rho\, \left| \mathbf\nabla \Phi \right|^2
+
\int_\eta^{+\infty} dz\; \rho'\, \left| \mathbf\nabla \Phi' \right|^2
\right].$

Use is made of the fluid being incompressible and its flow is irrotational (often, sensible approximations). As a result, both $\phi(x,y,z,t)$ and $\phi'(x,y,z,t)$ must satisfy the Laplace equation:Lamb (1994), §20

: $\nabla^2 \Phi = 0$ and $\nabla^2 \Phi' = 0.$

These equations can be solved with the proper boundary conditions: $\phi$ and $\phi'$ must vanish well away from the surface (in the "deep water" case, which is the one we consider).

Using Green's identity, and assuming the deviations of the surface elevation to be small (so the z–integrations may be approximated by integrating up to $z=0$ instead of $z = \eta$ ), the kinetic energy can be written as:

: $T \approx
\frac{1}{2} \iint dx\, dy\;
\left[
\rho\, \Phi\, \frac{\partial \Phi }{\partial z}\;
-\;
\rho'\, \Phi'\, \frac{\partial \Phi'}{\partial z}
\right]_{\text{at } z=0}.$

To find the dispersion relation, it is sufficient to consider a sinusoidal wave on the interface, propagating in the x–direction:

: $\eta = a\, \cos\, ( kx - \omega t) = a\, \cos\, \theta ,$

with amplitude $a$ and wave phase $\theta = kx - \omega t$ . The kinematic boundary condition at the interface, relating the potentials to the interface motion, is that the vertical velocity components must match the motion of the surface:

: $\frac{\partial\Phi}{\partial z} = \frac{\partial\eta}{\partial t}$ and $\frac{\partial\Phi'}{\partial z} = \frac{\partial\eta}{\partial t}$ at $z = 0$ .

To tackle the problem of finding the potentials, one may try separation of variables, when both fields can be expressed as:

: $\begin{align}
\Phi(x,y,z,t) & = + \frac{1}{|k$

\text{e}^{+|k|z}\, \omega a\, \sin\, \theta, \\ \Phi'(x,y,z,t)& = - \frac{1}

\text{e}^{-|k|z}\, \omega a\, \sin\, \theta. \end{align}

Then the contributions to the wave energy, horizontally integrated over one wavelength $\lambda = 2\pi/k$ in the x–direction, and over a unit width in the y–direction, become:Lamb (1994), §230.

: $\begin{align}
V_\text{g} &= \frac{1}{4} (\rho-\rho') g a^2 \lambda,
\\
V_\text{st} &= \frac{1}{4} \sigma k^2 a^2 \lambda,
\\
T &= \frac{1}{4} (\rho+\rho') \frac{\omega^2}$

a^2 \lambda. \end{align}

The dispersion relation can now be obtained from the Lagrangian $L = T - V$ , with $V$ the sum of the potential energies by gravity $V_{g}$ and surface tension $V_{st}$ :{{cite book | first=G. B. | last=Whitham | author-link=Gerald B. Whitham | title=Linear and nonlinear waves | publisher = Wiley-Interscience | year=1974 | isbn=0-471-94090-9 }} See section 11.7.

: $L = \frac{1}{4} \left[
(\rho+\rho') \frac{\omega^2}$

- (\rho-\rho') g - \sigma k^2 \right] a^2 \lambda.

For sinusoidal waves and linear wave theory, the phase–averaged Lagrangian is always of the form $L = D(\omega, k) a^{2}$ , so that variation with respect to the only free parameter, $a$ , gives the dispersion relation $D(\omega, k) = 0$ . In our case $D(\omega,k)$ is just the expression in the square brackets, so that the dispersion relation is:

: $\omega^2 = |k| \left( \frac{\rho-\rho'}{\rho+\rho'}\, g + \frac{\sigma}{\rho+\rho'}\, k^2 \right),$

the same as above.

As a result, the average wave energy per unit horizontal area, $(T + V)/\lambda$ , is:

: $\bar{E} = \frac{1}{2}\, \left[ (\rho-\rho')\, g + \sigma k^2 \right]\, a^2.$

As usual for linear wave motions, the potential and kinetic energy are equal (equipartition holds): $T = V$ .{{cite journal | title=On progressive waves | author=Lord Rayleigh (J. W. Strutt) | author-link=Lord Rayleigh | year=1877 | journal=Proceedings of the London Mathematical Society | volume=9 | pages=21–26 | doi=10.1112/plms/s1-9.1.21 | url=https://zenodo.org/record/1447762 }} Reprinted as Appendix in: Theory of Sound 1, MacMillan, 2nd revised edition, 1894.

Gallery

File:Surface waves and water striders.JPG|Ripples on water created by water striders

File:Narvijärvi ripples.jpg|Light breeze ripples in the surface water of a lake

Notes

References

{{cite journal | last=Longuet-Higgins | first=M. S. | author-link=Michael S. Longuet-Higgins | title=The generation of capillary waves by steep gravity waves | journal=Journal of Fluid Mechanics | volume=16 | issue=1 | year=1963 | issn=1469-7645 | pages= 138–159 | doi=10.1017/S0022112063000641 |bibcode = 1963JFM....16..138L | s2cid=119740891 }}
{{cite book | first=H. | last=Lamb | author-link=Horace Lamb | year=1994 | title=Hydrodynamics | publisher=Cambridge University Press | edition=6th| isbn=978-0-521-45868-9 }}
{{cite book| first=O. M. | last=Phillips | author-link=Owen Martin Phillips |title=The dynamics of the upper ocean |publisher=Cambridge University Press | year=1977 | edition=2nd | isbn=0-521-29801-6 }}
{{cite book | title=Water wave propagation over uneven bottoms | first=M. W. |last=Dingemans | year=1997 | series=Advanced Series on Ocean Engineering | volume=13 | publisher=World Scientific, Singapore | pages=2 Parts, 967 pages | isbn=981-02-0427-2 }}
{{cite book | first=Samuel | last=Safran | title=Statistical thermodynamics of surfaces, interfaces, and membranes | publisher=Addison-Wesley | year=1994 }}
{{Cite journal | first1=N. B. | last1=Tufillaro | first2=R. | last2=Ramshankar | first3=J. P. | last3=Gollub | title=Order-disorder transition in capillary ripples | journal=Physical Review Letters | volume=62 | issue=4 | pages=422–425 | year=1989 | doi=10.1103/PhysRevLett.62.422 | pmid=10040229 | bibcode=1989PhRvL..62..422T| url=http://scholarship.haverford.edu/cgi/viewcontent.cgi?article=1062&context=physics_facpubs }}