# Difference between revisions of "Capillary waves"

m (line break inserted) |
(Derivation made clearer) |
||

Line 53: | Line 53: | ||

(Again, we are neglecting the flow of the gas above for simplicity.) | (Again, we are neglecting the flow of the gas above for simplicity.) | ||

====Wave solutions==== | ====Wave solutions==== | ||

− | Let us | + | Let us suppose the surface of the liquid is described by a traveling plane wave: |

:<math>h(x,y,t)=\eta(t)e^{i\vec{q}\vec{r}},</math> | :<math>h(x,y,t)=\eta(t)e^{i\vec{q}\vec{r}},</math> | ||

− | where <math>\vec{q}=(q_x,q_y)</math> is a two dimensional wave number vector, | + | where <math>\eta(t)=\exp[i\omega t]</math> and <math>\vec{q}=(q_x,q_y)</math> is a two dimensional wave number vector, <math>\vec{r}=(x,y)</math> being the horizontal position. We may take <math>\vec{q}=(q,0)</math> without loss of generality: |

− | We may take <math>\vec{q}=( | ||

− | |||

− | :<math> | + | :<math>h(x,y,t)=\eta(t)e^{i q x}.</math> |

− | |||

− | + | In this case it is easy to perform the integrations involved in the expressions for the energies. The | |

− | To tackle the kinetic energy, suppose the fluid is | + | integration over <math>x</math> can taken over a period of oscillation <math>\lambda=2\pi/q</math>, then |

− | incompressible and its flow is irrotational (often, sensible approximations) - the flow will then be [[potential flow|potential]]: <math>\vec{v}=\nabla\psi</math>, and <math>\psi</math> is a potential (scalar field) which | + | multiplied by the number of oscillations in our very large (in principle, infinite) system: <math>L_x / \lambda</math>. The integration over <math>y</math> trivially yields <math>L_y</math>. Performing the integrations, keeping in mind that only the real part of complex numbers is to be taken as physical, one finds: |

− | must satisfy [[Laplace equation]] <math>\nabla^2\psi=0</math>. | + | |

+ | :<math>E_\mathrm{g} = \frac{A}{2} \frac{\rho g}{2} \eta^2,</math> | ||

+ | :<math>E_\mathrm{st} = \frac{A}{2} \frac{\sigma}{2} q^2 \eta^2,</math> | ||

+ | |||

+ | where <math>A=L_x\times L_y</math> is the area of the system. | ||

+ | |||

+ | To tackle the kinetic energy, suppose the fluid is incompressible and its flow is irrotational (often, sensible approximations) - the flow will then be [[potential flow|potential]]: <math>\vec{v}=\nabla\psi</math>, and <math>\psi</math> is a potential (scalar field) which must satisfy [[Laplace's equation]] <math>\nabla^2\psi=0</math>. | ||

If we try try separation of variables with the potential: | If we try try separation of variables with the potential: | ||

− | :<math>\psi(x,y,z,t)=\xi(t) f(z) e^{i | + | :<math>\psi(x,y,z,t)=\xi(t) f(z) e^{i q x},</math> |

− | with some function of time <math>\xi(t)</math>, and some function of vertical component (height) <math>f(z)</math> | + | with some function of time <math>\xi(t)</math>, and some function of vertical component (height) <math>f(z)</math>. |

− | Laplace equation then requires on the later | + | Laplace's equation then requires on the later |

:<math>f''(z)=-q^2 f(z) .</math> | :<math>f''(z)=-q^2 f(z) .</math> | ||

− | This equation can be solved with the proper boundary conditions: first, <math>\vec{v}</math> must vanish well below the surface (in the "deep water" case, which is the one we consider, otherwise a more general relation holds, which is also well known in | + | This equation can be solved with the proper boundary conditions: first, <math>\vec{v}</math> must vanish well below the surface (in the "deep water" case, which is the one we consider, otherwise a more general relation holds, which is also well known in oceanography). Therefore |

− | oceanography). Therefore | ||

:<math>f(z) = a \, \exp(|q| z) </math>, | :<math>f(z) = a \, \exp(|q| z) </math>, | ||

− | with some constant <math>a</math>. The less trivial condition is the matching between <math>\psi</math> and <math>h</math>: the potential field must correspond to a velocity | + | with some constant <math>a</math>. The less trivial condition is the proper matching between <math>\psi</math> and <math>h</math>: the potential field must correspond to a velocity field that is adjusted to the movement of the surface: |

− | field that is adjusted to the movement of the surface: | + | |

+ | :<math>v_z (z=h) =\partial h/\partial t</math>. | ||

− | + | This means that | |

− | <math>\xi(t)=\eta(t)'</math>, and | + | :<math>\xi(t)=\eta(t)'</math>, and <math>f'(z=h) = 1, </math> |

+ | |||

+ | so that | ||

:<math>f(z) = \exp( -|q|(h-z))/|q| </math>. | :<math>f(z) = \exp( -|q|(h-z))/|q| </math>. | ||

Line 107: | Line 112: | ||

the same dispersion as above if <math>\rho'</math> is neglected. | the same dispersion as above if <math>\rho'</math> is neglected. | ||

+ | |||

==External links== | ==External links== | ||

*[http://en.wikipedia.org/wiki/Capillary_wave Capillary wave entry in wikipedia] | *[http://en.wikipedia.org/wiki/Capillary_wave Capillary wave entry in wikipedia] |

## Revision as of 12:59, 29 October 2008

## Contents

## Thermal capillary waves

Thermal **capillary waves** are oscillations of an interface which are thermal in origin. These take place at the molecular level, where only the contribution due to surface tension is relevant.
Capillary wave theory is a classic account of how thermal fluctuations distort an interface (Ref. 1). It starts from some intrinsic surface that is distorted. In the Monge representation, the surface is given as . An increase in area of the surface causes a proportional increase of energy:

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

A Fourier analysis treatment begins by writing the intrinsic surface as an infinite sum of normal modes:

Since normal modes are orthogonal, the energy is easily expressible as a sum of terms . Each term of the sum is quadratic in the amplitude; hence equipartition holds, according to standard classical statistical mechanics, and the mean energy of each mode will be . Surprisingly, this result leads to a **divergent** surface (the width of the interface is bound to diverge with its area) (Ref 2). This divergence is nevertheless very mild; even for displacements on the order of meters, the deviation of the surface is comparable to the size of the molecules.
Moreover, the introduction of an external field removes this divergence: the action of gravity is sufficient to keep the width fluctuation on the order
of one molecular diameter for areas larger than about 1 mm^{2} (Ref. 2).
The action of gravity is taken into account by integrating the potential energy density due to gravity, from a reference height to the position of the surface, :

(For simplicity, one neglects the density of the gas above, which is often acceptable; otherwise, instead of the density the difference in densities appears).

Recently, a procedure has been proposed to obtain a molecular intrinsic
surface from simulation data (Ref. 3), the intrinsic sampling method. The density profiles obtained
from this surface are, in general, quite different from the usual
*mean density profiles*.

## Gravity-capillary waves

These are ordinary waves excited in an interface, such as ripples on a water surface. Their dispersion relation reads, for waves on the interface between two fluids of infinite depth:

where is the angular frequency, the acceleration due to gravity, the surface tension, and the mass density of the two fluids () and the is wavenumber.

### Derivation

This is a sketch of the derivation of the general dispersion relation, see Ref. 4 for a more detailed description. The problem is unfortunately a bit complex. As Richard Feynman put it (Ref. 6):

...[water waves], which are easily seen by everyone and which are used as an example of waves in elementary courses... are the worst possible example... they have all the complications that waves can have

#### Defining the problem

Three contributions to the energy are involved: the surface tension, gravity, and hydrodynamics. The parts due to surface tension (again the derivatives are taken to be small) and gravity are exactly as above. The new contribution involves the kinetic energy of the fluid:

where is the module of the velocity field . (Again, we are neglecting the flow of the gas above for simplicity.)

#### Wave solutions

Let us suppose the surface of the liquid is described by a traveling plane wave:

where and is a two dimensional wave number vector, being the horizontal position. We may take without loss of generality:

In this case it is easy to perform the integrations involved in the expressions for the energies. The integration over can taken over a period of oscillation , then multiplied by the number of oscillations in our very large (in principle, infinite) system: . The integration over trivially yields . Performing the integrations, keeping in mind that only the real part of complex numbers is to be taken as physical, one finds:

where is the area of the system.

To tackle the kinetic energy, suppose the fluid is incompressible and its flow is irrotational (often, sensible approximations) - the flow will then be potential: , and is a potential (scalar field) which must satisfy Laplace's equation . If we try try separation of variables with the potential:

with some function of time , and some function of vertical component (height) . Laplace's equation then requires on the later

This equation can be solved with the proper boundary conditions: first, must vanish well below the surface (in the "deep water" case, which is the one we consider, otherwise a more general relation holds, which is also well known in oceanography). Therefore

- ,

with some constant . The less trivial condition is the proper matching between and : the potential field must correspond to a velocity field that is adjusted to the movement of the surface:

- .

This means that

- , and

so that

- .

We may now find , which is . Performing the integration first we are left with

where we have dropped a factor of in the last step. The problem is thus specified by just a potential energy involving the square of and a kinetic energy involving the square of its time derivative: a regular harmonic oscillator. Its equation of motion will be

whose oscillatory solution is

the same dispersion as above if is neglected.

## External links

## References

- F. P. Buff, R. A. Lovett, and F. H. Stillinger, Jr. "Interfacial density profile for fluids in the critical region" Physical Review Letters
**15**pp. 621-623 (1965) - J. S. Rowlinson and B. Widom "Molecular Theory of Capillarity". Dover 2002 (originally: Oxford University Press 1982) ISBN 0486425444
- E. Chacón and P. Tarazona "Intrinsic profiles beyond the capillary wave theory: A Monte Carlo study", Physical Review Letters
**91**166103 (2003) - Samuel A. Safran "Statistical thermodynamics of surfaces, interfaces, and membranes" Addison-Wesley 1994 ISBN 9780813340791
- P. Tarazona, R. Checa, and E. Chacón "Critical Analysis of the Density Functional Theory Prediction of Enhanced Capillary Waves", Physical Review Letters
**99**196101 (2007) - R.P. Feynman, R.B. Leighton, and M. Sands "The Feynman lectures on physics" Addison-Wesley 1963. Section 51-4. ISBN 0201021153