Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Swimming, flying, and diving behaviors from a unified 2D potential model

## Abstract

Animals swim in water, fly in air, or dive into water to find mates, chase prey, or escape from predators. Even though these locomotion modes are phenomenologically distinct, we can rationalize the underlying hydrodynamic forces using a unified fluid potential model. First, we review the previously known complex potential of a moving thin plate to describe circulation and pressure around the body. Then, the impact force in diving or thrust force in swimming and flying are evaluated from the potential flow model. For the impact force, we show that the slamming or impact force of various ellipsoid-shaped bodies of animals increases with animal weight, however, the impact pressure does not vary much. For fliers, birds and bats follow a linear correlation between thrust lift force and animal weight. For swimming animals, we present a scaling of swimming speed as a balance of thrust force with drag, which is verified with biological data. Under this framework, three distinct animal behaviors (i.e., swimming, flying, and diving) are similar in that a thin appendage displaces and pressurizes a fluid, but different in regards to the surroundings, being either fully immersed in a fluid or at a fluid interface.

## Introduction

In nature, animals move in fluids with different locomotive modes: swimming, flying, jumping out of water, or diving into water. The animals typically gain their propulsive force by flapping fins or wings in a cyclic way. Animals’ flapping appendages are typically thin and wide, and effectively push and pressurize a fluid. Therefore, the motion of such a thin appendage in air or water can be simplified as a rigid thin plate moving in a fluid in order to help understand various animal motions (see Fig. 1).

Swimming or flying locomotion has been extensively studied in various aspects1,2,3,4,5,6,7. In water, most aquatic animals swim by flapping their fins or undulating part of the body3,8. The flapping motion displaces the surrounding fluid, which creates vortices and generates thrust force. Therefore, the aquatic animal can propel forward against drag. For flying animals, the flapping motion of the wings displaces the ambient air and also creates vortices, which produces thrust force. Hence, flying animals are able to lift their own weight and also move forward.

James Lighthill pioneered the small- or large-amplitude elongated body theory to understand the swimming speed through balancing the power generated by an animal with the rate of kinetic energy in a fluid1,9. This slender body approximation quantifies the efficiency of locomotion for aquatic animals analytically, and has been widely used. In another seminal work, Theodore Wu described animal locomotion using an inviscid potential flow2,10, which is an extension of the previously known potential flow of a thin plate. This calculation explains the pressure difference across a thin object while flapping, which is linked to the vortex generation and thrust force. However, these two studies are limited to swimming or flying animals while fully immersed in a fluid.

Animals not only locomote in a fluid, but also cross a air-water interface (e.g., diving). There are two types of diving animals. The first type are aquatic animals that jump out of and then re-enter the water11 (e.g., whales, dolphins, fishes, etc.). At the re-entering moment, animals experience huge impact force as they hit the water at high speeds. The other type are birds plunge-diving into water from air. Several bird species exhibit high-speed diving into water as a hunting mechanism12,13,14. These plunge-diving birds are not very common, but are widely spread in the phylogeny; the Sulidae family species (e.g., Northern Gannet, Brown Booby) and other species (e.g., Brown Pelican, Terns, and Kingfishers). Such a high-speed plunge-diving behavior allows the bird to gain momentum to dive faster and deeper underwater; however, it also induces great compressive force on the bird’s body at the moment of impact.

The impact (or slamming) dynamics of an object has been extensively studied by physicists or engineers in applications of marine craft hydrodynamics15,16,17,18,19. First, von Kármán15 and Wagner16 showed analytical solutions of the water-entry problem using a potential flow. A difference between these two models is whether the local uprise of the water on the impacting body is considered or not. Since then, there have been some advancements in modeling impact dynamics further17,18,19. The central idea of widely used theories for a water-entring body is to calculate the pressure and impact force from the velocity potential associated with a moving plate, which could be also useful to understand the impact force of diving animals.

In this study, we describe three different animal behaviors (i.e., swimming, flying, and diving) using one central mathematical framework based on a potential theory. From the mathematical point-of-view, the difference between swimming, flying, and diving is whether a body is moving while fully immersed in a fluid or at the interface. We are able to predict and quantify the thrust or impact force of animals from potential flow theory in order to understand the various locomotion modes. First, we review the previously known complex potential of a thin plate. Then, the impact force in diving and the thrust force in swimming or flying are estimated from the potential theory. These results are further verified with biological data. Finally, we conclude how this unified theory can explain diverse animal locomotion.

## Results

### Complex potential of a thin plate

We consider a system that a thin plate vertically moves at a speed of V as illustrated in Fig. 2. This canonical example has been already described in many books (e.g., pp. 336–372 in ref.20, pp. 136–139 in ref.21, pp. 304–309 in ref.18) and published articles22,23,24,25,26). However, we recap this classical potential model of the moving plate here, which will facilitate the move to the impact and thrust force calculations in the next section.

The complex potential, $$\Phi$$, for a moving plate is given as

\begin{aligned} \Phi = \phi + i \psi = -i Vz + i V (z^2 - c^2)^{1/2}, \end{aligned}
(1)

where $$\phi$$ is the velocity potential, i is the imaginary unit, $$\psi$$ is the streamfunction, $$z = x +i y$$ is the complex domain, and c is the half width of the plate. The complex velocity can be obtained by taking a derivative on $$\Phi$$ with respect to z. Here, the complex velocity is defined as $$u-iv$$ where u is the x-component velocity and v is the y-component velocity. The complex velocity from the above complex potential becomes

\begin{aligned} \frac{d \Phi }{dz} = u - i v = - iV + i V \frac{z}{(z^2-c^2)^{1/2}} \,. \end{aligned}
(2)

To further calculate other quantities in this canonical problem, we employ elliptical coordinates as

\begin{aligned} z-c = r_1 \, e^{i\theta _1}, ~~ z +c = r_2 \, e^{i\theta _2}. \end{aligned}
(3)

The first (or second) expression is based on the coordinate from the right end (or left end) of the plate. Then, the denominator in the second term of the complex velocity, Eq. (2), becomes

\begin{aligned} (z^2-c^2)^{1/2} = \sqrt{r_1 r_2} \, e^{i (\theta _1+\theta _2)/2}. \end{aligned}
(4)

This Cartesian-to-polar transformation as also illustrated in Fig. 2 is useful to check the boundary conditions in the following section.

#### Boundary conditions

Using the above complex potential and velocity, we can check whether this model satisfies the boundary conditions, i.e., $$v = V$$ on the plate and $$u=0$$ along the horizontal line outside the plate. First, we will evaluate the velocities at four different regions: the top and bottom sides of the plate, and the right and left horizon lines of the plate.

1. (i)

To evaluate quantities on the top of the plate, we choose the polar-coordinate parameters as ($$\theta _1=\pi$$, $$\theta _2=0^+$$; $$y = 0^+$$ & $$|x| < c$$). Then, part of the second term in Eqs. (1) or (2) becomes

\begin{aligned} (z^2-c^2)^{1/2} = \sqrt{r_1 r_2}\, e^{i \pi /2} = i (c^2-x^2)^{1/2}. \end{aligned}
(5)

Then, the velocity potential turns into

\begin{aligned} \phi _+= \mathrm {Re}[-iVz - V(c^2-x^2)^{1/2}] = V [y- (c^2-x^2)^{1/2}]. \end{aligned}
(6)
2. (ii)

On the bottom side of the plate ($$\theta _1=-\pi$$, $$\theta _2=0^-$$; $$y = 0^-$$& $$|x| < c$$), the same term is

\begin{aligned} (z^2-c^2)^{1/2} = \sqrt{r_1 r_2}\, e^{-i \pi /2} = - i (c^2-x^2)^{1/2}. \end{aligned}
(7)

Then, the velocity potential becomes

\begin{aligned} \phi _-= \mathrm {Re}[-iVz + V(c^2-x^2)^{1/2}] = V [ y+(c^2-x^2)^{1/2}]. \end{aligned}
(8)

This velocity potential is quite similar to the one on the top, Eq. (6), except for the sign of the second term.

In a similar way, we can calculate the complex velocity for both sides of the plate as

\begin{aligned} \left. \frac{d \Phi }{dz} \right| _\pm = \left. (u - i v)\right| _{\pm } = - iV \pm V \frac{x}{(c^2-x^2)^{1/2}}. \end{aligned}
(9)

Here, the subscript sign ($$+$$ or −) represents the top or bottom of the plate, respectively. The first term on the right hand side is a pure imaginary number, whereas the second term is a pure real number on the plate.

Equation (9) shows the vertical velocity of the plate as

\begin{aligned} v_\pm = - \mathrm {Im}\left[ \frac{d \Phi }{dz} \right] = V, \end{aligned}
(10)

which satisfies the kinematic boundary condition. The x component of velocity does not vanish as

\begin{aligned} u_{\pm } = \mathrm {Re}\left[ \frac{d \Phi }{dz} \right] = \pm V \frac{x}{(c^2-x^2)^{1/2}}. \end{aligned}
(11)

Like other potential flows, we do not expect no slip condition at a solid boundary (i.e., the plate), but allow a tangential slip. Equation (11) shows a diverging flow above and a converging flow beneath the plate when $$V>0$$, which indicates the fluid flowing around the plate.

3. (iii)

On the right horizon line of the plate ($$\theta _1=0$$, $$\theta _2=0$$; $$y = 0$$ & $$x > c$$),

\begin{aligned} (z^2-c^2)^{1/2} = \sqrt{r_1 r_2} = \sqrt{(x-c)(x+c)} = (x^2-c^2)^{1/2}. \end{aligned}
(12)
4. (iv)

On the left horizon line of the plate ($$\theta _1=-\pi$$, $$\theta _2=-\pi$$; $$y = 0$$ & $$x < -c$$),

\begin{aligned} (z^2-c^2)^{1/2} = (x^2-c^2)^{1/2}. \end{aligned}
(13)

Here, in both (iii) and (iv) regions, the complex velocity becomes

\begin{aligned} \left. \frac{d \Phi }{dz} \right| _{y=0, |x|>c} = u - i v = - iV + iV \frac{x}{(x^2-c^2)^{1/2}}. \end{aligned}
(14)

Now, let us check whether it satisfies boundary conditions outside the plate. The boundary conditions outside are only non-zero vertical velocity and an equipotential horizon, which means that both x-component velocity and velocity potential are zero.

\begin{aligned} \left. u \right| _{y=0, |x|>c}= & {} - \mathrm {Re}\left[ \frac{d \Phi }{dz} \right] = 0\, \nonumber \\ \left. \phi \right| _{y=0, |x|>c}= & {} \mathrm {Re}\left[ \Phi \right] =\mathrm {Re}\left[ i \left( -1 \pm \frac{x}{(x^2-c^2)^{1/2}} \right) \right] = 0. \nonumber \\ \end{aligned}
(15)

This condition of a constant velocity potential ($$\phi =0$$) along $$y=0$$ allows us to consider the horizontal surface as the free surface. Similarly, the free surface is modeled as an equipotential line of a velocity potential in many cases (e.g., p. 363 in27,28,29,30,31). This fact is useful to describe the case of “(I) Impacting plate” as in Fig. 1.

#### Circulation and vortex from a thin plate

Flow visualizations around locomoting animals have revealed that vortices are shed from the undulating body, especially near the tip of the fins or wings7,32,33,34,35,36. Moreover, the shed vortices in a fluid are connected with each other like a series of chains37. Hence, quantifying the vortices from the body might be useful to characterize fluid flows around the animals. In fluid mechanics analysis, circulation instead of vortices is widely used as a measure of rotation. The circulation is defined as an integral of the vorticity over an area:

\begin{aligned} \Gamma = \iint \varvec{\omega } \cdot {\hat{b}} \, dS = \iint \nabla \times {\mathbf{v}} dydx, \end{aligned}
(16)

where $$\varvec{\omega }$$ is the vorticity vector and $${\hat{b}}$$ is the unit vector $$(\equiv {\hat{x}} \times {\hat{y}})$$. Using Eqs. (10) and (11), the circulation around the plate can be further simplified as

\begin{aligned} \Gamma= & {} \iint (\partial _x v - \partial _y u) dydx = - \int ^c_{-c} (u_+ - u_-) dx \nonumber \\= & {} -2V \int ^c_{-c} \frac{x}{(c^2 - x^2)^{1/2}} dx. \end{aligned}
(17)

If it is integrated over the entire plate, the total circulation becomes zero (i.e., Kelvin’s theorem). However, it does not mean no vortex shed from the body. There would be equal and opposite signs of vortices shed from the edges. If we consider only the right edge of the plate, the circulation becomes

\begin{aligned} \Gamma ^{\rm {(right)}} = -2V \int _0^c \frac{x}{(c^2 - x^2)^{1/2}} dx = - 2 V c. \end{aligned}
(18)

Similarly, the circulation on the left side will be 2Vc. Hence, locomoting animals in a fluid (air or water) shed vortices with a circulation of 2Vc on each end-side of the appendage (i.e., wings for birds, and fins or flukes for aquatic animals).

#### Pressure difference across a thin plate

To calculate the force, we need to know the pressure on the plate first. The pressure from the unsteady Bernoulli equation is given as

\begin{aligned} p_\pm = p_0 - \rho \dfrac{\partial \phi _\pm }{\partial t} - \frac{\rho }{2} \left| \nabla \phi _\pm \right| ^2, \end{aligned}
(19)

where $$p_\pm$$ is the pressure above ($$+$$) and below (−) the plate and $$p_0$$ is the reference pressure (i.e., the atmospheric pressure for aerial fliers and hydrostatic pressure for aquatic swimmers).

(I) Using the velocity potential given in Eqs. (6) and (8), the pressure difference for a submerged plate moving in a fluid is given as

\begin{aligned} \frac{p_- -p_+}{\rho } = -\frac{\partial }{\partial t} (\phi _- - \phi _+) = -2 \frac{\partial }{\partial t} \left[ V(c^2-x^2)^{1/2} \right] . \end{aligned}
(20)

This pressure difference will be used to estimate the force generated by the undulating wings or fins in “Thrust force in swimming and flying”. It is worth noting that the last term in Eq. (19) does not contribute to the pressure on the plate at all due to the square of the velocity, which is the same on both the top and bottom sides of the plate.

(II) For the thin plate impacting a free surface ($$y= 0$$), the pressure on the upper side of the plate stays close to the atmopheric pressure (i.e., $$p_+ \simeq p_0$$) since the air density is so small compared to the water density. Hence, most pressure is built up on the water side not on the air side. Then, the pressure difference becomes

\begin{aligned} \frac{p_- -p_0}{\rho }= & {} -\frac{\partial }{\partial t} \left[ V(c^2-x^2)^{1/2} \right] - \frac{1}{2} \left[ V^2 + \frac{(Vx)^2 }{c^2-x^2} \right] \nonumber \\= & {} -\frac{\partial }{\partial t} \left[ V(c^2-x^2)^{1/2} \right] - \frac{1}{2} V^2 \left[ \frac{c^2 }{c^2-x^2} \right] . \end{aligned}
(21)

This second term is from the steady inertia term (i.e., the square of the velocity), which becomes singular at the ends of the plate. To avoid this singularity, there have been discussions in the previous literature29,38,39. In this present study, we will omit the last term for convenience.

### Impact force in diving

Some animals plunge-dive into water at high speeds: aquatic animals11 and aerial birds14,41,42,43. Most animals have a stream-lined body like a spheroidal head front for aquatic animals or a conical beak for birds, which might help reducing the likelihood of injury under high dynamic loadings while diving. To understand and quantify the impact force (i.e., dynamic loading) on the body, we will approximate the diving motion as a plate with its width increasing as the body penetrates the free surface and solve the potential flow as illustrated in Fig. 3. A similar trick has been used in the case of hull slamming problems15,16,17.

By integrating Eq. (21) over the plate, the impact force is obtained as

\begin{aligned} \frac{F}{\rho }= & {} \int _{-c}^c \frac{(p_- -p_0)}{\rho } dx \nonumber \\= & {} -V \frac{dc}{dt} \int _{-c}^c \frac{c}{(c^2-x^2)^{1/2}} dx - \frac{dV}{dt} \int _{-c}^c (c^2-x^2)^{1/2} dx \nonumber \\= & {} -\pi V \frac{d}{dt} \left( \frac{c^2}{2} \right) - {\pi } \frac{c^2}{2} \frac{dV}{dt}, \end{aligned}
(22)

where $$p_0$$ is the atmospheric pressure above the plate. Two integral identities ($$\int _{-c}^{c} c/\sqrt{c^2-x^2}\, dx = \pi c$$ and $$\int _{-c}^{c} \sqrt{c^2-x^2} dx = \pi c^2/2$$) are used to simplify the expression. The first term refers to the slamming force, which is positive for a downward plate ($$V<0$$). The second term represents the added-mass effect due to the body’s acceleration or deceleration.

By denoting the added mass of the 2D plate as $$M_{\rm {added}} =\rho {\pi } c^2$$, we can rewrite the above equation into a simple and generalized form as

\begin{aligned} F= & {} -V \frac{d}{dt} \left( \frac{M_{\rm {added}}}{2} \right) - \frac{M_{\rm {added}}}{2} \frac{d V}{dt} \nonumber \\= & {} -V \frac{d}{dt} \left( \frac{M_{\rm {added}}}{2} \right),~ \mathrm { with~constant~velocity}. \end{aligned}
(23)

This equation expressed in terms of the added mass is not limited to a flat plate anymore. Once an added mass value is known, we can calculate the impact force for any body shape. To briefly explain a formal way to evaluate the added mass, the velocity potential should be solved first and then integrated over the surface as $$M_{\rm {added}} = \rho \, U^{-2} \int \phi (\partial \phi /\partial n) d S$$ where n is the coordinate normal to the surface and S is the body surface (see p. 123 in27). Figure 4 summarizes the added mass of both 2D and 3D objects from previous literature.

Next, for simplicity, we approximate all aquatic animals as prolate spheroids; however the other shapes in Fig. 4 may be useful for future references. The added mass of a prolate spheroid is given as

\begin{aligned} M_{\rm {added}} = K_i \frac{4}{3} \pi \rho a b^2, \end{aligned}
(24)

where the coefficient $$K_i$$, has different values depending on the direction of the motion. We consider two moving directions: the axial direction along the major axis and the lateral direction along the minor axis. In terms of animal diving, the axial directional dive corresponds to a head-first dive, whereas the lateral dive corresponds to a belly-first dive. Coefficients, $$K_{\rm {Axial}}$$ and $$K_{\rm {Lateral}}$$, are functions of the eccentricity, $$e \equiv \sqrt{1-b^2/a^2}$$: the formulas are given in the caption of Fig. 4.

Equation (23) with Eq. (24) allows us to calculate the maximum impact force during diving. As the body penetrates the free surface, the water-contact depth along the diving direction increases over time. For simplicity, we can consider only “a” as a time-dependent variable for the axial dive or only “b” as a time-dependent variable for the lateral dive. Then, the diving speed, V, is approximated as da/dt for the axial dive and db/dt for the lateral dive. The maximum impact force occurs when the cross-sectional area on the free surface reaches its maximum. Therefore, the maximum impact force is estimated as

\begin{aligned} F_{\rm {Head-first}}\simeq & {} - K_{\rm {Axial}} \frac{2}{3} \pi \rho b^2 \, V^2, \nonumber \\ F_{\rm {Belly-first}}\simeq & {} - K_{\rm {Lateral}} \frac{4}{3} \pi \rho a b\, V^2. \end{aligned}
(25)

Next, we calculate the impact pressure acting on the body, which is defined as the maximum force divided by its wetted surface area. At the moment that an animal reaches its maximum impact force, only half of the body is in contact with water. So, the wetted surface area is approximated as half of the total surface area: $$\pi b^2 (1 + a/(b e) \cdot \arcsin {e} )$$ where the eccentricity is $$e\equiv \sqrt{1-b^2/a^2}$$. The impact pressure at the moment of reaching the maximum force is estimated as

\begin{aligned} P_{\rm {Head-first}}\simeq & {} - K_{\rm {Axial}} \frac{2}{3} \rho \left( 1 + \frac{a}{b e} \arcsin {e} \right) ^{-1} V^2, \nonumber \\ P_{\rm {Belly-first}}\simeq & {} - K_{\rm {Lateral}} \frac{4}{3} \rho \left( 1 + \frac{a}{b e} \arcsin {e} \right) ^{-1} \left( \frac{a}{b} \right) V^2. \end{aligned}
(26)

Figure 5 shows the maximum impact force and impact pressure versus animal mass. Green symbols represent the lateral impact (i.e., belly-first dive), whereas blue symbols represent the axial impact (i.e., head-first dive). In most cases, we do not have quantitative measurements of diving speed or postures. Thus, we estimate the diving speed from the jumping height as $$V^2 = 2 g H$$ as in11. In Fig. 5a, the belly-first dive produces more impact force than the head-first dive, which is quite intuitive due to the difference in the cross-sectional area depending on the diving direction. However, the impact pressure does not increase much with the animal mass as in Fig. 5b. This constant pressure would indicate that all animals feel a similar level of pressure and safely dive into water regardless of their body weight or length.

Furthermore, we develop a scaling argument for animal diving. From the previous study11, the jumping height is predicted as $$H \propto L^{2/3}$$. If the body area is assumed to be the square of the characteristic length $$L^2$$, the maximum impact force becomes proportional to $$F_{\rm {Impact}} \propto H L^2 = L^{8/3}$$. Using an allometric relation of the animal mass to the characteristic length as $$M \propto L^3$$, we get the impact force as $$F_{\rm {Impact}} \propto M^{8/9}$$. In terms of the scaling law of the impact pressure, the impact force ($$M^{8/9}$$) should be divided by its surface area ($$M^{2/3}$$). Hence, we anticipate that the impact pressure does not show any strong dependence on mass ($$M^{2/9} \simeq M^{0.22}$$). However, we observe a slightly decreasing trend in the impact pressure for large animals (see Fig. 5b). It might be due to two reasons. First, large diving animals have typically a more streamlined body shape than small aquatic animals as shown in the inset of Fig. 5a. If an animal has a streamlined body rather than a spherical shape, then the total surface area gets larger than that of a spherical body at a given volume or mass. In other word, the surface area does not follow $$M^{2/3}$$ strictly. Second, large jumping animals use a different jumping strategy called momentum jumping, whereas small animals use an impulsive jumping strategy11. Hence, our simple allometric scaling law disregarding the details of animal behavior or shape does not match with the data very well, especially for large animals.

### Thrust force in swimming and flying

Flying and swimming locomotions are induced by flapping motions, whose underlying mechanism is based on a similar fluid-mechanics principle with the diving motion. The pressure gradient developed across the thin appendage generates the thrust force as illustrated in Fig. 6, which is a key element to understand flying and swimming motions. There have been extensive studies44,45,46 to find a unifying scaling for swimming and flying motions. In contrast to the previous studies of scaling laws, our model roots in the potential theory to describe the animal’s diving and locomotion.

First, we assume that a plate is moving in a sinusoidal fashion as $$V = A (2 \pi f) e^{i 2 \pi f t}$$ where A and f are the amplitude and frequency of flapping. From Eq. (20), the thrust force per unit length can be calculated as

\begin{aligned} F= & {} \int _{-c}^c (p_- - p_+) dx = - \rho \int _{-c}^c \frac{\partial }{\partial t} (\phi _- - \phi _+) dx \nonumber \\= & {} - \rho \pi c^2 \frac{d V}{dt} = -i \rho \pi c^2 A (2 \pi f)^2 e^{i 2 \pi f t} \nonumber \\\propto & {} \rho \, c^2 A f^2, \text { for flapping appendage per unit length.} \end{aligned}
(27)

Strictly speaking, the time-averaged thrust force over a period will be zero if the pressure is purely periodic. However, real flying animals perform upstroke and downstroke in an asymmetric way by decreasing the angle of attack of the wing (or fin) and/or folding the wing (or fin) during the upstroke6,47,48. To take into account this effect, there should be an unknown non-zero prefactor less than one for the thrust force. Instead of finding or modeling details of the unknown prefactor, we approximate the total thrust force over the wing as $$\rho c^2 A f^2 L$$ where c is the cord half-length and L is the length of the wing span (or fin length). This thrust force becomes the lift force for fliers and the forward propulsive force for swimmers.

#### Flier: lift force balancing with weight

For fliers, the force generated by the wings is used to lift their own body as well as propel forward. However, most force is allocated to lifting the body since the animal weight is typically higher than the aerodynamic drag of the forward motion, especially for large animals. Hence, we assume that the force generated by the flapping wings balances with its own weight.

\begin{aligned} F_{\rm {weight}} (= \mathrm {Mass} \cdot g) \propto \rho \, c^2 A f^2 L. \end{aligned}
(28)

Figure 7a shows that the generated lift force based on our potential flow model is proportional to the animal weight quite well. Blue symbols are from bat species49 and green symbols are from birds50. Two solid lines represent our theoretical prediction of Eq. (28) with two different prefactors: one for bats and the other for birds. One possible reason of having the two prefactors is that bats and birds evolved flight independently51, which indicates that we do not expect one single curve to collapse all the data of bats and birds. We also observe that small fliers significantly deviate from our predicted linear lines, which indicates that small animals use more or less flapping-induced force to support their weight. This deviation is presumably due to some of the aerodynamic force spent for forward flight or the additional force gained from surrounding flows to compensate for its own weight.

#### Swimmer: forward-flying force balancing with drag

Aquatic animals do not need to support their body in water since their body density is close to water density. Instead, the thrust force generated from the fin is used to swim forward. There are two swimming regimes depending on the Reynolds number (i.e., a ratio of inertia to viscous force); defined as $$\mathrm {Re} = U_{\rm {swim}}L/\nu$$ where the characteristic velocity ($$U_{\rm {swim}}$$) is the swimming velocity, the characteristic length (L) is the body length, and $$\nu$$ is the kinematic viscosity of the fluid. The kinematic viscosity is about $$1.00\times 10^{-6}$$ $$\hbox {m}^2$$/s for freshwater and $$1.05\times 10^{-6}$$ $$\hbox {m}^2$$/s for seawater at the temperature of 20 $$^\circ$$C.

Most aquatic animals are bigger than a few centimeters and swim at about a few times its body length per second. So, their corresponding Reynolds number is more than a few thousands. At such high Reynolds numbers, the thrust force ($$\rho \, c^2 A f^2 L$$) balances with the form drag ($$F_{\rm {form\,drag}} \sim \frac{1}{2} \rho U_{\rm {swim}}^2L c$$). Then, the swimming speed for animals is given as

\begin{aligned} U_{\rm {swim}} \propto f \, (Ac)^{1/2}. \end{aligned}
(29)

This indicates that the swimming speed is proportional to the flapping frequency, f, and the geometric-mean stroke length, $$(Ac)^{1/2}$$. Gazoola et al.46 suggested a slightly different scaling as $$U_{\rm {swim}} \propto f A$$ based on scaling arguments. However, our prediction from the potential flow model results in the dependence on its cord length (c), which does not show up in the other model46. For small swimming animals (typically larvae smaller than a few centimeters), the skin drag might be dominant over the form drag as $$F_{\rm {skin\,drag}} \sim \rho U_{\rm {swim}}^2 \, c L \sqrt{ {\nu }/{U_{\rm {swim}} L}} = \rho (\nu L)^{1/2} U_{\rm {swim}}^{3/2}c$$. Then, balancing it with the thrust force, we get $$U_{\rm {swim}} \propto (Ac)^{2/3} f^{4/3} L^{1/3} \nu ^{-1/3}$$. However, there are not many aquatic animals belonging to this regime to confirm this prediction.

Our prediction shows that the swimming speed depends on the animal’s stroke amplitude (A), cord half-length (c), and frequency (f) as in Eq. (29). Figure 7b shows the relation between the swimming speed and the predicted speed of fishes ranging from tadpoles to whales (52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71; Many of the references were adapted from72). This linear relation can be associated with the Strouhal number, i.e., a ratio of unsteady to steady inertia. Then, we define the Strouhal number as

\begin{aligned} \mathrm {Strouhal~number} = \frac{f (Ac)^{1/2}}{U_{\rm {swim}} } = \mathrm {const.} \end{aligned}
(30)

Taylor et al.45 also showed the constant Strouhal number of locomoting animals using a slight different Strouhal number definition ($$f A/U_{\rm {swim}}$$). However, in terms of the order of magnitude, the Strouhal number in the previous studies45,72 is between 0.2 and 0.5, which is very close to what we observed in our study ($$0.31 \pm 0.19$$).

## Conclusions

In this paper, we reviewed the previously known potential model of a plate moving in a fluid using a complex potential and provided analogies to swimming, flying, and diving of animals. Additionally, using the unsteady Bernoulli equation, we calculated the circulation, pressure, and force on a locomotion body. The calculated force was decoded into the impact force for diving animals at the free surface or the thrust force for swimming or flying animals immersed in a fluid. Our prediction explained almost constant pressure on diving animals, the lift force balancing with weight for fliers, and the swimming speed as a result of thrust force balancing with drag for swimmers. Furthermore, measured kinematic data from various locomotion modes of both aquatic and flying animals support our theoretical predictions.

It is worth noting that there are three seminal works in analytical models for swimming animals; Wu’s model2 is based on a 2D potential (the same as presented in this paper) focusing on an undulating surface, whose results can be applied for animals swimming in a unbounded fluid. Lighthill’s model1 is based on a power balance of an elongated body; the total power is composed of the thrust-related work, the rate of the kinetic energy of wake at the trailing edge, and the rate of the kinetic energy ahead of the trailing edge. A more comparable model to our predicted swimming speed would be the work done by Gazzola et al.46. They developed a scaling model to describe both swimmers and fliers. Their scaling argument expression is quite similar to ours, but the main difference is that our model is based on the 2D potential flow and shows the importance of the width of flappers or the cord length of wings. Moreover, our study is unique as the first attempt to mathematically unify three distinct animal behaviors: swimming, flying, and diving. This calculation can be also useful in many examples of fluid-organism interactions in nature like a fluttering leaf73,74,75, spore/particle dispersal by a leaf’s motion76, a falling seed77,78, an animal lapping as a plate-like tongue moving out of the water79,80,81, and others.

## Material and methods

Animal data are obtained from previous publications49,50,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71 . For flying data, 23 bats in49 and 16 birds in50 are used. Other required data for Fig. 7a are the animal weight, wing span, cord length, and flapping amplitude. For bats, the flapping amplitude is not directly given in the paper with kinematic data, so we estimate it from the stroke angle and wing span. For the cord length, we approximately evaluate it as the wing area divided by the wing span. For fishes, we use the frequency, amplitude, cord length, and swimming speed of 32 fish species. Some data points are not explicitly given in text. In that case, we extract the value from the graph or best fitted lines. Cord lengths (i.e., fish or fluke width) of several species were not given in the same paper that described the kinematics. Then, we find the cord length from other papers of the same fish species. All these details are marked in Excel files uploaded in DOI:10.17605/OSF.IO/46SFV .

## Data availability

All matlab codes and data are freely available in https://doi.org/10.17605/OSF.IO/46SFV.

## References

1. 1.

Lighthill, J. Mathematical Biofluiddynamics (SIAM, 1973).

2. 2.

Wu, T. Y. Fish swimming and bird/insect flight. Annu. Rev. Fluid Mech. 43, 25 (2011).

3. 3.

Vogel, S. Life in Moving Fluids: The Physical Biology of Flow (Princeton University Press, 1994).

4. 4.

Biewener, A. & Patek, S. Animal Locomotion (Oxford University Press, 2018).

5. 5.

Alexander, R. M. N. Principles of Animal Locomotion (Princeton University Press, 2013).

6. 6.

Fish, F. E. & Lauder, G. V. Passive and active flow control by swimming fishes and mammals. Annu. Rev. Fluid Mech 38, 193 (2006).

7. 7.

Chin, D. D. & Lentink, D. Flapping wing aerodynamics: From insects to vertebrates. J. Exp. Biol. 219, 920 (2016).

8. 8.

Pedley, T. J. Pulmonary fluid dynamics. Annu. Rev. Fluid Mech. 9, 229 (1977).

9. 9.

Lighthill, M. J. Note on the swimming of slender fish. J. Fluid Mech. 9, 305 (1960).

10. 10.

Wu, T. Hydromechanics of swimming propulsion. part 1. Swimming of a two-dimensional flexible plate at variable forward speeds in an inviscid fluid. J. Fluid Mech. 46, 337 (1971).

11. 11.

Chang, B. et al. Jumping dynamics of aquatic animals. J. R. Soc. Interface 16, 1 (2019).

12. 12.

Huin, N. Diving depths of white-chinned petrels. Condor 96, 1111 (1994).

13. 13.

Prince, P. A., Huin, N. & Weimerskirch, H. Diving depths of albatrosses. Antarct. Sci. 6, 353 (1994).

14. 14.

Eliason, C. M., Straker, L., Jung, S. & Hackett, S. J. Morphological innovation and biomechanical diversity in plunge-diving birds. Evolution 74, 1514 (2020).

15. 15.

Von Kármán, T. The impact on seaplane floats during landing. in NACA Technical Note No. NACA-TN-32. Vol. 321, 2 (1929).

16. 16.

Wagner, H. Über Stoß- und Gleitvorgänge an der Oberfläche von Flüssigkeiten. ZAMM J. Appl. Math. Mech./Z. Angew. Math. Mech. 12, 193 (1932).

17. 17.

Abrate, S. Hull slamming. Appl. Mech. Rev. 64, 060803 (2011).

18. 18.

Faltinsen, O. M. Hydrodynamics of High-Speed Marine Vehicles (Cambridge University Press, 2006).

19. 19.

Seddon, C. M. & Moatamedi, M. Review of water entry with applications to aerospace structures. Int. J. Impact Eng. 32, 1045 (2006).

20. 20.

Kochin, N. K., Kibel, I. A., & Roze, N. V. Theoretical Hydromechanics. 577. (Interscience Publishers, 1964).

21. 21.

Currie, I. G. Fundamental Mechanics of Fluids, 4th edn. 603. (CRC Press, 2016).

22. 22.

Sarpkaya, T. Morison’s Equation and the Wave Forces on Offshore Structures. in Technical Report-ADA137438 (1981).

23. 23.

Berklite, R. B. Added mass of submerged objects of arbitrary shape. in Thesis in Naval Postgraduate School (1972)

24. 24.

Patton, K. T. Hydrodynamic Mass of Bodies in a Fluid. in DTIC Report- USL-TM-933-351-64 (1964).

25. 25.

Brennen, C. E. A Review of Added Mass and Fluid Inertial Forces. in Technical Report-ADA110190 (1982).

26. 26.

Faltinsen, O. M., and R. Zhao. "Water entry of ship sections and axisymmetric bodies. AGARD FDP and Ukraine Institute of Hydromechanics Workshop on High-Speed Body Motion in Water, Kiev, Ukraine." AGARD Rep 818 (1997).

27. 27.

Lamb, H. Hydrodynamics. 6th edition 768 (Dover, 1945).

28. 28.

Louf, J. F. et al. Cavity ripple dynamics after pinch-off. J. Fluid Mech. 850, 611 (2018).

29. 29.

Vincent, L., Xiao, T., Yohann, D., Jung, S. & Kanso, E. Dynamics of water entry. J. Fluid Mech. 846, 508 (2018).

30. 30.

Kim, S., Wu, Z., Esmaili, E., Dombroskie, J. J. & Jung, S. How a raindrop gets shattered on biological surfaces. Proc. Natl. Acad. Sci. USA 117, 13901 (2020).

31. 31.

Jung, S., Kim, E., Shelley, M. J. & Zhang, J. Surface waves on a semitoroidal water ring. Phys. Fluids 19, 58105 (2007).

32. 32.

Hedrick, T. L., Tobalske, B. W. & Biewener, A. A. Estimates of circulation and gait change based on a three-dimensional kinematic analysis of flight in cockatiels (Nymphicus hollandicus) and ringed turtle-doves (Streptopelia risoria). J. Exp. Biol. 205, 1389 (2002).

33. 33.

Müller, U. K., Smit, J., Stamhuis, E. J. & Videler, J. J. How the body contributes to the wake in undulatory fish swimming: Flow fields of a swimming eel (Anguilla anguilla). J. Exp. Biol. 204, 2751 (2001).

34. 34.

Liao, J. C. A review of fish swimming mechanics and behaviour in altered flows. Philos. Trans. R. Soc. B Biol. Sci. 362, 1973 (2007).

35. 35.

Tytell, E. D. & Lauder, G. V. The hydrodynamics of eel swimming: I. Wake structure. J. Exp. Biol. 207, 1825 (2004).

36. 36.

Müller, U. K. Fish ’n flag. Science 302, 1511 (2003).

37. 37.

Dickinson, M. H. How animals move: An integrative view. Science 288, 100 (2000).

38. 38.

Korobkin, A. Analytical models of water impact. Eur. J. Appl.Math. 15, 821 (2004).

39. 39.

Zhao, R. & Faltinsen, O. Water entry of two-dimensional bodies. J. Fluid Mech. 246, 593 (1993).

40. 40.

Baldwin, J. L. Vertical water entry of cones. in Technical report-AD723821 (1971).

41. 41.

Lee, D. N. & Reddish, P. E. Plummeting gannets: A paradigm of ecological optics. Nature 293, 293 (1981).

42. 42.

Ropert-Coudert, Y., Grémillet, D., Ryan, P., Kato, P., Naito, Y. & Le Maho, Y. Between air and water: The plunge dive of the Cape Gannet Morus capensis. Ibis 146, 281 (2004)

43. 43.

Chang, B. et al. How seabirds plunge-dive without injuries. Proc. Natl. Acad. Sci. 113, 12006 (2016).

44. 44.

Bejan, A. Unifying constructal theory for scale effects in running, swimming and flying. J. Exp. Biol. 209, 238 (2006).

45. 45.

Taylor, G. K., Nudds, R. L. & Thomas, A. L. Flying and swimming animals cruise at a Strouhal number tuned for high power efficiency. Nature 425, 707 (2003).

46. 46.

Gazzola, M., Argentina, M. & Mahadevan, L. Scaling macroscopic aquatic locomotion. Nat. Phys. 10, 758 (2014).

47. 47.

Tobalske, B. W. Biomechanics of bird flight. J. Exp. Biol. 210, 3135 (2007).

48. 48.

Liu, G. et al. Computational analysis of vortex dynamics and performance enhancement due to body-fin and fin-fin interactions in fish-like locomotion. J. Fluid Mech. 829, 65 (2017).

49. 49.

Bullen, R. D. & McKenzie, N. L. Scaling bat wingbeat frequency and amplitude. J. Exp. Biol. 205, 2615 (2002).

50. 50.

Pennycuick, C. J. Speeds and wingbeat frequencies of migrating birds compared with calculated benchmarks. J. Exp. Biol. 204, 3283 (2001).

51. 51.

Hedenström, A., Johansson, L. C. & Spedding, G. R. Bird or bat: Comparing airframe design and flight performance. Bioinspir. Biomimet. 4, 015001 (2009).

52. 52.

Fish, F. E. Comparative kinematics and hydrodynamics of odontocete cetaceans: Morphological and ecological correlates with swimming performance. J. Exp. Biol. 201, 2867 (1998).

53. 53.

Kojeszewski, T. & Fish, F. E. Swimming kinematics of the Florida manatee (Trichechus manatus latirostris): Hydrodynamic analysis of an undulatory mammalian swimmer. J. Exp. Biol. 210, 2411 (2007).

54. 54.

Haney, J. C., Stone, A. E. Seabird foraging tactics and water clarity: Are plunge divers really in the clear? Mar. Ecol. Prog. Ser. 49, 1 (1988).

55. 55.

Webb, P. W. & Keyes, R. S. Swimming kinematics of sharks. Fish. Bull. 80, 803 (1982).

56. 56.

Videler, B. Y. J. J. & Hess, F. Fast continuous swimming of two pelagic predators, Saithe (Pollachius virens) and Mackerel (Scomber scombrus): A kinematic analysis. J. Exp. Biol. 109, 209 (1984).

57. 57.

Donley, J. M. & Dickson, K. A. Swimming kinematics of juvenile kawakawa tuna (Euthynnus affinis) and chub mackerel (Scomber japonicus). J. Exp. Biol. 203, 3103 (2000).

58. 58.

Wardle, C. S., Videler, J. J., Arimoto, T., Franco, J. M., & He, P. The muscle twitch and the maximum swimming speed of giant bluefin tuna, Thunnus thynnus L. J. Fish Biol. 35, 129 (1989)

59. 59.

Dowis, H. J., Sepulveda, C. A., Graham, J. B. & Dickson, K. A. Swimming performance studies on the eastern Pacific bonito Sarda chiliensis, a close relative of the tunas (family Scombridae): II. Kinematics. J. Exp. Biol. 206, 2749 (2003).

60. 60.

Dewar, H., Graham, J. Studies of tropical tuna swimming performance in a large water tunnel. II. Thermoregulation. J. Exp. Biol. 192, 45 (1994).

61. 61.

Kayan, V. P., Kozlov, L. F. & Pyatetskii, V. E. Kinematic characteristics of the swimming of certain aquatic animals. Fluid Dyn. 13, 641 (1978).

62. 62.

Webb, P. W., Kostecki, P. T. & Stevens, E. D. The effect of size and swimming speed on locomotor kinematics of rainbow trout. J. Exp. Biol. 109, 77 (1984).

63. 63.

Webb, P. W. Steady swimming kinematics of tiger musky, an esociform accelerator, and rainbow trout, a generalist cruiser. J. Exp. Biol. 138, 51 (1988).

64. 64.

Webb, P. W. Effects of partial caudal fin amputation on the kinematics and metabolic rate of underyearling sockeye salmon (Oncorhynchus nerka) at steady swimming speeds. J. Exp. Biol. 59, 565 (1973).

65. 65.

Webb, P. W. Kinematics of plaice. Pleuronectes platessa, and cod, Gadus morhua, swimming near the bottom. J. Exp. Biol. 205, 2125 (2002).

66. 66.

Webb, P. W. Kinematics of lake sturgeon. Acipenser fulvescens, at cruising speeds. Can. J. Zool. 64, 2137 (1986).

67. 67.

Müller, U. K., Stamhuis, E. J. & Videler, J. J. Riding the waves: The role of the body wave in undulatory fish swimming. Integr. Comp. Biol. 42, 981 (2002).

68. 68.

Wassersug, R. J. & Von Seckendorf Hoff, K. The kinematics of swimming in anuran larvae. J. Exp. Biol. 119, 1 (1985).

69. 69.

Westneat, M. W. & Wainwright, S. A. Mechanical design for swimming: muscle, tendon, and bone. Fish Physiol. 19, 271 (2001).

70. 70.

Cengiz, Ö. On Maximum length record of the Chub Mackerel (Scomber japonicus Houttuyn, 1782) from Northern Aegean Sea (Turkey, eastern Mediterranean). Mar. Sci. Technol. Bull. 9, 173 (2020).

71. 71.

Tičina, V., Grubišić, L., Šegvić Bubić, T., & Katavić, I. Biometric characteristics of small Atlantic bluefin tuna (Thunnus thynnus, Linnaeus, 1758) of Mediterranean Sea origin. J. Appl. Ichthyol. 27, 971 (2011).

72. 72.

Eloy, C. Optimal Strouhal number for swimming animals. J. Fluids Struct. 30, 205 (2012) arXiv:1102.0223.

73. 73.

de Langre, E. Effects of wind on plants. Annu. Rev. Fluid Mech. 40, 141 (2008).

74. 74.

Gart, S., Mates, J. E., Megaridis, C. M. & Jung, S. Droplet impacting a cantilever: A leaf-raindrop system. Phys. Rev. Appl. 3, 044019 (2015).

75. 75.

Bhosale, Y., Esmaili, E., Bhar, K. & Jung, S. Bending, twisting and flapping leaf upon raindrop impact. Bioinspir. Biomimet. 15, 036007 (2020).

76. 76.

Kim, S., Park, H., Gruszewski, H. A., Schmale, D. G. & Jung, S. Vortex-induced dispersal of a plant pathogen by raindrop impact. Proc. Natl. Acad. Sci. 116, 4917 (2019).

77. 77.

Howe, F. & Smallwood, J. Ecology of seed dispersal. Annu. Rev. Ecol. Syst. 13 , 201 (1982).

78. 78.

Nave, G. K. et al. Wind dispersal of natural and biomimetic maple samaras. Biomimetics 6, 1 (2021) arXiv:2010.12553.

79. 79.

Reis, P. M., Jung, S., Aristoff, J. M. & Stocker, R. How cats lap: Water uptake by Felis catus. Science 330, 1231 (2010).

80. 80.

Gart, S., Socha, J. J., Vlachos, P. P. & Jung, S. Dogs lap using acceleration-driven open pumping. Proc. Natl. Acad. Sci. 112, 15798 (2015).

81. 81.

Jung, S. Pinch-off dynamics to elucidate animal lapping. Phys. Rev. Fluids. 6, 073102 (2021).

## Acknowledgements

The author thanks Ms. Jenna Rotheram for collecting data on swimming animals. This work was partially supported by the National Science Foundation Grant No. CBET-2002714 and CMMI-2042740.

## Author information

Authors

### Contributions

S.J. conceived the idea and wrote the manuscript.

### Corresponding author

Correspondence to Sunghwan Jung.

## Ethics declarations

### Competing interests

The author declares no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Jung, S. Swimming, flying, and diving behaviors from a unified 2D potential model. Sci Rep 11, 15984 (2021). https://doi.org/10.1038/s41598-021-94829-7