Difference between revisions of "Interaction Theory for Cylinders"
Line 156: | Line 156: | ||
-\infty}^{\infty} A_{\tau}^j H^{(1)}_{\tau - \nu} (k | -\infty}^{\infty} A_{\tau}^j H^{(1)}_{\tau - \nu} (k | ||
R_{jl}) \mathrm{e}^{\mathrm{i}(\tau -\nu) \varphi_{jl}} \Big] | R_{jl}) \mathrm{e}^{\mathrm{i}(\tau -\nu) \varphi_{jl}} \Big] | ||
− | + | + A_{\nu}^l = - | |
\frac{J'_\nu(k a_l)}{H^{(1)}_\nu{}'(k a_l)} | \frac{J'_\nu(k a_l)}{H^{(1)}_\nu{}'(k a_l)} | ||
\tilde{D}_{\nu}^{l}, | \tilde{D}_{\nu}^{l}, |
Revision as of 09:49, 14 May 2007
Introduction
We present an illustrative example of an interaction theory for the case of [math]\displaystyle{ n }[/math] Bottom Mounted Cylinders. This theory was presented in Linton and Evans 1990 and it can be derived from the Kagemoto and Yue Interaction Theory by simply assuming that each body is a cylinder.
Equations of Motion
After we have Removed the Depth Dependence the problem consists of [math]\displaystyle{ n }[/math] cylinders of radius [math]\displaystyle{ a_j }[/math] whose center is at [math]\displaystyle{ (x_j,y_j) }[/math] subject to Helmholtz's Equation
where [math]\displaystyle{ k }[/math] is the positive real root of the Dispersion Relation for a Free Surface
where [math]\displaystyle{ k_\infty }[/math] is the Infinite Depth wave number and [math]\displaystyle{ d }[/math] is the water depth.
Eigenfunction expansion of the potential
Each cylinder is subject to an incident potential and moves in response to this incident potential to produce a scattered potential. Each of these is expanded using the Cylindrical Eigenfunction Expansion The scattered potential of cylinder [math]\displaystyle{ j }[/math] can be expressed as
with discrete coefficients [math]\displaystyle{ A_{\mu}^j }[/math], where [math]\displaystyle{ (r_j,\theta_j) }[/math] are polar coordinates centered at center of the [math]\displaystyle{ j }[/math]th cylinder.
The incident potential upon cylinder [math]\displaystyle{ j }[/math] can be also be expanded in regular cylindrical eigenfunctions,
with discrete coefficients [math]\displaystyle{ D_{\nu}^j }[/math].
In these expansions, [math]\displaystyle{ J_\nu }[/math] and [math]\displaystyle{ H^{(1)}_\nu }[/math] denote Bessel and Hankel function respectively (Bessel functions) both of first kind and order [math]\displaystyle{ \nu }[/math]. For comparison with the Kagemoto and Yue Interaction Theory (which is written slightly differently), we remark that, for real [math]\displaystyle{ x }[/math],
with [math]\displaystyle{ I_\nu }[/math] and [math]\displaystyle{ K_\nu }[/math] denoting the modified Bessel functions of first and second kind, respectively, both of order [math]\displaystyle{ \nu }[/math].
Derivation of the system of equations
A system of equations for the unknown coefficients of the scattered wavefields of all cylinders is developed. This system of equations is based on transforming the scattered potential of cylinder [math]\displaystyle{ j }[/math] into an incident potential upon cylinder [math]\displaystyle{ l }[/math] ([math]\displaystyle{ j \neq l }[/math]). Doing this for all cylinders simultaneously, and relating the incident and scattered potential for each cylinder, a system of equations for the unknown coefficients is developed.
The scattered potential [math]\displaystyle{ \phi_j^{\mathrm{S}} }[/math] of cylinder [math]\displaystyle{ j }[/math] needs to be represented in terms of the incident potential [math]\displaystyle{ \phi_l^{\mathrm{I}} }[/math] upon cylinder [math]\displaystyle{ l }[/math], [math]\displaystyle{ j \neq l }[/math]. This can be accomplished by using Graf's Addition Theorem
where [math]\displaystyle{ (R_{jl},\varphi_{jl}) }[/math] are the polar coordinates of the mean centre position of cylinder [math]\displaystyle{ l }[/math] in the local coordinates of cylinder [math]\displaystyle{ j }[/math].
Making use of the eigenfunction expansion as well as Graf's Addition Theorem, the scattered potential of cylinder [math]\displaystyle{ j }[/math] can be expressed in terms of the incident potential upon cylinder [math]\displaystyle{ l }[/math] as
The ambient incident wavefield [math]\displaystyle{ \phi^{\mathrm{In}} }[/math] can also be expanded in the eigenfunctions corresponding to the incident wavefield upon cylinder [math]\displaystyle{ l }[/math] (cf. the example in Cylindrical Eigenfunction Expansion)
Let [math]\displaystyle{ \tilde{D}_{\nu}^{l} }[/math] denote the coefficients of this ambient incident wavefield in the incoming eigenfunction expansion for [math]\displaystyle{ \Delta_l }[/math] so that
The total
incident wavefield upon cylinder [math]\displaystyle{ j }[/math] can now be expressed as
which can be written as
Therefore it follows that
Final Equations
The scattered and incident potential can be related by the Diffraction Transfer Matrix for a Bottom Mounted Cylinder so that,
This gives the required equations to determine the coefficients of the scattered wavefields of all bodies,
[math]\displaystyle{ \mu \in \mathbb{Z} }[/math], [math]\displaystyle{ l=1,\dots,n }[/math].