View source for Cylindrical Eigenfunction Expansion

= Introduction =

There are any situations where we want to expand the three-dimensional linear water wave
solution in cylindrical co-ordinates. For example, scattering from a 
[[Bottom Mounted Cylinder]] or scattering from a [[Circular Floating Elastic Plate]]. In these cases it is easy to find
the solution by an expansion in the cylindrical eigenfunctions. If the depth dependence can be
removed the solution reduces to a two dimensional problem (see [[Removing The Depth Dependence]]). While
the theory here does apply in this two dimensional situtation, the theory is presented here
for the fully three dimensional (depth dependent) case. We begin by assuming the [[Frequency Domain Problem]].

= Outine of the theory = 

{{cylindrical equations}}
{{sommerfeld radiation condition three dimensions}}

{{separation of variables in cylindrical coordinates in finite depth}}

== The case of infinitely deep water == 

A solution will be developed for the same setting as before but under the
assumption of water of infinite depth. As in the previous section,
Laplace's equation must be solved in cylindrical coordinates
satisfying the free surface and the radiation condition. However,
instead of the bed condition, the water velocity potential is also required to
satisfy the depth condition. Therefore, <math>Z(z)</math> must be solved for satisfying the depth condition. It will turn out that in the case of
infinitely deep water an uncountable amount of separation constants
<math>\eta</math> is valid. 

As above, the general solution can be represented as
<center>
<math>
Z(z) = F \mathrm{e}^{\mathrm{i} \eta z} + G \mathrm{e}^{- \mathrm{i} \eta z}, \quad \eta \in \mathbb{C}
\backslash \{0\}.
</math>
</center>
Assuming <math>\eta</math> has got a positive
imaginary part, in order to satisfy the depth condition, <math>F<math> must be
zero. <math>Z(z)</math> then satisfies the free surface condition if <math>\eta</math> is a root of
<center>
<math>
-G \mathrm{i} \eta \mathrm{e}^{-\mathrm{i} \eta z} - \alpha G \mathrm{e}^{-\mathrm{i} \eta z} = 0, \quad  z=0,
</math>
</center>
which yields the dispersion relation 
<center>
<math>
\eta = - \mathrm{i} \alpha.
</math>
</center>
Therefore, <math>\eta</math> must even be purely imaginary. For <math>\Im \eta < 0</math>, 
this is also obtained, but with a minus sign in front of
<math>\eta</math>. However, this yields the same solution. One solution can
therefore be written as
<center>
<math>
Z(z) = G \mathrm{e}^{\alpha z}.
</math>
</center>
Now, <math>\eta</math> is assumed real. In this case, it is convenient to write
the general solution in terms of cosine and sine,
<center>
<math>
Z(z) = F \cos(\eta z) + G \sin(\eta z), \quad \eta \in \mathbb{R}
\backslash \{0\}.
</math>
</center>
This solution satisfies the depth condition automatically.
Making use of the free surface condition, it follows that
<center>
<math>
(-\eta F - \alpha G) \sin (\eta z) + (\eta G - \alpha F) \cos(\eta z)
= 0, \quad z=0,
</math>
</center>
which can be solved for <math>G</math>,
<center>
<math>
G = \frac{\alpha}{\eta} F.
</math>
</center>
Substituting this back gives
<center>
<math>
Z(z) = F \big( \cos(\eta z) + \frac{\alpha}{\eta} \sin(\eta z)
\big) , \quad \eta \in \mathbb{R} \backslash \{0\}. 
</math>
</center>
Obviously, a negative value of <math>\eta</math> produces the same
eigenfunction as the positive one. Therefore, only positive ones are
considered, leading to the definition
<center>
<math>
\psi(z,\eta) := \cos(\eta z) + \frac{\alpha}{\eta} \sin(\eta z), \quad
(z,\eta) \in \mathbb{R}_{\leq0} \times \mathbb{R}_{>0}.
</math>
</center>
This gives the vertical eigenfunctions in infinite depth.

For the radial and angular coordinate the same separation can be used
as in the finite depth case so that the general solution of problem
can be written as

<center><math>
\phi (r,\theta,z) = \mathrm{e}^{\alpha z} \sum_{\nu = -
\infty}^{\infty} \left[ E_\nu (-\mathrm{i} \alpha) I_\nu (-\mathrm{i} \alpha r) + 
F_{\nu} (-\mathrm{i} \alpha) K_\nu (-\mathrm{i} \alpha r) \right] \mathrm{e}^{\mathrm{i} \nu \theta} 
</math></center>
<center><math>
+ \int\limits_0^{\infty} \psi (z,\eta) \sum_{\nu = -
\infty}^{\infty} \left[ E_\nu I_\nu (\eta r) + F_{\nu} (\eta) K_\nu
(\eta r) \right] \mathrm{e}^{\mathrm{i} \nu \theta} \mathrm{d}\eta. 
</math></center>

Making use of the radiation condition as
well as the relations of the Bessel functions in the same way as in
the finite depth case, this can be rewritten as the eigenfunction 
expansion of the water velocity potential into cylindrical outgoing
waves in water of infinite depth, 
<center>
<math>
\phi (r,\theta,z) =  \mathrm{e}^{\alpha z} \sum_{\nu = -
\infty}^{\infty} A_{\nu} (\mathrm{i} \alpha) H_\nu^{(1)} (\alpha r) \mathrm{e}^{\mathrm{i} \nu
\theta} + \int\limits_0^{\infty} \psi (z,\eta) \sum_{\nu = -
\infty}^{\infty} A_{\nu} (\eta) K_\nu (\eta r) \mathrm{e}^{\mathrm{i} \nu
\theta} \mathrm{d}\eta.
</math>
</center>

= Example: Expansion of a plane wave = 

== Infinite depth ==

In Cartesian coordinates centred at the origin, the wavefield due to a plane incident wave travelling in the direction making an angle <math>\chi</math> with the <math>x</math>-axis is
given by

<center><math>
\phi^{\mathrm{In}}(x,y,z) = A \frac{g}{\omega} \, \mathrm{e}^{\mathrm{i}\alpha (x
\cos \chi + y \sin \chi)+ \alpha z},
</math></center>

where <math>A</math> is the amplitude (in displacement).
We want to express the ambient wavefield in the eigenfunction expansion of an
incoming wave in the local coordinates of a body whose mean-centre position is <math>O = (O_x,O_y)</math>. The ambient wave can be represented in an eigenfunction expansion centred at the origin as

<center><math>
\phi^{\mathrm{In}}(x,y,z) = A \frac{g}{\omega} \mathrm{e}^{ \alpha z}
 \sum_{\mu = -\infty}^{\infty} \mathrm{e}^{\mathrm{i}\mu (\pi/2 - \theta + \chi)}
 J_\mu(\alpha r)
</math></center>

(cf. [[Linton and McIver 2001]], p. 169).
Since the local coordinates of the body, that is <math>(r_l,\theta_l,z)</math>, are centred at its mean-centre position, a phase factor has to be defined which accounts for the position from the origin. Including this phase
factor, the ambient wavefield at body is given 
by

<center><math>
\phi^{\mathrm{In}}(r_l,\theta_l,z) = A \frac{g}{\omega} \, e^{\mathrm{i}\alpha (O_x^l
\cos \chi + O_x^l \sin \chi)} \, \mathrm{e}^{\alpha z}
\sum_{\mu = -\infty}^{\infty} \mathrm{e}^{\mathrm{i}\mu (\pi/2 -  \chi)}
J_\mu(\alpha r_l) \mathrm{e}^{\mathrm{i}\mu \theta_l}.
</math></center>


== Finite depth ==

It is easily seen that this is the same as in the infinite-depth case except that we need to replace <math>\mathrm{e}^{ \alpha z}</math> by <math>\frac{\cosh(k (z+d))}{\cosh kd}</math> and <math>\alpha</math> by <math>k</math>. Therefore, we have

<center><math>
\phi^{\mathrm{In}}(r_l,\theta_l,z) = A \frac{g}{\omega} \, e^{\mathrm{i} k (O_x^l
\cos \chi + O_x^l \sin \chi)} \, \frac{\cosh(k (z+d))}{\cosh kd}
\sum_{\mu = -\infty}^{\infty} \mathrm{e}^{\mathrm{i}\mu (\pi/2 -  \chi)}
J_\mu(k r_l) \mathrm{e}^{\mathrm{i}\mu \theta_l}.
</math></center>

[[Category:Problems with Cylindrical Symmetry]]
[[Category:Eigenfunction Matching Method]]