John Peach

Mind-sized STEM ideas and experiments, beyond the textbook.

Of-sailing-ships-hero.png

Home

keywords: hull design, Froude number, catamaran design, naval architecture, locomotion

In 18 hundred and 59, the engineer Brunel, Would build the greatest ship afloat, and rule the ocean's swell. Nineteen thousand tons of steel they used to shape the mighty keel, Forged inside the smelter where they made the gates of Hell... And the name upon the contract, Isambard Brunel.

—Ballad of the Great Eastern, by Sting

How are boats, velociraptors, and walking on the Moon related? We'll find out by designing a catamaran and making a mistake. Discover how Froude numbers unite naval engineering, paleontology, and space exploration. Learn how the same physics governs ship design, dinosaur movement, and astronaut locomotion.

The Great Eastern

Isambard Kingdom Brunel an English civil engineer described as "one of the most ingenious and prolific figures in engineering history" designed the SS Great Eastern which was launched in 1858. The Great Eastern was by far the largest ship built to date at 211 m in length and a displacement of 32,160 tons. Four steam engines with a combined power of 6 MW powered the Great Eastern.

great-eastern-1866.png

Brunel wanted a ship that could circumnavigate the globe without refueling and figured that it could carry a larger cargo with a smaller crew. But in the end, the Great Eastern couldn't maintain a speed much more than about 14 knots (1 knot = 1 nautical mile/hour = 1.15 mph) and wasn't the financial success Brunel hoped.

Engineers at the time didn't fully understand friction forces. Now, we can better predict the action of waves on a boat and the drag forces associated with the motion of a boat through the water. We'll use these equations to design a yacht, but in the process, we'll make a mistake that leads us back to The Great Eastern.

Design Parameters

In this post, I'll show you how to use the open-source program SMath Studio to write the equations required for designing a catamaran. The calculations below are from the Catamaransite which were derived from a design by Naval Architect Terho Halme.

Designing a sailing catamaran starts by selecting a few parameters to set the overall dimensions. Once you have the basic dimensions figured out you will know how much power it takes to drive the boat at different speeds, and this helps determine the size of the sails and the power requirements of the engine.

A consistent set of units helps the designer keep track of the parameters. If the calculations show your expected boat speed in meters \sqrt{\text{meters}} meters then obviously there's a problem somewhere with the units. We'll be using the metric system for units where lengths will be in meters and masses in kilograms.

The length overall L H L_H LH (often abbreviated as LOA) is the maximum length dimension. The waterline length L W L L_{WL} LWL is the maximum length touching the water if the water is still. In this example, they are nearly the same, L H = 12.2    m L_H = 12.2 \; m LH=12.2m and L W L = 12.0    m L_{WL} = 12.0 \; m LWL=12.0m. The canoe body draft (often just called draft), T C T_C TC, which is the maximum depth below the waterline, will be derived later.

side-view.png

In SMath Studio you can enter values for total length and waterline length like this:

smath-studio-length-parameters.png

SMath Studio lets you start equations anywhere on the page, although an equation calling on a previous equation needs to be below the first equation. To enter L H L_H LH type "L" followed by a period "." which starts a subscript, then "H". Type ":" to assign a value to your variable, and SMath Studio will fill in the equal sign after the colon. Now enter the variable value of 12.2 12.2 12.2 followed by "m" for meters. A units menu box will open showing available units beginning with the letter "m":

smath-studio-length-unit-entry.png

Press the TAB key to accept the choice of meters. Units are shown in blue on the worksheet. If you see a black "m" where you think a meter unit should be, go back and re-enter the unit.

Typing text somewhere on the page generates a comment. Click on the palette icon on the toolbar and you can change the background color of the comment.

smath-studio-text-color.png

That's almost everything you need to know about SMath Studio for this project.

Getting back to the design, everything else you need to input is dimensionless (that is, unitless) coefficients and ratios. The first one is the length-to-beam ratio which is the ratio of the waterline length to the maximum width of one hull at the waterline B W L B_{WL} BWL,

L B R = L W L B W L . L_{BR}=\frac{L_{WL}}{B_{WL}}. LBR=BWLLWL.

Select L B R L_{BR} LBR somewhere between 9 and 12 for a cruising catamaran. Lower values will make the hulls narrower, and decrease draft if all other parameters remain the same. For this example, set L B R = 11.0 L_{BR} = 11.0 LBR=11.0. In this figure, the inner curve is the waterline which is where the hull touches the surface of the water. The outer curve is the outline of the hull. single-hull-top-view.png

Reducing the width of the hull would result in a deeper draft, but doesn't here is because of the next ratio, the beam to draft ratio, B T R B_{TR} BTR, which is used to determine the canoe body draft

T C = B W L B T R . T_C = \frac{B_{WL}}{B_{TR}}. TC=BTRBWL.

The beam to draft ratio is generally between 1.5 and 2.8. Setting B T R ≈ 2 B_{TR} \approx 2 BTR2 minimizes friction resistance while wave resistance is minimized for values slightly below 2. We'll set B T R = 1.9 B_{TR} = 1.9 BTR=1.9 for now. Frictional resistance is proportional to the surface area of the hull in the water.

The boat's motion generates a wave at the bow. Wave resistance increases approximately as the cube of boat speed.

The midship coefficient, C m C_m Cm is a measure of the shape of the cross-section defined as the ratio of the cross-section area, A m A_m Am divided by the bounding rectangle formed by the waterline beam, B W L B_{WL} BWL, and the draft, T C T_C TC,

C m = A m B W L T C . C_m = \frac{A_m}{B_{WL}T_C}. Cm=BWLTCAm.

section-view.png

Since the area of a triangle is 1 2 w h \frac{1}{2} wh 21wh, a vee-shaped hull will have a minimum at C m = 0.5 C_m = 0.5 Cm=0.5, while elliptical cross-section results in C m = 0.785 C_m = 0.785 Cm=0.785, and a more rectangular cross-section might have C m = 0.9 C_m = 0.9 Cm=0.9. Large cargo ships have rectangular sections, but for sailboats, elliptical cross-sections are more typical.

James Wharram designs catamarans meant to be cut from sheets of marine plywood with no rounded cross-sections, such as the Ariki shown here. Construction is simpler, but the midship coefficient is quite low.

ariki.jpg

The area of any ellipse divided by the area of the enclosed bounding box may be calculated as follows. The parametric equation for an ellipse is ( a cos ⁡ t , b sin ⁡ t ) (a \cos t, b \sin t) (acost,bsint) with t ∈ [ 0 , 2 π ] t \in [0,2\pi] t[0,2π], ellipse.png so the area of a quarter ellipse is

A = ∫ 0 a y    d x = ∫ π 2 0 b sin ⁡ t ( − a sin ⁡ t )    d t = ∫ 0 π / 2 a b sin ⁡ 2 t    d t = a b 2 ( t − sin ⁡ ( 2 t ) 2 ) ∣ 0 π / 2 = a b π 4 . \begin{aligned} A &= \int_0^a y \; dx\\ &= \int_{\frac{\pi}{2}}^0 b \sin t (-a \sin t) \; dt \\ &= \int_0^{\pi/2} ab \sin^2 t \; dt \\ &= \left. \frac{ab}{2} \left( t - \frac{\sin(2t)}{2} \right) \right|_0^{\pi/2} = ab \frac{\pi}{4}. \end{aligned} A=0aydx=2π0bsint(asint)dt=0π/2absin2tdt=2ab(t2sin(2t)) 0π/2=ab4π.

Dividing by the area of the circumscribing rectangle A = a b A = ab A=ab gives the ratio π 4 ≈ 0.785 \frac{\pi}{4} \approx 0.785 4π0.785 which holds for any ellipse. Stronger hulls with less wetted surface area arise from elliptical cross-sections, so we'll assume for this design the hulls are elliptical, and set C m = 0.785 C_m = 0.785 Cm=0.785.

A boat's displacement is the volume of water with a mass equal to the mass of the boat. The prismatic coefficient is defined as

C p = D A m L W L C_p = \frac{D}{A_m L_{WL}} Cp=AmLWLD

where D D D is the displacement. For catamarans, 0.55 ≤ C p ≤ 0.64 0.55 \leq C_p \leq 0.64 0.55Cp0.64, and in this example, we've chosen C p = 0.59 C_p = 0.59 Cp=0.59. The prismatic coefficient for monohulls is usually in the range, 0.5 ≤ C p ≤ 0.6 0.5 \leq C_p \leq 0.6 0.5Cp0.6. The prismatic coefficient is the ratio of the actual displacement to the displacement that would be generated if the maximum cross-section were extended to the length of the waterline. Lower values of C p C_p Cp favor boats in light air while higher values are better for boats at or above hull speed (discussed below).

The waterplane coefficient is the ratio of the area of the boat at the waterline to the product of the waterline beam, B W L B_{WL} BWL, times the waterline length, L W L L_{WL} LWL,

C w = A w B W L L W L , C_w = \frac{A_w}{B_{WL}L_{WL}}, Cw=BWLLWLAw,

where A w A_w Aw is the waterplane area. For catamarans, the waterline beam for the aft 2 / 3 r d 2/3^{rd} 2/3rd of the hull is almost equal to B W L B_{WL} BWL, and the forward third narrows towards the bow. Looking at the sketch below I guessed that the waterplane area in the aft portion is about 90% of the rectangle and maybe 33% in the forward section. This gives

C w = . 9 × 2 3 + . 33 × 1 3 = 0.71. C_w = .9 \times \frac{2}{3} + .33 \times \frac{1}{3} = 0.71. Cw=.9×32+.33×31=0.71.

Lucky guess! This is exactly the value Terho Halme chose for the example. Values for C w C_w Cw are typically between 0.69 and 0.72 for catamarans.

Waterplane coefficient

The final design parameter is the length to beam ratio, L B R C LBRC LBRC, which is the ratio of the overall length to the distance between hull centers or centerline beam,

L B R C = L H B C B . LBRC = \frac{L_H}{B_{CB}}. LBRC=BCBLH.

two-hulls-top-view.png

A high value of L B R C LBRC LBRC means B C B B_{CB} BCB is small relative to the length and will reduce transverse stability, possibly leading to capsize. Low values mean that the boat will be wider and heavier, but able to carry more sail area. For this example, set L B R C = 2.2 LBRC = 2.2 LBRC=2.2.

Your SMath Studio worksheet should look something like this: smath-studio-design-parameters.png

I also included two constants, the number of hulls n H u l l = 2 n_{Hull} = 2 nHull=2 and the density of seawater, ρ S W = 1025 k g m 3 \rho_{SW} = 1025 \frac{kg}{m^3} ρSW=1025m3kg.

Derived Parameters

The remaining calculations are derived from these parameters and can be added to the worksheet somewhere below the design parameters. Many of them have been described already in terms of the design parameters.

When you add these equations to the worksheet include an equals sign "=" at the end to make SMath Studio immediately evaluate them. smath-studio-derived-parameters.png

The loaded displacement m L D C m_{LDC} mLDC is calculated from the volume enclosed by the waterline length and beam times the draft, then scaled by prismatic and midship coefficients and multiplied by the density of seawater. This is the mass of water displaced by the boat sitting in the water which must equal the mass of the boat. Empty and light displacements are also calculated as 70% and 80% of the loaded displacement, respectively.

Transversal and longitudinal stability are indicators of how well the boat can right itself. Stability is the length of the righting moment arm, so a unit of meters is correct.

Since the purpose was to show why keeping track of units is important, let's take a look at the last two calculations, engine power, and hull speed. Power should be in units of kilowatts and speed in knots, but instead, we get cubic meters and square root meters.

The first is easy to explain, and Halme says the engine power needed for a catamaran is about 4 KW / ton, which is the 4 in the equation. His equation for power is

P m = 4 m L D C 1025 P_m = 4 \frac{m{LDC}}{1025} Pm=41025mLDC

and I assumed 1025 1025 1025 was for the density of seawater with units of k g m 3 \frac{kg}{m^3} m3kg. Instead, 1025 1025 1025 should be replaced with the conversion from k g kg kg to metric tons, 1000    k g = 1    ton 1000 \; kg = 1 \; \text{ton} 1000kg=1ton. If we fix the equation we get

P m = 4 k W ton m L D C    k g 1000 k g ton = 4 ⋅ m L D C    k W = 28.5457    k W . P_m = 4\frac{kW}{\text{ton}} \frac{mLDC \; kg}{1000 \frac{kg}{\text{ton}}} = 4 \cdot mLDC \; kW = 28.5457 \; kW. Pm=4tonkW1000tonkgmLDCkg=4mLDCkW=28.5457kW.

If you want to calculate sail areas, daggerboard, and rudder dimensions or estimate the speed under sail, see Basics of Boat Design or How to dimension a sailing catamaran? written by Terho Halme.

This still leaves the problem of the hull speed units.

Froude Numbers

William Froude worked with Isambard Brunel on other engineering projects and had a minor role in the design of the Great Eastern. He realized that naval architects of the time had a poor understanding of engine power requirements and thought that by studying scale models of several ships he could get a better estimate of power requirements.

He built a towing tank at his home in Torquay, a seaside town in Devon, England, where he was able to measure the power required to pull the one, two, and four meter long models down the 100-meter tank. He found the ratio of the inertial force to the gravitational force remained constant as the scale dimensions changed, now called the Froude number,

F r = Inertial force Gravitational force . Fr = \frac{\text{Inertial force}}{\text{Gravitational force}}. Fr=Gravitational forceInertial force.

Force is mass times acceleration (Newton's Second Law of Motion). The units of force are in Newtons ( N N N), where 1    N = 1 k g − m s 2 1 \; N = 1 \frac{kg-m}{s^2} 1N=1s2kgm. The gravitational force acting on a ship is the acceleration due to gravity, 9.8 m s 2 9.8 \frac{m}{s^2} 9.8s2m times the mass of the ship in k g kg kg.

In fluid mechanics, the inertial force is often written as

ρ L 3 V t \rho L^3 \frac{V}{t} ρL3tV

where ρ \rho ρ is the density of the fluid and L L L is a characteristic length defining the scale of the system and used to generate dimensionless quantities such as the Froude number. In this case, the characteristic length is the length of the ship. Now, ρ L 3 \rho L^3 ρL3 is mass per unit volume times the volume of the characteristic length, k g m 3 × m 3 \frac{kg}{m^3}\times m^3 m3kg×m3, which becomes mass, and V t \frac{V}{t} tV is velocity over time, or acceleration, a = m s × 1 s = m s 2 a = \frac{m}{s} \times \frac{1}{s} = \frac{m}{s^2} a=sm×s1=s2m. So, we could write this as

ρ L 3 V t = m V t = m L V L t = m L V 2 . \rho L^3 \frac{V}{t} = m \frac{V}{t} = \frac{m}{L} V \frac{L}{t} = \frac{m}{L}V^2. ρL3tV=mtV=LmVtL=LmV2.

Froude's number squared is then

F r 2 = m L V 2 g m . Fr^2 = \frac{\frac{m}{L}V^2}{g m}. Fr2=gmLmV2.

Factoring out the mass m m m and taking square roots,

F r = V g L , Fr = \frac{V}{\sqrt{gL}}, Fr=gL V,

which is the usual representation. Froude also measured the drag, or viscous resistance, of a flat plate through the water. The dimensionless Reynolds number is the ratio of the inertial force to the viscous force and is used to estimate the viscous resistance. This plot shows the components forming the total resistance as a function of ship speed (from "Resistance and powering of ships"): resistance-and-powering-of-ships.png Now we can write the velocity in terms of the Froude number,

V = F r g L = 3.1 F r L    m s V = Fr \sqrt{gL} = 3.1 Fr \sqrt{L} \; \frac{m}{s} V=FrgL =3.1FrL sm

since g = 9.8 m s 2 g = 9.8 \frac{m}{s^2} g=9.8s2m. This means Halme's hull speed of 8.5    m s 8.5 \; \frac{m}{s} 8.5sm corresponds to a Froude number of 0.78 0.78 0.78.

From Thin or bulky: optimal aspect ratios for ship hulls, the Froude number for multihull sailboats seems to cluster around 1.25 1.25 1.25, a bit higher than Halme's estimate ( U U U is the velocity): optimal-aspect-ratios-for-ship-hulls.png What's so special about the hull speed? It's the speed where the boat length L L L matches the length of the wave L w Lw Lw generated by the bow. When this happens, the stern just catches up to the previous bow wave and can ride on it. Above this speed, the boat is always working uphill against the bow wave as can be seen here. This is why the characteristic length was chosen to be the length of the boat.

wavelength-vs-speed.png

Velociraptors

Mathematician and biologist D’Arcy Wentworth Thompson connected animal locomotion to Froude's number when he realized the strength of leg muscles is proportional to the cross-section area of the bones. He compared the similarity of animals to ship propulsion,

"In two similar and closely related animals, as also in two steam engines, the law is bound to hold that the rate of working must tend to vary with the square of the linear dimensions, according to Froude’s Law of steamship comparison."

Robert McNeill Alexander, Professor of Zoology at the University of Leeds compared the stride lengths of animals to hip heights to estimate bipedal running speeds. The hip height became the characteristic length in his equation,

λ L = 2.3 ( V 2 g L ) 0.3 \frac{\lambda}{L} = 2.3 \left( \frac{V^2}{gL} \right)^{0.3} Lλ=2.3(gLV2)0.3

where λ \lambda λ is the stride length and L L L is the height of the hip above the ground for bipedal locomotion. This can be solved for the speed,

V = 0.25 g λ 1.67 L 1.17 . V = 0.25 \sqrt{g} \frac{\lambda^{1.67}}{L^{1.17}}. V=0.25g L1.17λ1.67.

Alexander measured the stride length from dinosaur footprints and got the characteristic length from fossil bones and estimated the running speed of Tyrannosaurus Rex. Unlike what's portrayed in the movie "Jurassic Park" Alexander concluded, "Tyrannosaurus was not good at chasing Jeeps."

But, what about velociraptors? They will stalk you, hunt you down, and chew you to little bits. Is there any hope to outrun them?

velociraptors.jpg

Velociraptor stride lengths are usually about λ = 1 4    m \lambda = \frac{1}{4} \; m λ=41m, and the hip height L = 0.72    m L = 0.72 \; m L=0.72m giving a velocity estimate of

V = 0.25    9.8    0.2 5 1.67 0.7 2 1.17 = 0.11    m s ≈ 0.25    mph . V = 0.25 \; \sqrt{9.8} \; \frac{0.25^{1.67}}{0.72^{1.17}} = 0.11 \; \frac{m}{s} \approx 0.25 \; \text{mph}. V=0.259.8 0.721.170.251.67=0.11sm0.25mph.

The stride lengths were taken from footprints where the velociraptor may have been walking, so running speeds could be a bit higher. In the movie, velociraptors were scaled up quite a lot as the actual animal was probably about the size of a turkey and may have been more like this, um, vicious-looking dinosaur: vicious-velociraptor.jpg

Walking on the Moon

The velocity estimate depends on the gravitational acceleration constant g g g which is 9.8    m s 2 9.8 \; \frac{m}{s^2} 9.8s2m on Earth, but on the Moon, the force of gravity is about 16% of Earth's gravity, so the velocity becomes

V Moon = V Earth g Moon g Earth = 0.4 × V Earth . V_{\text{Moon}} = V_{\text{Earth}} \sqrt{\frac{g_{\text{Moon}}}{g_{\text{Earth}}}} = 0.4 \times V_{\text{Earth}}. VMoon=VEarthgEarthgMoon =0.4×VEarth.

The Apollo astronauts found that attempting to walk at normal terrestrial speeds on the Moon was quite difficult and wound up skipping instead. Even though gravity is reduced on the Moon, their walking speed and stride length didn't match what they were used to on Earth. A more complete description of the connection between the Froude number and walking and running speeds can be found in the paper by Vaughn and O'Malley, "Froude and the contribution of naval architecture to our understanding of bipedal locomotion."

Froude's Legacy: Connecting Earth, Sea, and Moon

What began as a simple exercise in catamaran design led us through an fascinating exploration of Froude numbers - a mathematical relationship that elegantly connects seemingly unrelated phenomena. William Froude's work on ship resistance not only helped naval architects understand why the mighty Great Eastern struggled to reach its intended speed, but also provided a fundamental framework that extends far beyond maritime engineering.

This same mathematical principle that governs how waves form around a ship's hull also helps explain why velociraptors weren't quite the terrifying sprinters portrayed in "Jurassic Park," and why Apollo astronauts found themselves bouncing rather than walking on the lunar surface. In each case, the ratio between inertial and gravitational forces - expressed through the Froude number - determines the natural dynamics of motion.

Whether we're calculating hull speeds for a catamaran, analyzing fossilized dinosaur footprints, or studying human locomotion in reduced gravity, the underlying physics remains the same. It's a beautiful reminder that mathematics often reveals unexpected connections in our universe, linking the waves of ancient sailing ships to the footsteps of both extinct predators and moonwalking astronauts.

Perhaps the most valuable lesson here isn't just about naval architecture or biomechanics - it's about the importance of understanding fundamental principles and keeping track of our units. As we saw with our initial power calculations, a simple unit error can lead us astray. But when we get the math right, it can illuminate connections across space, time, and even species, showing us how the same physical laws govern motion whether on Earth, on the Moon, or on the waves of the ocean.

William Froude's experiments show that a boat's hull speed, the length of a velociraptor's stride, and the effects of the gravitational field on the Moon are all related.

Image credits

Code and Software

SMath Studio: Tiny, powerful, free mathematical program with WYSIWYG editor and complete units of measurements support.

See all software used on Wild Peaches →

References and further reading

📬 Subscribe and stay in the loop. · Email · GitHub · Forum · Facebook

© 2020 – 2025 Jan De Wilde & John Peach