ShipCalculators.com

Propeller Theory: From Momentum to CFD

Contents

Why propeller theory matters

The propeller converts shaft torque into thrust with a typical open-water efficiency of 0.65 to 0.72 for modern merchant propellers. That figure sits at the centre of every fuel budget: a 3-percentage-point improvement in propulsive efficiency on a VLCC consuming 80 tonnes of fuel oil per day saves roughly 2.4 tonnes daily, or about 1,800at1,800 at 750/tonne bunker cost. Propeller theory is the quantitative basis for achieving those improvements. It also underpins prediction of cavitation erosion, hull-pressure vibration at blade rate, and the energy balance of the full propulsion chain from combustion to thrust.

The marine propeller article covers geometry, materials, and manufacturing. The marine propeller pitch and construction article covers pitch distributions, blade-area standards, and class society tolerance rules. This article works through the analytical hierarchy from first principles: momentum theory, blade-element theory, lifting-line theory, the Wageningen B-series open-water diagram, the propulsive coefficient chain, and cavitation prediction. Each tier adds fidelity and computational cost; each has a specific role in the design workflow.

Momentum theory and the actuator-disk ideal

The Rankine-Froude actuator disk

Rankine (1865) and Froude (1889) developed the simplest model of a propulsor: an infinitely thin disk of area A=πD2/4A = \pi D^2 / 4 that imparts uniform axial velocity to the through-passing flow. The disk is an idealisation that deliberately omits blade geometry, viscosity, and rotation; it delivers an upper bound on propeller efficiency.

The analysis applies the momentum theorem to a control volume enclosing the disk. Let VaV_a be the far-upstream axial velocity (the speed of advance), uu the induced axial velocity at the disk, and 2u2u the ultimate wake velocity increment. Conservation of mass, momentum, and energy gives the thrust:

T=ρA(Va+u)2uT = \rho A (V_a + u) \cdot 2u

The thrust loading coefficient is:

CT=T12ρVa2AC_T = \frac{T}{\frac{1}{2} \rho V_a^2 A}

Expressing uu in terms of CTC_T:

u=Va2(1+CT1)u = \frac{V_a}{2} \left( \sqrt{1 + C_T} - 1 \right)

The power delivered to the fluid equals the rate of kinetic-energy addition to the slipstream; the ideal efficiency is the ratio of useful thrust power TVaT V_a to total fluid power input:

ηideal=TVaTVa+slipstream KE loss=21+1+CT\eta_{\text{ideal}} = \frac{T V_a}{T V_a + \text{slipstream KE loss}} = \frac{2}{1 + \sqrt{1 + C_T}}

For a merchant tanker propeller with CT0.8C_T \approx 0.8, this gives ηideal0.87\eta_{\text{ideal}} \approx 0.87. A real four-bladed propeller of the same thrust loading will achieve ηO0.70\eta_O \approx 0.70 to 0.720.72, which is about 82% of the ideal bound. The 18% deficit comes from viscous losses on the blade sections, the tip-vortex trailing-vortex drag, hub vortex losses, and finite-blade-number effects that the actuator-disk model ignores entirely.

Slip

The apparent slip is the difference between the geometric pitch advance and the actual advance per revolution:

s=PVa/nP=1J(D/P)1s = \frac{P - V_a/n}{P} = 1 - \frac{J \cdot (D/P)}{1}

where JJ is the advance coefficient (next section) and P/DP/D is the pitch-diameter ratio. Typical design-point slip for merchant propellers is 0.05 to 0.15. Slip is not thermodynamic loss; it is the kinematic signature of momentum transfer. The fraction of slip that does represent energy loss is captured in the open-water efficiency rather than in the slip figure itself.

Limits of momentum theory

Momentum theory cannot predict blade shape, pitch, number of blades, or the radial distribution of loading. It gives no information about cavitation. It also assumes uniform axial induction and ignores swirl; in reality, the propeller imparts rotational kinetic energy to the wake, which is irreversible loss not captured by the basic actuator-disk. The extended Rankine-Froude theory corrects for rotational wake by treating the disk as also imparting angular momentum, but real blade geometry determines the swirl distribution. Momentum theory is best used as a design benchmark and for rapid parametric studies of disc loading limits.

Blade-element theory

Concept and setup

Blade-element theory (BET) divides each blade into thin radial strips of width drdr. Each strip is treated as a two-dimensional aerofoil section operating at the local resultant inflow velocity. This inflow is the vector sum of the axial component VaV_a and the rotational component ωr\omega r, where ω=2πn\omega = 2\pi n is the angular velocity and rr is the radius of the strip.

The effective angle of attack α\alpha of the strip depends on the geometric pitch angle ϕ\phi and the hydrodynamic inflow angle:

tan(ϕ+α)=Va(1+a)ωr(1a)\tan(\phi + \alpha) = \frac{V_a (1 + a)}{\omega r (1 - a')}

where aa is the axial induction factor and aa' is the rotational induction factor. At the design point, a=u/Vaa = u / V_a from the actuator-disk relation above; BET adds the rotational correction through aa'.

Lift and drag of the blade section

The local lift and drag on the strip follow from two-dimensional aerofoil theory:

dL=12ρW2cCLdrdL = \frac{1}{2} \rho W^2 c \cdot C_L \, dr

dD=12ρW2cCDdrdD = \frac{1}{2} \rho W^2 c \cdot C_D \, dr

where W=Va2+(ωr)2W = \sqrt{V_a^2 + (\omega r)^2} is the resultant inflow speed, cc is the local chord length, CLC_L is the lift coefficient, and CDC_D is the drag coefficient. The lift-to-drag ratio CL/CDC_L / C_D at the design pitch angle is the key section efficiency parameter: modern NACA 66-modified sections used in high-performance propellers have CL/CDC_L / C_D ratios of 50 to 80 at the design condition.

Resolving dLdL and dDdD in the axial (thrust) and tangential (torque) directions, then integrating over the blade radius from hub rhr_h to tip R=D/2R = D/2, and multiplying by the blade number ZZ:

T=ZrhR(dLcosϕdDsinϕ)drT = Z \int_{r_h}^{R} (dL \cos\phi - dD \sin\phi) \, dr

Q=ZrhRr(dLsinϕ+dDcosϕ)drQ = Z \int_{r_h}^{R} r (dL \sin\phi + dD \cos\phi) \, dr

Pure BET without any wake-induction correction overestimates thrust by 10 to 20% for typical merchant propellers, because it ignores the 3D effect of the trailing vortex sheet shed from each blade on the inflow seen by neighbouring strips. The correction is the principal motivation for lifting-line theory.

BET-momentum combination

The standard engineering approach (the combined BET-momentum or BEMT method) couples the momentum equations with the blade-element aerodynamics in an iterative loop: assume initial values of aa and aa', compute the local section loads from BET, equate those loads to the momentum theory predictions for the annular disk element, and iterate to convergence. BEMT remains the workhorse for rapid parametric studies and for propeller matching in commercial software such as NavCad (HydroComp) because it runs in fractions of a second per operating point.

Lifting-line theory

The vortex model

The trailing vortex sheet shed from a finite-span wing is the physical source of induced drag. For a propeller, each blade sheds a helical trailing vortex from every spanwise location where the bound circulation Γ(r)\Gamma(r) changes. Goldstein (1929) solved this problem exactly for a propeller with lightly loaded blades, obtaining the Goldstein circulation function κ(r)\kappa(r) that accounts for the mutual induction of the helical wake filaments.

The ideal bound circulation distribution that minimises the induced power loss for a given thrust was derived by Betz (1919) from an analogy with Prandtl’s minimum-induced-drag condition for a finite wing. The result is that the optimum propeller has a helical trailing vortex sheet that moves aft as a rigid helix (the Betz condition). Goldstein’s function κ(r)\kappa(r) is the correction to that ideal for a finite number of blades.

Lerbs (1952) extended lifting-line theory to moderately loaded propellers with an arbitrary prescribed circulation distribution, allowing the designer to specify the spanwise loading shape and compute the corresponding blade geometry. The Lerbs formulation remained the standard for early-design propeller calculation from the 1950s through the 1980s, and its physical structure underlies modern vortex-lattice propeller codes.

Practical lifting-line design workflow

In a lifting-line design run, the input is the required thrust TT, advance coefficient JJ, diameter DD, blade number ZZ, and the desired radial loading distribution. The output is the radial pitch distribution P(r)P(r) and the chord distribution c(r)c(r) that achieve the target performance. The steps are:

  1. Assume a circulation distribution Γ(r)\Gamma(r) (often starting with the ideal Betz distribution).
  2. Compute the induced velocities at the lifting line using Goldstein’s function or a numerical vortex-lattice integration.
  3. Compute the local hydrodynamic pitch angle from the total velocity (inflow plus induced).
  4. Check the section cavitation number and CLC_L against design limits; adjust the chord width if necessary.
  5. Recompute until the pitch and chord distributions converge.

Lifting-line theory gives accurate performance predictions (within 2 to 4% of model test) for conventional propellers with moderate pitch and no large spanwise variations in blade loading. It captures the tip-vortex-induced drag through the Goldstein function, which BET cannot do without a semi-empirical tip correction. The Prandtl tip-loss factor F=(2/π)arccos(exp(Z(Rr)/(2rsinϕ)))F = (2/\pi) \arccos(\exp(-Z(R-r)/(2r\sin\phi))) is a practical approximation to the Goldstein function used in BEMT codes.

The open-water diagram: J, KT, KQ, and eta_O

Non-dimensional parameters

The three primary non-dimensional parameters of propeller performance are:

Advance coefficient:

J=VanDJ = \frac{V_a}{n D}

where VaV_a is the speed of advance (m/s), nn is shaft speed (rev/s), DD is diameter (m). JJ is the kinematic similarity parameter: two geometrically similar propellers operating at the same JJ are in dynamically similar operating conditions, regardless of their absolute size or speed.

Thrust coefficient:

KT=Tρn2D4K_T = \frac{T}{\rho n^2 D^4}

Torque coefficient:

KQ=Qρn2D5K_Q = \frac{Q}{\rho n^2 D^5}

Open-water efficiency:

ηO=TVa2πnQ=J2πKTKQ\eta_O = \frac{T V_a}{2\pi n Q} = \frac{J}{2\pi} \cdot \frac{K_T}{K_Q}

These four quantities, plotted against JJ, form the open-water characteristic curves of the propeller. They are determined by model tests conducted to the ITTC Recommended Procedure 7.5-02-03-02: the propeller model (diameter typically 200 to 350 mm) is run in a towing tank or cavitation tunnel at a series of advance coefficients with the inflow velocity and shaft speed independently varied. The result is the KT(J)K_T(J), KQ(J)K_Q(J), ηO(J)\eta_O(J) diagram that is the central design and analysis tool.

Physical interpretation of the open-water diagram

At J=0J = 0 (zero advance, as in the bollard-pull condition), both KTK_T and KQK_Q are at their maximum; the propeller produces maximum thrust but no useful work and ηO=0\eta_O = 0. As JJ increases (higher ship speed relative to shaft speed), KTK_T and KQK_Q both fall; ηO\eta_O first rises to a maximum (the design JJ) then falls as the blade sections approach zero angle of attack. The propeller becomes unloaded and finally produces zero thrust at J0J_0 (the “run-away” point). The design point is chosen so that the maximum-efficiency JJ coincides with the operating advance coefficient at the service condition.

The propeller open-water efficiency calculator evaluates ηO\eta_O for any JJ, KTK_T, and KQK_Q values. The advance coefficient calculator converts between JJ, shaft speed, speed of advance, and diameter.

Wageningen B-series polynomial regression

The B-series is the most widely used systematic propeller series for merchant vessels. Systematic model tests at MARIN from the 1930s to 1960s covered blade numbers Z=2Z = 2 to 77, pitch-diameter ratios P/D=0.5P/D = 0.5 to 1.41.4, and expanded blade-area ratios AE/A0=0.30A_E/A_0 = 0.30 to 1.051.05. Oosterveld and Oossanen (1975) fitted these data to polynomial regressions of the form:

KT=iCT,iJsi(P/D)ti(AE/A0)uiZviK_T = \sum_{i} C_{T,i} \cdot J^{s_i} \cdot (P/D)^{t_i} \cdot (A_E/A_0)^{u_i} \cdot Z^{v_i}KQ=iCQ,iJsi(P/D)ti(AE/A0)uiZviK_Q = \sum_{i} C_{Q,i} \cdot J^{s_i} \cdot (P/D)^{t_i} \cdot (A_E/A_0)^{u_i} \cdot Z^{v_i}

The 1975 regression uses 47 terms for KTK_T and 39 terms for KQK_Q, fitted to model test data corrected to the ITTC 1978 friction line (Rn=2×106R_n = 2 \times 10^6). The regression is accurate to within 1 to 2% of model test for Z=3Z = 3 to 66, P/D=0.5P/D = 0.5 to 1.41.4, AE/A0=0.30A_E/A_0 = 0.30 to 1.051.05, and J=0J = 0 to J0J_0.

The Wageningen B-series calculator implements the full Oosterveld-Oossanen regression, returning KTK_T, KQK_Q, and ηO\eta_O for any point in the parameter space. The KT calculator and KQ calculator isolate the individual coefficient evaluations; the optimal pitch-diameter calculator finds the P/DP/D that maximises ηO\eta_O for a specified JJ or thrust requirement.

The principal limitation of the B-series is that it represents a specific family of propeller geometries from mid-twentieth-century design practice. Modern propellers with high skew, tip-fin modifications, or non-NACA section profiles may deviate from B-series predictions by 3 to 8% in efficiency, requiring proprietary regression data or panel-method analysis.

Lifting-surface theory and modern CFD

From lifting-line to lifting-surface

Lifting-line theory collapses the blade to a single concentrated vortex line. This is adequate for straight, moderately loaded blades but breaks down for highly skewed propellers, where the 3D curvature of the vortex sheet and the spanwise chord-loading distribution matter. Lifting-surface theory distributes bound and free vorticity over the actual blade surface and its helical wake sheet. The MIT codes PSF-2 and PUF-3 (Kerwin, 1978 onwards) established the academic standard; commercial implementations followed. The computational cost is higher than lifting-line by roughly one to two orders of magnitude, but the result is accurate blade-section pressure distributions, which are the direct input to cavitation analysis.

Viscous CFD

Modern propeller design uses RANS (Reynolds-averaged Navier-Stokes) CFD for final verification. The standard approach for a propeller in open water is a rotating-reference-frame simulation with the propeller blade geometry meshed at approximately 3 to 8 million cells. Steady MRF (multiple reference frame) or time-resolved sliding-mesh methods give KTK_T and KQK_Q to within 1 to 3% of model test for conventional propellers; for highly loaded or cavitating conditions, multiphase RANS with a cavitation model (Schnerr-Sauer or Kunz) is required.

CFD captures what no other method can at practical cost: tip-vortex formation, hub-vortex dynamics, propeller-rudder interaction in the behind-ship condition, and the fully unsteady loading in the non-uniform wake. A full open-water CFD campaign typically requires 10 to 30 operating points and 24 to 72 hours of compute time on a modern 64-core cluster; the cavitation tunnel model test it replaces costs 50,000to50,000 to 150,000 and takes 4 to 8 weeks.

The propulsive efficiency chain

Delivered-to-effective power

The chain from shaft power to effective towing power involves several efficiency ratios, each physically distinct. The overall quasi-propulsive coefficient (QPC or ηD\eta_D) relates the effective power PE=RTVP_E = R_T V (resistance times ship speed) to the delivered power PD=2πnQP_D = 2\pi n Q at the propeller:

ηD=PEPD=ηHηOηR\eta_D = \frac{P_E}{P_D} = \eta_H \cdot \eta_O \cdot \eta_R

where ηH\eta_H is the hull efficiency, ηO\eta_O is the propeller open-water efficiency, and ηR\eta_R is the relative-rotative efficiency. The shaft transmission efficiency ηS\eta_S (typically 0.97 to 0.99 for a conventional shafting arrangement with sterntube bearing) brings the chain from brake power PBP_B to PDP_D: PD=ηSPBP_D = \eta_S P_B. The quasi-propulsive coefficient calculator computes ηD\eta_D from the component efficiencies.

Wake fraction and hull efficiency

A propeller operating behind a ship hull does not see the free-stream ship speed VsV_s. The hull boundary layer reduces the mean axial velocity at the propeller plane. The Taylor wake fraction ww is defined so that:

Va=Vs(1w)V_a = V_s (1 - w)

For a single-screw full-form tanker or bulker, ww is in the range 0.25 to 0.40: the propeller plane sees only 60 to 75% of the ship speed as axial inflow. For a twin-screw container ship, ww is 0.10 to 0.18. Nominal wake (from a pitot-rake wake survey behind the hull model) and effective wake (accounting for propeller-hull interaction) differ by 5 to 15%.

The wake fraction also enters the design thrust calculation. The thrust the propeller must produce exceeds the bare-hull resistance because the propeller suction alters the stern pressure distribution, effectively increasing resistance. The thrust-deduction fraction tt captures this:

T(1t)=RTT (1 - t) = R_T

So T=RT/(1t)T = R_T / (1 - t). Typical values: t=0.12t = 0.12 to 0.220.22 for single-screw merchant ships, t=0.06t = 0.06 to 0.140.14 for twin-screw.

The hull efficiency combines these two:

ηH=1t1w\eta_H = \frac{1 - t}{1 - w}

When w>tw > t, as is typical for full-form single-screw ships, ηH>1.0\eta_H > 1.0. Values of 1.05 to 1.20 are standard for bulkers and tankers. This does not violate energy conservation: the hull boundary layer delivers pre-decelerated water to the propeller, which then re-accelerates it. The propeller does less work than it would in free stream because the water is already in relative motion with the hull. The hull efficiency calculator and wake fraction & thrust deduction calculator handle these parameters; the thrust deduction calculator and wake fraction calculator isolate each factor.

Relative-rotative efficiency

ηR\eta_R is the ratio of propeller torque in the actual non-uniform wake to the torque the same propeller would absorb in the uniform open-water test at the same thrust. The non-uniform wake means the blade sections see varying angles of attack as they rotate, which slightly alters the integrated torque at the same mean thrust. For most merchant propellers, ηR\eta_R lies in the range 0.97 to 1.02. Single-screw full-form ships with a strong wake non-uniformity tend toward 1.02 to 1.05; twin-screw ships with more uniform inflow tend toward 0.97 to 0.99.

Shaft efficiency and the full chain

The shaft efficiency ηS\eta_S covers friction in the sterntube bearing, intermediate shaft bearings, and thrust block (excluding the main engine). Modern self-lubricated composite bearings in sterntube arrangements reduce sterntube bearing losses; values of ηS=0.97\eta_S = 0.97 to 0.990.99 are typical. The full propulsion chain is then:

PB=RTVsηHηOηRηSP_B = \frac{R_T V_s}{\eta_H \cdot \eta_O \cdot \eta_R \cdot \eta_S}

For a Panamax bulker with ηH=1.12\eta_H = 1.12, ηO=0.68\eta_O = 0.68, ηR=1.02\eta_R = 1.02, ηS=0.98\eta_S = 0.98, the QPC is 1.12×0.68×1.02×0.98=0.7611.12 \times 0.68 \times 1.02 \times 0.98 = 0.761. So about 24% of brake power is lost between the engine output flange and the ship’s forward motion through water.

Comparison of propeller design theories

MethodFidelityTypical accuracy (KT)CPU costPrimary use
Actuator-disk momentum theoryLowN/A (no blade)NegligibleEfficiency bound, disc-loading scaling
BET-momentum (BEMT)Low-medium5-15%<1 sRapid parametric design, matching
Lifting-line (Lerbs-type)Medium2-5%1-10 sEarly-stage optimisation, systematic series
Wageningen B-series polynomialMedium1-2% (within series)<1 sCommercial design, all stages
Lifting-surface (panel method)High1-3%1-10 minDetailed blade design, cavitation analysis
RANS CFD (MRF, steady)High1-3%HoursVerification, behind-ship condition
RANS CFD (sliding mesh, unsteady)Very high0.5-2%DaysPressure pulses, cavitation, ESD interaction

The Wageningen B-series is anomalously efficient in terms of accuracy-per-CPU-second because it is a regression of real model test data rather than a theoretical approximation: within the tested parameter space it is a proxy for the model basin rather than a calculation.

Propeller geometry and its effect on performance

Pitch-diameter ratio

The pitch-diameter ratio P/DP/D is the primary performance-tuning parameter in the Wageningen B-series. Increasing P/DP/D at fixed JJ increases both KTK_T and KQK_Q; the maximum-efficiency JJ also shifts. For a given ship speed and shaft speed (fixed JJ), the P/DP/D that maximises ηO\eta_O can be found by differentiating the B-series polynomials and setting dηO/d(P/D)=0d\eta_O / d(P/D) = 0. Typical design P/DP/D values: 0.65 to 0.75 for slow-speed two-stroke direct-drive propellers (low shaft RPM, large diameter), 0.95 to 1.05 for medium-speed four-stroke geared propellers, and 1.10 to 1.25 for high-speed naval propellers.

Expanded blade-area ratio

Increasing AE/A0A_E/A_0 (expanded area ratio) at constant P/DP/D and ZZ reduces KTK_T and KQK_Q slightly (the blade sections operate at lower lift coefficient for the same total thrust) but, more importantly, reduces blade loading and delays cavitation inception. The Keller criterion gives a practical minimum AE/A0A_E/A_0 to avoid cavitation:

(AEA0)min=(1.3+0.3Z)T(p0pv)D2+k\left(\frac{A_E}{A_0}\right)_{\min} = \frac{(1.3 + 0.3 Z) T}{(p_0 - p_v) D^2} + k

where p0p_0 is the static pressure at the shaft centreline, pvp_v is vapour pressure, and k=0k = 0 for twin-screw ships and k=0.2k = 0.2 for single-screw ships. The Keller cavitation criterion calculator evaluates this limit. Exceeding AE/A0=0.80A_E/A_0 = 0.80 on a merchant propeller incurs manufacturing complexity and efficiency loss; above 1.0 the blades become so broad that wake interaction between them degrades performance.

Blade number and skew

Increasing blade number ZZ reduces the loading per blade and the amplitude of periodic thrust variation (important for ship vibration). Five-bladed propellers are now the norm for container ships; six-bladed propellers are used on high-speed ferries and warships. Odd blade numbers are preferred when even numbers would produce a harmonic coincidence with the engine firing frequency. Highly skewed designs (skew angle 25 to 45 degrees) reduce hull pressure-pulse amplitudes by distributing the blade’s passage through the wake deficit over a longer arc; the penalty is a more complex casting and higher root-bending moments under IACS UR M55.

Cavitation: physics, prediction, and design

The cavitation number

Cavitation occurs when the local static pressure on the blade suction surface falls below the water vapour pressure pvp_v (2,340 Pa at 20°C). The dimensional parameter governing inception is the cavitation number referred to the section at 0.7R:

σ0.7R=p0+ρghpv12ρW0.7R2\sigma_{0.7R} = \frac{p_0 + \rho g h - p_v}{\frac{1}{2} \rho W_{0.7R}^2}

where p0p_0 is atmospheric pressure (101,325 Pa), hh is the depth of the shaft centreline below the waterline, and W0.7RW_{0.7R} is the resultant inflow speed at 0.7R: W0.7R=Va2+(0.7πnD)2W_{0.7R} = \sqrt{V_a^2 + (0.7\pi n D)^2}. A higher σ\sigma means greater resistance to cavitation. Typical values for merchant propellers at design condition: σ0.7R=1.5\sigma_{0.7R} = 1.5 to 4.04.0. The cavitation inception calculator evaluates σ0.7R\sigma_{0.7R}.

Types of propeller cavitation

Sheet cavitation forms as a stable cavity on the suction face of the blade, attached at the leading edge and closing mid-chord. Moderate sheet cavitation (< 30% of blade span) is acceptable on most merchant propellers; it degrades performance but is not immediately erosive. When the sheet cavity closes unsteadily and sheds vortex structures, it becomes cloud cavitation, which is highly erosive: the collapse of bubble clouds generates local pressure spikes measured in hundreds of MPa.

Tip-vortex cavitation forms in the low-pressure core of the trailing vortex shed from the blade tip. It is the first type to appear as speed increases and is visible as a continuous helical thread of white bubbles. By itself it is relatively benign, but it broadens into an area of intense bubble collapse when it impinges on the downstream hull, rudder, or energy-saving devices.

Face cavitation occurs on the pressure face during off-design operation at high-pitch or very low ship speed (backing manoeuvres). It is typically less of a design concern than suction-face sheet cavitation.

Root cavitation forms near the blade root where local flow velocities are low but the geometry creates a suction region. It interacts with the hub vortex and can cause steady erosion of the blade root fillet.

The Burrill chart

Burrill and Emerson (1963) presented an empirical chart relating the suction-side thrust-loading coefficient:

τc=T12ρW0.7R2AP\tau_c = \frac{T}{\frac{1}{2} \rho W_{0.7R}^2 \cdot A_P}

where APA_P is the projected blade area, to the local cavitation number σ0.7R\sigma_{0.7R}. Contour lines on the chart indicate the percentage of blade back area in cavitation at inception. For merchant ship single-screw propellers, the design point should fall at or below the 2.5% back cavitation contour. For twin-screw warships, the 10% contour may be acceptable. The chart converts directly into a minimum AE/A0A_E/A_0 requirement for the given thrust and speed, and it is the first check after a preliminary B-series selection. The Burrill chart calculator automates this check.

Controlling cavitation through design

The four principal design levers are: (1) increasing DD to reduce the disc loading T/AT/A and lower CTC_T; (2) increasing AE/A0A_E/A_0 to reduce blade-section lift coefficients; (3) increasing ZZ to spread the thrust over more blades; and (4) optimising the radial loading distribution (via pitch or circulation) to reduce the peak suction-surface pressure coefficient. From a hydrodynamic standpoint, a larger, slower propeller is always better for cavitation; the constraint is the hull clearance and the engine power-speed curve. The ship resistance and powering article covers how hull resistance determines the required thrust and hence the disc loading.

Fixed-pitch versus controllable-pitch propellers

The design choice between a fixed-pitch propeller (FPP) and a controllable-pitch propeller (CPP) is fundamentally about operational flexibility versus mechanical simplicity.

CharacteristicFixed-pitch (FPP)Controllable-pitch (CPP)
Efficiency at design pointBaselineWithin 1-2% of FPP
Efficiency at off-designDegraded (fixed P/DP/D)Maintained (adjust P/DP/D)
Reversing manoeuvreEngine reversal (slow)Pitch reversal (fast)
Hub complexitySimple, solid hubHydraulic hub with oil passages
Hub diameter ratio0.16 to 0.220.24 to 0.32
MaintenanceLow; no moving parts in hubOil seals, hydraulic servo require periodic maintenance
Shaft speed controlVariable via governorCan run at constant shaft speed
Typical applicationVLCC, bulk carrier, container shipFerry, naval vessel, offshore supply, cruise ship
Cavitation behaviourOptimised for one design P/DP/DCan unload blade to reduce cavitation in shallow water

For most deep-sea merchant ships with slow-speed two-stroke engines, the FPP is standard: the engine reverses for manoeuvring, the design speed is narrow, and hub simplicity reduces maintenance. CPPs are standard where fast astern power is needed without engine reversal (ferries, RoRo) or where the vessel operates across a wide speed range (offshore vessels on DP, frigates varying from 5 to 30 knots). The CPP adds roughly 4 to 6% to propeller procurement cost and requires an oil distribution box (OD box) and hydraulic power unit.

Modern high-efficiency propeller designs

Kappel tip-fin propeller

The Kappel propeller replaces the standard blade tip with a smooth winglet-like curve that bends smoothly toward the suction side of the blade. This reduces the intensity of the tip vortex (and hence tip-vortex-induced drag) without the blade-tip loading that would accelerate cavitation. Full-scale sea trial measurements on vessels fitted with Kappel propellers in place of conventional propellers of the same diameter have shown propulsive power reductions of 3 to 6% at the same thrust, corroborated by MARIN towing tank tests. MAN Energy Solutions has held the commercial licence for Kappel propeller manufacture since 2007.

Ducted propellers and the Kort nozzle

A Kort nozzle is an accelerating duct (internal diameter narrowing from inlet to outlet) surrounding the propeller. The duct accelerates the inflow, increasing the effective mass flow through the disc. For highly loaded propellers (CT>2C_T > 2), the added thrust from the duct more than compensates for the duct’s own drag, improving total bollard-pull efficiency by 20 to 30% over the open propeller. The efficiency advantage diminishes as speed increases: above about 14 knots for typical merchant applications, the duct drag exceeds its thrust augmentation and an open propeller is more efficient. Kort nozzles are standard on oceangoing tugs, push-barges, and trawlers operating at high bollard-pull fractions; they are rare on conventional merchant ships at design speeds above 12 knots.

Energy-saving devices at the propeller

Several energy-saving devices (ESDs) target the propeller-local flow:

  • Propeller boss cap fins (PBCF): small radially arranged fins on the aft face of the propeller boss cap that disrupt the hub vortex and recover 1 to 2% of delivered power.
  • Pre-swirl stators: fixed fins ahead of the propeller that impart a counter-rotation to the inflow, reducing the rotational kinetic energy deposited in the propeller race by 1 to 4%.
  • Mewis duct: a combination of a partial duct and pre-swirl fins, targeting both the accelerating-duct thrust augmentation and the pre-swirl recovery; claimed gains of 3 to 8% in the behind-ship condition.
  • Contra-rotating propellers (CRP): the forward propeller delivers a swirling race; the aft propeller, rotating in the opposite direction, converts that swirl to additional thrust. Theoretical recovery efficiency approaches 10 to 15%, limited in practice by the mechanical complexity of the coaxial shaft arrangement.

The energy-saving devices article covers these in detail. The bulbous bow retrofits article covers hull-side changes that alter the wake field seen by the propeller.

Propeller-engine matching and the load diagram

The propeller curve and engine envelope

A propeller operating in a given condition absorbs a torque that, at constant P/DP/D and resistance, varies as Qn2Q \propto n^2 (since KQK_Q is roughly constant near the design JJ). Expressed in delivered power: PD=2πnQn3P_D = 2\pi n Q \propto n^3. This cubic relationship is the propeller curve. The main engine must develop at least this power at the corresponding RPM. The engine manufacturer supplies a load diagram (also called the power envelope) bounded by the maximum continuous rating (MCR) at 100% RPM and torque limits, the minimum combustion stability limit, and the speed limits for critical torsional resonance avoidance.

Matching is correct when the propeller curve at the clean-hull, deep-sea design condition intersects the engine load diagram at or near the “M-point” (MCR with a 5 to 15% engine margin). Over the vessel’s life, resistance increases due to hull fouling, trim changes, and eventually dry-docking roughness, shifting the propeller curve to higher torque at the same RPM. The propeller is therefore typically designed 3 to 5% “light” (lower pitch or larger diameter than the exact optimum) at delivery, to run heavy as the hull ages.

Off-design operation

Two conditions regularly push the propeller away from its design point:

Heavy weather: wave-added resistance of 20 to 50% of calm-water resistance in significant wave heights of 4 to 6 m is common on North Atlantic routes. The engine load increases at the same shaft speed; if the propeller curve already runs close to MCR, the master slows down to stay within the torque limit. This is the normal “sea margin” scenario.

Ballast condition: a VLCC in ballast carries roughly one-quarter of the loaded displacement, cutting resistance by 40 to 55%. At the same shaft speed, the propeller is severely underloaded and operates at very low KTK_T. Shaft speed typically increases toward the engine maximum-RPM limit. CPP vessels manage this by reducing pitch; FPP vessels accept the reduced fuel efficiency in ballast.

The interaction with slow steaming is direct: a propeller designed for 15 knots at MCR operating at 12 knots (typical slow-steaming condition) sees a JJ reduction of roughly 20% if shaft speed is also reduced proportionally. Provided the engine output curve allows operation in that region, the propeller can remain close to its maximum-efficiency JJ by reducing both speed and RPM simultaneously; this is the basis of the “optimal slow steaming” condition described in the ship resistance and powering and slow steaming articles. CII compliance depends on exactly this matching: the what-is-cii article discusses the regulatory framework.

The ITTC 1978 performance prediction method

The standard method for predicting ship performance from model tests is the ITTC 1978 method (ITTC Procedure 7.5-02-03-01.4). It extrapolates propeller open-water characteristics from model to full scale by correcting ΔKT\Delta K_T and ΔKQ\Delta K_Q for the Reynolds number difference between model tests (Re typically 2×1062 \times 10^6 at 0.7R) and full-scale operation (Re typically 7×1087 \times 10^8 to 2×1092 \times 10^9 at 0.7R). The correction terms are:

ΔKT=0.3ΔCDc/D(P/D)\Delta K_T = -0.3 \Delta C_D \cdot c / D \cdot (P/D)

ΔKQ=+0.25ΔCDc/D\Delta K_Q = +0.25 \Delta C_D \cdot c / D

where ΔCD\Delta C_D is the change in section drag coefficient from model to ship scale, computed from the ITTC 1957 skin-friction line applied to the blade section at 0.75R with an appropriate roughness allowance (ISO 484 class I finish: mean roughness Ra1.6R_a \leq 1.6 µm; class S: Ra0.8R_a \leq 0.8 µm). A higher surface roughness (from fouling, erosion, or poor polish) reduces ηO\eta_O by 0.5 to 2%, which directly maps to increased fuel consumption. The trim optimisation article covers how trim interacts with the effective wake and hence the propulsion efficiency.

Limitations

Several assumptions and practical boundaries constrain the predictive accuracy of propeller theory at full scale:

Wake non-uniformity: all performance prediction methods, including the ITTC 1978 method, use a mean effective wake fraction derived from the self-propulsion test. The actual propeller plane wake is strongly non-uniform (particularly behind single-screw ships), with local axial velocity variations of ±\pm40% relative to the mean. The blade sections pass through this varying field once per revolution, generating unsteady thrust and torque fluctuations that lifting-line and even CFD RANS steady methods do not fully capture.

Scale effects in cavitation: cavitation tunnel tests at model scale operate at Reynolds numbers 100 to 1,000 times lower than full scale. The nuclei content of the tunnel water affects inception; dissolved-gas content affects bubble dynamics. Nuclei seeding at measured full-scale values is standard practice in modern cavitation tunnel tests but is imperfectly controlled, so inception cavitation-number predictions carry ±\pm10 to 20% uncertainty.

B-series applicability limits: the Wageningen B-series regression is valid for Z=2Z = 2 to 77, P/D=0.5P/D = 0.5 to 1.41.4, AE/A0=0.30A_E/A_0 = 0.30 to 1.051.05, and J=0J = 0 to J0J_0. Outside these bounds the polynomial extrapolates, and errors grow rapidly. Kappel tip-fin propellers, contra-rotating configurations, and ducted propellers all lie outside the B-series parameter space; dedicated regression data or higher-fidelity methods are required.

Viscous interaction with the hull: the propeller affects the pressure distribution on the hull stern, which in turn changes the effective wake fraction. This propulsion-hull interaction is captured in the self-propulsion test thrust-deduction factor tt, but tt is load-dependent. As the propeller loading changes with speed or fouling condition, tt changes, and the simple constant-tt model used in most performance predictions introduces an error of 1 to 3% in the predicted effective power at off-design conditions.

Blade-root and hub-vortex losses: lifting-line theory and most panel codes treat the hub as a solid cylinder and do not compute hub-vortex losses. In reality, the hub vortex extracts 1 to 2% of delivered power on a conventional hub; PBCF energy-saving devices target exactly this source.

Material deformation and blade-gap effects for CPP: controllable-pitch propeller blades flex slightly under load, changing their effective pitch by 0.5 to 2% relative to the set geometric pitch. This load-dependent deflection is not captured in rigid-blade propeller theory and requires structural finite-element analysis coupled to the hydrodynamic solver.

See also

Calculators used in propeller analysis:

Related wiki articles:

Frequently asked questions

What is the advance coefficient J in propeller theory?
J = Va / (n D), where Va is the speed of advance at the propeller plane (m/s), n is shaft speed in revolutions per second, and D is propeller diameter in metres. J is the primary similarity parameter for propeller performance: all propellers that are geometrically similar and operate at the same J produce identical non-dimensional thrust and torque coefficients KT and KQ.
What is the ideal efficiency of an actuator disk?
The Rankine-Froude momentum theory gives the ideal efficiency as eta_ideal = 2 / (1 + sqrt(1 + CT)), where CT = T / (0.5 rho Va^2 A) is the thrust loading coefficient. For a merchant ship propeller with CT in the range 0.5 to 1.0, the ideal efficiency lies between 0.85 and 0.90. Real propellers reach roughly 75 to 85 percent of this bound because of viscous, tip-vortex, and finite-blade losses.
What is hull efficiency and why can it exceed unity?
Hull efficiency eta_H = (1 - t) / (1 - w), where t is the thrust-deduction fraction and w is the Taylor wake fraction. For typical single-screw merchant ships, w > t, so the denominator is smaller than the numerator and eta_H > 1.0 (commonly 1.05 to 1.20). This does not violate thermodynamics: the hull boundary layer delivers pre-accelerated water to the propeller at a lower energy cost than the propeller would pay in open water, effectively subsidising propulsive work.
What is the Burrill chart used for?
The Burrill chart (1943) plots the suction-side thrust-loading coefficient tau_c against the cavitation number sigma_0.7R at the 0.7R blade radius, with contour lines of allowable percentage back-cavitation (typically 2.5 to 10 percent for merchant ships). It is the standard practical tool for selecting the minimum expanded blade-area ratio AE/A0 that avoids erosive cavitation at the design thrust and speed.
How does the Wageningen B-series predict propeller performance?
The B-series regression (Oosterveld and Oossanen, 1975) expresses KT and KQ as polynomial functions of J, P/D, AE/A0, and blade number Z, fitted to systematic open-water model tests at MARIN. Evaluating the polynomials at any operating point gives KT and KQ directly, from which open-water efficiency eta_0 = (J / 2pi) * (KT / KQ) follows. The calculator at /calculators/tech-propeller-wageningen covers the full B-series parameter space.