Abstract

The best swimmers have a streamlined shape that ensures an attached flow pattern and a laminar boundary layer at rather large values of the Reynolds number. Simple expressions may be obtained for the volumetric drag coefficient of an ideal body of revolution under laminar unseparated flow conditions together with estimations of a critical value of the Reynolds number. A measure, the capacity-efficiency factor, calculated for different organisms and underwater vehicles, shows that information about animal shapes and locomotion is of utmost biological interest and could be useful to improve robot fish and underwater vehicles as well.

1. Introduction

From the hydromechanical point of view fish swimming is a very complicated unsteady phenomenon. A description of the diversity of fish locomotion and a classification of swimming modes, categories, and styles can be found in Blake [1]. Moreover, recent interest in robot fish requires answering questions about optimal shape and the power necessary for fish locomotion. This study focuses on the most simple estimations of fish drag and power requirements during quasisteady motion when changes in body shape can be neglected. This approach would not be very reliable in the case of anguilliform propulsion, but it is acceptable for the carangi- and thunniforms of the best swimmers (see [24]).

The best swimmers must have a small drag coefficient: to achieve high underwater velocity. In this case a body of volume moving at a constant speed in a fluid with density exhibits low drag . Owing to the huge difference in water and air densities, an underwater body has approximately 800 times larger drag than the same body moving in air at the same velocity, since, at equal Reynolds numbers, the drag coefficients are approximately equal. For large Reynolds number the boundary layer thickness can be neglected and fluid outside the thin layer can be treated as ideal and incompressible. For underwater swimming the wave drag, important for ships or animals moving on the water surface, is negligible.

The total drag can be divided into pressure and friction components. The pressure drag can be reduced almost to zero when the boundary layer does not separate from the body surface. In this case the d’Alembert paradox applies; that is, a closed rigid body in unbounded flow of ideal incompressible fluid has zero drag (see, e.g., [5]). Experiments with animals that are good swimmers (like dolphins) have shown that during gliding (inertial movement without manoeuvring and shape change) they exhibit unseparated flow pattern (e.g., [6, 7]). Attached flow has also been observed around the special-shaped body UA-2 [8]. Its UA-2c version is shown in Figure 1.

To reduce friction, a laminar boundary layer must be maintained as large as possible over the wetted surface. Therefore, an ideal good swimmer must have an unseparated laminar boundary layer over its whole surface. In this paper the drag coefficient and the capacity characteristics of a corresponding slender body of revolution (with a small ratio of maximum diameter to length ) have been calculated for different animals, human athletes, and submarines in order to compare them in terms of swimming efficiency.

2. Materials and Methods

2.1. Frictional Drag on a Slender Body of Revolution

To estimate laminar frictional drag on a slender body of revolution, the Mangler-Stepanov transformations (see, e.g., [5]), which reduce the rotationally symmetric boundary layer equations to a two-dimensional case, can be used. The following relations between coordinates , for the rotationally symmetric boundary layer (shown in Figure 2) and the corresponding two-dimensional coordinates , are valid (coordinates are dimensionless based on the body length): where is the dimensionless radius of the rotationally symmetric body based on its length . The flow velocity at the outer edge of the boundary layer, the displacement thickness, and the skin-friction coefficient are related as follows (see, e.g., [5]): All the values in (3) are dimensionless, based on ambient flow velocity , body length , and , respectively. These equations are valid for an arbitrary rotationally symmetric body provided that the thickness of the boundary layer is small in comparison with the radius; that is, the flow is unseparated. For a slender body, the coordinate can be calculated along the body’s axis and the velocity can be supposed to be equal to unity, neglecting the thickness of the boundary layer and the pressure distribution peculiarities (see, e.g., [10]). From the first equation of (3) the value of will also be equal to unity; that is, the rotationally symmetric boundary layer on a slender body can be reduced to the flat plate one [11]. According to the Blasius expression, the equality holds for a laminar flow, where the Reynolds number and is the kinematic viscosity (see, e.g., [5]). Introducing the variable , using (2) and (3), the laminar skin-friction drag coefficients of a slender rotationally symmetric body may be obtained as and using (1), where the volumetric Reynolds number .

Note that the volumetric frictional drag coefficient does not depend on the slender body shape, provided its volume remains constant (see also [11]). Anyway, this is valid for laminar attached boundary layer and at limited only. In the next section the critical value of the Reynolds number will be calculated. But in any case (5) is a reliable estimate for the minimum possible drag on a rigid body of revolution.

2.2. Critical Values of the Reynolds Number

The laminar-to-turbulent flow transition in the boundary layer influences the skin-friction drag and depends on many parameters such as pressure gradient, surface roughness, and pulsations in the ambient flow (e.g., [5]). Nevertheless, according to the Tollmien-Schlichting-Lin theory (e.g., [12]), the boundary layer on a flat plate remains laminar for any frequencies of disturbances if This inequality, taking into account the Blasius expression for displacement thickness (e.g., [5]) , can be rewritten as follows: In the previous section it has been shown that the boundary layer around a slender body of revolution can be reduced to that on a plate with the use of the Mangler-Stepanov transformation. Given (2) and (7), the condition for the axisymmetrical boundary layer to remain laminar can be written as follows [13]: If the boundary layer remains laminar over the entire surface, the integral in (8) can be substituted by Thus, To calculate the critical values of the Reynolds numbers and , corresponding to the right-hand parts of the inequalities (10), information about the body shape is necessary. Supposing the optimal body shape to be close to that of the dolphin or UA-2c, the ratio can be determined directly after calculating the corresponding shape with the use of the method proposed in Nesteruk [8]. Typical shapes for various ratios are shown in Figure 3. The ratio can be calculated from the following approximate expression: with the values of parameter ; 0.2998; 0.3335 for ; 0.21; 0.3, respectively. To simplify the calculations, the average value is used.

Both the forebody and the tail of the shape corresponding to the smallest thickness ratio are concave, while for less slender bodies ( and ) only the tail is concave (see Figure 3). Some fast-swimming fish have a concave forebody too (e.g., the Mediterranean spearfish Tetrapturus belone, Indo-Pacific sailfish Istiophorus platypterus, black marlin Makaira indica, or swordfish Xiphias gladius).

With the use of (10) and (11) the critical Reynolds number can be estimated as follows: Equations (12) show that the boundary layer remains laminar on slender bodies of revolution at rather large Reynolds numbers and the critical value of the Reynolds number increases with the diminishing of the thickness ratio .

2.3. Power Requirements for Swimming

The power balance for steady swimming at maximum velocity , when the thrust is equal to the drag , can be written as follows: where is the physiological maximum of the available animal power per unit mass, is that mass, and is the propulsion efficiency, which takes into account the fin drag. Assuming a neutrally buoyant animal (, in the case of negatively buoyant animals the difference between weight and buoyancy does not exceed 10%, e.g., [14], and can be neglected), the capacity-efficiency factor can be estimated as If the animal body shape is close to the slender rotationally symmetric one and ensures an unseparated flow pattern, then for a purely laminar boundary layer, (5) and (14) yield In dimensionless form (15) can be written as follows:

3. Results

3.1. Which Shape Is Better?

When designing robot fish or other underwater vehicles, the goal is to minimize the drag for a given volume and a prescribed velocity range. To solve this problem, the minimum volumetric drag coefficient for the hull (the main part of the body) must be achieved at the given volumetric Reynolds number . Thus, the question arises: what kind of shape must be chosen?

Since for laminar attached flow on a slender body of revolution is independent of the shape (see (5)), the hull shape may be arbitrary provided it ensures a laminar attached flow. For example, any shape shown in Figure 3 can be used at the corresponding subcritical Reynolds number range. If , the tail part of the hull has a turbulent boundary layer and the optimal shape depends on the thickness ratio and tail shape peculiarities.

It is rather difficult to calculate the frictional drag for supercritical Reynolds numbers, especially in the laminar-to-turbulent transition region. Expression (5) has been used for the laminar forebody and the flat plate concept (e.g., [15]) for the turbulent tail. The results are represented in Figure 3. For subcritical Reynolds numbers the universal straight line corresponds to (5) in logarithmic coordinates. For as shown in Figure 3 this line bifurcates into different colour lines corresponding to the shapes with different values of the thickness ratio ; 0.21; 0.3 (green, yellow, and red lines, resp.).

For different animals the typical Reynolds numbers are shown in Figure 3. For comparison, the cases of a human sportsman (during underwater dolphin-kick swimming at velocity 2.7 m/s) and two different submarines are also presented in Figure 3. All names shown in green correspond to the values of the thickness ratio , in yellow , and in red .

To emphasize the fact that the theoretical curves in Figure 3 show the minimal possible values of the drag, all names are written over the corresponding lines. For example, the real total drag coefficients for submarines (based on their mass, maximum underwater velocity, capacity, propeller efficiency 0.85, and ; also shown in Figure 3 by green markers “+”) are 3–5 times greater than the minimum theoretical values shown by the green line. This fact can be explained by the presence of separation on the submarines’ hulls. Small organisms (e.g., mosquitofish) have greater values of than the submarines (see Figure 3).

The minimum possible values of correspond to animals with very slender bodies (e.g., narrow-barred Spanish mackerel Scomberomorus commerson, wahoo Acanthocybium solandri, Mediterranean spearfish Tetrapturus belone, Indo-Pacific sailfish Istiophorus platypterus, and saltwater crocodile Crocodylus porosus) and are located at the subcritical Reynolds number range . For plumper animals (e.g., dolphins, sharks) the same range is already supercritical and greater values of are expected. Whales and some fast-swimming fish (e.g., swordfish Xiphias gladius, black marlin Makaira indica) swim at supercritical Reynolds numbers.

3.2. Capacity-Efficiency Factor

Equation (15) allows estimating the capacity-efficiency factor for different animals and vehicles with the use of their data for maximum speed, volume, and viscosity of water. Large values of correspond to the better swimmers (i.e., the animals that can produce more energy per unit time and per unit mass and use it effectively for locomotion). Smaller values of this factor correspond to animals or vehicles that have a large drag (greater than the theoretical minimum (5)) and do not need to (or cannot) use a lot of their energy for fast-swimming and/or cannot use their energy in proper way (e.g., due to the small efficiency ).

Examples of calculation of the capacity-efficiency factor are presented in Figure 4 and in Tables 1 and 2. In Table 1 the data from Aleyev [6] (crosses in Figure 4) have been used. In Table 2 the different data about the animals’ mass and maximal velocity have been used for calculations (circles in Figure 4). Different colours of points correspond to the different values of the thickness ratio , as mentioned in the previous section.

Unfortunately, for the majority of the data points, it is very difficult to determine the relevant value of the viscosity of water, which varies from at 0°C to at °C. In Aleyev [6] only data about animal length, the ratio , and are available. To calculate the velocities the average value was used. The volume was estimated with the use of (11). Inaccuracy in the maximum velocity, the mass, and viscosity data is a reason for discrepancies in values obtained for the same animals in Tables 1 and 2.

4. Discussion

4.1. Comparison with the Literature

Since the experimental drag data for live animals are very limited, expression (5) is here related only to rigid bodies of revolution. The body Dolphin was manufactured and tested by North American Aviation in 1967-1968 (see [16]). The profile NACA-66 was chosen for the shape of this body of revolution; its parameters were ;  m;  m; . Tests revealed the minimal value of at ; (5) yields . The more slender body of Hansen and Hoyt [17] (;  m;  m) has the minimum experimental drag coefficient at that is closer to the theoretical value , which can be calculated from (5). The optimal shape X-25 () for an unclosed (tail-boom) body of revolution calculated in Parsons et al. [18] has at . For the same Reynolds number range the theoretical values of the minimum drag coefficient are (see (5)). The higher vales of the drag can be explained by the presence of separation and turbulence on all the above-mentioned bodies of revolution.

The theoretical drag values of different bodies of revolution calculated in Parsons et al. [18], Dodbele et al. [19], Zedan et al. [20], and Lutz and Wagner [21] are rather different. This can probably be explained as due to different semiempirical criteria for the laminar-turbulent transition in the boundary layer. The values of calculated in Parsons et al. [18], Dodbele et al. [19], Zedan et al. [20], and Lutz and Wagner [21] exceed estimation (5). For example, in Zedan et al. [20] the theoretical value at was obtained, whereas the shape calculated by Parsons et al. [18] has 2.4 times less drag and (5) yields a ten times smaller value.

A similar comparison can be performed for the critical values of the Reynolds number calculations (12); for the bodies Dolphin, Hansen and Hoyt, and X-35 the critical Reynolds numbers are , , and , respectively. At the predictions of the laminar-to-turbulent transition coordinate presented by Parsons et al. [18] are approximately and for bodies Dolphin and X-35, respectively. This very slight dependence on the Reynolds number can be explained by the use of the cross-section of the laminar separation as the transition point in Parsons et al. [18], whereas (12) are obtained for shapes without any separation.

Expressions (12) yield or for the bottlenose dolphin with body ratio . These estimations resolve the well-known Gray paradox, since the Reynolds number taken for estimations in Gray [22] corresponds to laminar flow on the dolphin (see also [9]).

Thus slender bodies of revolution can delay laminar-turbulent transitions on their surfaces and reduced skin-friction drag. It must be stressed that relations (12) are valid only for a flow pattern without separation. That is why the effect of the turbulization delay has not been achieved on standard (separated) slender bodies of revolution. The difference in shape can be hardly perceptible (see, e.g., Figure 5), but the pressure distribution is very sensitive to small changes in shape and similar shapes can have very different pressure gradients and separation behaviour. In particular, the body shape UA-2c has , a negative-pressure-gradient forebody ending at the minimum pressure point , a long positive-pressure-gradient region (approximately 45% of the total body length), and a negative-pressure-gradient near its tail. Separation on Goldschmied’s body  m, , with a long negative-pressure-gradient forebody (appox. 76% of the total hull length), a short zone of pressure increase (its length %), and a negative-pressure-gradient region near the tail, was removed only with the use of boundary layer suction; see Goldschmied [23].

4.2. Ranking of Different Animals and Vehicles

It can be seen from Figure 4 and Tables 1 and 2 that the best swimmers are fish whose shape corresponds to the minimal possible values of . These are the Indo-Pacific sailfish, Mediterranean spearfish, narrow-barred Spanish mackerel, and wahoo. Some flying fish and some molluscs (e.g., southern shortfin squid Illex coindetii, European squid Loligo vulgaris) have similar values of . These animals are both perfect swimmers and rather good fliers. In particular, the squid change their shape during flight to create lift forces (see, e.g., [24]).

The capacity-efficiency factor of the best swimmers is approximately 100 times greater than that of common good swimmers. For example, burst swimming corresponds to values of circa 10 body lengths per second for subcarangiform fish of between 10 and 20 cm in length [25]. Applying expressions (11) and (15) and the average value yields where must be taken in metres. For the largest fish of  m and the value of is 3.1 W/kg. According to the data of Azuma [26] the capacity per unit of body mass varies from 10 to 160 W/kg.

The smallest values of are associated with nonstreamlined animals (e.g., ocean sunfish Mola mola, yellow boxfish Ostracion tuberculatus, bowhead whale Balaena mysticetus, or sperm whale Physeter catodon) and submarines. The shapes of these animals and vehicles obviously cannot ensure any attached flow pattern. The low values of for submarines are both the result of the large supercritical Reynolds numbers at which they move (there are huge differences in the theoretical values of both the laminar and the turbulent friction, shown in Figure 3) and of the separation that increases the drag 3- to 5-fold in comparison to the value possible for an attached flow pattern.

Whales have a rather wide range of values of , larger for the hydrodynamically “better shaped” animals; for example, the sei whale Balaenoptera borealis has approximately 10 times greater capacity-efficiency factor than the bowhead whale Balaena mysticetus. The main predator of whales—the killer whale Orcinus orca—has approximately twice as large a value of as the sei whale and comes close to the characteristics of its relative, the bottlenose dolphin Tursiops truncatus.

Humans are not the best swimmers. For example, the world records men have a value similar to some turtles, sturgeons, and the blue whale Balaenoptera musculus, 2000 times smaller than the capacity-efficiency of the best swimmers.

The sharks have a very broad range of capacity-efficiency factor that decreases with increasing Reynolds number. For example, the juvenile blue shark belongs to the best swimmers, while barracuda and adult blue shark have 7–9 times smaller values of . The largest great white sharks that swim at supercritical Reynolds numbers have the smallest value of . The same large difference can be seen in the case of birds. For example, the small and fast chinstrap penguin Pygoscelis antarcticus has approximately 100 times greater value of than the large and slow emperor penguin Aptenodytes forsteri.

The capacity-efficiency factor can be sometimes very close for juvenile and adult animals (e.g., sturgeon) and sometimes very different (e.g., blue shark). In the case of the blue shark the large difference can be explained by the fact that juvenile animals swim at subcritical Reynolds number, whereas the adults swim at transitional and supercritical values of .

4.3. Final Remarks

It can be concluded that the best swimmers have a streamlined shape that ensures an attached flow pattern and a laminar boundary layer at rather large values of the Reynolds number. As a result the hydrodynamical drag can be much smaller in comparison with the “poorly” shaped animals and vehicles where separation and/or turbulence occur.

The large difference in corresponding values (see Figure 5 and Tables 1 and 2) shows that information about animal shapes and locomotion is not only of biological interest but very useful to improve the capabilities of robot fish and underwater vehicles as well. Better measurements of the maximum velocity, mass, and water temperature are necessary to determine the top swimmers among the animals.

The volumetric drag coefficient of an ideal laminar unseparated body of revolution with a prescribed volume is independent of the shape and can be calculated from expression (5), which can be a rather good estimate of the drag on the best-shaped fish. The critical value of the Reynolds number depends on the shape peculiarities and increases as the slenderness ratio decreases (see (12)). At supercritical Reynolds numbers a turbulent boundary layer develops along the body’s tail and the drag increases drastically.

The capacity-efficiency factor (15) can be used to estimate the swimming efficiency of both animals and underwater vehicles. Animals that swim at supercritical Reynolds numbers have much smaller value of in comparison with the best swimmers at subcritical values of .

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Julyan Cartwright for very useful discussions of the results and his assistance in preparing the paper. The study was supported by the EU-financed Project EUMLS (EU-Ukrainian Mathematicians for Life Sciences) Grant agreement PIRSES-GA-2011-295164-EUMLS.