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.

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

2. Three-dimensional case<br />

The potential of a three-dimensional source is:<br />

S<br />

1<br />

ϕ D<br />

⊲1⊳<br />

4 jEx Exqj<br />

Figure 6.12 shows a triangular patch ABC and a source S. Quadrilateral<br />

patches may be created by combining two triangles. The zero-flow condition<br />

for this patch is<br />

V ⊲Ea ð Eb⊳1<br />

2<br />

A<br />

b<br />

c<br />

C �<br />

C<br />

i<br />

a<br />

B<br />

iMi D 0<br />

Figure 6.12 Source point S and patch ABC<br />

The first term is the volume flow through ABC due to the uniform flow; the<br />

index 1 indicates the x component (of the vector product of two sides of<br />

the triangle). The flow M through a patch ABC induced by a point source<br />

of unit strength is ˛/⊲4 ⊳. ˛ is the solid angle in which ABC is seen<br />

from S. The rules of spherical geometry give ˛ as the sum of the angles<br />

between each pair of planes SAB, SBC, and SCA minus :<br />

˛ D ˇSAB,SBC C ˇSBC,SCA C ˇSCA,SAB<br />

where, e.g.,<br />

ˇSAB,SBC D arctan<br />

[⊲EA ð EB⊳ ð ⊲EB ð EC⊳] Ð EB<br />

⊲EA ð EB⊳ Ð ⊲EB ð EC⊳jEBj<br />

Here EA, EB, EC are the vectors pointing from the source point S to the panel<br />

corners A, B, C. The solid angle may be approximated by AŁ /d2 if the<br />

distance d between patch centre and source point exceeds a given limit.<br />

AŁ is the patch area projected on a plane normal to the direction from the<br />

source to the patch centre:<br />

Ed D 1<br />

3⊲EA C EB C EC⊳<br />

A Ł D 1<br />

2 ⊲Ea ð Eb⊳ Ed<br />

d<br />

With known source strengths i, one can determine the potential and its<br />

derivatives r at all patch corners. From the values at the corners A, B,<br />

C, the average velocity within the triangle is found as:<br />

Ev D r D<br />

A C<br />

En 2 EnAB C<br />

AB<br />

B A<br />

En 2 EnAC<br />

AC

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

Saved successfully!

Ooh no, something went wrong!