The aim of this article is to provide a theoretical basis upon which to advance and deploy novel tandem flapping foil systems for efficient marine propulsion. We put forth three key insights into tandem flapping foil hydrodynamics related to their choreography, propulsive efficiency, and unsteady loading. in particular, we propose that the performance of the aft foil depends on a new nondimensional number, s/Uτ, which is the interfoil separation s normalized by the distance that the freestream U advects in one flapping period, τ. Additionally, we show how unsteady loading can be mitigated through choice of phase lag.

The extinct ocean-going plesiosaurs were unique within vertebrates because they used two flipper pairs identical in morphology for propulsion. Although fossils of these Mesozoic marine reptiles have been known for more than two centuries, the function and dynamics of their tandem-flipper propulsion system has always been unclear and controversial. We address this question quantitatively for the first time in this study, reporting a series of precisely controlled water tank experiments that use reconstructed plesiosaur flippers scaled from well-preserved fossils. Our aim was to determine which limb movements would have resulted in the most efficient and effective propulsion. We show that plesiosaur hind flippers generated up to 60% more thrust and 40% higher efficiency when operating in harmony with their forward counterparts, when compared with operating alone, and the spacing and relative motion between the flippers was critical in governing these increases. The results of our analyses show that this phenomenon was probably present across the whole range of plesiosaur flipper motion and resolves the centuries-old debate about the propulsion style of these marine reptiles, as well as indicating why they retained two pairs of flippers for more than 100 million years.

The propulsive performance of a pair of tandem flapping foils is sensitively dependent on the spacing and phasing between them. Large increases in thrust and efficiency of the hind foil are possible, but the mechanisms governing these enhancements remain largely unresolved. Two-dimensional numerical simulations of tandem and single foils oscillating in heave and pitch at a Reynolds number of 7000 are performed over a broad and dense parameter space, allowing the effects of inter-foil spacing and phasing to be investigated over a range of non-dimensional frequencies (or Strouhal number, ). Results indicate that the hind foil can produce from no thrust to twice the thrust of a single foil depending on its spacing and phasing with respect to the fore foil, which is consistent with previous studies that were carried out over a limited parameter space. Examination of instantaneous flow fields indicates that high thrust occurs when the hind foil weaves between the vortices that have been shed by the fore foil, and low thrust occurs when the hind foil intercepts these vortices. Contours of high thrust and minimal thrust appear as inclined bands in the parameter space and this behaviour is apparent over the entire range of Strouhal numbers considered . A novel quasi-steady model that utilises kinematics of a virtual hind foil together with data obtained from simulations of a single flapping foil shows that performance augmentation is primarily determined through modification of the instantaneous angle of attack of the hind foil by the vortex street established by the fore foil. This simple model provides estimates of thrust and efficiency for the hind foil, which is consistent with data obtained through full simulations. The limitations of the virtual hind foil method and its physical significance is also discussed.