28.02.2013 Views

Practical Ship Hydrodynamics

Practical Ship Hydrodynamics

Practical Ship Hydrodynamics

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

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

similarly employ a Kutta condition, but this is often omitted. If a Kutta condition<br />

is employed in frequency-domain computations, the wake will oscillate<br />

harmonically in strength. This can be modelled by discrete dipole elements<br />

of constant strength, but for high frequencies this approach requires many<br />

elements. The use of special elements which consider the oscillating strength<br />

analytically is more efficient and accurate, but also more complicated. Such a<br />

‘Thiart element’ has been developed by Professor Gerhard Thiart of Stellenbosch<br />

University and is described in detail in Bertram (1998b), and Bertram<br />

and Thiart (1998). The oscillating ship creates a vorticity. The problem is<br />

similar to that of an oscillating airfoil. The circulation is assumed constant<br />

within the ship. Behind the ship, vorticity is shed downstream with ship speed<br />

V. Then:<br />

� ∂<br />

∂t<br />

�<br />

∂<br />

V<br />

∂x<br />

⊲x, z, t⊳ D 0<br />

is the vortex density, i.e. the strength distribution for a continuous vortex<br />

sheet. The following distribution fulfils the above condition:<br />

⊲x, z, t⊳ D Re⊲ Oa⊲z⊳ Ð e i⊲ωe/V⊳⊲x xa⊳ Ð e iωet ⊳ for x xa<br />

Here Oa is the vorticity density at the trailing edge xa (stern of the ship).<br />

We continue the vortex sheet inside the ship at the symmetry plane y D 0,<br />

assuming a constant vorticity density:<br />

⊲x, z, t⊳ D Re⊲ Oa⊲z⊳ Ð e iωet ⊳ for xa x xf<br />

xf is the leading edge (forward stem of the ship). This vorticity density is<br />

spatially constant within the ship.<br />

A vortex distribution is equivalent to a dipole distribution if the vortex<br />

density and the dipole density m are coupled by:<br />

D dm<br />

dx<br />

The potential of an equivalent semi-infinite strip of dipoles is then obtained<br />

by integration. This potential is given (except for a so far arbitrary ‘strength’<br />

constant) by:<br />

�� � xf zo<br />

8⊲x, y, z⊳ D Re<br />

1 zu<br />

with r D p ⊲x ⊳ 2 C y 2 C ⊲z ⊳ 2 and:<br />

Om⊲ ⊳ y<br />

�<br />

d d eiωet D Re⊲ϕ⊲x, y, z⊳ Ð e 3<br />

r iωet<br />

⊳<br />

�<br />

xf for xa xf<br />

Om⊲ ⊳ D<br />

V �<br />

i⊲ωe/V⊳⊲ xa⊳ 1 e<br />

iωe<br />

� C ⊲xf xa⊳ for 1 xa

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

Saved successfully!

Ooh no, something went wrong!