Module 1: Hydrodynamics & Swimming

In 1936, the Cambridge zoologist James Gray estimated the drag on a fast-swimming dolphin from the then-current rigid-body hydrodynamics. His calculation showed that a dolphin should need roughly seven times more muscle powerthan it actually possesses to sustain the observed speeds of 10–20 m/s. This discrepancy, known ever since as Gray's paradox, forced biologists to look for unusual drag-reduction mechanisms. Eighty years later we understand cetacean locomotion as a near-optimal solution to a well-defined fluid dynamics problem, with the dolphin's compliant skin, the lunate-fluke geometry, and the universal Strouhal optimum all contributing to an animal that moves through water with extraordinary efficiency.

1. Reynolds Number and Drag at Cetacean Scale

The relevant dimensionless number for a dolphin-sized swimmer is the Reynolds number\(\,Re = \rho U L / \mu\,\). For a 2.5 m bottlenose dolphin cruising at\(\,U = 5\,\text{m/s}\,\) in seawater (\(\rho = 1025\,\text{kg/m}^3\),\(\mu = 1.08\times10^{-3}\,\text{Pa s}\)):

\[ Re = \frac{1025 \cdot 5 \cdot 2.5}{1.08\times10^{-3}} \approx 1.2\times10^7 \]

This is firmly in the turbulent regime for a rigid body; Blasius flat-plate transition occurs at \(Re \sim 5\times10^5\). Gray assumed a rigid body in fully turbulent flow and used the Prandtl 1/7-power drag law:

\[ C_{D,\text{turb}} = \frac{0.074}{Re^{1/5}} \approx 0.0028 \quad\text{at}\quad Re = 10^7 \]

With \(A_{wet} \approx 1.8\,\text{m}^2\) the drag force is \(D = \tfrac{1}{2}\rho U^2 C_D A_{wet} \approx 65\,\text{N}\)

The required sustained power is \(P = D\cdot U \approx 325\,\text{W}\). Gray's analysis gave a higher effective drag coefficient (~0.008) because he included form drag from a rigid bluff body; he concluded a dolphin needed ~2 kW of muscle power, far above the ~280 W that dolphin muscle mass should deliver.

1.1 Resolution: Laminar Flow via Compliant Skin

Max Kramer (1957–60) and subsequent workers (Carpenter, Bushnell, Fish) have shown that dolphin skin is not a rigid surface. The epidermis contains a fine network of dermal papillae and an underlying layer of compliant, fluid-filled ridges. These features damp incipient Tollmien-Schlichting waves β€” the precursors to boundary-layer transition β€” by extracting energy from them elastically. The result is that a substantial portion of the dolphin's boundary layer remains laminar even at Re β‰ˆ 107. For laminar flow, the Blasius drag coefficient is much lower:

\[ C_{D,\text{lam}} = \frac{1.328}{\sqrt{Re}} \approx 0.0004 \quad\text{at}\quad Re = 10^7 \]

Real dolphin drag coefficients (measured in tow-tank experiments and inferred from glide decelerations) lie around \(C_D \approx 0.0036\), substantially below fully turbulent and above fully laminar predictions β€” consistent with a mixed laminar/transitional boundary layer.

2. The Lunate Fluke and the Strouhal Optimum

Cetaceans generate thrust by oscillating a lunate tail fluke(high-aspect-ratio, crescent-shaped) in the vertical plane. Fluke amplitude, frequency, and forward speed are related via the Strouhal number:

\[ St = \frac{f A}{U} \]

where \(f\) is stroke frequency (Hz), \(A\) is peak-to-peak tip amplitude (m), and \(U\) is forward speed (m/s)

The remarkable empirical observation, first emphasised by Triantafyllou et al. (1991, 1993) and Taylor et al. (2003), is that virtually all efficient swimming and flying animals cruise in the narrow band \(St^* \approx 0.25\text{--}0.35\) β€” a universal biomechanical optimum that transcends phylogeny, body plan, and medium. Bottlenose dolphins, bluefin tuna, salmon, great white sharks, pigeons, hummingbirds, and fruit flies all converge on this narrow range. Theoretical work on two-dimensional oscillating foils confirms that thrust efficiency peaks at\(St \approx 0.3\) because this is where the shed vortex street (a β€œreverse KΓ‘rmΓ‘n street”) is most coherent.

2.1 Lighthill's Elongated-Body Theory

Sir James Lighthill (1960, 1970) derived the mean thrust of an oscillating caudal fin using slender-body potential-flow theory. For an airfoil of chord \(c\)oscillating in heave with amplitude \(A\) at angular frequency\(\omega\), the mean thrust per unit span is:

\[ \bar T = \tfrac{\pi}{4} \rho c^2 \omega^2 A^2 \cdot [\text{reduction factor}(k)] \]

with reduced frequency \(k = \omega c /(2U)\); the β€œvirtual mass per unit span” is \(m_{\text{virt}} = \tfrac{1}{4}\pi\rho c^2\).

The Theodorsen reduction factor accounts for wake vorticity that reduces thrust below the quasi-steady limit. At \(k \approx 0.3\) (typical cetacean cruising) the reduction factor is ~0.6. Thus thrust scales as \(\bar T \propto \rho c^2 \omega^2 A^2\). Matching thrust to drag yields the swimmer's equilibrium cruising speed.

2.2 Why Lunate?

The crescent (lunate) fluke shape minimizes induced drag for a given thrust by maximizing span-to-chord (aspect) ratio. Just as albatross wings and the horizontal tails of fast pelagic fish (tuna, marlin, swordfish) are high-aspect-ratio and often crescent-shaped, the cetacean fluke has evolved under the same hydrodynamic pressure. The key is to keep the tip vortices small. For an oscillating thin wing, the induced (vortex) drag scales as\(D_{ind} \propto L^2/(\rho U^2 b^2)\) where \(b\) is the span. High aspect ratio (large \(b\) for a given area) minimizes this.

3. Flow Visualization and Fluke Kinematics

Dolphin Flow Field: Laminar Forward, Turbulent Aft, Reverse-KΓ‘rmΓ‘n WakeLaminar boundary layer(damped by compliant skin)Turbulent transition(adverse pressure gradient)Reverse KΓ‘rmΓ‘n street(thrust-producing wake)U β‰ˆ 5 m/sdorsal fin (stabilizer)lunate fluke (high AR)pectoral flipper
Strouhal Number Optimum: Efficient Swimming Converges on St β‰ˆ 0.3Strouhal number St = fA/UThrust efficiency Ξ·0.00.10.20.30.40.50.60.7Optimal band St ∈ [0.25, 0.35]DolphinTunaSalmonSharkBlue whaleKiller whaleΞ·(St) peaks β‰ˆ 0.3(Triantafyllou et al. 1993)

4. Leaping, Bow-Riding, and Biomimetic Surfaces

Cetaceans frequently leap clear of the water β€” breaching in whales, porpoising in dolphins. One long-standing hypothesis (Au & Weihs 1980) is that porpoising is energetically favorable because ballistic flight through air (where drag is 800Γ— lower than in water) is cheaper than swimming near the surface where wave drag adds to form and friction drag. Surface wave drag scales with the Froude number:

\[ Fr = \frac{U}{\sqrt{gL}} \]

For \(L = 2.5\,\text{m}\), surface wave drag peaks at \(U \approx 3\,\text{m/s}\) (\(Fr \approx 0.6\)).

Above this crossover speed, the drag penalty of staying just beneath the surface becomes large enough that leaping is cheaper. Bow-riding β€” where dolphins align with the pressure field at the bow of a moving vessel and catch a β€œfree ride” from the ship's displacement wave β€” is another exploitation of pressure gradients for drag reduction (Williams 1989).

4.1 Biomimetic Skin Surfaces

Two radically different drag-reduction strategies have evolved in swift aquatic predators:

  • Shark denticles (riblets): tiny V-shaped scales of dentine that align with the flow, break up turbulent vortices near the wall, and reduce skin friction by 5–10% at high Re. These work by restricting lateral movement of turbulent structures.
  • Dolphin compliant dermis: an elastic dermal layer that absorbs energy from incipient Tollmien-Schlichting waves and delays laminar-turbulent transition. These work by postponing turbulence altogether.

Both strategies have inspired engineering applications: aircraft manufacturers have experimented with riblet films that yield 2–8% fuel savings, and compliant coatings for ship hulls have been investigated since Kramer's original patent.

5. Power, Cost of Transport, and Optimal Cruising Speed

The power required to swim at speed \(U\) is\(P_{\text{hydro}} = D\,U = \tfrac{1}{2}\rho U^3 C_D A_{wet}\). However the metabolic power \(P_{\text{met}}\) the animal must expend to produce this hydrodynamic power must include two additional terms: the basal metabolic rate (BMR) and the mechanical-to-chemical efficiency \(\eta_{\text{musc}} \approx 0.20\text{--}0.25\) of vertebrate muscle:

\[ P_{\text{met}}(U) \approx \text{BMR} + \frac{1}{\eta_{\text{musc}}}\left( \tfrac{1}{2}\rho U^3 C_D A_{wet} \right) \]

The cost of transport \(COT = P_{\text{met}}/(M_{body} U)\) is the energy required to move 1 kg of body mass 1 m. BMR contributes a term \(\propto 1/U\); the drag term contributes \(\propto U^2\). The sum has a minimum at the optimal cruising speed:

\[ U^* = \left(\frac{2\,\eta_{\text{musc}}\,\text{BMR}}{\rho C_D A_{wet}}\right)^{1/3} \]

For a bottlenose dolphin this predicts \(U^* \approx 2\text{--}3\,\text{m/s}\), in good agreement with observed cruising speeds.

Empirically measured minimum COT for a bottlenose dolphin is roughly 1.3 J/(kg m), about half the mass-specific COT of a human swimmer and about 3Γ— better than a running terrestrial mammal of comparable size. Among cetaceans, blue whales achieve the lowest mass-specific COT known for any animal: ~0.4 J/(kg m) at their preferred speed of ~3 m/s. This energetic efficiency enables the enormous migrations discussed in Module 7.

5.1 Wave Drag and the Optimum Submergence Depth

Surface-wave drag, as noted, peaks near Froude number Fr β‰ˆ 0.5. For a body swimming at depth \(h\) beneath the surface, the wave resistance decays roughly as\(e^{-4\pi h/\lambda_w}\) where \(\lambda_w = 2\pi U^2/g\) is the wavelength of the Kelvin-wave pattern. A dolphin swimming at \(U = 5\,\text{m/s}\)generates \(\lambda_w \approx 16\,\text{m}\); to avoid wave drag it should stay at depths \(h \gtrsim \lambda_w/(4\pi) \approx 1.3\,\text{m}\). Observations show cetaceans routinely cruise at exactly this depth or below.

6. Swim Muscles, Tail Mechanics, and Gait Transitions

Cetacean locomotion is powered by the epaxial (dorsal) and hypaxial (ventral) lumbar muscles, which derive directly from the spinal flexors of their terrestrial ancestors. These muscles, aggregated into massive epaxial and hypaxial blocks along the spine, contract alternately to generate the dorsoventral bending of the tail stock.

6.1 Tendon Energy Storage

The cetacean tail stock contains large elastic tendons and collagenous septa that store elastic energy during each half-stroke and release it during the reversal. An elastic spring of stiffness \(k\) and natural frequency \(\omega_0=\sqrt{k/m}\)resonates with driven oscillations at \(\omega \approx \omega_0\). Observed tail beat frequencies match the estimated natural frequency of the muscle-tendon system, suggesting resonant energy recovery analogous to the Achilles tendon in running mammals. A cetacean cruising at its preferred gait is thus moving in synchrony with the resonance of its own tail.

6.2 Burst vs Cruise Mechanics

Dall's porpoise (Phocoenoides dalli) reaches top bursts of ~55 km/h (15 m/s) and the orca ~56 km/h. Burst speeds require a shift from pure aerobic, slow, oxidative muscle fibers (type I) to anaerobic, fast, glycolytic fibers (type IIb). The cost of this mode is lactate accumulation and restricted duration; all cetaceans show a sharp separation between cruising (aerobic) and sprinting (anaerobic) regimes.

6.3 Allometric Scaling of Swimming Speed

Large cetaceans are not, in the absolute sense, the fastest swimmers. Despite a blue whale's immense mass, its cruise speed is only ~3–4 m/s. Speed scales with body length roughly as \(U \propto L^{1/2}\) across swimming animals (Bainbridge relation for fish extended to cetaceans). The physical basis is that muscle cross-section scales as \(L^2\), force as \(L^2\), and drag as\(L^2 U^2\); setting force equal to drag gives \(U \propto L^0\)β€” constant. Empirically the scaling is closer to \(L^{0.4}\) due to shape and Reynolds corrections.

5. Simulation: Hydrodynamics of Cetacean Swimming

This four-panel simulation (i) plots Strouhal number across swimming and flying animals to reveal the universal optimum; (ii) derives the dolphin's drag coefficient as a function of speed, comparing laminar, turbulent, and compliant-skin predictions to Gray's rigid-body estimate; (iii) situates cetaceans on a Reynolds-number vs body-length diagram of aquatic life; and (iv) computes the Lighthill thrust of an oscillating fluke as a function of amplitude.

Python
script.py166 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Key Observations

  • Panel 1: Cetaceans, fish, pinnipeds, and flyers all cluster in the narrow band St β‰ˆ 0.25–0.35. This is a biological universal.
  • Panel 2: A dolphin with Gray's rigid-body assumption would require C_D β‰ˆ 0.008, giving the famous paradox. Real dolphins operate much closer to the laminar (Blasius) prediction.
  • Panel 3: Large cetaceans (blue whale, sperm whale) reach Re β‰ˆ 108, the highest steady-state Re of any animal.
  • Panel 4: Lighthill thrust goes as A2. The observed dolphin tail amplitude (~20–25% of body length) produces thrust matching drag at cruising speed.

Module Summary

Gray's Paradox (1936)

Rigid-body drag estimates predict dolphins should need 7Γ— more muscle power than available

Compliant Skin

Dermal papillae damp Tollmien-Schlichting waves, keeping much of boundary layer laminar

Dolphin Drag

Measured C_D β‰ˆ 0.0036, between Blasius laminar (0.0004) and Prandtl turbulent (0.0028)

Reynolds Number

Bottlenose at 5 m/s: Re β‰ˆ 1.2Γ—10⁷. Blue whale: Re β‰ˆ 10⁸ (highest of any animal)

Strouhal Optimum

St = fA/U β‰ˆ 0.25–0.35 across cetaceans, fish, birds, bats β€” universal biomechanical limit

Lighthill Thrust

T ∝ ρc²ω²AΒ², virtual mass πρcΒ²/4, Theodorsen reduction factor ~0.6

Lunate Fluke

High aspect ratio minimizes induced (vortex) drag β€” convergent with tuna tails, albatross wings

Froude Crossover

Above U β‰ˆ 3 m/s porpoising saves energy by escaping surface wave drag

References

  1. Gray, J. (1936). Studies in animal locomotion. VI. The propulsive powers of the dolphin. Journal of Experimental Biology, 13, 192–199.
  2. Kramer, M.O. (1960). Boundary layer stabilization by distributed damping. Journal of the American Society of Naval Engineers, 72, 25–34.
  3. Carpenter, P.W. & Garrad, A.D. (1985). The hydrodynamic stability of flow over Kramer-type compliant surfaces. JFM, 155, 465–510.
  4. Lighthill, M.J. (1960). Note on the swimming of slender fish. Journal of Fluid Mechanics, 9(2), 305–317.
  5. Lighthill, M.J. (1970). Aquatic animal propulsion of high hydromechanical efficiency. JFM, 44, 265–301.
  6. Triantafyllou, M.S., Triantafyllou, G.S. & Gopalkrishnan, R. (1991). Wake mechanics for thrust generation in oscillating foils. Physics of Fluids A, 3(12), 2835–2837.
  7. Taylor, G.K., Nudds, R.L. & Thomas, A.L.R. (2003). Flying and swimming animals cruise at a Strouhal number tuned for high power efficiency. Nature, 425, 707–711.
  8. Fish, F.E. (1993). Power output and propulsive efficiency of swimming bottlenose dolphins. Journal of Experimental Biology, 185, 179–193.
  9. Fish, F.E. & Lauder, G.V. (2006). Passive and active flow control by swimming fishes and mammals. Annual Review of Fluid Mechanics, 38, 193–224.
  10. Au, D. & Weihs, D. (1980). At high speeds dolphins save energy by leaping. Nature, 284, 548–550.
  11. Williams, T.M. et al. (1992). Travel at low energetic cost by swimming and wave-riding bottlenose dolphins. Nature, 355, 821–823.