Difference between revisions of "Dispersion Relation for a Free Surface"

From WikiWaves
Jump to navigationJump to search
Line 1: Line 1:
 
The dispersion equation for a free surface is one of the most important equations
 
The dispersion equation for a free surface is one of the most important equations
 
in linear water wave theory.  
 
in linear water wave theory.  
 
 
It arises when separating
 
It arises when separating
 
variables subject to the boundary conditions for a free surface.
 
variables subject to the boundary conditions for a free surface.
Line 7: Line 6:
 
and we present here the two-dimensional version. We denote the vertical coordinate
 
and we present here the two-dimensional version. We denote the vertical coordinate
 
by <math>z</math> which is point vertically up and the free surface is at
 
by <math>z</math> which is point vertically up and the free surface is at
<math>z=0.</math>.
+
<center><math>z=0.</math></center>.
 
The equations for the
 
The equations for the
 
[[Frequency Domain Problem]] with radial frequency <math>\,\omega</math> in terms of
 
[[Frequency Domain Problem]] with radial frequency <math>\,\omega</math> in terms of
 
the potential alone which are
 
the potential alone which are
 
+
<center>
<math>\frac{\partial \phi}{\partial z} - k_{\infty} \phi, \, z=0</math>
+
<math>\frac{\partial \phi}{\partial z} - \alpha \phi, \, z=0</math>
 
+
</center>
where <math>k_{\infty}</math> is the wavenumber in [[Infinite Depth]] which is given by  
+
where <math>\alpha</math> is the wavenumber in [[Infinite Depth]] which is given by  
<math>k_{\infty}=\omega^2/g</math> where <math>g</math> is gravity. We also have
+
<math>\alpha=\omega^2/g</math> where <math>g</math> is gravity. We also have
 
the equations within  the fluid  
 
the equations within  the fluid  
 
+
<center>
 
<math>\nabla^2\phi =0 </math>
 
<math>\nabla^2\phi =0 </math>
 
+
</center>
<math>\frac{\partial \phi}{\partial z} = 0, \, z=-H.</math>
+
<center>
 
+
<math>\frac{\partial \phi}{\partial z} = 0, \, z=-h.</math>
 +
</center>
 
We then find a separation of variables solution to [[Laplace's Equation]] and  
 
We then find a separation of variables solution to [[Laplace's Equation]] and  
apply the boundary condition at <math>z=-H</math> and we obtain the
+
apply the boundary condition at <math>z=-h</math> and we obtain the
 
following expression for the velocity potential
 
following expression for the velocity potential
 
+
<center>
<math>\phi(x,z) = e^{ikx} \cosh k(z+H) \,</math>
+
<math>\phi(x,z) = e^{kx} \cos k(z+H) \,</math>
 
+
</center>
 
If we then apply the condition at <math>z=0</math> we see that the constant <math>\,k</math>
 
If we then apply the condition at <math>z=0</math> we see that the constant <math>\,k</math>
 
(which corresponds to the wavenumber) is given by
 
(which corresponds to the wavenumber) is given by
 
+
<center>
<math>  k \sinh(kH) = k_{\infty} \cosh(kH)  </math>
+
<math>  k \sin(kH) = \alpha \cos(kh)  </math>
 
+
</center>
 
or
 
or
 
+
<center>
 
<math>  k \tanh(kH) = k_{\infty}\,\,\,(1)</math>
 
<math>  k \tanh(kH) = k_{\infty}\,\,\,(1)</math>
 
+
</center>
 
 
 
This is the dispersion equation for a free surface.
 
This is the dispersion equation for a free surface.
  
 
We can also write the separation of variables as  
 
We can also write the separation of variables as  
 
+
<center>
<math>\phi(x,z) = e^{kx} \cos k(z+H) \,</math>   
+
<math>\phi(x,z) = e^{ikx} \cosh k(z+H) \,</math>   
 
+
</center>
 
in which case the dispersion equation becomes
 
in which case the dispersion equation becomes
 
+
<center>
<math>  k \tan(kH) = -k_{\infty}\,\,\,(2)</math>
+
<math>  k \tanh(kh) = \alpha\,\,\,(2)</math>
 
+
</center>
Equation (1) has one real positive solution (plus imaginary solutions)
 
and equation (2) has an infinite number of positive real solutions (plus
 
imaginary solutions).  Sometimes (especially in older works) both equations
 
are used so that only real solutions need to be considered. This separation
 
certainly makes sense in numerical solutions but does adds unnecessarily
 
to the notation.
 
  
 
= Solution of the dispersion equation =
 
= Solution of the dispersion equation =
  
The solution of equation (1) consists of one real and infinite
+
Equation (1) has two purely imaginary solutions plus a countable number
number of imaginary roots with positive part plus their negatives.  
+
of positive and negative real solutions. Note that the equation is
These solutions multiplied by <math>i</math> are the solutions to
+
even in <math>k</math> so that for every solution the negative is also a solution.  
equation (2). We denote the solutions to (2) by <math>k_n</math>
+
The solutions of equation(2) are just <math>i</math> times the solutions
 +
to equation (1). Sometimes (especially in older works) both equations
 +
are used so that each equation needs only to be solved for real solutions.
 +
We denote the solutions to (1) by <math>k_n</math>
 
where <math>k_0</math> is the imaginary solution with
 
where <math>k_0</math> is the imaginary solution with
 
positive imaginary part and <math>k_n</math>
 
positive imaginary part and <math>k_n</math>
Line 68: Line 64:
 
The dispersion equation is a classical Sturm-Liouville equation.
 
The dispersion equation is a classical Sturm-Liouville equation.
 
The vertical eigenfunctions <math>\cos k_n (z-H)</math>
 
The vertical eigenfunctions <math>\cos k_n (z-H)</math>
form complete set for <math>L_2[-H,0]\,</math> and they are orthogonal.
+
form complete set for <math>L_2[-h,0]\,</math> and they are orthogonal.
 
Also, as <math>n\to\infty</math> <math>k_n \to \ n\pi/H</math>
 
Also, as <math>n\to\infty</math> <math>k_n \to \ n\pi/H</math>
 
so that in the limit the vertical eigenfunctions become the same as
 
so that in the limit the vertical eigenfunctions become the same as
the Fourier cosine series for <math>L_2[-H,0]\,</math> (remembering
+
the Fourier cosine series for <math>L_2[-h,0]\,</math> (remembering
 
that the eigenfunctions satisfy the boundary conditions of zero
 
that the eigenfunctions satisfy the boundary conditions of zero
normal derivative at <math>z=H</math> which is why we have the  
+
normal derivative at <math>z=h</math> which is why we have the  
 
cosine series).
 
cosine series).
  
 
[[Category:Linear Water-Wave Theory]]
 
[[Category:Linear Water-Wave Theory]]

Revision as of 05:38, 6 March 2008

The dispersion equation for a free surface is one of the most important equations in linear water wave theory. It arises when separating variables subject to the boundary conditions for a free surface. The same equation arises when separating variables in two or three dimensions and we present here the two-dimensional version. We denote the vertical coordinate by [math]\displaystyle{ z }[/math] which is point vertically up and the free surface is at

[math]\displaystyle{ z=0. }[/math]

.

The equations for the Frequency Domain Problem with radial frequency [math]\displaystyle{ \,\omega }[/math] in terms of the potential alone which are

[math]\displaystyle{ \frac{\partial \phi}{\partial z} - \alpha \phi, \, z=0 }[/math]

where [math]\displaystyle{ \alpha }[/math] is the wavenumber in Infinite Depth which is given by [math]\displaystyle{ \alpha=\omega^2/g }[/math] where [math]\displaystyle{ g }[/math] is gravity. We also have the equations within the fluid

[math]\displaystyle{ \nabla^2\phi =0 }[/math]

[math]\displaystyle{ \frac{\partial \phi}{\partial z} = 0, \, z=-h. }[/math]

We then find a separation of variables solution to Laplace's Equation and apply the boundary condition at [math]\displaystyle{ z=-h }[/math] and we obtain the following expression for the velocity potential

[math]\displaystyle{ \phi(x,z) = e^{kx} \cos k(z+H) \, }[/math]

If we then apply the condition at [math]\displaystyle{ z=0 }[/math] we see that the constant [math]\displaystyle{ \,k }[/math] (which corresponds to the wavenumber) is given by

[math]\displaystyle{ k \sin(kH) = \alpha \cos(kh) }[/math]

or

[math]\displaystyle{ k \tanh(kH) = k_{\infty}\,\,\,(1) }[/math]

This is the dispersion equation for a free surface.

We can also write the separation of variables as

[math]\displaystyle{ \phi(x,z) = e^{ikx} \cosh k(z+H) \, }[/math]

in which case the dispersion equation becomes

[math]\displaystyle{ k \tanh(kh) = \alpha\,\,\,(2) }[/math]

Solution of the dispersion equation

Equation (1) has two purely imaginary solutions plus a countable number of positive and negative real solutions. Note that the equation is even in [math]\displaystyle{ k }[/math] so that for every solution the negative is also a solution. The solutions of equation(2) are just [math]\displaystyle{ i }[/math] times the solutions to equation (1). Sometimes (especially in older works) both equations are used so that each equation needs only to be solved for real solutions. We denote the solutions to (1) by [math]\displaystyle{ k_n }[/math] where [math]\displaystyle{ k_0 }[/math] is the imaginary solution with positive imaginary part and [math]\displaystyle{ k_n }[/math] are the real solutions positive solutions ordered so that they are increasing.

The dispersion equation is a classical Sturm-Liouville equation. The vertical eigenfunctions [math]\displaystyle{ \cos k_n (z-H) }[/math] form complete set for [math]\displaystyle{ L_2[-h,0]\, }[/math] and they are orthogonal. Also, as [math]\displaystyle{ n\to\infty }[/math] [math]\displaystyle{ k_n \to \ n\pi/H }[/math] so that in the limit the vertical eigenfunctions become the same as the Fourier cosine series for [math]\displaystyle{ L_2[-h,0]\, }[/math] (remembering that the eigenfunctions satisfy the boundary conditions of zero normal derivative at [math]\displaystyle{ z=h }[/math] which is why we have the cosine series).