MN Ch. 1 — Ship Resistance and Flow: Introduction

1. Overview — why resistance prediction matters

Chapter 1 sets the stage for the whole book. It answers three questions: why a naval architect must predict resistance accurately, how resistance can be predicted (three families of methods), and how the book is organised. The central idea is that the resistance of a ship — together with the propulsive efficiency — fixes the engine power needed to reach a given speed, and that the resistance itself can only be understood by studying the flow of water around the hull.

Reading guide. This is an introductory, mostly historical chapter. No equations are derived here (they begin in Section 2). The one quantitative statement is Froude's law of comparison, given below as a relation in words; the maths that justify it come in Section 3 (Similarity).

Fonte: LARSSON, Lars; RAVEN, Hoyte C. Ship Resistance and Flow. The Principles of Naval Architecture Series. J. R. Paulling (ed.). Jersey City: SNAME, 2010. Chapter 1 — Introduction.

Edital: Anexo 2-A, Área I, itens 1 e 5 (Resistência do Navio e métodos para sua obtenção) · Anexo 2-B, Área I, item 3 (Larsson & Raven), Chapter 1.

2. The importance of accurate resistance predictions (§1.1)

A central problem for the practicing naval architect is predicting the resistance of a new design at an early stage of the project — before the hull exists. When a ship is ordered, the owner and the shipyard sign a contract that specifies, among the strictest clauses, the contract speed: the speed reached at a stated power consumption during a trial run before delivery.

That trial is meant to take place under ideal conditions — no wind, no seaway, no influence from restricted water or currents. In reality, corrections almost always have to be applied for those factors. If the corrected speed falls below the contract speed, the yard pays the owner a penalty depending on the difference between the achieved speed and the contract speed; if the gap is too large, the owner may even refuse the ship.

The designer's dilemma. Two opposite pressures squeeze the prediction:
  • Compete or lose the order. Because shipyards compete fiercely on the global market, the offer must be at least as good as the rivals'. A few percent more power for a given speed can cost the contract.
  • Over-optimism is expensive. If the prediction is too optimistic and the ship misses the specification, the penalty falls on the yard.

The power needed to drive the ship at a given speed depends on more than resistance. The propulsive efficiency — how well the propeller performs and interacts with the hull — and the losses in the power train both matter. Even so, resistance is the single most important factor determining the required power.

Resistance, like the other forces on the hull, results from shear and normal stresses (pressures) exerted on the hull surface by the water flow. Knowing the flow around the ship is therefore essential both to understand the individual resistance components and to design the hull well. The flow around the stern also fixes the operating conditions of the propeller — which is why the book puts heavy emphasis on describing the hull flow.

As in every design project, conflicting demands compete. The hydrodynamic qualities are only one set among many; they also include seakeeping and manoeuvring, which — together with propulsive efficiency — are covered in other volumes of the Principles of Naval Architecture series, not here.

3. Different ways to predict resistance (§1.2)

There are three families of methods for predicting resistance: model testing, empirical methods, and computational techniques. The next three sub-sections trace the history of each; the fourth explains which method is used at which stage of design.

3.1 Model testing (§1.2.1)

Because ship resistance is complicated, naval architects turned to experiments very early. Leonardo da Vinci (1452–1519) is recorded testing three models with different fore-and-aft displacement distributions. Samuel Fortrey (1622–1681) towed small wooden models in a tank with falling weights. Colonel Beaufoy, under the Society for the Improvement of Naval Architecture (founded London, 1791), ran between 9,000 and 10,000 towing experiments from 1791 to 1798 in the Greenland Dock. Fredrik af Chapman in Sweden published resistance tests on simple shapes in 1795, and Benjamin Franklin (1764) was probably the first American to test models — to confirm that resistance rises in shallow water.

The big obstacle was scaling: how to extrapolate the model's towing force to full scale, and at what speed to tow the model to represent a given full-scale speed. The French scientist Ferdinand Reech (1844) first solved the problem but never pursued it, so the solution is credited to the Englishman William Froude, who stated his law of comparison in 1868.

Froude's law of comparison (1868). "The (residuary) resistance of geometrically similar ships is in the ratio of the cube of their linear dimensions if their speeds are in the ratio of the square roots of their linear dimensions."

Put as relations between model (subscript m) and ship (subscript s) with linear scale ratio $\lambda = L_s/L_m$, the corresponding-speed condition and the residuary-resistance ratio are:

$$ \frac{V_s}{\sqrt{L_s}} = \frac{V_m}{\sqrt{L_m}} \qquad\Longrightarrow\qquad \frac{R_{r,s}}{R_{r,m}} = \lambda^3 $$

The same idea, written as a single speed parameter, is that the ratio $V/\sqrt{L}$ must be equal at both scales — the origin of the speed–length ratio (which is itself not dimensionless; the dimensionless form is the Froude number $V/\sqrt{gL}$). The residuary resistance is the total resistance minus that of an equivalent flat plate (a rectangular plank of the same area, length and speed as the hull). Froude's plan: split total resistance into friction (measured from planks, which make no waves, at both scales) and residuary resistance from the waves (found by subtracting friction from the model test, then scaled by displacement and added to full-scale plank friction).

Froude ran his first model experiments in 1863 in a rainwater tank with a falling weight, then persuaded the British Admiralty (1868) to fund a proper tank, completed near Torquay in 1871: 85 m long, 11 m wide at the surface, 3 m deep, with a mechanically propelled towing carriage — the forerunner of today's towing tanks.

The friction problem and the ITTC lines. Froude's weak point was his friction formula. Correct friction scaling waited for Reynolds (1883), who found the scaling parameter is a dimensionless number — later the Reynolds number. Schoenherr (1932) introduced it into model testing with a plank-friction formula, but only in 1957 did the ITTC recommend Reynolds-based friction scaling, through a different formula — the "ITTC-57" friction line. The modified scheme (ITTC-57 line replacing Froude's friction) is still called "Froude scaling."

By the early 1960s, very high block-coefficient ships made a finer split necessary: not all viscous effects fit into plank friction, so a second viscous component — tied to the three-dimensional (3D) shape of the hull — had to be added. This "3D extrapolation" technique was proposed for general use by the ITTC in 1978, giving the "ITTC-78" procedure, used by most tanks today for normal displacement hulls.

3.2 Empirical methods (§1.2.2)

Model tests are time-consuming, especially when many alternative designs must be screened very early. Hence the need for fast, if less accurate, estimates. These come in two types: systematic series and statistical formulas from unsystematic data.

Milestones of empirical resistance prediction (§1.2.2).
Method / seriesAuthor(s)TypeKey point
Taylor standard seriesAdm. Taylor (1933); reanalysed by Gertler (1954)Systematic seriesFirst comprehensive systematic tests (early 1900s, Washington EMB); based on the old ship Leviathan (designed 1900); Gertler re-corrected the data using Schoenherr friction
Series 60SNAME + ATTC, from 1948; reported by Todd (1963)Systematic seriesSeveral parent models (one per block coefficient) → realistic hulls for all variations; also gives self-propelled (delivered-power) results
Delft seriesKeuning & Sonnenberg (1998)Systematic seriesSailing-yacht models since the mid-1970s; > 50 models, continuously extended; rare modern systematic effort
Doust & O'Brien regressionDoust & O'Brien (1959)Statistical formulaFirst statistical formula on unsystematic data (150 fishing vessels); polynomial in six shape parameters — but no physics, just a polynomial
Holtrop–Mennen methodHoltrop & Mennen (1978)Statistical formulaUses a theoretical wave-resistance expression (two travelling pressure disturbances: bow and stern); coefficients fitted to 334 hulls; the most widely used rapid-estimate method today

The crucial contrast is in the last column: Doust & O'Brien's formula involves no physics — it is merely a polynomial in the tested parameters — whereas Holtrop & Mennen built a physical wave-resistance expression into the regression and split resistance by the 3D extrapolation procedure. That physical grounding is why Holtrop–Mennen remains the standard rapid method.

3.3 Computational techniques (§1.2.3)

Thanks to fast computer development over the past ~50 years, computational ship hydrodynamics matured over a shorter span than the experimental methods. The history runs along two tracks — first inviscid (potential-flow) methods, then viscous (boundary-layer → RANS) methods.

Inviscid track. The first method that counts as computational hydrodynamics came in a landmark paper by the Australian mathematician Michell (1898), over a century ago. Neglecting viscosity, he derived the inviscid flow around a "slender," narrow-beam ship in a uniform stream and integrated the fore-and-aft pressure components to get total wave resistance. To make it solvable he linearised the boundary conditions — applying the hull condition to the centreplane (so the result held for a vanishingly thin ship) and the free-surface condition to the flat undisturbed surface. Havelock and coworkers (early 20th c.) instead measured wave-making resistance from the energy of the wave system, and introduced sources and sinks distributed on the centreplane (positive forward, negative aft, proportional to local waterline angle). Inui and coworkers (Japan, 1960s–70s) advanced hull optimisation for wave resistance; in the same period, methods for experimentally determining wave resistance from wave cuts near the hull were developed — the landmark paper being Eggers, Sharma & Ward (1967).

From analytical to numerical panel methods. With computers (1960s), numerical methods began. Hess & Smith (1962) at Douglas Aircraft built an inviscid panel method — body surface discretised by flat quadrilateral panels — applicable to arbitrary 3D bodies, useful at once in aerodynamics. The decisive step for ships was Dawson (1977), who added the free surface, linearising the free-surface condition about a "double-model" solution (surface treated as a plane of symmetry). Because for a bluff hull the double-model flow is much closer to reality than Michell's linearisation about the undisturbed flow, Dawson's approach is less approximate, and — like all panel methods — it satisfies the exact hull boundary condition on the actual hull surface. Its drawback (exact on the hull, linearised on the free surface) was later removed: Larsson, Kim & Zhang (1989), refined and validated by Janson (1997), with parallel work by Jensen (Germany, 1988) and Raven (Holland, 1996).

Panel methods are now used extensively, but the assumption of zero viscosity is an inherent weakness, so agreement with measurement is unlikely to improve much further — to go beyond, viscosity must be included.

Viscous track. Computational viscous-flow methods also appeared with the computer. The 1960s brought several 2D boundary-layer methods; 3D boundary-layer research ran through most of the 1970s. A 1980 Gothenburg workshop concluded that the boundary layer was well computed over the forward and middle hull, but the stern flow could not be predicted at all by the boundary-layer approximation — that needed RANS (Reynolds-Averaged Navier-Stokes) methods. RANS had just been applied to ship flows by Spalding's group at Imperial College (1978).

The RANS workshops — progress of viscous CFD (§1.2.3).
WorkshopYearState of the art
1st Gothenburg1980Boundary-layer methods OK forward/midship; stern flow not predictable
2nd Gothenburg199017 of 19 methods now RANS; much better stern flow; but wake contours in the propeller plane too smooth (under-predicted vortex)
Tokyo199410 methods with free-surface capability; but compute power too small for fine free-surface resolution — potential-flow panels still gave better waves
4th (Gothenburg)2000Big accuracy gains in wake and waves; flow details still imperfect; turbulence modelling remains the hard limit

Better resolution of RANS solutions can keep improving accuracy, but the turbulence-modelling problem cannot be avoided. To beat it, far more expensive methods — Large Eddy Simulation (LES) or Direct Numerical Simulation (DNS) — are required, which demands very large increases in computer power.

3.4 Use of the methods (§1.2.4)

The three methods serve different stages of design:

Which method, which design stage (§1.2.4).
StageMethodWhy
Early basic design — explore the design space (length, beam, draft, block coefficient, LCB)EmpiricalTime is short; a fast, reasonable estimate is enough. Most CAD ship-design packages embed a module — usually Holtrop–Mennen
Forebody optimisation — refine main dimensions and local shapeNumerical (potential flow)Now standard; tune bulb size/shape and fore-shoulder radius, mainly to minimise wave resistance
Afterbody optimisation — stern sections (V/U/bulb), bilge radius, rudder effectNumerical (viscous / RANS)Boundary layer is large at the stern; boundary-layer theory too crude for the wake → RANS needed; aim is to minimise delivered power, not just resistance, so hull–propeller interaction must be estimated
Final hull decision — accurate resistance & powerModel testingStill used for most new ships; numerical results not yet as reliable, so the few best candidates are model-tested before the final decision

Optimisation may be done manually (systematic shape variation) or formally, linking a CFD code to a hull-shape program (often a CAD tool). With constraints, single-objective optimisation (e.g. minimise delivered power) or multi-objective optimisation (also pressure fluctuations, or different qualities like seakeeping) is possible.

The trajectory. The author expects regular model testing to be replaced by numerical prediction sooner or later. Tanks and other facilities would then serve mainly for advanced investigations and for validating new computational techniques — but numerical predictions have not yet reached model-test reliability.

4. The structure of this book (§1.3)

The volume aims to provide three things: a basic understanding of the resistance problem for ships and other marine vehicles; insight into the three prediction methods; and practical design guidelines.

The first objective is covered by the next six sections, then the prediction methods, then the guidelines:

Roadmap of the book (§1.3).
SectionTopic
2 Governing EquationsEquations governing the flows of interest, with boundary conditions
3 SimilarityUses Section 2's equations to prove the similarity laws for extrapolating model data to full scale
4 Decomposition of ResistanceSplits the total resistance of four widely different ships into components, briefly described
5 Wave resistanceThe wave-resistance component in detail
6 Viscous resistanceThe viscous-resistance component in detail
7 "Other components"The remaining resistance components in detail
8 Experimental techniquesFirst of the three prediction techniques
9 Numerical methodsSecond prediction technique
10 Empirical predictionsThird prediction technique
11 Hull design guidelinesPractical guidelines for designing a ship with good resistance properties
How Sections 5–7 map to the three components. After Section 4 divides resistance into components, Sections 5, 6 and 7 discuss them in detail and in order: wave resistance, viscous resistance, and "other components," respectively. Then Sections 8–10 cover the three prediction techniques (experimental, numerical, empirical), and Section 11 closes with design guidelines.