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 ellipsoidshaped 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 aspects^{1,2,3,4,5,6,7}. In water, most aquatic animals swim by flapping their fins or undulating part of the body^{3,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 largeamplitude elongated body theory to understand the swimming speed through balancing the power generated by an animal with the rate of kinetic energy in a fluid^{1,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 flow^{2,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 airwater interface (e.g., diving). There are two types of diving animals. The first type are aquatic animals that jump out of and then reenter the water^{11} (e.g., whales, dolphins, fishes, etc.). At the reentering moment, animals experience huge impact force as they hit the water at high speeds. The other type are birds plungediving into water from air. Several bird species exhibit highspeed diving into water as a hunting mechanism^{12,13,14}. These plungediving 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 highspeed plungediving 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 hydrodynamics^{15,16,17,18,19}. First, von Kármán^{15} and Wagner^{16} showed analytical solutions of the waterentry 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 further^{17,18,19}. The central idea of widely used theories for a waterentring 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 pointofview, 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 articles^{22,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
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 \(uiv\) where u is the xcomponent velocity and v is the ycomponent velocity. The complex velocity from the above complex potential becomes
To further calculate other quantities in this canonical problem, we employ elliptical coordinates as
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
This Cartesiantopolar 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.

(i)
To evaluate quantities on the top of the plate, we choose the polarcoordinate 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^2c^2)^{1/2} = \sqrt{r_1 r_2}\, e^{i \pi /2} = i (c^2x^2)^{1/2}. \end{aligned}$$(5)Then, the velocity potential turns into
$$\begin{aligned} \phi _+= \mathrm {Re}[iVz  V(c^2x^2)^{1/2}] = V [y (c^2x^2)^{1/2}]. \end{aligned}$$(6) 
(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^2c^2)^{1/2} = \sqrt{r_1 r_2}\, e^{i \pi /2} =  i (c^2x^2)^{1/2}. \end{aligned}$$(7)Then, the velocity potential becomes
$$\begin{aligned} \phi _= \mathrm {Re}[iVz + V(c^2x^2)^{1/2}] = V [ y+(c^2x^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^2x^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^2x^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.

(iii)
On the right horizon line of the plate (\(\theta _1=0\), \(\theta _2=0\); \(y = 0\) & \(x > c\)),
$$\begin{aligned} (z^2c^2)^{1/2} = \sqrt{r_1 r_2} = \sqrt{(xc)(x+c)} = (x^2c^2)^{1/2}. \end{aligned}$$(12) 
(iv)
On the left horizon line of the plate (\(\theta _1=\pi\), \(\theta _2=\pi\); \(y = 0\) & \(x < c\)),
$$\begin{aligned} (z^2c^2)^{1/2} = (x^2c^2)^{1/2}. \end{aligned}$$(13)
Here, in both (iii) and (iv) regions, the complex velocity becomes
Now, let us check whether it satisfies boundary conditions outside the plate. The boundary conditions outside are only nonzero vertical velocity and an equipotential horizon, which means that both xcomponent velocity and velocity potential are zero.
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 in^{27,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 wings^{7,32,33,34,35,36}. Moreover, the shed vortices in a fluid are connected with each other like a series of chains^{37}. 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:
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
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
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 endside 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
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
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
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 literature^{29,38,39}. In this present study, we will omit the last term for convenience.
Impact force in diving
Some animals plungedive into water at high speeds: aquatic animals^{11} and aerial birds^{14,41,42,43}. Most animals have a streamlined 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 problems^{15,16,17}.
By integrating Eq. (21) over the plate, the impact force is obtained as
where \(p_0\) is the atmospheric pressure above the plate. Two integral identities (\(\int _{c}^{c} c/\sqrt{c^2x^2}\, dx = \pi c\) and \(\int _{c}^{c} \sqrt{c^2x^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 addedmass 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
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 in^{27}). 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
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 headfirst dive, whereas the lateral dive corresponds to a bellyfirst dive. Coefficients, \(K_{\rm {Axial}}\) and \(K_{\rm {Lateral}}\), are functions of the eccentricity, \(e \equiv \sqrt{1b^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 watercontact depth along the diving direction increases over time. For simplicity, we can consider only “a” as a timedependent variable for the axial dive or only “b” as a timedependent 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 crosssectional area on the free surface reaches its maximum. Therefore, the maximum impact force is estimated as
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{1b^2/a^2}\). The impact pressure at the moment of reaching the maximum force is estimated as
Figure 5 shows the maximum impact force and impact pressure versus animal mass. Green symbols represent the lateral impact (i.e., bellyfirst dive), whereas blue symbols represent the axial impact (i.e., headfirst 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 in^{11}. In Fig. 5a, the bellyfirst dive produces more impact force than the headfirst dive, which is quite intuitive due to the difference in the crosssectional 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 study^{11}, 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 strategy^{11}. 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 fluidmechanics 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 studies^{44,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
Strictly speaking, the timeaveraged 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 upstroke^{6,47,48}. To take into account this effect, there should be an unknown nonzero 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 halflength 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.
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 species^{49} and green symbols are from birds^{50}. 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 independently^{51}, 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 flappinginduced 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: forwardflying 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
This indicates that the swimming speed is proportional to the flapping frequency, f, and the geometricmean 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 model^{46}. 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 halflength (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 from^{72}). 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
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 studies^{45,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 model^{2} 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 model^{1} is based on a power balance of an elongated body; the total power is composed of the thrustrelated 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 fluidorganism interactions in nature like a fluttering leaf^{73,74,75}, spore/particle dispersal by a leaf’s motion^{76}, a falling seed^{77,78}, an animal lapping as a platelike tongue moving out of the water^{79,80,81}, and others.
Material and methods
Animal data are obtained from previous publications^{49,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 in^{49} and 16 birds in^{50} 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.
Lighthill, J. Mathematical Biofluiddynamics (SIAM, 1973).
 2.
Wu, T. Y. Fish swimming and bird/insect flight. Annu. Rev. Fluid Mech. 43, 25 (2011).
 3.
Vogel, S. Life in Moving Fluids: The Physical Biology of Flow (Princeton University Press, 1994).
 4.
Biewener, A. & Patek, S. Animal Locomotion (Oxford University Press, 2018).
 5.
Alexander, R. M. N. Principles of Animal Locomotion (Princeton University Press, 2013).
 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.
Chin, D. D. & Lentink, D. Flapping wing aerodynamics: From insects to vertebrates. J. Exp. Biol. 219, 920 (2016).
 8.
Pedley, T. J. Pulmonary fluid dynamics. Annu. Rev. Fluid Mech. 9, 229 (1977).
 9.
Lighthill, M. J. Note on the swimming of slender fish. J. Fluid Mech. 9, 305 (1960).
 10.
Wu, T. Hydromechanics of swimming propulsion. part 1. Swimming of a twodimensional flexible plate at variable forward speeds in an inviscid fluid. J. Fluid Mech. 46, 337 (1971).
 11.
Chang, B. et al. Jumping dynamics of aquatic animals. J. R. Soc. Interface 16, 1 (2019).
 12.
Huin, N. Diving depths of whitechinned petrels. Condor 96, 1111 (1994).
 13.
Prince, P. A., Huin, N. & Weimerskirch, H. Diving depths of albatrosses. Antarct. Sci. 6, 353 (1994).
 14.
Eliason, C. M., Straker, L., Jung, S. & Hackett, S. J. Morphological innovation and biomechanical diversity in plungediving birds. Evolution 74, 1514 (2020).
 15.
Von Kármán, T. The impact on seaplane floats during landing. in NACA Technical Note No. NACATN32. Vol. 321, 2 (1929).
 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.
Abrate, S. Hull slamming. Appl. Mech. Rev. 64, 060803 (2011).
 18.
Faltinsen, O. M. Hydrodynamics of HighSpeed Marine Vehicles (Cambridge University Press, 2006).
 19.
Seddon, C. M. & Moatamedi, M. Review of water entry with applications to aerospace structures. Int. J. Impact Eng. 32, 1045 (2006).
 20.
Kochin, N. K., Kibel, I. A., & Roze, N. V. Theoretical Hydromechanics. 577. (Interscience Publishers, 1964).
 21.
Currie, I. G. Fundamental Mechanics of Fluids, 4th edn. 603. (CRC Press, 2016).
 22.
Sarpkaya, T. Morison’s Equation and the Wave Forces on Offshore Structures. in Technical ReportADA137438 (1981).
 23.
Berklite, R. B. Added mass of submerged objects of arbitrary shape. in Thesis in Naval Postgraduate School (1972)
 24.
Patton, K. T. Hydrodynamic Mass of Bodies in a Fluid. in DTIC Report USLTM93335164 (1964).
 25.
Brennen, C. E. A Review of Added Mass and Fluid Inertial Forces. in Technical ReportADA110190 (1982).
 26.
Faltinsen, O. M., and R. Zhao. "Water entry of ship sections and axisymmetric bodies. AGARD FDP and Ukraine Institute of Hydromechanics Workshop on HighSpeed Body Motion in Water, Kiev, Ukraine." AGARD Rep 818 (1997).
 27.
Lamb, H. Hydrodynamics. 6th edition 768 (Dover, 1945).
 28.
Louf, J. F. et al. Cavity ripple dynamics after pinchoff. J. Fluid Mech. 850, 611 (2018).
 29.
Vincent, L., Xiao, T., Yohann, D., Jung, S. & Kanso, E. Dynamics of water entry. J. Fluid Mech. 846, 508 (2018).
 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.
Jung, S., Kim, E., Shelley, M. J. & Zhang, J. Surface waves on a semitoroidal water ring. Phys. Fluids 19, 58105 (2007).
 32.
Hedrick, T. L., Tobalske, B. W. & Biewener, A. A. Estimates of circulation and gait change based on a threedimensional kinematic analysis of flight in cockatiels (Nymphicus hollandicus) and ringed turtledoves (Streptopelia risoria). J. Exp. Biol. 205, 1389 (2002).
 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.
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.
Tytell, E. D. & Lauder, G. V. The hydrodynamics of eel swimming: I. Wake structure. J. Exp. Biol. 207, 1825 (2004).
 36.
Müller, U. K. Fish ’n flag. Science 302, 1511 (2003).
 37.
Dickinson, M. H. How animals move: An integrative view. Science 288, 100 (2000).
 38.
Korobkin, A. Analytical models of water impact. Eur. J. Appl.Math. 15, 821 (2004).
 39.
Zhao, R. & Faltinsen, O. Water entry of twodimensional bodies. J. Fluid Mech. 246, 593 (1993).
 40.
Baldwin, J. L. Vertical water entry of cones. in Technical reportAD723821 (1971).
 41.
Lee, D. N. & Reddish, P. E. Plummeting gannets: A paradigm of ecological optics. Nature 293, 293 (1981).
 42.
RopertCoudert, 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.
Chang, B. et al. How seabirds plungedive without injuries. Proc. Natl. Acad. Sci. 113, 12006 (2016).
 44.
Bejan, A. Unifying constructal theory for scale effects in running, swimming and flying. J. Exp. Biol. 209, 238 (2006).
 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.
Gazzola, M., Argentina, M. & Mahadevan, L. Scaling macroscopic aquatic locomotion. Nat. Phys. 10, 758 (2014).
 47.
Tobalske, B. W. Biomechanics of bird flight. J. Exp. Biol. 210, 3135 (2007).
 48.
Liu, G. et al. Computational analysis of vortex dynamics and performance enhancement due to bodyfin and finfin interactions in fishlike locomotion. J. Fluid Mech. 829, 65 (2017).
 49.
Bullen, R. D. & McKenzie, N. L. Scaling bat wingbeat frequency and amplitude. J. Exp. Biol. 205, 2615 (2002).
 50.
Pennycuick, C. J. Speeds and wingbeat frequencies of migrating birds compared with calculated benchmarks. J. Exp. Biol. 204, 3283 (2001).
 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.
Fish, F. E. Comparative kinematics and hydrodynamics of odontocete cetaceans: Morphological and ecological correlates with swimming performance. J. Exp. Biol. 201, 2867 (1998).
 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.
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.
Webb, P. W. & Keyes, R. S. Swimming kinematics of sharks. Fish. Bull. 80, 803 (1982).
 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.
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.
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.
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.
Dewar, H., Graham, J. Studies of tropical tuna swimming performance in a large water tunnel. II. Thermoregulation. J. Exp. Biol. 192, 45 (1994).
 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.
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.
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.
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.
Webb, P. W. Kinematics of plaice. Pleuronectes platessa, and cod, Gadus morhua, swimming near the bottom. J. Exp. Biol. 205, 2125 (2002).
 66.
Webb, P. W. Kinematics of lake sturgeon. Acipenser fulvescens, at cruising speeds. Can. J. Zool. 64, 2137 (1986).
 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.
Wassersug, R. J. & Von Seckendorf Hoff, K. The kinematics of swimming in anuran larvae. J. Exp. Biol. 119, 1 (1985).
 69.
Westneat, M. W. & Wainwright, S. A. Mechanical design for swimming: muscle, tendon, and bone. Fish Physiol. 19, 271 (2001).
 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.
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.
Eloy, C. Optimal Strouhal number for swimming animals. J. Fluids Struct. 30, 205 (2012) arXiv:1102.0223.
 73.
de Langre, E. Effects of wind on plants. Annu. Rev. Fluid Mech. 40, 141 (2008).
 74.
Gart, S., Mates, J. E., Megaridis, C. M. & Jung, S. Droplet impacting a cantilever: A leafraindrop system. Phys. Rev. Appl. 3, 044019 (2015).
 75.
Bhosale, Y., Esmaili, E., Bhar, K. & Jung, S. Bending, twisting and flapping leaf upon raindrop impact. Bioinspir. Biomimet. 15, 036007 (2020).
 76.
Kim, S., Park, H., Gruszewski, H. A., Schmale, D. G. & Jung, S. Vortexinduced dispersal of a plant pathogen by raindrop impact. Proc. Natl. Acad. Sci. 116, 4917 (2019).
 77.
Howe, F. & Smallwood, J. Ecology of seed dispersal. Annu. Rev. Ecol. Syst. 13 , 201 (1982).
 78.
Nave, G. K. et al. Wind dispersal of natural and biomimetic maple samaras. Biomimetics 6, 1 (2021) arXiv:2010.12553.
 79.
Reis, P. M., Jung, S., Aristoff, J. M. & Stocker, R. How cats lap: Water uptake by Felis catus. Science 330, 1231 (2010).
 80.
Gart, S., Socha, J. J., Vlachos, P. P. & Jung, S. Dogs lap using accelerationdriven open pumping. Proc. Natl. Acad. Sci. 112, 15798 (2015).
 81.
Jung, S. Pinchoff 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. CBET2002714 and CMMI2042740.
Author information
Affiliations
Contributions
S.J. conceived the idea and wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The author declares no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Jung, S. Swimming, flying, and diving behaviors from a unified 2D potential model. Sci Rep 11, 15984 (2021). https://doi.org/10.1038/s41598021948297
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598021948297
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.