Green Function Method for a Floating Body on the Surface

From WikiWaves
Jump to: navigation, search


Introduction

The problem of a two-dimensional floating body which has negligible submergence is solved using a green function. The problem of a dock is solved in Green Function Method for a Finite Dock and for a floating elastic plate is solved in Green Function Methods for Floating Elastic Plates

Equations for a Finite Plate in Frequency Domain

We consider the problem of small-amplitude waves which are incident on finite floating body occupying water surface for [math]-L\ltx\ltL[/math]. The submergence of the body is considered negligible. We assume that the problem is invariant in the [math]y[/math] direction.

[math] \Delta \phi = 0, \;\;\; -h \lt z \leq 0, [/math]
[math] \partial_z \phi = 0, \;\;\; z = - h, [/math]
[math] \partial_z\phi=\alpha\phi, \,\, z=0,\,x\lt-L,\,\,{\rm or}\,\,x\gtL [/math]

where [math]\alpha = \omega^2[/math]. The equation under the body consists of the kinematic condition

[math] \mathrm{i}\omega w = \partial_z \phi,\,\,\, z=0,\,\,-L\leq x\leq L [/math]

plus the kinematic condition. The body motion is expanded using the modes for heave and pitch. Using the expression [math]\partial_n \phi =\partial_t w[/math], we can form

[math] \frac{\partial \phi}{\partial z} = i\omega \sum_{n=0,1} \xi_n X_n(x) [/math]

where [math]\xi_n \,[/math] are coefficients to be evaluated. The functions [math]X_n(x)[/math] are given by

[math] X_0 = \frac{1}{\sqrt{2L}} [/math]

and

[math] X_1 = \sqrt{\frac{3}{2L^3}} x [/math]

Note that this numbering is non-standard for a floating body and comes from Eigenfunctions for a Uniform Free Beam.


Equation in Terms of the Modes

The equations are

[math] \Delta \phi = 0, \;\;\; -h \lt z \leq 0, [/math]
[math] \partial_z \phi = 0, \;\;\; z = - h, [/math]
[math] \partial_z\phi=\alpha\phi, \,\, z=0,\,x\lt-L,\,\,{\rm or}\,\,x\gtL [/math]
[math] \mathrm{i}\omega\sum_{n=0,1}\zeta_{n}X_{n} =\partial_{z}\phi,\,\,x\in (-L,L),\,\, z=0, [/math]
[math] \left(- \alpha\gamma + 1\right) \zeta_{n}=-i\omega \int_{-L}^{L}\phi X_{n}\mathrm{d}x, \,\,x\in (-L,L),\,\, z=0, [/math]

We solve for the potential (and displacement) as the sum of the diffracted and radiation potentials in the standard way, as for a rigid body.

[math] \phi=\phi^{\mathrm{D}}+\phi^{\mathrm{R}} ,\, [/math]

We begin with the diffraction potential [math]\phi^{\mathrm{D}}[/math] which satisfies the following equations

[math] \Delta\phi^{\mathrm{D}} =0,\,\,-h\ltz\lt0, [/math]
[math] \partial_{z}\phi^{\mathrm{D}} =0,\,\,z=-h, [/math]
[math] \partial_{z}\phi^{\mathrm{D}} =\alpha \phi^{\mathrm{D}},\,\,x\notin(-L,L),\,\, z=0, [/math]
[math] \partial_{z}\phi^{\mathrm{D}} =0,\,\,x\in(-L,L),\,\,z=0. [/math]

[math]\phi^{\mathrm{D}}[/math] satisfies the Sommerfeld Radiation Condition

[math] \frac{\partial}{\partial x} \left(\phi^{\mathrm{D}}-\phi^{\rm I} \right) \pm k_0\left( \phi^{\mathrm{D}}-\phi^{\rm I}\right) = 0 ,\,\,\mathrm{as} \,\,x\rightarrow\infty. [/math]

[math]\phi^{\mathrm{I}}\,[/math] is a plane wave travelling in the [math]x[/math] direction,

[math] \phi^{\mathrm{I}}(x,z)=A \phi_0(z) e^{\mathrm{i} k x} \, [/math]

where [math]A [/math] is the wave amplitude (in potential) [math]\mathrm{i} k [/math] is the positive imaginary solution of the Dispersion Relation for a Free Surface (note we are assuming that the time dependence is of the form [math]\exp(-\mathrm{i}\omega t) [/math]) and

[math] \phi_0(z) =\frac{\cosh k(z+h)}{\cosh k h} [/math]

We now consider the radiation potentials [math]\phi^{\mathrm{R}}[/math]. We can express the radiation potential as:

[math] \phi^{\mathrm{R}}=\sum_{n=0,1}\zeta_n \phi_n^{\mathrm{R}} [/math]

which satisfy the following equations

[math] \Delta\phi_n^{\mathrm{R}} =0,\,\,-h\ltz\lt0, [/math]
[math] \partial_{z}\phi_n^{\mathrm{R}} =0,\,\,z=-h, [/math]
[math] \partial_{z}\phi_n^{\mathrm{R}} =\alpha\phi_n^{\mathrm{R}},\,\,x\notin(-L,L),\, \, z=0 [/math]
[math] \partial_{z}\phi_n^{\mathrm{R}} = i\omega X_{n},\,\,x\in(-L,L),\,\,z=0. [/math]

The radiation condition for the radiation potential is

[math] \frac{\partial\phi_n^{\mathrm{R}}}{\partial x}\pm ik\phi_n^{\mathrm{R}}=0,\,\,\mathrm{as} \,\,x\rightarrow\pm\infty. [/math]

Therefore we find the potential as

[math] \left( - \alpha\gamma + 1\right) \zeta_{n}=-i\omega \int_{-L}^{L}\phi^{\mathrm{D}} X_{n}\mathrm{d}x + \sum_{m=0,1}\left(\omega^2 a_{mn}(\omega) - i\omega b_{mn}(\omega)\right) \zeta_{m}, [/math]

where the functions [math]a_{mn}(\omega)[/math] and [math]b_{mn}(\omega)[/math] are given by

[math] \omega^2 a_{mn}(\omega) -i\omega b_{mn}(\omega) = - i\omega\int_{-L}^{L}\phi_m^{\mathrm{R}}X_{n}\mathrm{d}x, [/math]

and they are referred to as the added mass and damping coefficients (see Added-Mass, Damping Coefficients And Exciting Forces) respectively.

Note that for this simple example the added mass and damping matrices are diagonal.

Solution for the Radiation and Diffracted Potential

We use the Free-Surface Green Function for two-dimensional waves, with singularity at the water surface since we are only interested in its value at [math]z=0[/math] (details about this method can be found in Integral Equation for the Finite Depth Green Function at Surface). Using this we can transform the system of equations to

[math] \phi^{\mathrm{D}}(x) = \phi^{\mathrm{I}}(x) + \int_{-L}^{L}G(x,\xi) \alpha\phi^{\mathrm{D}}(\xi) \mathrm{d} \xi [/math]

and

[math] \phi_n^{\mathrm{R}}(x) = \int_{-L}^{L}G(x,\xi) \left( \alpha\phi_n^{\mathrm{R}}(\xi) - i\omega X_n(\xi) \right)\mathrm{d} \xi [/math]

Reflection and Transmission Coefficients

The Reflection and Transmission Coefficients represent the ratio of the amplitude of the reflected or transmitted wave to the amplitude of the incident wave. Conservation of energy means that [math] |R|^2+|T|^2=1\,[/math].

A diagram depicting the area [math]\Omega\,[/math] which is bounded by the rectangle [math]\partial \Omega \,[/math]. The rectangle [math]\partial \Omega \,[/math] is bounded by [math] -h \leq z \leq 0 \,[/math] and [math]-\infty \leq x \leq \infty \,[/math] or [math]-N \leq x \leq N\,[/math]

We can calculate the Reflection and Transmission coefficients by applying Green's theorem to [math]\phi\,[/math] and [math]\phi^{\mathrm{I}}\,[/math] [math]\phi^{\mathrm{I}}\,[/math] is a plane wave travelling in the [math]x[/math] direction,

[math] \phi^{\mathrm{I}}(x,z)=A \phi_0(z) e^{\mathrm{i} k x} \, [/math]

where [math]A [/math] is the wave amplitude (in potential) [math]\mathrm{i} k [/math] is the positive imaginary solution of the Dispersion Relation for a Free Surface (note we are assuming that the time dependence is of the form [math]\exp(-\mathrm{i}\omega t) [/math]) and

[math] \phi_0(z) =\frac{\cosh k(z+h)}{\cosh k h} [/math]

We assume that [math]A=1[/math]. This gives us

[math] \iint_{\Omega}(\phi\Delta\phi^{\mathrm{I}} - \phi^{\mathrm{I}}\Delta\phi)\mathrm{d}x\mathrm{d}z = \int_{\partial\Omega}(\phi \partial_n \phi^{\rm I} - \phi^{\rm I}\partial_n\phi)\mathrm{d}s = 0, [/math]

This means that (using the far field behaviour of the potential [math]\phi[/math])

[math] \int_{\partial\Omega_{B}} (\phi \partial_n \phi^{\rm I} - \phi^{\rm I}\partial_n\phi)\mathrm{d}s + 2k_0 R \int_{-h}^{0} \left(\phi_0(z)\right)^2 \mathrm{d}z = 0, [/math]

For the present case the body is present only on the surface and we therefore have

[math] \int_{-L}^{L} e^{-k_0 x} \left(\alpha \phi(x) - \partial_n \phi(x)\right)\mathrm{d}x + 2k_0 R \int_{-h}^{0} \left(\phi_0(z)\right)^2 \mathrm{d}z = 0 [/math]

Therefore

[math] R = -\frac{\int_{-L}^{L} e^{-k_0 x} \left(\alpha \phi(x) - \partial_n \phi(x)\right)\mathrm{d}x } {2 k_0 \int_{-h}^{0} \left(\phi_0(z)\right)^2 \mathrm{d}z}. [/math]

and using a wave incident from the right we obtain

[math] T = 1 - \frac{\int_{-L}^{L} e^{k_0 x} \left(\alpha \phi(x) - \partial_n \phi(x)\right)\mathrm{d}x } {2 k_0 \int_{-h}^{0} \left(\phi_0(z)\right)^2 \mathrm{d}z}. [/math]

Note that an expression for the integral in the denominator can be found in Eigenfunction Matching for a Semi-Infinite Dock

Matlab Code

A program to calculate the solution in elastic modes can be found here

rigid_plate_modes.m

Additional code

This program requires