28.02.2013 Views

Practical Ship Hydrodynamics

Practical Ship Hydrodynamics

Practical Ship Hydrodynamics

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

216 <strong>Practical</strong> <strong>Ship</strong> <strong>Hydrodynamics</strong><br />

computed by simple numerical integration if the integrands are transformed<br />

analytically to remove singularities. In the formulae for this element, En is the<br />

unit normal pointing outward from the body into the fluid, c � the integral over<br />

S excluding the immediate neighbourhood of Exq, andr the Nabla operator<br />

with respect to Ex.<br />

1. Two-dimensional case<br />

A Rankine source distribution on a closed body induces the following potential<br />

at a field point Ex:<br />

�<br />

⊲Ex⊳ D ⊲Exq⊳G⊲Ex, Exq⊳ dS<br />

S<br />

S is the surface contour of the body, the source strength, G⊲Ex, Exq⊳ D<br />

⊲1/2 ⊳ ln jEx Exqj is the Green function (potential) of a unit point source.<br />

Then the induced normal velocity component is:<br />

�<br />

vn⊲Ex⊳ DEn⊲Ex⊳r ⊲Ex⊳ D C ⊲Exq⊳En⊲Ex⊳rG⊲Ex, Exq⊳ dS C<br />

S<br />

1<br />

2 ⊲Exq⊳<br />

Usually the normal velocity is given as boundary condition. Then the important<br />

part of the solution is the tangential velocity on the body:<br />

�<br />

vt⊲Ex⊳ D Et⊲Ex⊳r ⊲Ex⊳ D C ⊲Exq⊳Et⊲Ex⊳rG⊲Ex, Exq⊳ dS<br />

S<br />

Without further proof, the tangential velocity (circulation) induced by a<br />

distribution of point sources of the same strength at point Exq vanishes:<br />

�<br />

C rG⊲Ex, Exq⊳Et⊲Ex⊳ dS D 0<br />

S⊲Ex⊳<br />

Exchanging the designations Ex and Exq and using rG⊲Ex, Exq⊳ D rG⊲Exq, Ex⊳,<br />

we obtain:<br />

�<br />

C rG⊲Ex, Exq⊳Et⊲Exq⊳ dS D 0<br />

S<br />

We can multiply the integrand by ⊲Ex⊳ – which is a constant as the integration<br />

variable is Exq – and subtract this zero expression from our initial<br />

integral expression for the tangential velocity:<br />

�<br />

�<br />

vt⊲Ex⊳ D C ⊲Exq⊳Et⊲Ex⊳rG⊲Ex, Exq⊳ dS C ⊲Ex⊳rG⊲Ex, Exq⊳Et⊲Exq⊳ dS<br />

S<br />

S<br />

� �� �<br />

�<br />

D0<br />

D C [ ⊲Exq⊳Et⊲Ex⊳<br />

S<br />

⊲Ex⊳Et⊲Exq⊳]rG⊲Ex, Exq⊳ dS<br />

For panels of constant source strength, the integrand in this formula tends to<br />

zero as Ex !Exq, i.e. at the previously singular point of the integral. Therefore<br />

this expression for vt can be evaluated numerically. Only the length 1S of<br />

the contour panels and the first derivatives of the source potential for each<br />

Ex, Exq combination are required.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!