13'o Congress

IMAM 2009, istanbul, Turkey,

of Intl. Maritime Assoc. of Mediterranea

12_15 Oct.200

Computerized method for propeller design of optimum diameter and rpm

M.M. GAAFARY

AssociateProress"?r:,;!:::;";'#!g:j

j'"*'::-:#:;TJ#""""Engineering; ".:

ABSTRACT: The process of marine propeller dosign, using theoretical or experimental charts of methodica series' consumes time and effort. This work proposes a computerized desigrr process for optimization o propeller diameter and number of revolutions for better characteristics and higher performance. In thir process, for known propeller rpm, thrust, and blad e area ratio, anoth.r *ol.ston of propeller thrusl coefficient is developed such that the unknown diameter is eliminated. Hence, a new developed curve ol thrust coefficient can be plofted on any arbitrary propeller performance curyes such as ssp; ;;;;r-"1'rn*.0 B-Series charts of prdpeller performance. using- Newton's divided differgnce method of, interpolation, the required optimum propeller diameter can be determined' A similar process has been appried to determine the optimum propeller number of revolutions o1 known propeller diameter. Flow charts and computer programs have been developed to represent the proposed procedure. Applications have been performed on the design of three different propellers of the existing ships. complete agreement is achieved between results of the present work and the real ships propeller
diameters.

NOMENCLATI.IRE

diameter,

Ka

torque coeffi cient of open propeller,

(*)
^o

A^

Kr
blade area ratio, projected area, Wageningen series coefficient of Rn torque coefficient, thrust coefficient,
Po

AP

Pv

Cijt
Cq

Rn
n

effect,

tlirust coefficient of open propeller, static pressure at propeller CL, psi water-vapor pressure, psi Reynolds number, propeller number of revolutions per
second.

Cr D
i,
J

PID

j,

k,l

propeller diameter, integers,

advance ratio,

a
T
t

Kd

thrust coefficient with eliminated

number of revolutions, thrust coefficient with eliminated

Vu

V.

propeller pitch/diameter ratio propeller torque, propeller thrust, thrust deduction factor. propeller advance speed, service ship speed,

K"

211

____.{
.a

,:9

:.i

t,;

,;:.

w Z q p
1.

ship's wake fraction,

number of blades,

propeller efficiency,
water density.

INTRODUCTION

To perform a successful propeller design, a naval

architect is faced with sorne

requirements, such as:

irnportant

1- High efficiency,

2- Sufficient thrust 3- Proper strenglh, 4- No cavitation erosion,

5- Acceptable low vibration and
noise excitation. habitable

Generally, the most efficient propeller is that having the largest possible diameter and lowest rpm. Usually there are two ways to develop a
successful propeller design, namely, the theoretical and experimental methods. There are many moderrr theoretical methods, such as lifting line and lifting

Turnock, et al (2006). It is the main target of this paper to develop a direct computerized design procedure for optimization of propeller diameter and rpm that guarantees the highest possible propeller performance. In this process, for known propeller rpm, thrust, and blade area ratio, another expression of propeller thrust coefficient is developed such that the unknown diameter is eliminated. Hence, a new developed curve of thrust coefficient is plotted on Oosterveld's form of Wageningen B-Series charls of propeller performance, (1975). Using Newton's divided difference method of interpolation, Carnahan, et al ,(1976), the required optimum propeller diameter has been determined from the advance ratio of the highest propeller efficiency. A similar process can be applied to determine the optimum propeller number of revolutions when the propeller diameter is known. Also, this proposed procedure is valid to apply when the propeller torque is known, instead of the propeller thrust, hence, the optimum diameter and rpm can be determined.

surface theories, Breslin, et aI (1994), Kerwin, (1981), Bahgat, (1966), Morgan, (1979, Andersen, et al (1919), Gaafary, 0987 and 1995), and Carlton, (1994). Also, there are many experimental propeller series test results, such as the NSMB Wageningen A- and B-Series, Troost, (1951), SSPA-Series, Lindgren, (196I), and many others. The Wageningen B-Series is still up to now applicable, and it is very efficient method of conventional propeller design, but the area ratios are limited, and a designer usually interpolates between chafts, which makes it a difficult process.

Flow charts and computer programs have been developed to represent both of the Oosterveld's
mathematical polynomial of Wageningen series and

the proposed procedure. Applications have been performed on three different propellers of existing
ships, where complete agreement has been achieved

for propeller particulars deduced method and those of the real ship.
2.

by the proposed

Oosterveld,

et al (1915) has

developed

mathematical polynomial to represent the propeller performance curves of B-Series propeller type. This

useful and simple polynomial form is expressed by its multi-terms of arbttrary propeller's number of

blades, blade area ratio, pitch/diameter ratio,

advance ratio and Reynolds number.

Either theoretical or experimental method is applied the main task of propeller design is the optimization process of propeller particulars, to achieve the highest possible efficiency, Berlram, (2000), and

FORMULATION OF THE METHOD This work presents a proposed procedure for propeller design, where it deals with charts of propeller performance curves of (Kr, Kq, and q) versus (J). These curves are based on theoretical or experimental results of methodical series of marine propellers. Aside from the well known useful charts of Bp - E, of propeller series of B-type of NSMB, Troost, (1951) and Lewis, (1989), the present work introduces some kind of substitute of these design charts. The main target of this work is to develop a simple computerized procedure of propeller design that achieves optimum propeller diameter and number of revolutions. This is to develop the proper theory, then computerize that method, apply the

212

computer program on existing ships of ih" ,u-" propeller t1'pe, and finaliy, eompare the resuits. This might be a direct verification of this computerized method.
2.1 Optimum Propeller Diameter

Using this equation one can get the blad

'' 'A area rario. ( + ), where,

(Ao: no, ).
The propeller thrust can be determined on basis of accurate determination of the total ship resistance, Rr, and the thrust deduction

_^o

Consider a specific systematic series of

conventional marine propellers with available performance curves, such as those of Troost,

factor, (t), as:

. o

"r ! - (1-r)
For the number of propeller blades, (Z), as' it increases, the propeller efficiency and
perf,ormance wi I I i ncrease.

P-

(1951), and Lindgren, (1961), the optimum propeller diameter can be accomplished by eliminating the diameter between the relations defining thrust coefficient, (Kr) and advance ralio,
(J). To reach there a naval architect need to know the following propeller particulars: propeller rpm, 2. speed ofadvance, 3.. blade arearatio, 4. required propeller thrust, and 5. propeller number of blades. These particulars can be determined according to the following:

Now, we might return to our main track of finding the optimum propeller diameter, where

1.

Kr:

Pn2Pa

(1)

tI:vo nD Introduce a new variable, K-"

Breslin, et al
(

(2)
as

developed by

Kn:
get:

To avoid the harm of cavitation, the lower propeller rpm, the less chance of cavitation
occurrence could be,

Manipulating equations

fr

994), n2T
1

(3)

(1,2,

and

3) together

to

. Knowing the ship's

Va:

wake,

(w) at

ship

(4) KnJn This relation of equation (4) can be graphed on a

Kr:

service speed, the speed ofadvance is:

Vs (1 -w) As the blade area ratio increases, there will be less chance of cavitation occwrence, and less propeller efficiency. However, to avoid cavitation or reduce it, an expression of projected blade area has been developed by NSMB, as in Lewis, (1989). When this expression is plotted, it gives a line just

specific Wageningen B-Series chart of marine propellers performance curves of (Kr, Kq, and q) versus (J), of specific number of blades, and blade
possible manual solution can follows,

area ratio, as shown and depicted on Figure (1). A

be secured as

1.

Determine each intersection of (Kr; of (4),

with those (Krt culves which

dependent on varying (PiD) ratio,

are

limit of Burrill's cavitation for merchant diagram ships, Lewis, (1989),

above the upper

2.
a

Determine the different values of (r1) for each intersection. and then draw a curve
connecting them,

this expression is,

(Ar)' :
AI1F -

I-

L360(Po-Pv)I.sVo
Ap

;sqft

J.

where, according to Burrill,

Consider the point at which the curve o{ (q) osculates the envelope of maximum efficiency curyes, The optimum diameter is secured at a specific manually interpolated (P/D) ratio

213

----|r--

that coresponds to the point where the fwo curyes of (n) osculate each other,

5.

The advance ratio at this osculation point determines the optimum diameter,

Dopti.*u

- #

(5)

l. propeller diameter, 2. speed ofadvance, 3. blade area ratio, 4. required propeller thrust, and 5. propeller number of blades. In the same manner,
these parameters can be
determined as previously described.

Until now, this procedr,r" is performed manually for propeller design according to on. of the available methodical series charts of propeller
performance curves.

Now, after eliminating the propeller number of revolutions, (n), between the expressions of (Kr)
and (J). one may get,

It is the target of .this paper to perform this procedure using a newly developed computer
program. During the development of the program, some special steps are made to treat the differences befween applying the method numerically and manually.

Kr
where,

Ka

'J'

Q)
(8)

V r\d and

Tc.
o

wT-

Dffirences when the procedure ls

apptied nttmerically rctther than manually : l. Not to find intersections between Kr-curve of (4), with (Kr; curves of varying (p/D) ratio, but to determine numerically, .by Newton's interpolation, the (P,4D) ratio of Kr- curve at every value of (J) with equal small steps, 2. At every specific value of (J) and (P/D), determine the (Kq) value, hence, determine the propeller efficiency for each value of
(J),

o.s p

1(1v]

(e)

To determine the optimum propeller rpm for a given diameter, one may apply a very similar
procedure to -that previously described. A new curve of Kr, based on equation (7), canbe graphed on a specific Wageningen B-Series chart of marine propellers performance curves of (Kr, Kq, and 11) versus (J), of specific number of blades, and blade arearatio, this process is depicted on Figure (1). A possible manual solution can be secured when following the same first four steps mentioned in the case of optimum diameter, while the fifth step is given instead, advance ratio at this osculation point determines the optimum rpm,

3.

Find the highest efficiency at every specific (J), and hence, the optimum diameter can be determined.

5. The

The applied expression of Newton's

divided difference method of interpolation as given in Carnahan, et al (1976), is, F (xn, xn-1-, ... ,

RPMopumu*

xo:

ILo

-"P_, ttj=s\^i-xj)

The present developed computer program is applicable to determine the optimum propeller
rpm, just as in the case of optimum diameter.

- #

(10)

Hence, this is the formulation of the optimum diameter, now we fum to the optimum rpm formulation.
2.2 Optimum Propeller RPM

3.

PROPOSED COMPUTEzuZED METHOD

3.1 Computer Programfor Wageningen B-Series

For propeller design pulposes, the

is

In a very similar fashion, a similar procedure

developed to yield the optimum propeller rpm for given:

useful experimental test results of Wageningin methodical series of B-type propeller performance of NSMB,

Troost, (1951), were plotted in charts form known

214

as Bp - 6 diagrams. It is true that these B-type propeller chafts are still applicable and yields accurate propeller design. Howevet, aside from Bp -6 diagrams, it is the target of this work to apply the curves of propeller performance of (Kr, Kq, and q)
versus (J) at different (P/D) ratios. Oosterveld, et al

3.2 Computer Programfor Optimum Diameter

(1975), has developed a mathematical polynomial to represent these propeller perfonnance curves of B-Series propeller type. In the present wotk, a flow chart has been designed to represent the logical steps of processing numerically the Oosterveld's mathematical polynomial of Wageningen series, as shown on Figure (2). Based on that flow char1, a computer program has been designed and developed in this work. As an application of the developed computer program to determine the propeller performance values has been performed for a specific propeller of the following particulars:

Considering the formulation steps of the optimum diameter introduced in2.7, another flow chart has been designed to describe the procedure's logic and its executing steps, as shown on Figure ( ) This include5 representation of the mathematical polynornial of Wageningen series of propeller characteristics, calculation of (Kr) as shown in equation (4), interpolation fot (P/D) value at every step of (J) which is done according to Newton's method of interpolat ion. When Newton's equation of interpolation, (6) is applied numerically to determine the required (P/D), using the "pre-" and "post-" values of (Kr; and their corresponding (P/D) values, it yields,

:5 t Z . P/D :0.70 :0.75 . ou no

, ,

/D),"qurred : PD(j) + IPD(j + 1) (12) PD(i\1 v / ' * 1--!I:!IJi)-l 'Kr(j+1)-Kr(j)' At that specific (P/D), calculations of (Kq and q)
(P
are performed. Repeating this process -is required at

o Rn

':2

x1,06

In this numerical

approach, the mathematical

polynomial representation of Wageningen BSeries type of propellers, developed by Oosterveld, et al (1975), is considered for
applications:

every stepped (J). Hence, a decision must be taken by choosing the highest efficiency and its corresponding advance ratio, (J), to determine the optimum diameter as in equation (5). Based on that flow chart, a computer program has been designed and developed specially to perform this proposed procedure numerically, which will save time and effort as a contribution of this present work.
3.3 Computer Programfor Optimum kPM In a very similar manner, the formulation steps of the optimum rpm introduced in sub-section 2.2 are

Kr:

Ito If =o ZJ=oE?=o c,iur. ,t . (X)t . {f,)u 'l'

....(11)
The propeller performance cur-ves of (Kr, Kq, and q) are plotted versus (J) as shown on Figure (3).

Important Feature I: the propeller performance curves shown on Figure (1) of Wageningen series can be reproduced by using this developed program. Almost for any arbitrary B-type propeller, the performance cur-ves can be developed at any
intermediate blade area ratio, any number of blades, (Z), and also at any pitch/diameter ratio, (P/D).

considered, and minor modifications on the flow chart of finding optimum diameter, shown on Figure (4), should be made to represent the procedure's logic and its executing steps. This includes representation of the mathematical polynomial of Wageningen series of propeller characteristics, calculation of (Kr) as shown in equation (7), interpolation for (P/D) value at every step of (J) which is done according to Newton's method. At that specific (P/D) of equation (12),

calculations

of (Kq and n) are performed.

Repeating this process is required at every stepped

ztc

---

(J). Hence, a decision must be taken by choosing

Doptimum

the highest efficiency and its

corresponding advance ratio, (J), to determine the optimum lpm, as in equation (10).

:8.26 m and P/D :0.792. Comment: There is agreement betw-een

the

diameters and P/D ratios.

Minor modifications need to be done to the computer program of finding the optimum
diameter numerically, which effort in finding oplimum rpm.

will

save time and

Important Feature II' A naval alchitect is not obligated to use charls such as Bp - 6 which has fixed blade ratios and then interpolate between them to get the optimum propeller particulars, a process which will consume time and effort. This work provides a quick way to design a propeller with optimum diameter and rpm.

Vs : 16 knots; w : 0.37; t: 0.22 EHP :23,200; Z :5: P/D:0.74 N :85'p*; (ff) :0.733;D :8.8m
Present Method Results:
Doptimum

4.3 Optimum Propeller Diameter of 250,000 Ton Dwt Tanker This tanker has the following data:

'diameters

:8.92 m and P,1D :0.74. Comment: There is agreement befween

and P/D ratios.

the

4.

APPLICATIONS OF THE PROPOSED

present

Important Feature III: The computer program gives accurate output results that almost agree with the real ship propeller diameters.

METHOD

For the purpose of applying the

computerized method on some existing ships, where each set of data has been prepared as an input data to the developed computer program. In the following, there are three different examples

5. ANALYSIS OF THE RESULTS From the present work,.the following remarks and comments can be drawn: l- The proposed computerized method has been developed to be applied almost on any
propeller performance curves, such
as

provided herein together

program output data.

with the

computer

24.1 Optimum Propeller Diameter of SL-7 Container Ship The SL-7 container ship has the following data: Vs :33 knots; w : 0.125; t:0.098

Wageningen B-series propellers, as shown and depicted on Figure (1). Figure (2) shows a flow chart that has been developed in this present work, it gives the logical sequence of executing steps, in order to determine the propeller

performance curyes

EHP : 81,028; Z :6; P/D: 1.13 N : 134 tp*; (Yo) :0.65;D :7.0 m
Present Method Results:
Doptimum

for almost

any

3-

propeller methodical series. A computer program has been developed by following the designed flow chart. As an

:7.2 m and P/D : I.12. Comment: There is agreement between

the

application of the program, performance curves of one of

the the

diameters and PiD ratios.

Wageningen B-series propellers, as shown on Figure (3).

4.2 Optimum Propeller Diameter of a Merchant

Ship

4-

This merchant ship has the following data: Vs : 19 knots; w :0.25;t- 0.25 EHP :26,606; Z : 5: PID : 0.794

Figure (4) shows another flow charl that has been developed to represent the proposed computerized method that might determine the optimum propeller diameter
propeller methodical series. Applications of the present method on three different existing ship propellers, yields

and rym for almost any

10s

tp*;

f#)

:0.75;D

8.2 m

5-

Present Method Results:

216

almost a complete agreement of the results with those of reai ships.

methodical series, such as Gawn-Burrill series, Ktype propeller series ofducted propellers.
REFERENCES
Andersen, P. and Boil, P., Nov. 1979. Propeller Design by Lifting Line Theory,ISH Design Basic Program, Yol.2" Institute Skibs-og Havteknik, DTH, Denrnark. Bahgat. F.. 19b6. Marine Propellets. Al-Maaref Establishment, Alexandria, Egypt. Bertram, V., 2000. Practical Ship Hydrodynamics, Buttelwofth - Heineman-r.r, New York. Breslin, J. P., and AnderserrP., 1994. Hydrodynamics of Ship Propellers, First Edirion, Cambridge

6. A

CONCLUSIONS

proposed computerized method of propeller design has been developed and checked for

applications

on Wageningen B-series

of

conventional marine propellers. The following conclusions can be drawn from the present work:

- 1- For known propeller

determined

rpm, blade area ratio, ship required thrust, and number of blades,
be

the optimum propeller diameter can

Universify Press.
Carlton, J.S., 1994. Marine Propellers & Propuluon, Butterworth and Heinemann, Oxford. Carnahan, B., Luther, H.A., and Wiikes, J.O., 1916. Applied Numerical Methods. John Wiley, New

2-

3-

4-

5-

6-

78-

On the curve of required KT, the Newton's divided difference method of interpolation has been applied to determine intermediate values of pitch/diameter ratios at different advance coelflcients. A computer program has been designed to represent l"he developed method. Optimizing the different values of Kr and Kq for the highest efficiency, yields the determination of optimum propeller diameter at specific (J), and (P/D). The method could be considered as a quick alternative tool that saves time and effort that are consumed in propeller design by the available charts. Applications of the present method on three different existing ship propellers, yields almost a complete agreement of the results with those of real ships. The method is also applicable to determine the optimum propeller rpm. This computerized method is also valid to be applied when the propeller torque is known instead of propeller thrust, by following the same procedure to determine

York. N.Y..
Gaafary, M. M., l.4ay 1987. Forces on Ship Propeller

Blades of Low Aspect-Ratio, Ph.D. Dissertation, Department of Ocean Engineering, Stevens Institut of Tech., Hoboken, N.J., USA. Gaafary, M. M., Younis, G. M., Mosaad, M. A., and Hamdi, T. A., 1995. Hydrodynamics of High Speed

Marine Propellers -- New Approach, IMAM'95,

Dubrovnik, Croatia. Kerwin, J.E., Feb. 1981. Hydrodynamic T'heoryfor Propeller Design and Analysls, Cambridge, Mass., MIT Deparlment of Ocean Engineering. Lewis, Edward V., 1989. Principles of Naval Architecture, Volume Ii, The Sociefy of Naval Architects and Marine Engineers, SNAME. Lindgren, H., 1961 . Model Test with a Family of Three and Five Bladed Propel/ers, SSPA, the Swedish State Shipbuilding Experimental Tank, Nr. 47. Morgan, W"B., Silovic, V., and Denny, S.8., 1979. Propeller Lifting Surface Corrections, Trans., SNAME, Vol. 86. Oosterueld, M.W.C., and Oossanen, P.V., July 1975. Further Computer-Analyzed Data of the
Wageningen B-Screw Series, Intemational Shipbuilding Progress, Y o1. 22,No. 25 1. Troost, L., 1951. Open Water Test Series with Modern Propeller Fornts, NSMB, Wageningen, Part II, and Part III.

Tumock, S.R., Pashias, C. and Rogers,E.,2006. Flow .feattrre identi.fcation for capture of propeller tip

vortex evolution, Proceedings of the

the optimum propeller diameter and

optimum rpm.

26th Symposium on Naval Hydrodynamics. Rome, Italy, INSEAN Italian Ship Model Basin / Office of Naval
Research.

Future Work:

The

present

method might have more

development for

the application on

other

a4a ztt

Kr

to

Ko
rl

K- lor

,\ \ \ \,
/\-to*o \
EHVLOPE

Ha q q

dlch/diem=r.4 / =p,/oL

il !ryigl_]i:Pf qqgg $lpdc.rrcJ RaLioJO,65

i
= t.Lt

, -:
l

OF

Nr: 2 !'

xo ro, nit.r,/arJ),\

clEs

^*.--:(-l.o

Figure

(l): A

scheme for the procedure of finding an optimum propeller diameter using Wageningen 8-6-65 performance curves at a *.nge of (p/D) values, breslin, et al (1994\.

Wageningen B-series of Z=5, Ae/Ad=O.75 at PID=O.7

0,7 0,6 0,5
F tu
od

0,4

***0.7KT at PD= -{-0.7KQ at PD= 10

:a o Fl F

ct

n?
0,2
0,1 0
I I

:z 0.7ETA at PD=

0,1 0,2 0,3 0.4

Advance Ratio,

0.5
J

0,6 0,7

0,8

Figure (3); Present proposed method computer program output of wageningen 8-5-75 propeller perfofinance curves.

218

lnput Data : Blade Area Ratio, Number of blades, n, P/D

Do k= 1.86
Read : Coefficients of Wageningen Sdries oolvnomials. C{k).s(k)- t{k}. u{kl.

EKT{I) =0.

Do K= 1,39 Estimate : thrust coeff. (KT)

EKT(l)=

Estimate

Torque coeff. (KQ)

EKQ(t) =EKQ(t)+Ka(t)

Estimate : Propeller Efficiency (ETA)

J{l)=111;+0.02

OUtput Data : KT, KQ,

',i |:.I
.1 .::

,! ;l

''l
,
tq!-.::!:

Figure (2): Flowchart to determine the propeller performance curves of Wageningen B-series.
,-

:i

.,,;

:i -l:..r r'

::'

iit,,
't:i,:''.
j. r:.

'

'"i'j,.. -,:1'.;

iLif
',i:*ii,,

219

Input Data : Blade Area Ratio, Number of blades, n, P/D

Read : Coefficients of Wageningen .series polynomials, ctk). slk). tlk). ufk). vfk)

Estimate : thrust coeff. (KT)

.q()
-,@ 9!

EKT(l)= EKT(l)+KT(l)
PD(i+1)= PD(i) + 0.075

o6
+.

'tsk r6

9l o9 Eo
F5 Ya
aO

EStiMAtC : KN, KT1, D

xo * q()
F?
t{ c)d NN

Do J= 1, 15

OH
p.6
t;

<0
=Q

Estimate PD1 by Newton's

Method of lnterpolation

9o oo go
\JH

AH

o\< I o.

=*. Fa
Do K= 40, 86

oti 9H

(H*

ov

Estimate

Toroue coeff. (KQ)

-0) >f,

EKQ(I)= rKQ(l)+KQ(l)
vo

Estimate : Propeller Efficiencv (ETA)

Yes

9= boo

ETA < ETA (r-1)

J(l)=1111 +O.02

Output Data: J,KT.KQ, ETA,PD