F Series WAGENINGEN
F Series WAGENINGEN
The hub diameter and chord length at the hub are both a
polynomial resulting from the analysis of the MARIN
database.
The most commonly applied blade profiles by the industry
are used for the F-series as well. For the camber, the NACA
mean line is selected with a=0.6 at the hub to delay flow
separation and a=0.8 at 0.6R towards the tip. For thickness,
Figure 1: Thrust loading coefficient 𝑪𝑻 and cavitation the NACA66 DTMB mod distribution is applied perpen-
number 𝝈𝑵 as function of PD and BAR for the 5-bladed dicular to the nose-tail line.
propellers
The tip speed ranges from 30 m/s for bulkers and LNG 3 F-SERIES DESIGN
carriers towards 50 m/s for fast displacement ships and
fast container ships. The tip speed is important to specify 3.1 Design Approach
the relation between the hydrodynamic pressure and the Nowadays, a competitive propeller can be designed using
hydrostatic pressure. optimization techniques. Therefore, the F-series design w-
as approached using optimization techniques and pol-
The clearance between the propeller tip and the hull
ynomial fitting of the optimized parameters that describe
surface above the propeller was defined as well to
the geometry. The optimization and smoothing was itera-
evaluate hull excitation. The clearance ranges from 20%
tive between the supplied design input, the description of
to 50% of the diameter for a fast displacement ships and
the design variables and the optimization results. This des-
bulkers, respectively.
ign iteration loop is schematically presented by Figure 2.
To estimate the global response levels due to propeller
induced hull excitation, DNV (as a participant to this JIP)
proposed to use an empirical transfer function to relate
the excitation force to the vertical vibration in the aft ship.
The first harmonic excitation force on the square virtual
plane of size 1.5 times the propeller diameter is mult-
iplied by a correlation factor. Hence, vibration is assumed
to be linear with the virtual force. The correlation factor
was determined based on the ship type and its displ-
acement.
The prescription of the astern torque is important for the
propeller blade thickness in the area around 0.8R,
following from the DNV rules. Despite the ambition to
allow full torque astern for each propeller in the F-series Figure 2: Design iteration loop
Performance optimization was performed for a subset of 64 For the design and smoothing towards the geometry of the
propellers which are well distributed within the F-series F-series, the winner from each individual optimization was
matrix of 150 propellers and represent most of the typical hand-selected, for which Table 2 provides a scheme. The
ships. Pareto fronts for the subset of 64 combinations of the main
parameters Z, PD and BAR. From all fronts a winner was
3.2 Definition of Optimization Problem
selected and polynomial regression determined the relation
Multi-objective optimisation techniques provide the of all geometric parameters as a function of these main
various trade-offs between conflicting objectives and parameters, yielding a series description of the geometry.
indicate the influence of design constraints on the
achievable objectives. Table 1: Objectives and constraints within the pro-
To allow large design freedom, the propeller is fully peller optimization
parameterized in the radial pitch, camber and outline (both Objectives Constraints
chord and skew) distributions. These are defined as Bézier Maximize behind Avoid bubble cavitation with at least
curves of which the parameters (control points) form the efficiency 5% margin on cavitation number 𝜎𝑁
optimization variables. Within MARIN and the
Cooperative Research Ships (CRS), the multi-objective Minimize hull Comply to the cavity-planform
optimisation method PropArt has been developed for excitation, criterion for suction-side sheet-
propeller design optimisation studies. calculated as a cavitation such that the sheet
virtual force on a cavitation merges smoothly with tip
A genetic optimization algorithm was used to control the
flat plate above the vortex cavitation.
computational process and to provide input to the geometry
propeller
generator. Typically, 96 geometries are generated within Comply to the time ratio and aspect
each generation. Through subsequent generations, the Minimize ratio criteria for pressure-side
genetic algorithms modifies the design parameters to cavitation volume cavitation as determined in the EU
converge towards the best propellers within the design on all blades EROCAV project as reported by
space. Moulijn et al. (2006).
Each propeller geometry is scaled to the correct PD and
BAR, after which each geometry is subjected to a series of Table 2: Selection process for the F-series survivor
three PROCAL computations: from each particular optimization
1. Steady computation of open water characteristics
2. Unsteady computation in the specified wake field,
iterated to the specified thrust coefficient 𝑪𝑻 using the
inflow speed as iteration variable.
3. Cavitation computation for the same condition as the
unsteady fully-wetted computation.
Propeller computations were done with PROCAL
(Bosschers et al, 2008), a boundary element method (BEM)
which computes the inviscid flow around the propeller
operating in the ship’s wake field. Besides efficiency, the The entire series of 150 propellers would be analysed or its
calculations give the pressure distribution and the extent performance showing goal and constraint behaviour over
and dynamics of the sheet cavitation developing on the the wider design space.
propeller blades. PROCAL has been developed by MARIN In each design iteration (see Figure 2) the best geometries
within the framework of CRS and has been extensively from the optimizations were selected and smoothed as
validated. function of the main parameters PD, BAR and Z.
The results of the PROCAL computations are evaluated Subsequently, the input was further adjusted towards better
and converted into objectives and constraints. Three convergence and smoother polynomial-based description
objectives are defined for the design of the F-series, as of the F-series geometry. High curvature in the polynomial
given in Table 1. Further, in addition to the objectives, surfaces was avoided for the purposes of consistency and
three constraints are used to avoid cavitation erosion risk smoothness, which was regarded as more important than
in the prescribed design condition. the recording of every detail from the individual
optimizations.
3.3 Optimization results
For a normal optimization the Pareto front shows the trade- 3.4 Geometry and Smoothing
offs of all solutions. An example is given in Figure 11 and The following summarizes the iterative path towards the
Figure 12 at the end of this paper. The radial distributions geometry description of the F-series:
are shown as well as the Pareto Fronts (converging from The blade thickness of the F-series follows the DNV
blue to red) and the open water diagrams. Individual rules. The final polynomial description defines the radial
G73P35 was chosen as a winner from this optimization. distribution. For simplicity reasons, it was chosen to
describe the thickness with a simple cubic Bezier curve,
such that both the thickness at 0.6R and 0.8R fulfil the at 0.7R is mainly driven by cavitation constraints, either
rules. The thickness at the tip is equal (0.0035D) for all by the bubble cavitation criterion, the cavity planform or
propellers. the pressure-side cavitation. After smoothing, clear
A fit was made on the input parameters of the blade trends were found. The camber at 0.7R increases with PD
outline function from the selected propellers from the and decreases with BAR and Z.
optimizations. It was found that the chord distribution
was captured by a 4th order Bezier with control points at
the hub, 0.8R and the tip. The majority of the propellers
is relatively well captured by this simple and smooth
description of the chord length.
For the skew a similar procedure was followed. A
parabolic Bezier curve with quadratic polynomial input
captures the majority of the propellers sufficiently
accurate and provides consistency within the F-series. In
general, the wider blades (higher BAR or lower Z) have
higher skew. Figure 4: Camber polynomials at root and 0.7R
The rake of the F-series balances the skew-induced rake
and keeps the blade as tight as possible above the hub so
4 COMPUTED PERFORMANCE
that the propeller can lay down on the pressure side for
storage. The rake was not optimized for hydrodynamic The performance of the 150 F-series propellers was
performance. No special tip-rake was applied. computed by PROCAL.
For the pitch, initial optimizations allowed much freedom 4.1 Efficiency
with a 5th order Bezier curve. Nonetheless, it was found The behind efficiency of the F-series as calculated by
that a simple parabolic Bezier curve, always forced PROCAL will be higher than the B-series as shown by
towards maximum pitch at 0.7R, would describe the Figure 5 (left).
optimization results sufficiently accurate. Only the pitch
reduction at hub and tip are required as design
parameters. Bilinear polynomials were found to capture
the trends from the optimizations well. An impression is
given in Figure 3. The low BAR propellers have more
pitch reduction, compared to higher BAR propellers,
both at the hub and at the tip. According to the
optimizations, a pitch distribution like the F-series
provides an optimum balance between efficiency, hull
excitation and the stability and non-erosiveness of the
sheet cavitation.
The efficiency level is at a level that may be expected from Sketches of the computed cavitation for some of the
contemporary propeller designs. PROCAL is close propellers are presented in Figure 7 for the F5-69 and in
concerning efficiency when close to the design point at J = Figure 8 for two other propellers with extreme BARs.
0.75, and when compared with the scale effect corrected
results. Nonetheless, PROCAL gives overestimation of KT
Figure 7: Cavitation outline for the F5-69 (see table 3
and significant overestimation of KQ at lower J values.
for explanation of colors)
These kind of differences between PROCAL and model
tests are observed regularly for propellers in the MARIN However, there are some outliers at the high PD, high
database. BAR, high Z that have some potentially erosive root
cavitation at both suction side and pressure side, as shown
4.3 Cavitation
in Figure 8 for the F5-120 PD = 1.6 in Figure 8.
The maximum cavitation area in percentage of the blade
area of one blade is presented by Figure 5 (right). In the
wake peak, typically the cavitation area is 20 to 40 percent
of the blade area, except some of the high BAR propellers
with high pitch.
According to the computations, the sheet-cavitation for
most propellers in the F-series behaves exemplary. The
cavitation is stable, increasing in length with radii and
retracts towards the tip, smoothly ending in the tip vortex.
The cavitation behaves in similar consistent manner over
the pitch as shown in Figure 7. The outline of the propeller
blade is presented in black, the cavitation is presented in
different colours according to Table 3. Figure 8: Cavitation outline for two other F-series
Table 3: Legend for the colours in the cavitation figures propellers (see table 3 for explanation of colors)
Black Blue Green Red Purple No clear difference in cavitation behaviour as function of
Non-erosive Growing Complies to Erosion Pressure-
blade number was seen, when compared for equal area of
region phase, outline risk side
one blade, although more super cavitation is observed for
non- criterion, cavitation
the propellers with more blades.
erosive non-erosive
For reference, Figure 9 shows an example of the cavitation Figure 10 (middle) considers the relation between the
behaviour of the B-series, computed in the same maximum pressure fluctuation and the integrated vertical
conditions. In general, the B-series show less stable sheet force which shows a poor correlation.
cavitation. The cavitation of the F-series at similar The empirical approach as proposed by DNV was used to
condition is much closer to contemporary designs. relate the first harmonic virtual force to vertical vibrations
Due to the focus on high efficiency, the cavitation volume levels. Figure 10 (right) presents the predicted vibration
of the F-series is not always less compared to the B-series, levels. Levels up to 4 and 10 mm/s are classified as average
but the behaviour is significantly more favourable. For for ships with aft and forward accommodation
higher BAR, the differences between the behaviour of the respectively.
B-series and F-series become smaller. The approach by DNV only uses the first harmonic
contribution to the excitation force due to the substantial
uncertainties in the measurement or computational
methods for the higher harmonics. Therefore, the higher
harmonics are correlated within the empiricism of the
approach. Nevertheless, during the optimization, the
choice was made to minimize the total virtual force, instead
of just the first harmonic force in order to avoid high levels
for the higher harmonics.
Although both the computations with PROCAL and the
DNV approach are engineering tools, they were considered
capable to steer the characteristics of the F-series
propellers. This analysis indicates that the F-Series are
Figure 9: B-series propellers computed in similar cond- designed with a focus on high efficiency, while staying at
itions as the F-series the average comfort levels.