Difference between revisions of "Dispersion Relation for a Floating Elastic Plate"

From WikiWaves
Jump to navigationJump to search
 
(22 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
{{incomplete pages}}
 +
 +
This is an generalisation of the [[Dispersion Relation for a Free Surface]], to the case where the
 +
surface condition is given by a [[:Category:Floating Elastic Plate|Floating Elastic Plate]]
 +
 +
 
== Separation of Variables ==
 
== Separation of Variables ==
  
Line 11: Line 17:
 
[[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 is  
 
the potential alone is  
 
+
<center>
 
<math>D\frac{\partial^4}{\partial x^4} \frac{\partial \phi}{\partial z}
 
<math>D\frac{\partial^4}{\partial x^4} \frac{\partial \phi}{\partial z}
 
+ \left(\rho g
 
+ \left(\rho g
- \omega^2 \rho_i h \right) \frac{\partial \phi}{\partial z} =  
+
- \omega^2 \rho_i d \right) \frac{\partial \phi}{\partial z} =  
 
  - \rho \omega^2 \phi, \, z=0</math>
 
  - \rho \omega^2 \phi, \, z=0</math>
 
+
</center>
 
plus the equations within  the fluid  
 
plus 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>
 
where <math> \,g </math> is the acceleration due to gravity,  <math> \,\rho_i </math> and <math> \,\rho </math>  
 
where <math> \,g </math> is the acceleration due to gravity,  <math> \,\rho_i </math> and <math> \,\rho </math>  
are the densities of the plate and the water respectively, <math> \,h </math> and <math> \,D </math> are
+
are the densities of the plate and the water respectively, <math> \,d </math> and <math> \,D </math> are
 
the thickness and flexural rigidity of the plate.
 
the thickness and flexural rigidity of the plate.
  
 
We then look for a separation of variables solution to [[Laplace's Equation]] and obtain the
 
We then look for a separation of variables solution to [[Laplace's Equation]] and obtain the
 
following expression for the velocity potential
 
following expression for the velocity potential
 
+
<center><math>
<math>\phi(x,z) = e^{ikx} \cosh k(z+H) \,</math>
+
\phi(x,z) = e^{ikx} \cosh 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>
 +
-D k^5 \sinh(kh) - k \left(\rho g - \omega^2 \rho_i d \right) \sinh(kh) =
 +
-\rho \omega^2 \cosh(kh)  \,\,\,(1)
 +
</math></center>
 +
This is the dispersion equation for a floating elastic plate.
  
<math> -D k^5 \sinh(kH) - k \left(\rho g - \omega^2 \rho_i h \right) \sinh(kH) =  
+
We can also separate variables as
-\rho \omega^2 \cosh(kH) \, </math>
+
<center><math>
 +
\phi(x,z) = e^{kx} \cos k(z+h) \,
 +
</math></center>
 +
in which case the dispersion relations becomes
 +
<center><math>
 +
D k^5 \sin(kh) + k \left(\rho g - \omega^2 \rho_i h \right) \sin(kh) =  
 +
-\rho \omega^2 \cos(kh)
 +
</math></center>
  
 
== Solution of the dispersion equation ==
 
== Solution of the dispersion equation ==
Line 50: Line 68:
 
== Non-dimensional form ==
 
== Non-dimensional form ==
  
The dispersion equation is often given in non-dimensional form. The form used
+
The dispersion equation (1) is often given in non-dimensional form. The form used
 
by [[Tayler_1986a|Tayler 1986]] is to scale length with respect
 
by [[Tayler_1986a|Tayler 1986]] is to scale length with respect
 
to <math>L\,</math> and time with respect to <math>\sqrt{L/g}\,</math>. The  
 
to <math>L\,</math> and time with respect to <math>\sqrt{L/g}\,</math>. The  
 
non-dimensional equations then become
 
non-dimensional equations then become
 
+
<center>
<math> -\beta k^5 \sinh(kH) - k \left(1 - \alpha \gamma \right) \sinh(kH) =  
+
<math> -\beta k^5 \sinh(kh) - k \left(1 - \alpha \gamma \right) \sinh(kh) =  
-\alpha \cosh(kH)  \, </math>
+
-\alpha \cosh(kh)  \, </math>
 
+
</center>
where <math>H</math> is the nondimensional water depth, <math>\alpha = \omega^2,\, </math> <math>\beta = D/(\rho g L^4)\,</math>  
+
where <math>h</math> is the nondimensional water depth, <math>\alpha = L\omega^2/g,\, </math> <math>\beta = D/(\rho g L^4)\,</math>  
and <math>\gamma = \rho_i h/(\rho L)\,</math>.
+
and <math>\gamma = \rho_i d/(\rho L)\,</math>. This equation has four parameters
 
+
<math>\alpha,\beta,\gamma,h\,</math>. We can trivially remove one of these by setting the
This reduces to the following form if the length parameter is set to the
+
<math>L=h\,</math>. Alternatively we can
 +
set the length parameter to the
 
[[Characteristic Length]] <math>L = (D/(\rho g))^{1/4}\,</math> (so that in
 
[[Characteristic Length]] <math>L = (D/(\rho g))^{1/4}\,</math> (so that in
the notation above <math>\beta = 1\,</math>):
+
the notation above <math>\beta = 1\,</math> and we also have three parameters). It should
 +
be noted that in many cases, especially for ice floes, we have infinitely deep water and
 +
the parameter <math>\gamma</math> is small so that the [[Characteristic Length]]
 +
can be see as a single parameter which define the ice properties. It is then very useful
 +
to compare the value of the [[Characteristic Length]] for real ice and ice in the model
 +
basin (which much be deliberately made weaker to keep the [[Characteristic Length]] similar).
  
<math>-k^5 \sinh(kH) - k \left(1 - \alpha \gamma \right) \sinh(kH) =-\alpha \cosh(kH).</math>
+
It turns out that by choosing <math>L^5=D/(\rho\omega^2)\,</math> we obtain a dispersion
 +
equation with only two parameters and an alternative (and more straightforward)
 +
derivation of this is as follow.
 +
[[Tim Williams]] noticed that the dispersion equation can be written as
 +
<center>
 +
<math> -\frac{D}{\rho \omega^2} k^5 \sinh(kh) - k \frac{\left(\rho g - \omega^2 \rho_i d \right)}
 +
{\rho \omega^2} \sinh(kh) =  
 +
- \cosh(kh) </math></center>
 +
and if we simply non-dimensionalise with respect to a length <math>L^5=D/(\rho\omega^2)\,</math>
 +
we obtain the following dispersion equation which depends on only two parameters
 +
<center>
 +
<math> -k^5 \sinh(kh) - k \frac{\rho g - \omega^2 \rho_i d}
 +
{\rho \omega^2 L} \sinh(kh) =
 +
- \cosh(kh). </math></center>
 +
This can be written as
 +
<center>
 +
<math>\cosh(k h)-(k^4+\varpi)k\sinh(k h)=0,</math></center>
 +
where
 +
<center><math>\varpi=\frac{\rho g - \omega^2 \rho_i d}{\rho \omega^2 L}
 +
=(1-k_\infty\sigma)/(k_\infty L),\quad k_\infty=\omega^2/g,\quad\sigma=\rho_ih/\rho,\quad L^5=D/(\rho\omega^2).</math></center>
 +
<math> k_\infty </math> is the wave number for a wave traveling in open water of infinite depth, and
 +
<math>\sigma\,</math> is the amount that the plate would be submerged relative to a region of open water. Note that if <math>\sigma\,</math> is greater than <math> 1/k_\infty, </math> or <math>2\pi\,</math> times the infinite depth open water wavelength, then the parameter <math>\varpi</math> becomes negative. This is a highly non-physical situation
 +
since for the plate to float <math>\rho_i<\rho</math> and we are therefore requiring the thickness to
 +
be similar in size to the wavelength which violates our assumption of [[Negligible Submergence]]. However,
 +
this situation does have important theoretical applications, for example connecting the
 +
[[Free-Surface Green function for a Floating Elastic Plate]] with the Green function for an Elastic Plate
 +
in a vaccum.
  
The most sophisticated form, which is due to [[Tim Williams]] and has certain theoretical and practical advantages, is the following,
+
The dispersion relation relates the wavenumber <math> k/L\, </math> and thus wave speed <math> L\omega/k\, </math> to the above
  
<math>\cosh(k H)-(k^4+\varpi)k\sinh(k H)=0,</math>
+
= Matlab Code =
  
where <math>\varpi=(1-k_\infty\sigma)/(k_\infty L),\quad k_\infty=\omega^2/g,\quad\sigma=\rho_ih/\rho,\quad L^5=D/(\rho\omega^2).</math>
+
A program to calculate solutions to the dispersion relation for a free surface
 +
[http://www.math.auckland.ac.nz/~meylan/code/dispersion/dispersion_elastic_surface.m dispersion_elastic_surface.m]
  
<math> k_\infty </math> is the wave number for a wave traveling in open water of infinite depth, and
 
<math>\sigma\,</math> is the amount that the plate would be submerged relative to a region of open water. Note that if <math>\sigma\,</math> is greater than <math> 1/k_\infty, </math> or <math>2\pi\,</math> times the infinite depth open water wavelength, then the parameter <math>\varpi</math> becomes negative.
 
  
The dispersion relation relates the wavenumber <math> k/L\, </math> and thus wave speed to the above parameters.
+
[[Category:Floating Elastic Plate]]

Latest revision as of 09:36, 20 October 2009


This is an generalisation of the Dispersion Relation for a Free Surface, to the case where the surface condition is given by a Floating Elastic Plate


Separation of Variables

The dispersion equation arises when separating variables subject to the boundary conditions for a Floating Elastic Plate of infinite extent. The same equation arises when separating variables in two or three dimensions and we present here the two-dimensional version. The equations are described in detail in the Floating Elastic Plate page and we begin with the equations The equations of motion for the Frequency Domain Problem with radial frequency [math]\displaystyle{ \,\omega }[/math] in terms of the potential alone is

[math]\displaystyle{ D\frac{\partial^4}{\partial x^4} \frac{\partial \phi}{\partial z} + \left(\rho g - \omega^2 \rho_i d \right) \frac{\partial \phi}{\partial z} = - \rho \omega^2 \phi, \, z=0 }[/math]

plus the equations within the fluid

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

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

where [math]\displaystyle{ \,g }[/math] is the acceleration due to gravity, [math]\displaystyle{ \,\rho_i }[/math] and [math]\displaystyle{ \,\rho }[/math] are the densities of the plate and the water respectively, [math]\displaystyle{ \,d }[/math] and [math]\displaystyle{ \,D }[/math] are the thickness and flexural rigidity of the plate.

We then look for a separation of variables solution to Laplace's Equation and obtain the following expression for the velocity potential

[math]\displaystyle{ \phi(x,z) = e^{ikx} \cosh 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{ -D k^5 \sinh(kh) - k \left(\rho g - \omega^2 \rho_i d \right) \sinh(kh) = -\rho \omega^2 \cosh(kh) \,\,\,(1) }[/math]

This is the dispersion equation for a floating elastic plate.

We can also separate variables as

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

in which case the dispersion relations becomes

[math]\displaystyle{ D k^5 \sin(kh) + k \left(\rho g - \omega^2 \rho_i h \right) \sin(kh) = -\rho \omega^2 \cos(kh) }[/math]

Solution of the dispersion equation

The dispersion equation was first solved by Fox and Squire 1994 and they determined that the solution consists of one real, two complex, and infinite number of imaginary roots plus their negatives. Interestingly the eigenfunctions form an over complete set for [math]\displaystyle{ L_2[-H,0]\, }[/math]. Also, there are some circumstances (non-physical) in which the complex roots become purely imaginary. The solution of this dispersion equation is far from trivial and the optimal solution method has been developed by Tim Williams and is described below.

Non-dimensional form

The dispersion equation (1) is often given in non-dimensional form. The form used by Tayler 1986 is to scale length with respect to [math]\displaystyle{ L\, }[/math] and time with respect to [math]\displaystyle{ \sqrt{L/g}\, }[/math]. The non-dimensional equations then become

[math]\displaystyle{ -\beta k^5 \sinh(kh) - k \left(1 - \alpha \gamma \right) \sinh(kh) = -\alpha \cosh(kh) \, }[/math]

where [math]\displaystyle{ h }[/math] is the nondimensional water depth, [math]\displaystyle{ \alpha = L\omega^2/g,\, }[/math] [math]\displaystyle{ \beta = D/(\rho g L^4)\, }[/math] and [math]\displaystyle{ \gamma = \rho_i d/(\rho L)\, }[/math]. This equation has four parameters [math]\displaystyle{ \alpha,\beta,\gamma,h\, }[/math]. We can trivially remove one of these by setting the [math]\displaystyle{ L=h\, }[/math]. Alternatively we can set the length parameter to the Characteristic Length [math]\displaystyle{ L = (D/(\rho g))^{1/4}\, }[/math] (so that in the notation above [math]\displaystyle{ \beta = 1\, }[/math] and we also have three parameters). It should be noted that in many cases, especially for ice floes, we have infinitely deep water and the parameter [math]\displaystyle{ \gamma }[/math] is small so that the Characteristic Length can be see as a single parameter which define the ice properties. It is then very useful to compare the value of the Characteristic Length for real ice and ice in the model basin (which much be deliberately made weaker to keep the Characteristic Length similar).

It turns out that by choosing [math]\displaystyle{ L^5=D/(\rho\omega^2)\, }[/math] we obtain a dispersion equation with only two parameters and an alternative (and more straightforward) derivation of this is as follow. Tim Williams noticed that the dispersion equation can be written as

[math]\displaystyle{ -\frac{D}{\rho \omega^2} k^5 \sinh(kh) - k \frac{\left(\rho g - \omega^2 \rho_i d \right)} {\rho \omega^2} \sinh(kh) = - \cosh(kh) }[/math]

and if we simply non-dimensionalise with respect to a length [math]\displaystyle{ L^5=D/(\rho\omega^2)\, }[/math] we obtain the following dispersion equation which depends on only two parameters

[math]\displaystyle{ -k^5 \sinh(kh) - k \frac{\rho g - \omega^2 \rho_i d} {\rho \omega^2 L} \sinh(kh) = - \cosh(kh). }[/math]

This can be written as

[math]\displaystyle{ \cosh(k h)-(k^4+\varpi)k\sinh(k h)=0, }[/math]

where

[math]\displaystyle{ \varpi=\frac{\rho g - \omega^2 \rho_i d}{\rho \omega^2 L} =(1-k_\infty\sigma)/(k_\infty L),\quad k_\infty=\omega^2/g,\quad\sigma=\rho_ih/\rho,\quad L^5=D/(\rho\omega^2). }[/math]

[math]\displaystyle{ k_\infty }[/math] is the wave number for a wave traveling in open water of infinite depth, and [math]\displaystyle{ \sigma\, }[/math] is the amount that the plate would be submerged relative to a region of open water. Note that if [math]\displaystyle{ \sigma\, }[/math] is greater than [math]\displaystyle{ 1/k_\infty, }[/math] or [math]\displaystyle{ 2\pi\, }[/math] times the infinite depth open water wavelength, then the parameter [math]\displaystyle{ \varpi }[/math] becomes negative. This is a highly non-physical situation since for the plate to float [math]\displaystyle{ \rho_i\lt \rho }[/math] and we are therefore requiring the thickness to be similar in size to the wavelength which violates our assumption of Negligible Submergence. However, this situation does have important theoretical applications, for example connecting the Free-Surface Green function for a Floating Elastic Plate with the Green function for an Elastic Plate in a vaccum.

The dispersion relation relates the wavenumber [math]\displaystyle{ k/L\, }[/math] and thus wave speed [math]\displaystyle{ L\omega/k\, }[/math] to the above

Matlab Code

A program to calculate solutions to the dispersion relation for a free surface dispersion_elastic_surface.m