Capillary wave
A capillary wave is a wave traveling along the phase boundary of a fluid, whose dynamics are dominated by the effects of surface tension.
Capillary waves are common in nature and the home, and are often referred to as ripples. The wavelength of capillary waves in water is typically less than a few centimeters.
When generated by light wind in open water, a nautical name for them is "cat's paw" waves, since they may resemble paw prints. Light breezes 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.
A gravity–capillary wave on a fluid interface is influenced by both the effects of surface tension and gravity, as well as by fluid inertia.
Contents
Capillary waves, proper
The dispersion relation for capillary waves is
\[ \omega^2=\frac{\sigma}{\rho+\rho'}\, |k|^3,\] where ω is the angular frequency, σ the surface tension, ρ the density of the heavier fluid, ρ' the density of the lighter fluid and k the wavenumber. The wavelength is \( \lambda=\frac{2 \pi}{k}.\)
Gravity–capillary waves
In general, waves are also affected by gravity and are then called gravity–capillary waves. Their dispersion relation reads, for waves on the interface between two fluids of infinite depth:[1][2]
\[ \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, ρ and ρ‘ are the mass density of the two fluids (ρ > ρ‘). Notice the factor \((\rho-\rho')/(\rho+\rho')\) in the first term is the Atwood number.
Gravity wave regime
For large wavelengths (small k = 2π/λ), 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), 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, an interesting and common situation occurs when 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 λc are dominated by surface tension, and much above by gravity. The value of this wavelength is[1]:
- \[\lambda_c = 2 \pi \sqrt{ \frac{\sigma}{(\rho-\rho') g}}.\]
For the air–water interface, λc is found to be 1.7 cm.[1]
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 caustics which corresponds to the minimal group velocity.[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."[4] The derivation of the general dispersion relation is therefore quite involved.[5]
Therefore, first the assumptions involved are pointed out. 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 σ. For gravity, an assumption is made of the density of the fluids being constant (i.e., incompressibility), and likewise g (waves are not high for gravitation to change appreciably). For surface tension, the deviations from planarity (as measured by derivatives of the surface) are supposed to be small. Both approximations are excellent for common waves.
The last contribution involves the kinetic energies of the fluids, and is the most involved. One must use a hydrodynamic framework to tackle this problem. 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; again, these are typically 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.)
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, ρ and ρ’ 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 ∇Φ and ∇Φ’, respectively. Here Φ(x,y,z,t) and Φ’(x,y,z,t) are velocity potentials.
Three contributions to the energy are involved: the potential energy Vg due to gravity, the potential energy Vst due to the surface tension and the kinetic energy T of the flow. The part Vg due to gravity is the simplest: integrating the potential energy density due to gravity, ρ g z (or ρ’ g z) from a reference height to the position of the surface, z = η(x,y,t)[6]: \[ 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[7]: \[ 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[8]: \[ T= \frac{1}{2} \iint dx\, dy\; \left[ \int_{-\infty}^\eta dz\; \rho\, \left| \bold\nabla \Phi \right|^2 + \int_\eta^{+\infty} dz\; \rho'\, \left| \bold\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 Φ(x,y,z,t) and Φ’(x,y,z,t) must satisfy the Laplace equation[9]: \[\nabla^2 \Phi = 0\] and \(\nabla^2 \Phi' = 0.\) These equations can be solved with the proper boundary conditions: Φ and Φ’ 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=η), the kinetic energy can be written as[8]: \[ 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[7]: \[\eta = a\, \cos\, ( kx - \omega t) = a\, \cos\, \theta ,\] with amplitude a and wave phase θ = kx - ω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[7]: \[\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[7]: \[ \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}{|k|} \text{e}^{-|k|z}\, \omega a\, \sin\, \theta. \end{align} \] Then the contributions to the wave energy, horizontally integrated over one wavelength λ = 2π/k in the x–direction, and over a unit width in the y–direction, become[7][10]: \[ \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}{|k|} 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 Vg and surface tension Vst[11]: \[ L = \frac{1}{4} \left[ (\rho+\rho') \frac{\omega^2}{|k|} - (\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(ω,k) a², so that variation with respect to the only free parameter, a, gives the dispersion relation D(ω,k) = 0.[11] In our case D(ω,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 ) / λ, 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.[12] |
See also
- Capillary action
- Dispersion (water waves)
- Fluid pipe
- Thermal capillary wave
- Two-phase flow
- Ocean surface wave
- Wave-formed ripple
Gallery
- 2006-01-14 Surface waves.jpg
Ripples on water
- Surface waves and water striders.JPG
Ripples on water created by water striders
- Plughole.JPG
Ripples of tapwater over a plughole
- Ripple - in rail.jpg
Notes
- ↑ 1.0 1.1 1.2 Lamb (1994), §267, page 458–460.
- ↑ Dingemans (1997), Section 2.1.1, p. 45.</br>Phillips (1977), Section 3.2, p. 37.
- ↑ Script error
- ↑ R.P. Feynman, R.B. Leighton, and M. Sands (1963). The Feynman Lectures on Physics. Addison-Wesley. Volume I, Chapter 51-4.
- ↑ See e.g. Safran (1994) for a more detailed description.
- ↑ Lamb (1994), §174 and §230.
- ↑ 7.0 7.1 7.2 7.3 7.4 Lamb (1994), §266.
- ↑ 8.0 8.1 Lamb (1994), §61.
- ↑ Lamb (1994), §20
- ↑ Lamb (1994), §230.
- ↑ 11.0 11.1 Script error See section 11.7.
- ↑ Script error Reprinted as Appendix in: Theory of Sound 1, MacMillan, 2nd revised edition, 1894.
References
- Script error
- Script error
- Script error
- Script error
- Script error
External links