Difference between revisions of "Rankine Intergral Equations For Ship Flows"

From WikiWaves
Jump to navigationJump to search
m
Line 32: Line 32:
 
We obtain:
 
We obtain:
  
<center><math> \frac{1}{2} \phi (\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} dS_X + \iint_{S_F(Z=0)} \left[ \phi( \vec{X} \frac{\partial G}{\partial Z} (\vec{X}; \vec{\xi}) - G(\vec{X};\vec{\xi}) \frac{\partial\phi}{\partial Z} \right] dX dY </math></center>
+
<center><math> \frac{1}{2} \phi (\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} \mathrm{d}S_X + \iint_{S_F(Z=0)} \left[ \phi( \vec{X} \frac{\partial G}{\partial Z} (\vec{X}; \vec{\xi}) - G(\vec{X};\vec{\xi}) \frac{\partial\phi}{\partial Z} \right] \mathrm{d}X \mathrm{d}Y </math></center>
<center><math> = \iint_{S_B} G(\vec{X};\vec{\xi}) B(\vec{X}) d S_X \, </math></center>
+
<center><math> = \iint_{S_B} G(\vec{X};\vec{\xi}) B(\vec{X}) \mathrm{d} S_X \, </math></center>
  
 
* Note that the integral over the free surface (Z=0) does not vanish since we have not use the relevant wave Green function.
 
* Note that the integral over the free surface (Z=0) does not vanish since we have not use the relevant wave Green function.
Line 45: Line 45:
 
Upon substitution, the second integral over <math> S_F\,</math> becomes:
 
Upon substitution, the second integral over <math> S_F\,</math> becomes:
  
<center><math> I_F = \iint_{Z=0} \left[ \phi(\vec{X}) \frac{\partial G(\vec{G};\vec{\xi})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] dS </math></center>
+
<center><math> I_F = \iint_{Z=0} \left[ \phi(\vec{X}) \frac{\partial G(\vec{G};\vec{\xi})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] \mathrm{d}S </math></center>
  
 
It follows that over <math> Z=0\,</math>, only values and tangential gradients of <math>\phi(\vec{X})\,</math> are now present leading to an integro-differential equation:
 
It follows that over <math> Z=0\,</math>, only values and tangential gradients of <math>\phi(\vec{X})\,</math> are now present leading to an integro-differential equation:
  
<center><math> \frac{1}{2} \phi(\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} dS_X </math></center>
+
<center><math> \frac{1}{2} \phi(\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} \mathrm{d}S_X </math></center>
<center><math> + \iint_{S_F(Z=0)} \left[ \phi(\vec{X}) \frac{\partial G(\vec{X};{\vec{\xi}})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t}- U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] dX dY = \iint_{S_B} V(\vec{X}) G(\vec{X};\vec{\xi}) dS_X \, </math></center>
+
<center><math> + \iint_{S_F(Z=0)} \left[ \phi(\vec{X}) \frac{\partial G(\vec{X};{\vec{\xi}})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t}- U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] \mathrm{d}X \mathrm{d}Y = \iint_{S_B} V(\vec{X}) G(\vec{X};\vec{\xi}) \mathrm{d}S_X \, </math></center>
  
 
* Unknown is <math>\phi(\vec{X})\,</math> over <math>S_B\,</math> & <math>S_F\,</math>. Its X-derivatives may be approximated by carefully selected numerical differentiation schemes forming a core part of Rankine panel methods, discussed below.
 
* Unknown is <math>\phi(\vec{X})\,</math> over <math>S_B\,</math> & <math>S_F\,</math>. Its X-derivatives may be approximated by carefully selected numerical differentiation schemes forming a core part of Rankine panel methods, discussed below.

Revision as of 05:43, 24 December 2008

Rankine integral equations for ship flow problems with forward speed

  • The Green integral equation extends easily to flows past ships in calm water and in waves when the free surface condition is more complex than that of the [math]\displaystyle{ U=0\, }[/math] frequency domain problem.
  • Neumann-Kelvin problem in time domain

Consider a vessel which starts from rest at [math]\displaystyle{ t=0\, }[/math] and translates forward with constant velocity [math]\displaystyle{ U\, }[/math] and also possibly oscillating with amplitudes [math]\displaystyle{ \xi_i(t)\, }[/math] if ambient waves are present.

It was shown earlier that the simplest forward speed free surface condition for the forward problem takes the form:

[math]\displaystyle{ \begin{cases} \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \phi + g \frac{\partial\phi}{\partial Z} = 0, \quad Z=0 \\ \frac{\partial\phi}{\partial n} = V, \quad \mbox{on} \ \bar{S}_B \end{cases} }[/math]

Relative to the ship frame. The normal velocity [math]\displaystyle{ V\, }[/math] on [math]\displaystyle{ \bar{S}_B\, }[/math] can be of three forms:

[math]\displaystyle{ V(\vec{X}) = \begin{cases} U n_1, \quad t\gt 0: \ \mbox{forward translation} \\ n_i \xi_i (t), \quad t\gt 0: \ \mbox{radiation} \\ - \frac{\partial\phi_I}{\partial n}, \quad t\gt 0: \ \mbox{diffraction} \end{cases} }[/math]
  • More general free-surface conditions with space dependent coefficients arising from gradients of the double body flow exist and are described in the literature the steps in deriving the relevant integral equations are very similar to the ones that follows:
  • Wave Green functions that satisfy analytically the time-domain free surface condition stated above exist and are derived in W&L. Their evaluation is how ever time-consuming and they apply only to the Neumann-Kelvin formulation.
  • Proceeding with the derivation of the Green integral equation as above and using the Rankine source as the Green function:
[math]\displaystyle{ \phi_2(\vec{X}) = - \frac{1}{4\pi} \left| \vec{X} - \vec{\xi} \right|^{-1} \equiv G(\vec{X};\vec{\xi}) \, }[/math]

We obtain:

[math]\displaystyle{ \frac{1}{2} \phi (\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} \mathrm{d}S_X + \iint_{S_F(Z=0)} \left[ \phi( \vec{X} \frac{\partial G}{\partial Z} (\vec{X}; \vec{\xi}) - G(\vec{X};\vec{\xi}) \frac{\partial\phi}{\partial Z} \right] \mathrm{d}X \mathrm{d}Y }[/math]
[math]\displaystyle{ = \iint_{S_B} G(\vec{X};\vec{\xi}) B(\vec{X}) \mathrm{d} S_X \, }[/math]
  • Note that the integral over the free surface (Z=0) does not vanish since we have not use the relevant wave Green function.
  • Otherwise the remaining integral over [math]\displaystyle{ S_B\, }[/math] retains its form. The integral over [math]\displaystyle{ S_\infty\, }[/math] can be shown to vanish. The proof if non-trivial and may be found in references.

Over [math]\displaystyle{ Z=0\, }[/math], it follows from the free-surface condition:

[math]\displaystyle{ \frac{\partial\phi}{\partial Z} = - \frac{1}{g} \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right) ^2 \phi, \quad Z=0 }[/math]

Upon substitution, the second integral over [math]\displaystyle{ S_F\, }[/math] becomes:

[math]\displaystyle{ I_F = \iint_{Z=0} \left[ \phi(\vec{X}) \frac{\partial G(\vec{G};\vec{\xi})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] \mathrm{d}S }[/math]

It follows that over [math]\displaystyle{ Z=0\, }[/math], only values and tangential gradients of [math]\displaystyle{ \phi(\vec{X})\, }[/math] are now present leading to an integro-differential equation:

[math]\displaystyle{ \frac{1}{2} \phi(\vec{\xi}) + \iint_{S_B} \phi(\vec{X}) \frac{\partial G(\vec{X};\vec{\xi})}{\partial n_X} \mathrm{d}S_X }[/math]
[math]\displaystyle{ + \iint_{S_F(Z=0)} \left[ \phi(\vec{X}) \frac{\partial G(\vec{X};{\vec{\xi}})}{\partial Z} + \frac{1}{g} G(\vec{X};\vec{\xi}) \left( \frac{\partial}{\partial t}- U \frac{\partial}{\partial X} \right)^2 \phi(\vec{X}) \right] \mathrm{d}X \mathrm{d}Y = \iint_{S_B} V(\vec{X}) G(\vec{X};\vec{\xi}) \mathrm{d}S_X \, }[/math]
  • Unknown is [math]\displaystyle{ \phi(\vec{X})\, }[/math] over [math]\displaystyle{ S_B\, }[/math] & [math]\displaystyle{ S_F\, }[/math]. Its X-derivatives may be approximated by carefully selected numerical differentiation schemes forming a core part of Rankine panel methods, discussed below.
  • The integral over the infinite free surface [math]\displaystyle{ S_F (Z=0) \, }[/math] is truncated at some finite distance from the ship as drawn below
  • A domain denoted by the shaded area is also introduced defined as the "beach". This is located as the outer boundary of [math]\displaystyle{ S_F\, }[/math] and selected so that over its surface the following free surface condition is enforced:
[math]\displaystyle{ \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right)^2 \phi + g \frac{\partial\phi}{\partial Z} + 2 V \left( \frac{\partial}{\partial t} - U \frac{\partial}{\partial X} \right) \phi + V^2 \phi = 0, \quad Z=0 }[/math]
  • This condition differs from the Neumann-Kelvin condition by the addition of the terms that are multiplied by the dissipative parameter [math]\displaystyle{ V(\vec{X})\, }[/math] which varies from [math]\displaystyle{ V=0\, }[/math] at the inner boundary to a finite value at the outer boundary of the beach.
  • It can be shown that converting from the time to the frequency domain via [math]\displaystyle{ \frac{\partial}{\partial t} \to i\omega, V \, }[/math] is the familiar Rayleigh viscosity that plays a key role in the enforcement of the radiation conditions (See W&L).


Ocean Wave Interaction with Ships and Offshore Energy Systems