Difference between revisions of "Free-Surface Green Function for a Floating Elastic Plate"

From WikiWaves
Jump to navigationJump to search
 
(46 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
{{incomplete pages}}
 +
 
This is a special version of the free-surface Green function which applied when the [[Floating Elastic Plate]]
 
This is a special version of the free-surface Green function which applied when the [[Floating Elastic Plate]]
boundary condition applies at the free-surface
+
boundary condition applies at the free-surface. It reduced to the [[Free-Surface Green Function]] in the limit as the plate terms tend to zero and to the Green function for an infinite plate in the limit as the water terms tend to zero.
  
 
= Two Dimensions =
 
= Two Dimensions =
  
The Green function <math>G(x,z)</math> representing  
+
== Green function with singularity on the surface ==
outgoing waves as <math>|x|\rightarrow \infty</math> satisfies
+
For the case of a floating thin plate we almost always make the assumption of shallow draft and the Green function satisfies the following equations
 +
The Green function <math>G(x,x^{\prime},z)</math> representing outgoing waves as <math>|x|\rightarrow \infty</math> satisfies
 
<center><math>
 
<center><math>
\nabla^2 G = 0, -h<z<0,  
+
\nabla^2 G = 0, -H<z<0,  
 
</math></center>
 
</math></center>
 
<center><math>
 
<center><math>
\frac{\partial G}{\partial z} =0, z=-h,
+
\frac{\partial G}{\partial z} =0, z=-H,
 
</math></center>
 
</math></center>
 
<center><math>
 
<center><math>
 
{\left( \beta \frac{\partial^4}{\partial x^4} -  
 
{\left( \beta \frac{\partial^4}{\partial x^4} -  
 
\gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G}
 
\gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G}
  = \delta(x), z=0,
+
  = \delta(x-x^{\prime}), z=0,
 
</math></center>
 
</math></center>
(note we are only considering singularities at the free surface).  
+
(note we are only considering singularities at the free surface). It is important to realise that these equations are based on the [[Frequency Domain Problem]] and that the exact form of the Green function is dependent on whether we have <math>e^{i\omega t}</math> or <math>e^{-i\omega t}</math> dependence. In what follows we assume <math>e^{i\omega t}</math> dependence.  
  
This problem can be solved (a solution is given below) to give
+
This problem can be solved (a solution is given below and also in [[Evans and Porter 2006]]) to give
 
<center><math>
 
<center><math>
G(x,z) = -i\sum_{n=-2}^\infty\frac{\sin{(k_n h)}\cos{(k_n(z+h))}}{2\alpha C(k_n)}e^{-k_n|x|},
+
G(x,x^{\prime},z) = -\sum_{n=-2}^\infty\frac{\sin{(k_n H)}\cos{(k_n(z+H))}}{2\alpha C(k_n)}e^{-k_n|x-x^{\prime}|},
 
</math></center>
 
</math></center>
 
where
 
where
 
<center><math>
 
<center><math>
C(k_n)=\frac{1}{2}\left(h + \frac{(5\beta k(n)^4 + 1 - \alpha\gamma)\sin^2{(k(n)h)}}{\alpha}\right),
+
C(k_n)=\frac{1}{2}\left(h - \frac{(5\beta k_n^4 + 1 - \alpha\gamma)\sin^2{(k_n H)}}{\alpha}\right),
 
</math></center>
 
</math></center>
 
and <math>k_n</math> are the solutions of the [[Dispersion Relation for a Floating Elastic Plate]],
 
and <math>k_n</math> are the solutions of the [[Dispersion Relation for a Floating Elastic Plate]],
Line 33: Line 36:
 
</math></center>
 
</math></center>
 
with <math>n=-1,-2</math> corresponding to the complex solutions with positive real part,  
 
with <math>n=-1,-2</math> corresponding to the complex solutions with positive real part,  
<math>n=0</math> corresponding to the imaginary solution with positive imaginary part and
+
<math>n=0</math> corresponding to the imaginary solution with positive imaginary part, since we have assumed <math>e^{i\omega t}</math> dependence (it would be of negative imaginary part if we had chosen <math>e^{-i\omega t}</math> dependence)
 
<math>n>0</math> corresponding to the real solutions with positive real part.
 
<math>n>0</math> corresponding to the real solutions with positive real part.
 
We can also write the Green function as
 
We can also write the Green function as
Line 41: Line 44:
 
\cosh q_n\left(  z+H\right)  .
 
\cosh q_n\left(  z+H\right)  .
 
</math></center>
 
</math></center>
where <math>r=|x|</math>, <math>q_n=ik_n</math> and we have used a different non-dimensionalisation so that
+
where <math>r=|x-x^{\prime}|</math>, <math>q_n=ik_n</math> and we have used a different non-dimensionalisation so that
 
the [[Dispersion Relation for a Floating Elastic Plate]] is written as
 
the [[Dispersion Relation for a Floating Elastic Plate]] is written as
 
<center><math>  
 
<center><math>  
Line 51: Line 54:
 
R\left(  q_n\right)  =\frac{\omega^{2}q_n}{\omega^{2}\left(  5q_n^{4}+u\right)
 
R\left(  q_n\right)  =\frac{\omega^{2}q_n}{\omega^{2}\left(  5q_n^{4}+u\right)
 
+H\left[  \left(  q_n^{5}+uq\right)  ^{2}-\omega^{4}\right]  }.  
 
+H\left[  \left(  q_n^{5}+uq\right)  ^{2}-\omega^{4}\right]  }.  
 +
</math></center>
 +
 +
== Green function with singularity in the water ==
 +
 +
We can define the Green function <math>G(x,x^{\prime},z)</math>, representing
 +
outgoing waves as <math>|x|\rightarrow \infty</math>, which satisfies
 +
<center><math>
 +
\nabla^2 G = \delta(x-x^{\prime})\delta(z-z^{\prime}), -H<z<0,
 +
</math></center>
 +
<center><math>
 +
\frac{\partial G}{\partial z} =0, z=-H,
 +
</math></center>
 +
<center><math>
 +
{\left( \beta \frac{\partial^4}{\partial x^4} -
 +
\gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G}
 +
= 0, z=0.
 +
</math></center>
 +
Note that this Green function reduces the [[Free-Surface Green Function]] in the limit as <math>\beta</math> and <math>\gamma</math> tend to zero.
 +
 +
This problem can be solved to give [[Evans and Porter 2003]]
 +
<center><math>
 +
G(x,x^{\prime},z) = -\sum_{n=-2}^\infty\frac{\cos{(k_n(z^{\prime}+H))}\cos{(k_n(z+H))}}{2k_n C(k_n)}e^{-k_n|x-x^{\prime}|},
 +
</math></center>
 +
where <math>
 +
C(k_n)</math>
 +
and <math>k_n</math> are as before.
 +
 +
== Relationship between the two Green functions ==
 +
 +
If we denote the Green function with the singularity at the free-surface by <math>G_s</math> and the Green function with
 +
the singularity in the water by <math>G_w</math>then we have the following relationship between them
 +
<center><math>
 +
G_s = -\frac{1}{\alpha}\left. \frac{\partial G_w}{\partial z^{\prime}} \right|_{z^{\prime} =0}
 
</math></center>
 
</math></center>
  
Line 56: Line 92:
  
 
<center><math>
 
<center><math>
\phi_{=P= }\left(  r,z\right)  =-\frac{\omega}{2}\sum_{q\in K^{
+
G\left(  r,z\right)  =-\frac{\omega}{2}\sum_{n=-2}^{\infty}
}}\frac{R\left(  q\right)  }{\sinh qH}H_{0}^{\left(  1\right)  }\left(
+
\frac{R\left(  q_n\right)  }{\sinh q_nH}H_{0}^{\left(  1\right)  }\left(
qr\right)  \cosh q\left(  z+H\right)
+
q_n r\right)  \cosh q_n\left(  z+H\right)
 
</math></center>
 
</math></center>
  
Line 69: Line 105:
 
==Previous Work ==
 
==Previous Work ==
  
Radially symmetric waves in a thin plate without
+
In the following sections, the
the hydrodynamic support, have studied by [[Sneddon 1945]] using the two
 
dimensional radial Fourier transform. This integral transform is commonly
 
called the Hankel transform
 
<center><math>
 
\hat{w}\left(  \gamma\right)  =\frac{1}{2\pi}\int_{0}^{\infty}w\left(
 
r\right)  rJ_{0}\left(  \gamma r\right)  dr.
 
</math></center>
 
A complete analytical solution for deflection of an infinite plate due to a
 
localized static load was derived by [[Wyman 1950]] by choosing the
 
bounded solutions of
 
<center><math>
 
w_{rr}+\frac{1}{r}w_{r}+\mathrm{i}w=0 (3-81)
 
</math></center>
 
whose solutions satisfy the radial form of the (non-dimensional) thin plate
 
equation
 
<center><math>
 
\left[  \frac{\partial^{2}}{\partial r^{2}}+\frac{1}{r}\frac{\partial
 
}{\partial r}\right]  ^{2}w+w=0.
 
</math></center>
 
 
 
If we consider a plate floating on water we have to
 
to solve Laplace's equation, and thus the application of the Fourier
 
transform is not as straightforward as in the two examples cited (Sneddon
 
[[Sneddon45]] and Wyman [[Wyman50]]). Response of an infinite thin
 
elastic plate under fluid loading is studied in a series of papers by Crighton
 
[[Crighton72, Crighton77, Crighton79]]. In the following sections, the
 
 
deflection of a plate is expressed by the infinite number of modes that exist
 
deflection of a plate is expressed by the infinite number of modes that exist
 
in the vibrating plate floating on the water, instead of a finite number of
 
in the vibrating plate floating on the water, instead of a finite number of
 
modes that represent the free oscillation of an elastic plate without
 
modes that represent the free oscillation of an elastic plate without
 
restraints. The fact that waves in an incompressible fluid can be expressed by
 
restraints. The fact that waves in an incompressible fluid can be expressed by
an infinite series of natural modes is shown by John [[fritz, fritz2]],
+
an infinite series of natural modes is shown by [[John 1949]] and [[John 1950]]
each mode being a Hankel function of the first kind. It is known that the
+
and discussed in [[Cylindrical Eigenfunction Expansion]]
oscillation of an infinite plate can also be expressed in terms of Hankel
+
Each mode being a Hankel function of the first kind in the horizontal
functions. Kheisin, in 1967 ([[Kheisin67]] chapter IV), studied the same
+
(in the three-dimensional case) and a hyperbolic function
problem solved in this chapter. Kheisin derived the same dispersion equation
+
in the vertical. [[Kheisin 1967]] chapter IV, studied the same
for the physical variables, and then studied the properties of the inverse
+
problem solved in this chapter and derived the same dispersion equation
transform for the simple shallow water and static load cases, which are
+
for the physical variables. He studied the properties of the inverse
approximated versions of the solutions reported here.
+
transform only for the simple shallow water and static load cases, which are
 +
approximated versions of the full solutions which is presented here.
  
 
==Mathematical model==
 
==Mathematical model==
Line 114: Line 125:
 
<center><math>
 
<center><math>
 
\left[  \nabla^{4}-m\omega^{2}+1\right]  w+\mathrm{i}\omega
 
\left[  \nabla^{4}-m\omega^{2}+1\right]  w+\mathrm{i}\omega
\phi=p_{=a= }.  
+
\phi=p_n.  
 
</math></center>
 
</math></center>
plus Laplace's equation for the velocity potential and the bottom condition
+
where <math>p_n</math> is the pressure and <math>n</math> is the dimension. We assume
 +
that <math>p_2 = \delta(x) </math> and that <math>p_3 = \delta(x) \delta(y)</math>.
 +
<math>w</math> is the displacement of the surface which is related to the velocity
 +
potential by the equation
 +
<math>\frac{\partial\phi}{\partial z} = \frac{\partial w}{\partial t}</math>
 +
We also have Laplace's equation for the velocity potential and the bottom condition
 
<center><math>
 
<center><math>
 
\nabla^{2}\phi=0
 
\nabla^{2}\phi=0
Line 131: Line 147:
  
 
We solve the system of equations using the [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform] in <math>\left(  x,y\right)  </math>-plane for three-dimensions and with
 
We solve the system of equations using the [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform] in <math>\left(  x,y\right)  </math>-plane for three-dimensions and with
a Fourier Transform in
+
a [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform] in
 
<math>x</math> for two dimensions. We choose the Fourier transform with respect to <math>x</math>
 
<math>x</math> for two dimensions. We choose the Fourier transform with respect to <math>x</math>
 
and <math>y</math> defined as
 
and <math>y</math> defined as
Line 137: Line 153:
 
\hat{\phi}\left(  \alpha,k,z\right)  =\int_{-\infty}^{\infty}\int_{-\infty
 
\hat{\phi}\left(  \alpha,k,z\right)  =\int_{-\infty}^{\infty}\int_{-\infty
 
}^{\infty}\phi\left(  x,y,z\right)  e^{\mathrm{i}\left(  \alpha
 
}^{\infty}\phi\left(  x,y,z\right)  e^{\mathrm{i}\left(  \alpha
x+ky\right)  }dxdy
+
x+ky\right)  }\mathrm{d}x\mathrm{d}y
 
</math></center>
 
</math></center>
and the inverse Fourier transform defined as
+
and the inverse [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform] defined as
 
<center><math>
 
<center><math>
 
\phi\left(  x,y,z\right)  =\frac{1}{4\pi^{2}}\int_{-\infty}^{\infty}
 
\phi\left(  x,y,z\right)  =\frac{1}{4\pi^{2}}\int_{-\infty}^{\infty}
 
\int_{-\infty}^{\infty}\hat{\phi}\left(  \alpha,k,z\right)
 
\int_{-\infty}^{\infty}\hat{\phi}\left(  \alpha,k,z\right)
e^{-\mathrm{i}\left(  \alpha x+ky\right)  }d\alpha dk.
+
e^{-\mathrm{i}\left(  \alpha x+ky\right)  }\mathrm{d}\alpha \mathrm{d}k.
 
</math></center>
 
</math></center>
 
For the one-dimensional case, the definitions are
 
For the one-dimensional case, the definitions are
 
<center><math>
 
<center><math>
 
\hat{\phi}\left(  \alpha,z\right)    =\int_{-\infty}^{\infty}\phi\left(
 
\hat{\phi}\left(  \alpha,z\right)    =\int_{-\infty}^{\infty}\phi\left(
x,z\right)  e^{\mathrm{i}\alpha x}dx,  
+
x,z\right)  e^{\mathrm{i}\alpha x}\mathrm{d}x,  
 
</math></center>
 
</math></center>
 
<center><math>
 
<center><math>
 
\phi\left(  x,z\right)    =\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat{\phi
 
\phi\left(  x,z\right)    =\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat{\phi
}\left(  \alpha,z\right)  e^{-\mathrm{i}\alpha x}d\alpha.  
+
}\left(  \alpha,z\right)  e^{-\mathrm{i}\alpha x}\mathrm{d}\alpha.  
 
</math></center>
 
</math></center>
We denote the spatial Fourier transform by using a hat over <math>w</math> and <math>\phi</math>.
+
We denote the spatial [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform]
 +
by using a hat over <math>w</math> and <math>\phi</math>.
  
The Fourier transform of both sides of the Laplace's equation  
+
The [http://en.wikipedia.org/wiki/Fourier_Transform Fourier Transform] of both sides of the Laplace's equation  
 
becomes an ordinary differential equation with respect to <math>z</math>,
 
becomes an ordinary differential equation with respect to <math>z</math>,
 
<center><math>
 
<center><math>
Line 218: Line 235:
 
}{\gamma\tanh\gamma H}.  
 
}{\gamma\tanh\gamma H}.  
 
</math></center>
 
</math></center>
Where <math>d\left(  \gamma,\omega\right)  </math>\ in the  
+
Where <math>d\left(  \gamma,\omega\right)  </math> in the  
 
[[Dispersion Relation for a Floating Elastic Plate]] function, and the associated algebraic equation <math>d\left(  \gamma
 
[[Dispersion Relation for a Floating Elastic Plate]] function, and the associated algebraic equation <math>d\left(  \gamma
 
,\omega\right)  =0</math> the dispersion equation for waves propagating through an
 
,\omega\right)  =0</math> the dispersion equation for waves propagating through an
 
ice sheet. This dispersion relation (and the Fourier
 
ice sheet. This dispersion relation (and the Fourier
transform was derived by  
+
transform of it) was derived by  
 
[[Kheisin 1967]] chapter IV, and [[Fox and Squire 1994]].
 
[[Kheisin 1967]] chapter IV, and [[Fox and Squire 1994]].
  
Line 228: Line 245:
  
 
Our task now is to derive the inverse Fourier transform). We notice that the roots of the  
 
Our task now is to derive the inverse Fourier transform). We notice that the roots of the  
[[Dispersion Relation for a Floating Elastic Plate]] for a fixed <math>\omega</math>, are the poles of the function <math>\hat{w}\left(  \gamma\right)  </math>, which are necessary for calculation of integrals involved in the inverse Fourier transform. There are a pair of two real roots <math>\pm q_0</math>, an infinite number of pure imaginary roots <math>\left\{ \pm q_{n}\right\} {n=1,2,...}</math>, (<math>\Im q_{n}>0</math>)plus four complex roots  
+
[[Dispersion Relation for a Floating Elastic Plate]] for a fixed <math>\omega</math>, are the poles of the function <math>\hat{w}\left(  \gamma\right)  </math>, which are necessary for calculation of integrals involved in the inverse Fourier transform. There are a pair of two real roots <math>\pm q_0</math>, an infinite number of pure imaginary roots  
 +
<math>  \pm q_{n},\,n\in\mathbb{N}</math>, (<math>\mathrm{Im} q_{n}>0</math>) plus four complex roots  
 
(in all physical situations) occur at plus and minus complex-conjugate pairs <math>\pm q_{-1}</math> and  
 
(in all physical situations) occur at plus and minus complex-conjugate pairs <math>\pm q_{-1}</math> and  
<math>\pm q_{-2}^{}</math> with <math>Im \left(  q_{-1,-2}\right)  >0</math>).
+
<math>\pm q_{-2}</math> with <math>\mathrm{Im} \left(  q_{-1,-2}\right)  >0</math>).
  
 
We note that the dispersion equation represents a relationship between the spatial wavenumbers <math>\gamma</math> and the radial frequency <math>\omega</math>, which is how the name "dispersion" came about.
 
We note that the dispersion equation represents a relationship between the spatial wavenumbers <math>\gamma</math> and the radial frequency <math>\omega</math>, which is how the name "dispersion" came about.
Line 241: Line 259:
 
which is also a function of <math>\gamma</math> only and has the same poles as <math>\hat{w}\left(  \gamma\right)  </math> does.
 
which is also a function of <math>\gamma</math> only and has the same poles as <math>\hat{w}\left(  \gamma\right)  </math> does.
  
The Fourier transform of the differential equations can be interpreted as wave-like forcing of the ice sheet, i.e.,
+
We calculate the displacement <math>w\left(  x,y\right)  </math> using the
<center><math>
+
inverse Fourier transform of <math>\hat{w}\left(  \gamma\right)  </math>. Since <math>\hat{w}\left(  \gamma\right)  </math> is radially symmetric, the inverse transform may be written ([[Bracewell 1965]])
p_{=a= }\left(  x,y,t\right)  =\hat{p}_{=a= }\left(  \alpha
 
,k,\omega\right)  e^{\mathrm{i}\left(  \omega t-\alpha x-ky\right)  }
 
</math></center>
 
where <math>\hat{p}_{=a= }\left(  \alpha,k,\omega\right)  </math> is the amplitude of the wave-like forcing or a Fourier transform of the wave-like force which we set to be one. Therefore, superposition or the inverse Fourier transform of the solution under this wave-like force represents the response of a plate to a localized force expressed by the delta function in Eqn.~((3-23)). Since the waves are radially symmetric, <math>p_{=a= }</math> can be written as a function of <math>r</math>,
 
<center><math>
 
p_{=a= }\left(  r,t\right)  =\hat{p}_{=a= }\left(  \gamma
 
,\omega\right)  J_{0}\left(  \gamma r\right)  e^{\mathrm{i}\omega t}
 
</math></center>
 
where <math>\hat{p}_{=a= }\left(  \gamma,\omega\right)  </math> is the Fourier component of radially symmetric wave-like function.
 
 
 
====The inverse Fourier transform (sec:3)====
 
 
 
We calculate the displacement <math>w\left(  x,y\right)  </math> by performing the two dimensional inverse Fourier transform of <math>\hat{w}\left(  \gamma\right)  </math> in Eqn.~((3-8)) and, since <math>\hat{w}\left(  \gamma\right)  </math> is radially symmetric, the inverse transform may be written (Bracewell [[Bracewell65]])
 
 
<center><math>
 
<center><math>
w_{=P= }\left(  r\right)  =\frac{1}{2\pi}\int_{0}^{\infty}\hat{w}\left(
+
w\left(  r\right)  =\frac{1}{2\pi}\int_{0}^{\infty}\hat{w}\left(
\gamma\right)  \gamma J_{0}\left(  \gamma r\right)  d\gamma (3-11)
+
\gamma\right)  \gamma J_{0}\left(  \gamma r\right)  \mathrm{d}\gamma  
 
</math></center>
 
</math></center>
 +
for the three dimensional problem,
 
where <math>J_{0}</math> is Bessel function and <math>r</math> is the distance from the point of
 
where <math>J_{0}</math> is Bessel function and <math>r</math> is the distance from the point of
forcing. The response to line forcing is given by the inverse Fourier
+
forcing. For the two-dimensional problem the inverse Fourier
 
transform of <math>\hat{w}\left(  \gamma\right)  </math> in the <math>x</math>-axis and since
 
transform of <math>\hat{w}\left(  \gamma\right)  </math> in the <math>x</math>-axis and since
 
<math>\hat{w}\left(  \gamma\right)  </math> is an even function, this is
 
<math>\hat{w}\left(  \gamma\right)  </math> is an even function, this is
 
<center><math>
 
<center><math>
w_{=L= }\left(  r\right)  =\frac{1}{\pi}\int_{0}^{\infty}\hat{w}\left(
+
w\left(  r\right)  =\frac{1}{\pi}\int_{0}^{\infty}\hat{w}\left(
\gamma\right)  \cos\left(  \gamma r\right)  d\gamma (3-12)
+
\gamma\right)  \cos\left(  \gamma r\right)  \mathrm{d}\gamma  
 
</math></center>
 
</math></center>
where again <math>r=\left|  x\right|  </math> is the distance from the line of forcing. Note that the factors <math>1/2\pi</math> and <math>1/\pi</math> result from the form of the Fourier transform in two and one dimensional spaces that we defined in Eqn.~((3-22)) and Eqn.~((3-36)).
+
where again <math>r=\left|  x\right|  </math>. Note that the factors <math>1/2\pi</math> and <math>1/\pi</math> result from the form of the Fourier transform in two and one dimensional spaces.
  
The integrals of Eqn.~((3-11)) and Eqn.~((3-12)) are calculated using the singularities of <math>\hat{w}</math>, i.e., the roots of the dispersion equation. First, since <math>\hat{w}\left(  \gamma\right)  </math> is an even fractional function that equals zero when <math>\gamma=0</math> and is bounded in the whole plane except in regions around its poles, we find that <math>\hat{w}\left(  \gamma\right)  </math> can be expressed using the Mittag-Leffler expansion in section 3.2. Pairing each pole <math>q</math> with its negative counterpart <math>-q</math>, gives
+
The integrals are calculated using the [[Mittag-Leffler Expansion for the Floating Elastic Plate Dispersion Relation]]
 +
where we show that
 
<center><math>
 
<center><math>
\hat{w}\left(  \gamma\right)  =\sum_{q\in K^{}}\frac{2qR\left(
+
\hat{w}\left(  \gamma\right)  =\sum_{n=-2}^{\infty}\frac{2q_nR\left(
q\right)  }{\gamma^{2}-q^{2}} (3-50)
+
q_n\right)  }{\gamma^{2}-q_n^{2}}  
 
</math></center>
 
</math></center>
where <math>R\left(  q\right)  </math> is the residue of <math>\hat{w}\left(  \gamma\right)  </math> at <math>\gamma=q</math>. We denoted the set of roots of the dispersion equation with positive imaginary part and a positive real root by <math>K^{}</math>. Thus, this set is <math>K^{}=\left\{  q_{=T= },q_{=D= },-q_{=D=  }^{},\mathrm{i}q_{1},\mathrm{i}q_{2},\mathrm{i} q_{3},\cdots\right\}  </math>. Note that the rest of the roots of the dispersion equation are the negative of the values in <math>K^{}</math>. By substituting this expansion into the integrals in Eqn.~((3-11)) and Eqn.~((3-12)), we are able to perform the integration and write each result as a summation over the roots <math>q\in K^{}</math>.
+
where <math>R\left(  q_n\right)  </math> is the residue of <math>\hat{w}\left(  \gamma\right)  </math> at <math>\gamma=q_n</math>
 
+
given by
\begin{figure}[tbh]
 
\begin{center}
 
\includegraphics[height=5.0918cm,width=7.5037cm]{polyres.eps}
 
\caption{Graphs of <math>\left|  R\left(  \mathrm{i}q\right)  \right|  </math> as function of real number <math>q</math>, with <math>\tanh</math>-functions (solid line) and polynomial expression (dashed line).}
 
(fig:3-31)
 
\end{center}
 
\end{figure}
 
 
 
The residues <math>R\left(  q\right)  </math> can be calculated using the usual formula.
 
Since each of the poles of <math>\hat{w}\left(  \gamma\right)  </math> is simple, the
 
residue <math>R\left(  q\right)  </math> at a pole <math>q</math> can be found using the expression
 
<center><math>\begin{matrix}
 
R\left(  q\right)  & =\left[  \left.  \frac{d}{d\gamma}d\left(  \gamma
 
,\omega\right)  \right|  _{\gamma=q}\right]  ^{-1}\\
 
& =\left[  4q^{3}+\omega^{2}\left(  \frac{qH+\tanh qH-qH\tanh^{2}qH}
 
{q^{2}\tanh^{2}qH}\right)  \right]  ^{-1}. (3-16)
 
\end{matrix}</math></center>
 
As each pole <math>q</math> is a root of the dispersion equation, we may substitute
 
<math>\tanh qH=\omega^{2}/\left(  q^{5}+uq\right)  </math>, where for brevity we have
 
defined <math>u=\left(  1-m\omega^{2}\right)  </math>. The residue may then be given as
 
the rational function of the pole
 
 
<center><math>
 
<center><math>
R\left(  q\right)  =\frac{\omega^{2}q}{\omega^{2}\left(  5q^{4}+u\right)
+
R\left(  q_n\right)  =\frac{\omega^{2}q_n}{\omega^{2}\left(  5q_n^{4}+u\right)
+H\left[  \left(  q^{5}+uq\right)  ^{2}-\omega^{4}\right]  }. (3-17)
+
+H\left[  \left(  q_n^{5}+uq_n\right)  ^{2}-\omega^{4}\right]  }.
</math></center>
+
</math></center>.
This form avoids calculation of the hyperbolic tangent which becomes small at the imaginary roots. This causes numerical round-off problems and rapid variation in computing Eqn.~((3-16)) since <math>q_{n}H</math> tends to <math>n\pi</math> as <math>n</math> becomes large, which makes <math>\tan q_{n}H</math> become small. Fig.~((fig:3-31)) shows the graphs of the two expressions, Eqn.~((3-16)) and Eqn.~((3-17)) of the residue for imaginary argument. The imaginary roots <math>\mathrm{i}q_{1},\mathrm{i}q_{2},\mathrm{i}q_{3},...</math> are located where the two curves in Fig.~((fig:3-31)) coincide (the spiky parts of the solid curve). There are two such points at a spike and the root is the one on the left. Eqn.~((3-16)), from the direct calculation, is a rapidly varying function near the roots <math>\left\{  \mathrm{i} q_{n}\right\}  _{n=1}^{\infty}</math>, hence a small numerical error in the values of the roots will result in a large error in the residue. In contrast, Eqn.~((3-17)) gives us a smooth function, and the resulting calculation of the residue is stable.
 
  
Using the identities (Abramowitz and Stegun [[Stegun]] formula 11.4.44 with <math>\nu=0</math>, <math>\mu=0</math>, <math>z=-\mathrm{i}q</math> and <math>a=r</math>)
+
Using the identities [[Abramowitz and Stegun 1964]] formula 11.4.44 with <math>\nu=0</math>, <math>\mu=0</math>, <math>z=-\mathrm{i}q</math> and <math>a=r</math>)
 
<center><math>
 
<center><math>
 
\int_{0}^{\infty}\frac{\gamma}{\gamma^{2}-q^{2}}J_{0}\left(  \gamma r\right)
 
\int_{0}^{\infty}\frac{\gamma}{\gamma^{2}-q^{2}}J_{0}\left(  \gamma r\right)
d\gamma=K_{0}\left(  -\mathrm{i}qr\right)  (3-24)
+
\mathrm{d}\gamma=K_{0}\left(  -\mathrm{i}qr\right)   
</math></center>
 
for <math>Im q>0</math>, <math>r>0</math>, where <math>K_{0}</math> is a modified Bessel function and an identity between the modified Bessel function, we notice that <math>q_{=T= }</math> term of Eqn.~((3-50)) may pose a problem. However, considering that <math>q_{=T= }=\lim_{\varepsilon\searrow0}\left(  q_{=T=  }+\mathrm{i}\varepsilon\right)  </math>, we are able to apply Eqn.~((3-24)) to all terms of Eqn.~((3-50)). Alternatively, we may put an additional imaginary term <math>\mathrm{i}\beta\omega</math>, <math>\beta>0</math>, representing damping, to the dispersion equation
 
<center><math>
 
\gamma^{4}+\mathrm{i}\beta\omega+1-m\omega^{2}-\frac{\omega^{2}
 
}{\gamma\tanh\gamma H}.
 
 
</math></center>
 
</math></center>
Every element in <math>K^{}</math> has positive imaginary part. Addition of the damping term proportional to <math>w_{t}</math> ensures that the inverse transform satisfies the radiation condition. We omit the damping term here for the sake of algebraic simplicity.
+
for <math>Im q>0</math>, <math>r>0</math>, where <math>K_{0}</math> is a modified Bessel function and an identity between the modified Bessel function.
  
We can calculate the integral of the inverse Fourier transform of <math>\hat{w}</math> and, using the identity between <math>K_{0}</math> and Hankel function of the first kind (Abramowitz and Stegun [[Stegun]] formula 9.6.4), <math>K_{0}\left( \zeta\right)  =\frac{\mathrm{i}\pi}{2}H_{0}^{\left(  1\right)  }\left( \mathrm{i}\zeta\right)  </math> for <math>Re \zeta\geq0</math>, the displacement for point forcing may be written in the equivalent forms
+
We can calculate the integral of the inverse Fourier transform of <math>\hat{w}</math> and, using the identity between <math>K_{0}</math> and Hankel function of the first kind ([[Abramowitz and Stegun 1964]] formula 9.6.4),  
<center><math>\begin{matrix}
 
w_{=P= }\left(  r\right)  & =\frac{1}{\pi}\sum_{q\in K^{}
 
}qR\left(  q\right)  K_{0}\left(  -\mathrm{i}qr\right)
 
(3-13)\\
 
& =\frac{\mathrm{i}}{2}\sum_{q\in K^{}}qR\left(  q\right)
 
H_{0}^{\left(  1\right)  }\left(  qr\right)  (3-14)
 
\end{matrix}</math></center>
 
where the subscript P of <math>w_{=P= }</math> indicates the response to a point load.
 
 
 
The identity (Erdelyi [[bateman]] formula 1.2 (11) with <math>x=\gamma</math>, <math>y=r</math>, and <math>a=-\mathrm{i}q</math>)
 
 
<center><math>
 
<center><math>
\int_{0}^{\infty}\frac{\cos\left( \gamma r\right)  }{\gamma^{2}-q^{2}}
+
K_{0}\left( \zeta\right)  =\frac{\mathrm{i}\pi}{2}H_{0}^{\left(  1\right)  }\left( \mathrm{i}\zeta\right)   
d\gamma=-\frac{\pi}{\mathrm{i}2q}\exp\left( \mathrm{i}qr\right)
 
  (3-25)
 
 
</math></center>
 
</math></center>
for <math>Im q>0</math>, <math>r>0</math>. The above integration gives the surface displacement for line forcing
+
for <math>Re \zeta\geq0</math>, the displacement for three-dimensions may be written in the equivalent forms
 
<center><math>
 
<center><math>
w_{=L= }\left(  r\right) =\mathrm{i}\sum_{q\in K^{}
+
w\left(  r\right)   =\frac{1}{\pi}\sum_{n=-2}^{\infty}
}R\left(  q\right)  \exp\left(  \mathrm{i}qr\right)  (3-15)
+
q_n R\left(  q_n\right)  K_{0}\left(  -\mathrm{i}q_n r\right)
 
</math></center>
 
</math></center>
where the subscript L of <math>w_{=L= }</math> indicates the response to a line load.
 
 
====Modal expansion of the solutions====
 
 
We have shown that the fundamental solution for finite water depth <math>H</math> is
 
found by first finding the roots in the upper-half plane, <math>K^{}</math>,
 
of the dispersion Eqn.~((3-9)). Numerical calculation is achieved by
 
truncating the sum after some finite number of
 
roots. Straightforward computer code (in MatLab) to find the roots has been
 
given by Fox and Chung [[report]]. After calculating the residue for each
 
root given in Eqn.~((3-17)), the sum in Eqn.~((3-14)) can be
 
calculated by separating it into
 
 
<center><math>
 
<center><math>
w_{=P= }\left(  r\right)  =\frac{\mathrm{i}}{2}q_{=T=
+
=\frac{\mathrm{i}}{2}\sum_{n=-2}^{\infty}q_n R\left(  q_n\right)
}R_{=T= }H_{0}^{\left(  1\right)  }\left(  q_{=T= }r\right)
+
H_{0}^{\left(  1\right)  }\left(  q_n r\right).
-Im \left[  q_{=D= }R_{=D= }H_{0}^{\left(  1\right)
 
}\left(  q_{=D= }r\right)  \right] +\frac{1}{\pi}\sum_{n=1}^{\infty
 
}\mathrm{i}q_{n}R_{n}K_{0}\left(  q_{n}r\right) (3-18)
 
 
</math></center>
 
</math></center>
We used the identities <math>-\mathrm{i}\left(  -q_{=D= }\right)
 
=\left(  -\mathrm{i}q_{=D= }\right)  ^{}<math>, </math>R\left(  q^{
 
}\right)  =\left(  R\left(  q\right)  \right)  ^{}<math> and </math>H_{0}^{\left(
 
1\right)  }\left(  -\zeta^{}\right)  =-\left(  H_{0}^{\left(  1\right)
 
}\left(  \zeta\right)  \right)  ^{}</math>, and denote the residues for the
 
poles in <math>K^{}</math> by <math>R_{=T= }=R\left(  q_{=T= }\right)  </math>,
 
<math>R_{=D= }=R\left(  q_{=D= }\right)  </math>, and <math>R_{n}=R\left(
 
\mathrm{i}q_{n}\right)  <math>, </math>n\in\mathbf{N}</math>, respectively. Note that
 
<math>R\left(  -\gamma_{=D= }^{}\right)  \rightarrow-R_{=D= }^{}</math>
 
as <math>\beta\rightarrow0</math>. Since, from Eqn.~((3-17)), <math>\left|  qR\left(
 
q\right)  \right|  \propto\left|  q\right|  ^{-8}<math> as </math>\left|  q\right|  </math>
 
increases, the sum may be terminated after a relatively small number of terms
 
without significant error, as shown in Fig.~((fig:rsd)). We also notice
 
that the smaller the depth <math>H</math>, the fewer number of roots we need to achieve
 
the desired accuracy, since for a given <math>n</math>, <math>q_{n}</math> increases as <math>H</math> becomes small.
 
  
\begin{figure}[tbh]
+
In two-dimensions the identity ([[Erdelyi et. al.1954]] formula 1.2 (11) with <math>x=\gamma</math>, <math>y=r</math>, and <math>a=-\mathrm{i}q</math>)
\begin{center}
 
\includegraphics[height=7.15cm,width=7.3126cm]{rsd.eps}
 
\caption{Log-log plot of the residues corresponding to the evanescent modes. +, * and <math>\diamondsuit</math> indicate for <math>\omega=0.1,1.0,10</math> respectively. The water depth is <math>H=20\pi</math> (deep water).}
 
(fig:rsd)
 
\end{center}
 
\end{figure}
 
 
 
The surface displacement for line forcing is given by the sum
 
 
<center><math>
 
<center><math>
w_{=L= }\left(  r\right)  =\mathrm{i}R_{=T= }\exp\left(
+
\int_{0}^{\infty}\frac{\cos\left(  \gamma r\right)  }{\gamma^{2}-q^{2}}
\mathrm{i}q_{=T= }r\right)  -2Im \left[  R_{=D=
+
\mathrm{d}\gamma=-\frac{\pi}{\mathrm{i}2q}\exp\left(  \mathrm{i}qr\right)
}\exp\left(  \mathrm{i}q_{=D= }r\right)  \right]
 
+\mathrm{i}\sum_{n=1}^{\infty}R_{n}\exp\left(  -q_{n}r\right)
 
. (3-19)
 
</math></center>
 
We have written the solutions ((3-18)) and ((3-19)) in terms of the
 
travelling, damped travelling and evanescent mode in order to emphasize the
 
behavior of each mode.
 
 
 
We consider here the energy flux due to the wave produced by the force at
 
<math>r=0</math>. From the equation of motion ((2-41)) and <math>\mathbf{v}=\nabla\phi</math>,
 
the energy crossing a surface <math>S</math> in a time period <math>\left[  t,t+T\right]  </math>
 
is
 
<center><math>
 
\rho\int_{t}^{t+T}\int_{S}\phi_{t}\left(  x,y,z,t\right)  \phi_{n}\left(
 
x,y,z,t\right)  d\sigma dt.
 
</math></center>
 
For time-harmonic waves, i.e., <math>\phi=Re \left[  \phi\exp\left(
 
i\omega t\right) \right]  <math>, and </math>T=2\pi/\omega</math>, the above equation becomes
 
<center><math>
 
-i2\rho\pi\int_{S}\left[  \phi^{}\phi_{n}-\phi\phi_{n}^{}\right]
 
d\sigma=2\rho\piIm \int_{S}\phi^{}\phi_{n}d\sigma
 
. (3-52)
 
 
</math></center>
 
</math></center>
We notice that if <math>S</math> is the closed surface of domain <math>\mathcal{V}</math>, then from Green's theorem, we have
+
for <math>Im q>0</math>, <math>r>0</math> can be used. This gives us
 
<center><math>
 
<center><math>
\int_{S}\left[ \phi^{}\phi_{n}-\phi\phi_{n}^{}\right] d\sigma
+
w\left(  r\right) =\mathrm{i}\sum_{n=-2}^{\infty}
=\int_{\mathcal{V}}\left[ \phi^{}\nabla^{2}\phi-\phi\nabla^{2}\phi^{
+
R\left(  q\right) \exp\left( \mathrm{i}qr\right).
}\right]  dV=0, (3-82)
 
 
</math></center>
 
</math></center>
which states that there is no net energy propagation to infinity, i.e., the
 
law of energy conservation. Hence, only the travelling mode <math>H_{0}^{\left(
 
1\right)  }\left(  q_{=T= }r\right)  <math> for the point load and </math>\exp\left(
 
iq_{=T= }r\right)  </math> for the line load carry the energy to infinity since
 
the rest of the solution decays exponentially. Note that the decaying rate of
 
the Hankel function of a real variable is <math>1/\sqrt{r}</math> and the integral on the cylindrical surface gives <math>d\sigma=rdrd\theta dz</math>, thus we have non-zero energy flux for the travelling mode of the point load response.
 
  
The evanescent modes in Eqn.~((3-19)) and ((3-18)) are always real since <math>\mathrm{i}R_{n}</math> in Eqn.~((3-18)) and Eqn.~((3-19)) is real for each <math>n=1,2,\cdots</math>. Hence, the only imaginary term in the responses due to point or line load that give non-zero in Eqn.~((3-52)) is the coefficient of the travelling wave. This corresponds to the travelling waves being the only modes that carry energy away from the load. The damped-travelling and evanescent modes contribute motion that is in phase with the forcing, while the travelling mode has a component in quadrature to the forcing.
+
==Solution for the Velocity Potential==
  
The velocity potential in the water can be found using Eqn.~((3-5)) to
+
The velocity potential in the water can be found in two-dimensions
Eqn.~((3-6)), giving
 
 
<center><math>
 
<center><math>
\phi_{=P= }\left(  r,z\right)  =-\frac{\omega}{2}\sum_{q\in K^{
+
\phi\left(  r,z\right)  =-\frac{\omega}{2}\sum_{n=-2}^{\infty}
}}\frac{R\left(  q\right)  }{\sinh qH}H_{0}^{\left(  1\right)  }\left(
+
\frac{R\left(  q_n\right)  }{\sinh q_nH}H_{0}^{\left(  1\right)  }\left(
qr\right)  \cosh q\left(  z+H\right)
+
q_nr\right)  \cosh q_n\left(  z+H\right)
 
</math></center>
 
</math></center>
for the point load and
+
and for three-dimensions
 
<center><math>
 
<center><math>
\phi_{=L= }\left(  r,z\right)  =-\omega\sum_{q\in K^{}}
+
\phi\left(  r,z\right)  =-\omega\sum_{n=-2}^{\infty}
\frac{R\left(  q\right)  }{q\sinh qH}\exp\left(  \mathrm{i}qr\right)
+
\frac{R\left(  q_n\right)  }{q_n\sinh q_nH}\exp\left(  \mathrm{i}q_nr\right)
\cosh q\left(  z+H\right) .
+
\cosh q_n\left(  z+H\right).
 
</math></center>
 
</math></center>
for the line load.
 
 
Note that the zeros of <math>d\left(  \gamma\right)  \gamma\tanh\gamma H</math>\ are the
 
same as those of the dispersion equation. Thus the singularities of <math>\hat
 
{\phi}<math> and </math>\hat{w}<math> are the same. The term </math>\sinh qH</math> is close to zero for
 
most imaginary roots so these expressions are not directly suitable for
 
numerical computation. The substitution following Eqn.~((3-16)) may be
 
used to give computationally stable expressions.
 
 
Each term of <math>w_{=P= }</math> and <math>w_{=L= }</math> in Eqn.~((3-18)) and
 
Eqn.~((3-19)) is a natural mode of a floating ice sheet whose wavenumber
 
is a root of the dispersion equation. We showed that the mode of <math>q_{=T=
 
}</math>, the travelling mode carries the wave energy outwards, and the modes of
 
<math>q_{=D= }</math> and <math>\mathrm{i}q_{n}</math>, <math>n=1,2,...</math>, the damped
 
travelling and evanescent modes, are exponentially decaying.
 
 
====Summary of the analytic structure of <math>\hat{w====\left(
 
\gamma\right)  </math>}
 
  
Here we summarize and clarify the analytic properties of the complex valued
 
function <math>\hat{w}\left(  \gamma\right)  </math> that have been discussed so far.
 
 
It may be worth reminding ourselves that the Mittag-Leffler expansion of <math>\hat{w}</math> given by Eqn.~((ap-20)) has an extra term <math>\hat{w}\left( 0\right)  =0</math> and is unique. Hence, even though a function such as <math>\hat{w}\left(  \gamma\right)  +1</math> has the same singularities and residues, the resulting series expansion will be different. This uniqueness of the expansion guaranties the uniqueness of the solution derived using the inverse Fourier transform of the series expansion form.
 
 
We consider the question of whether we can reconstruct the Fourier transform <math>\hat{w}</math> from the modified formula of residues given by Eqn.~((3-17)) and the positions of the roots given by the dispersion equation <math>d\left( \gamma\right)  =0</math>. In addition to the residues and poles, we require that <math>\hat{w}\left(  0\right)  =0</math>, to uniquely reconstruct <math>\hat{w}</math>. The answer to the question is yes and no depending on how the formulae are used. It follows that the function <math>\hat{w}\left(  \gamma\right)  </math> for <math>\gamma\in C</math> can be reconstructed from <math>q\in K^{}</math> and <math>R\left(  q\right)  </math> (either Eqn.~((3-16)) or Eqn.~((3-17))) as
 
<center><math>
 
\hat{w}\left(  \gamma\right)  =\sum_{q\in K^{}}\frac{2qR\left(
 
q\right)  }{\gamma^{2}-q^{2}}.
 
</math></center>
 
However, we find the following
 
<center><math>
 
\hat{w}\left(  \gamma\right)  \neq R\left(  \gamma\right)  \sum_{q\in
 
K^{}}\frac{2q}{\gamma^{2}-q^{2}}. (3-80)
 
</math></center>
 
We note that near <math>\gamma=q\in K^{}</math> the left and the
 
right hand sides share the same analytic property, i.e.,
 
<center><math>
 
\int_{C_{\varepsilon}\left(  q\right)  }\hat{w}\left(  \gamma\right)
 
d\gamma=\int_{C_{\varepsilon}\left(  q\right)  }R\left(  \gamma\right)
 
\sum_{q\in K^{}}\frac{2q}{\gamma^{2}-q^{2}}\,d\gamma
 
</math></center>
 
where <math>C_{\varepsilon}\left(  q\right)  </math> is a circular contour of small radius <math>\varepsilon</math> around a pole <math>q</math>. The reason for ((3-80)) is because the function <math>R\left(  \gamma\right)  </math> defined by Eqn.~((3-16)) and Eqn.~((3-17)) introduce zeros and singularities of their own, which change the analytic properties of the right hand side of Eqn.~((3-80)).
 
  
It may seem trivial that the formulae of the residues Eqn.~((3-16)) and Eqn.~((3-17)) are valid only at the pole <math>\gamma=q</math>, nevertheless we have confirmed Eqn.~((3-80)).
+
[[Category:Floating Elastic Plate]]

Latest revision as of 16:40, 8 December 2009


This is a special version of the free-surface Green function which applied when the Floating Elastic Plate boundary condition applies at the free-surface. It reduced to the Free-Surface Green Function in the limit as the plate terms tend to zero and to the Green function for an infinite plate in the limit as the water terms tend to zero.

Two Dimensions

Green function with singularity on the surface

For the case of a floating thin plate we almost always make the assumption of shallow draft and the Green function satisfies the following equations The Green function [math]\displaystyle{ G(x,x^{\prime},z) }[/math] representing outgoing waves as [math]\displaystyle{ |x|\rightarrow \infty }[/math] satisfies

[math]\displaystyle{ \nabla^2 G = 0, -H\lt z\lt 0, }[/math]
[math]\displaystyle{ \frac{\partial G}{\partial z} =0, z=-H, }[/math]
[math]\displaystyle{ {\left( \beta \frac{\partial^4}{\partial x^4} - \gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G} = \delta(x-x^{\prime}), z=0, }[/math]

(note we are only considering singularities at the free surface). It is important to realise that these equations are based on the Frequency Domain Problem and that the exact form of the Green function is dependent on whether we have [math]\displaystyle{ e^{i\omega t} }[/math] or [math]\displaystyle{ e^{-i\omega t} }[/math] dependence. In what follows we assume [math]\displaystyle{ e^{i\omega t} }[/math] dependence.

This problem can be solved (a solution is given below and also in Evans and Porter 2006) to give

[math]\displaystyle{ G(x,x^{\prime},z) = -\sum_{n=-2}^\infty\frac{\sin{(k_n H)}\cos{(k_n(z+H))}}{2\alpha C(k_n)}e^{-k_n|x-x^{\prime}|}, }[/math]

where

[math]\displaystyle{ C(k_n)=\frac{1}{2}\left(h - \frac{(5\beta k_n^4 + 1 - \alpha\gamma)\sin^2{(k_n H)}}{\alpha}\right), }[/math]

and [math]\displaystyle{ k_n }[/math] are the solutions of the Dispersion Relation for a Floating Elastic Plate,

[math]\displaystyle{ \beta k^5 \sin(kH) - k \left(1 - \alpha \gamma \right) \sin(kH) = -\alpha \cos(kH) \, }[/math]

with [math]\displaystyle{ n=-1,-2 }[/math] corresponding to the complex solutions with positive real part, [math]\displaystyle{ n=0 }[/math] corresponding to the imaginary solution with positive imaginary part, since we have assumed [math]\displaystyle{ e^{i\omega t} }[/math] dependence (it would be of negative imaginary part if we had chosen [math]\displaystyle{ e^{-i\omega t} }[/math] dependence) [math]\displaystyle{ n\gt 0 }[/math] corresponding to the real solutions with positive real part. We can also write the Green function as

[math]\displaystyle{ G\left( r,z\right) = -\omega\sum_{n=-2}^{\infty} \frac{R\left( q_n\right) }{q_n\sinh q_nH}\exp\left( \mathrm{i}q_nr\right) \cosh q_n\left( z+H\right) . }[/math]

where [math]\displaystyle{ r=|x-x^{\prime}| }[/math], [math]\displaystyle{ q_n=ik_n }[/math] and we have used a different non-dimensionalisation so that the Dispersion Relation for a Floating Elastic Plate is written as

[math]\displaystyle{ -q^5 \sinh(kH) - k \left(1 - m\omega^2 \right) \sinh(kH) = -\omega^2 \cosh(kH) \, }[/math]

and

[math]\displaystyle{ R\left( q_n\right) =\frac{\omega^{2}q_n}{\omega^{2}\left( 5q_n^{4}+u\right) +H\left[ \left( q_n^{5}+uq\right) ^{2}-\omega^{4}\right] }. }[/math]

Green function with singularity in the water

We can define the Green function [math]\displaystyle{ G(x,x^{\prime},z) }[/math], representing outgoing waves as [math]\displaystyle{ |x|\rightarrow \infty }[/math], which satisfies

[math]\displaystyle{ \nabla^2 G = \delta(x-x^{\prime})\delta(z-z^{\prime}), -H\lt z\lt 0, }[/math]
[math]\displaystyle{ \frac{\partial G}{\partial z} =0, z=-H, }[/math]
[math]\displaystyle{ {\left( \beta \frac{\partial^4}{\partial x^4} - \gamma\alpha + 1\right)\frac{\partial G}{\partial z} - \alpha G} = 0, z=0. }[/math]

Note that this Green function reduces the Free-Surface Green Function in the limit as [math]\displaystyle{ \beta }[/math] and [math]\displaystyle{ \gamma }[/math] tend to zero.

This problem can be solved to give Evans and Porter 2003

[math]\displaystyle{ G(x,x^{\prime},z) = -\sum_{n=-2}^\infty\frac{\cos{(k_n(z^{\prime}+H))}\cos{(k_n(z+H))}}{2k_n C(k_n)}e^{-k_n|x-x^{\prime}|}, }[/math]

where [math]\displaystyle{ C(k_n) }[/math] and [math]\displaystyle{ k_n }[/math] are as before.

Relationship between the two Green functions

If we denote the Green function with the singularity at the free-surface by [math]\displaystyle{ G_s }[/math] and the Green function with the singularity in the water by [math]\displaystyle{ G_w }[/math]then we have the following relationship between them

[math]\displaystyle{ G_s = -\frac{1}{\alpha}\left. \frac{\partial G_w}{\partial z^{\prime}} \right|_{z^{\prime} =0} }[/math]

Three Dimensions

[math]\displaystyle{ G\left( r,z\right) =-\frac{\omega}{2}\sum_{n=-2}^{\infty} \frac{R\left( q_n\right) }{\sinh q_nH}H_{0}^{\left( 1\right) }\left( q_n r\right) \cosh q_n\left( z+H\right) }[/math]

Derivation of the Green function

We present here a derivation of the Green function for both the two and three dimensional problem. The devivation uses the Mittag-Leffler Expansion for the Floating Elastic Plate Dispersion Equation. This derivation is based on the PhD thesis of Hyuck Chung.

Previous Work

In the following sections, the deflection of a plate is expressed by the infinite number of modes that exist in the vibrating plate floating on the water, instead of a finite number of modes that represent the free oscillation of an elastic plate without restraints. The fact that waves in an incompressible fluid can be expressed by an infinite series of natural modes is shown by John 1949 and John 1950 and discussed in Cylindrical Eigenfunction Expansion Each mode being a Hankel function of the first kind in the horizontal (in the three-dimensional case) and a hyperbolic function in the vertical. Kheisin 1967 chapter IV, studied the same problem solved in this chapter and derived the same dispersion equation for the physical variables. He studied the properties of the inverse transform only for the simple shallow water and static load cases, which are approximated versions of the full solutions which is presented here.

Mathematical model

In terms of non-dimensional variables the Floating Elastic Plate equation becomes

[math]\displaystyle{ \left[ \nabla^{4}-m\omega^{2}+1\right] w+\mathrm{i}\omega \phi=p_n. }[/math]

where [math]\displaystyle{ p_n }[/math] is the pressure and [math]\displaystyle{ n }[/math] is the dimension. We assume that [math]\displaystyle{ p_2 = \delta(x) }[/math] and that [math]\displaystyle{ p_3 = \delta(x) \delta(y) }[/math]. [math]\displaystyle{ w }[/math] is the displacement of the surface which is related to the velocity potential by the equation [math]\displaystyle{ \frac{\partial\phi}{\partial z} = \frac{\partial w}{\partial t} }[/math] We also have Laplace's equation for the velocity potential and the bottom condition

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

for [math]\displaystyle{ -\infty\lt x,y\lt \infty, }[/math] [math]\displaystyle{ -H\lt z\lt 0 }[/math] and

[math]\displaystyle{ \phi_{z}=0\quad \mathrm{at} \quad \;z=-H. }[/math]

The system of equations from together with the Sommerfeld Radiation Condition form the BVP, which we will solve for [math]\displaystyle{ w }[/math] and [math]\displaystyle{ \phi }[/math] for a given [math]\displaystyle{ \omega }[/math].

Spatial Fourier transform

We solve the system of equations using the Fourier Transform in [math]\displaystyle{ \left( x,y\right) }[/math]-plane for three-dimensions and with a Fourier Transform in [math]\displaystyle{ x }[/math] for two dimensions. We choose the Fourier transform with respect to [math]\displaystyle{ x }[/math] and [math]\displaystyle{ y }[/math] defined as

[math]\displaystyle{ \hat{\phi}\left( \alpha,k,z\right) =\int_{-\infty}^{\infty}\int_{-\infty }^{\infty}\phi\left( x,y,z\right) e^{\mathrm{i}\left( \alpha x+ky\right) }\mathrm{d}x\mathrm{d}y }[/math]

and the inverse Fourier Transform defined as

[math]\displaystyle{ \phi\left( x,y,z\right) =\frac{1}{4\pi^{2}}\int_{-\infty}^{\infty} \int_{-\infty}^{\infty}\hat{\phi}\left( \alpha,k,z\right) e^{-\mathrm{i}\left( \alpha x+ky\right) }\mathrm{d}\alpha \mathrm{d}k. }[/math]

For the one-dimensional case, the definitions are

[math]\displaystyle{ \hat{\phi}\left( \alpha,z\right) =\int_{-\infty}^{\infty}\phi\left( x,z\right) e^{\mathrm{i}\alpha x}\mathrm{d}x, }[/math]
[math]\displaystyle{ \phi\left( x,z\right) =\frac{1}{2\pi}\int_{-\infty}^{\infty}\hat{\phi }\left( \alpha,z\right) e^{-\mathrm{i}\alpha x}\mathrm{d}\alpha. }[/math]

We denote the spatial Fourier Transform by using a hat over [math]\displaystyle{ w }[/math] and [math]\displaystyle{ \phi }[/math].

The Fourier Transform of both sides of the Laplace's equation becomes an ordinary differential equation with respect to [math]\displaystyle{ z }[/math],

[math]\displaystyle{ \frac{\partial^{2}\hat{\phi}}{\partial z^{2}}\left( \alpha,k,z\right) -\left( \alpha^{2}+k^{2}\right) \hat{\phi}\left( \alpha,k,z\right) =0 }[/math]

with a solution

[math]\displaystyle{ \hat{\phi}\left( \alpha,k,z\right) =A\left( \gamma\right) e^{\gamma z}+B\left( \gamma\right) e^{-\gamma z} }[/math]

where [math]\displaystyle{ \gamma^{2}=\alpha^{2}+k^{2} }[/math] and [math]\displaystyle{ A }[/math], [math]\displaystyle{ B }[/math] are unknown coefficients to be determined by the boundary condition. We find that [math]\displaystyle{ \hat{\phi} }[/math] is a function only of the magnitude of the Fourier transform variables, [math]\displaystyle{ \gamma=\left\| \mathbf{\gamma}\right\| =\left| \left| \left( \alpha,k\right) \right| \right| }[/math], hence we may now denote [math]\displaystyle{ \hat{\phi }\left( \alpha,k,z\right) }[/math] by [math]\displaystyle{ \hat{\phi}\left( \gamma,z\right) }[/math]. We can reach the same conclusion using the fact that [math]\displaystyle{ w\left( x,y\right) }[/math] and [math]\displaystyle{ \phi\left( x,y,z\right) }[/math] are functions of [math]\displaystyle{ r=\sqrt{x^{2}+y^{2}} }[/math] thus the Fourier transforms must be functions of [math]\displaystyle{ \gamma=\sqrt{\alpha^{2}+k^{2}} }[/math].

We find the unknown coefficients [math]\displaystyle{ A }[/math] and [math]\displaystyle{ B }[/math] from the Fourier transformed ocean floor condition that [math]\displaystyle{ \left. \hat{\phi}_{z}\right| _{z=-H}=0 }[/math] to be [math]\displaystyle{ A\left( \gamma\right) =Ce^{\gamma H} }[/math], and [math]\displaystyle{ B\left( \gamma\right) =Ce^{-\gamma H} }[/math]. Thus, we obtain the depth dependence of the Fourier transform of the potential

[math]\displaystyle{ \hat{\phi}\left( \gamma,z\right) =\hat{\phi}\left( \gamma,0\right) \frac{\cosh\gamma\left( z+H\right) }{\cosh\gamma H}. }[/math]

At the surface, [math]\displaystyle{ z=0 }[/math], differentiating both sides with respect to [math]\displaystyle{ z }[/math] the vertical component of the velocity is

[math]\displaystyle{ \hat{\phi}_{z}\left( \gamma,0\right) =\hat{\phi}\left( \gamma,0\right) \gamma\tanh\gamma H. }[/math]

Using this relationship to substitute for [math]\displaystyle{ \hat{\phi}_{z} }[/math] in the non-dimensional form of the kinematic condition, we find that

[math]\displaystyle{ \hat{\phi}\left( \gamma,0\right) =\frac{\mathrm{i}\omega\hat {w}\left( \gamma\right) }{\gamma\tanh\gamma H}. }[/math]

Thus, once we find [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] then we can find [math]\displaystyle{ \hat{\phi}\left( \gamma,z\right) }[/math].

The Fourier transform of the plate equation is also an algebraic equation in the parameter [math]\displaystyle{ \gamma }[/math]

[math]\displaystyle{ \left( \gamma^{4}+1-m\omega^{2}\right) \hat{w}+\mathrm{i}\omega \hat{\phi}=1. }[/math]

Hence, we have the single algebraic equation for each spatial Fourier component of [math]\displaystyle{ w }[/math]

[math]\displaystyle{ \left( \gamma^{4}+1-m\omega^{2}-\frac{\omega^{2}}{\gamma\tanh\gamma H}\right) \hat{w}=1. }[/math]

Notice that we have used the fact that the Fourier transform of the delta function is [math]\displaystyle{ 1 }[/math].

We find that the spatial Fourier transform of the displacement of the ice sheet for the localized forcing, both point and line, is

[math]\displaystyle{ \hat{w}\left( \gamma\right) =\frac{1}{d\left( \gamma,\omega\right) } }[/math]

where

[math]\displaystyle{ d\left( \gamma,\omega\right) =\gamma^{4}+1-m\omega^{2}-\frac{\omega^{2} }{\gamma\tanh\gamma H}. }[/math]

Where [math]\displaystyle{ d\left( \gamma,\omega\right) }[/math] in the Dispersion Relation for a Floating Elastic Plate function, and the associated algebraic equation [math]\displaystyle{ d\left( \gamma ,\omega\right) =0 }[/math] the dispersion equation for waves propagating through an ice sheet. This dispersion relation (and the Fourier transform of it) was derived by Kheisin 1967 chapter IV, and Fox and Squire 1994.

Inverse Fourier Transform

Our task now is to derive the inverse Fourier transform). We notice that the roots of the Dispersion Relation for a Floating Elastic Plate for a fixed [math]\displaystyle{ \omega }[/math], are the poles of the function [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math], which are necessary for calculation of integrals involved in the inverse Fourier transform. There are a pair of two real roots [math]\displaystyle{ \pm q_0 }[/math], an infinite number of pure imaginary roots [math]\displaystyle{ \pm q_{n},\,n\in\mathbb{N} }[/math], ([math]\displaystyle{ \mathrm{Im} q_{n}\gt 0 }[/math]) plus four complex roots (in all physical situations) occur at plus and minus complex-conjugate pairs [math]\displaystyle{ \pm q_{-1} }[/math] and [math]\displaystyle{ \pm q_{-2} }[/math] with [math]\displaystyle{ \mathrm{Im} \left( q_{-1,-2}\right) \gt 0 }[/math]).

We note that the dispersion equation represents a relationship between the spatial wavenumbers [math]\displaystyle{ \gamma }[/math] and the radial frequency [math]\displaystyle{ \omega }[/math], which is how the name "dispersion" came about.

We may also solve for the velocity potential at the surface of the water

[math]\displaystyle{ \hat{\phi}\left( \gamma,0\right) =\frac{\mathrm{i}\omega w} {\gamma\tanh\gamma H}=\frac{\mathrm{i}\omega}{\gamma\tanh\left( \gamma H\right) d\left( \gamma,\omega\right) }, (3-10) }[/math]

which is also a function of [math]\displaystyle{ \gamma }[/math] only and has the same poles as [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] does.

We calculate the displacement [math]\displaystyle{ w\left( x,y\right) }[/math] using the inverse Fourier transform of [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math]. Since [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] is radially symmetric, the inverse transform may be written (Bracewell 1965)

[math]\displaystyle{ w\left( r\right) =\frac{1}{2\pi}\int_{0}^{\infty}\hat{w}\left( \gamma\right) \gamma J_{0}\left( \gamma r\right) \mathrm{d}\gamma }[/math]

for the three dimensional problem, where [math]\displaystyle{ J_{0} }[/math] is Bessel function and [math]\displaystyle{ r }[/math] is the distance from the point of forcing. For the two-dimensional problem the inverse Fourier transform of [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] in the [math]\displaystyle{ x }[/math]-axis and since [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] is an even function, this is

[math]\displaystyle{ w\left( r\right) =\frac{1}{\pi}\int_{0}^{\infty}\hat{w}\left( \gamma\right) \cos\left( \gamma r\right) \mathrm{d}\gamma }[/math]

where again [math]\displaystyle{ r=\left| x\right| }[/math]. Note that the factors [math]\displaystyle{ 1/2\pi }[/math] and [math]\displaystyle{ 1/\pi }[/math] result from the form of the Fourier transform in two and one dimensional spaces.

The integrals are calculated using the Mittag-Leffler Expansion for the Floating Elastic Plate Dispersion Relation where we show that

[math]\displaystyle{ \hat{w}\left( \gamma\right) =\sum_{n=-2}^{\infty}\frac{2q_nR\left( q_n\right) }{\gamma^{2}-q_n^{2}} }[/math]

where [math]\displaystyle{ R\left( q_n\right) }[/math] is the residue of [math]\displaystyle{ \hat{w}\left( \gamma\right) }[/math] at [math]\displaystyle{ \gamma=q_n }[/math] given by

[math]\displaystyle{ R\left( q_n\right) =\frac{\omega^{2}q_n}{\omega^{2}\left( 5q_n^{4}+u\right) +H\left[ \left( q_n^{5}+uq_n\right) ^{2}-\omega^{4}\right] }. }[/math]

.

Using the identities Abramowitz and Stegun 1964 formula 11.4.44 with [math]\displaystyle{ \nu=0 }[/math], [math]\displaystyle{ \mu=0 }[/math], [math]\displaystyle{ z=-\mathrm{i}q }[/math] and [math]\displaystyle{ a=r }[/math])

[math]\displaystyle{ \int_{0}^{\infty}\frac{\gamma}{\gamma^{2}-q^{2}}J_{0}\left( \gamma r\right) \mathrm{d}\gamma=K_{0}\left( -\mathrm{i}qr\right) }[/math]

for [math]\displaystyle{ Im q\gt 0 }[/math], [math]\displaystyle{ r\gt 0 }[/math], where [math]\displaystyle{ K_{0} }[/math] is a modified Bessel function and an identity between the modified Bessel function.

We can calculate the integral of the inverse Fourier transform of [math]\displaystyle{ \hat{w} }[/math] and, using the identity between [math]\displaystyle{ K_{0} }[/math] and Hankel function of the first kind (Abramowitz and Stegun 1964 formula 9.6.4),

[math]\displaystyle{ K_{0}\left( \zeta\right) =\frac{\mathrm{i}\pi}{2}H_{0}^{\left( 1\right) }\left( \mathrm{i}\zeta\right) }[/math]

for [math]\displaystyle{ Re \zeta\geq0 }[/math], the displacement for three-dimensions may be written in the equivalent forms

[math]\displaystyle{ w\left( r\right) =\frac{1}{\pi}\sum_{n=-2}^{\infty} q_n R\left( q_n\right) K_{0}\left( -\mathrm{i}q_n r\right) }[/math]
[math]\displaystyle{ =\frac{\mathrm{i}}{2}\sum_{n=-2}^{\infty}q_n R\left( q_n\right) H_{0}^{\left( 1\right) }\left( q_n r\right). }[/math]

In two-dimensions the identity (Erdelyi et. al.1954 formula 1.2 (11) with [math]\displaystyle{ x=\gamma }[/math], [math]\displaystyle{ y=r }[/math], and [math]\displaystyle{ a=-\mathrm{i}q }[/math])

[math]\displaystyle{ \int_{0}^{\infty}\frac{\cos\left( \gamma r\right) }{\gamma^{2}-q^{2}} \mathrm{d}\gamma=-\frac{\pi}{\mathrm{i}2q}\exp\left( \mathrm{i}qr\right) }[/math]

for [math]\displaystyle{ Im q\gt 0 }[/math], [math]\displaystyle{ r\gt 0 }[/math] can be used. This gives us

[math]\displaystyle{ w\left( r\right) =\mathrm{i}\sum_{n=-2}^{\infty} R\left( q\right) \exp\left( \mathrm{i}qr\right). }[/math]

Solution for the Velocity Potential

The velocity potential in the water can be found in two-dimensions

[math]\displaystyle{ \phi\left( r,z\right) =-\frac{\omega}{2}\sum_{n=-2}^{\infty} \frac{R\left( q_n\right) }{\sinh q_nH}H_{0}^{\left( 1\right) }\left( q_nr\right) \cosh q_n\left( z+H\right) }[/math]

and for three-dimensions

[math]\displaystyle{ \phi\left( r,z\right) =-\omega\sum_{n=-2}^{\infty} \frac{R\left( q_n\right) }{q_n\sinh q_nH}\exp\left( \mathrm{i}q_nr\right) \cosh q_n\left( z+H\right). }[/math]