Conformal Mapping in Hydrodynamic: C. Ciobanu, M. Caţă, A. R. Anghel
Conformal Mapping in Hydrodynamic: C. Ciobanu, M. Caţă, A. R. Anghel
Conformal Mapping in Hydrodynamic: C. Ciobanu, M. Caţă, A. R. Anghel
1. Introduction
During the last 20 years, in the area of naval hydrodynamics as well as in other
domains of mechanical engineering, a growing interest has occurred towards the
algorithmic methods of solving some definite problems.
The shapes of the ships have been described around a point as it is difficult to describe
them according to analytical mappings.
It is necessary to approximate the ship's shape by continuous functions, in order to get
some practical results. A method which has imposed itself during the last few years is that
of multi-parameter conformal mapping, with good results also in the case of extreme
bulbous forms.
The advantage of conformal mapping is that the velocity potential of the fluid around
an arbitrary shape of a cross section in a complex plane can be derived from the more
convenient circular section in another complex plane. In this manner, hydrodynamic
problems can be solved directly by using the coefficients of the mapping function.
The general transformation formula is given by:
n
f ( Z ) s a 2 k 1 Z 2 k 1 ,
k 0
with f ( Z ) z , z x iy is the plane of the ship's cross section.
Z ie e i - is the plane of the unit circle,
s is the scale factor,
a 1 1,
a 2 k 1 are the conformal mapping coefficients ( k 1, ... , n ),
n is the number of parameters.
Therefore we can write in turn:
*
Department of Mathematics-Informatics, Naval Academy “Mircea cel Batran”, Constanta
**
IT Department, Naval Academy “Mircea cel Batran”, Constanta
***
Nautical Sciences Department, Naval Academy “Mircea cel Batran”, Constanta
BULLETIN OF THE TRANSILVANIA UNIVERSITY OF BRAŞOV
n
x iy s a 2 k 1 (ie e i ) ( 2 k 1) ,
k 0
n
x iy s a 2 k 1 (i ) ( 2 k 1) e ( 2 k 1) [cos( ) i sin( )] 2 x 1 ,
k 0
n
x iy s ( 1) k a 2 k 1e ( 2 k 1) [i cos(2k 1) sin(2k 1) ] .
k 0
From the relation between the coordinates in the z - plane (the ship's cross section) and
the variables in the Z - plane (the circular cross section), it follows:
n
x s ( 1) k a 2 k 1e ( 2 k 1) sin(2k 1) ,
k 0
n
y s (1) k a 2 k 1e ( 2 k 1) cos(2k 1) .
k 0
Now by using conformal mapping approximations, the contour of the ship's cross
section, follows from putting 0 in previous sections. We get:
n
x 0 s (1) k a 2 k 1 sin( 2k 1) ,
k 0
n
y 0 s ( 1) k a 2 k 1 cos(2k 1) .
k 0
The breadth on the waterline of the approximate ship's cross section is defined by
n
B0 2 s , with a
k 0
2 k 1 ,
d0
with s .
A simple transformation of the cross sectional hull form will be obtained with n = 2 in
the Lewis transformation. An extended description of the representation of ship hull
forms by Lewis two-parameter conformal mapping is given by Kerczek and Tuck [5] in
1969.
The two-parameter Lewis transformation of a cross section is defined by:
z f ( Z ) s a 1 Z s a1 Z 1 s a3 Z 3 ,
where a 1 1, s is the scale factor and the conformal mapping coefficients a1 and
a 3 are called Lewis coefficients.
Then, for z x iy and Z ie e i , that is Z ie cos i sin , we
have:
x s [e sin( ) a1e sin a3 e 3 sin 3 ], hence
x s (e sin a1e sin a3 e 3 sin 3 )
and
y s [e cos( ) a1e cos a3 e 3 cos 3 ], hence
y s (e cos a1e cos a3 e 3 cos 3 ) .
For 0 we obtain the contour of the so- called Lewis form expressed as:
x0 s (sin a1 sin a 3 sin 3 ) ,
y 0 s (cos a1 cos a3 cos 3 )
or as:
x0 s [(1 a1 ) sin a 3 sin 3 ] ,
y 0 s ((1 a1 ) cos a3 cos 3 ) ,
where the scale factor s is:
Bs Ds
s or s ,
2(1 a1 a3 ) 2(1 a1 a 3 )
in which Bs is the sectional breadth on the load waterline and Ds is the sectional draft.
The half breadth to draft ratio H 0 is given by:
Bs 1 a1 a3
H0 .
2 Ds 1 a1 a 3
An integration of the Lewis form delivers the sectional area coefficient s :
As 1 a12 3a32
s
Bs Ds 4 (1 a 3 ) 2 a12
in which As is the area of the cross section,
d
As 2 xdy s2 (1 a12 3a32 )
0 2
and
Bs Ds 2[(1 a3 ) 2 a12 ] .
Now the coefficients a1 , a 3 and the scale factor s will be determined in such a
manner that the sectional breadth, the draft and the area of the approximate cross section
and of the actual cross section are identical.
BULLETIN OF THE TRANSILVANIA UNIVERSITY OF BRAŞOV
( H 0 1)(1 a3 )
From H 0 expression, we obtain a1 and, by putting that into the
1 H0
expression of s , we get a quadratic equation in a 3 :
2 2
H 1 H 1 3 2
s (1 a3 ) 0
2
(1 a3 ) 2 s 0 (1 a 3 ) 2 a3 ,
H 0 1 4 4 H0 1 4
therefore
4 s 4 s H 0 1 2 H 0 1 8 H 0 1 8 s
2 2 2
3 1 a 2 a3
H 0 1 4 H 0 1 H 0 1
3 s
4
4 H 1 4 s
2
s 0 1 1 0
4 H 0 1
Now, putting
2
H 0 1
m 3 4 s 1 4 s ,
H 0 1
2
H 0 1
n 8 s 2 8 s 2m 6 ,
H 0 1
2
H 0 1
p 4 s 1 4 s 1 m 4 ,
H 0 1
we obtain the equation
ma32 na3 p 0,
that means:
m 3 9 2m
a3 .
m
The Lewis form with the other solution of a 3 in the quadratic equation, namely:
m 3 9 2m
a3
m
are looped. As ships are „better behaving” thus solution is not considered.
It is obvious that a transformation of a half immersed circle with radius r will result for
s r , a1 0 and a 3 0 .
A typical and realistic form is presented below. More precisely we have considered a
dry bulk carrier of 55.000 tone deadweight capacity.
Bulletin of the Transilvania University of Braşov • Vol. 13(48) - 2006 89
3. Java application
The following application was made in Java language and creates both a text file and a
graphical chart.
c1[i]=3+4*ss[i]/Math.PI+(1-4*ss[i]/Math.PI)*((h0[i]-1)/(h0[i]+1))*((h0[i]-1)/
(h0[i]+1));
a3[i]=(-c1[i]+3+Math.sqrt(9-2*c1[i]))/c1[i];
a1[i]=(h0[i]-1)*(1+a3[i])/(1+h0[i]);
if((1+a1[i]+a3[i])==0) ms[i]=0;
else ms[i]=bs[i]/(2*(1+a1[i]+a3[i]));
for(int j=0;j<=18;j++)
{ x0[i][j]=ms[i]*((1+a1[i])*Math.sin(Math.PI*j/36)-
a3[i]*Math.sin(3*Math.PI*j/36));
y0[i][j]=ms[i]*((1-
a1[i])*Math.cos(Math.PI*j/36)+a3[i]*Math.cos(3*Math.PI*j/36));
}
BULLETIN OF THE TRANSILVANIA UNIVERSITY OF BRAŞOV
The pairs of coordinates are obtained from the following initial values of the
parameters: the half breadth (x), the draft (y), and the area of the cross section (as), for a
dry bulk carrier. The half breadth (x), and the draft (y) have been measured from the
midship cross section, and the area of the cross section (as) has been calculated with
Trapezium Method.
x={0, 7.2, 11.6, 14.4, 16, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4, 16.4,
16.4, 16, 13.6, 8, 0};
y={0, 4.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4, 12.4,
12.4, 12.4, 12.4, 12.4, 0};
as={0, 35.092, 147.86, 243.61, 318.06, 368.77, 386.73, 395.16, 397.3, 397.3, 397.3,
397.3, 397.3, 397.3, 397.3, 396.74, 391.72, 368.86, 297, 184.04, 26.66};
The graphical representation of the points shows the contour of the ship’s cross section
of the dry bulk carrier.
This is a mathematical solution in order to obtain the contour of the ship’s cross section
of different types of ships, using conformal mapping approximations.
The advantage of conformal mapping is that the velocity potential of the fluid around
an arbitrary shape of a cross section in a complex plane can be derived from the more
convenient circular section in another complex plane.
References