Universal Scaling Laws for a Generic Swimmer Model

Abstract

We have developed a minimal model of a swimmer without body deformation based on force and torque dipoles which allows accurate 3D Navier-Stokes calculations. Our model can reproduce swimmer propulsion for a large range of Reynolds numbers, and generate wake vortices in the inertial regime, reminiscent of the flow generated by the flapping tails of real fish. We performed a numerical exploration of the model from low to high Reynolds numbers and obtained universal laws using scaling arguments. We collected data from a wide variety of micro-organisms, thereby extending the experimental data presented in (M. Gazzola et al., Nature Physics 10, 758, 2014). Our theoretical scaling laws compare very well with experimental data across the different regimes, from Stokes to turbulent flows. We believe that this model, due to its relatively simple design, will be very useful for obtaining numerical simulations of collective effects within fish schools composed of hundreds of individuals.

Introduction The wide variety of means employed by living creatures to move in aquatic environments is fascinating Childress (1981). Motion generally involves a complex interplay between the deformation of the body and the surrounding fluid. From the smallest organisms, like bacteria, to colossal blue whales Berg (2004); S.Powar et al. (2022); Garcia et al. (2011); Smits (2019); Gray (1936); Wolfgang et al. (1999); Dabiri et al. (2006); Müller et al. (2001), the differences in length scale $L$, velocity $v$ and mode of locomotion are so vast that the elaboration of a universal model to describe swimming across these scales might seem impossible.

The importance of inertia with respect to viscous dissipation is quantified by the Reynolds number, which expresses the ratio of stress due to inertia to stress due to viscosity: $Re=\rho vL/\eta$, where $\rho$ and $\eta$ represent fluid density and viscosity, respectively. The Reynolds numbers associated with aquatic living species span several decades, typically ranging from $Re^{\rm fish}\sim 10^{3}-10^{6}$ for fish to $Re^{\rm micro-org}\lesssim 10^{-3}$ for micro-organisms like spermatozoa. Consequently, the drag force exerted by the fluid on the swimmer is largely dependent on the species considered, originating from either fully viscous dissipation at the small scale, or turbulent inertia at the large scale.

As a result, each swimming model Taylor (1951); Spagnolie and Lauga (2010); Lauga (2020); Lighthill (1960); Liu and Kawachi (1999); Kern and Koumoutsakos (2006) tends to offer a tailored approach specific to the flow regime considered and the corresponding body deformation. Most models focus on periodic body deformation, coupled with the surrounding fluid, and resolve the full swimming cycle. This provides specific approaches that address swimming at small scales, where viscous forces dominate (e.g. for micro-organisms Lauga and Powers (2009)), differently from the macroscopic strategies of large fish or mammals, which can leverage the inertia of the surrounding fluid to break time reversibility Purcell (1977). The highly diverse physical origins of particular swimming patterns represent an obstacle to the exploration of a more comprehensive and universal viewpoint.

In this letter, we propose a different approach. Deformation kinematics, such as undulations, oscillations and pulsations, are ignored, and locomotion is described using force and torque dipoles applied by a solid body of finite size $L$ on a fluid. While a similar description has already been used for micro-swimmers at low $Re$ Hernandez et al. (2005); Mehandia and Nott (2008); Jibuti et al. (2014), to the best of our knowledge, it has never yet been employed for high $Re$, when inertia starts to dominate viscous forces. In this work, we study the motion described by our model over $8$ decades of Reynolds numbers ($10^{-5}\lesssim Re\lesssim 10^{4}$), theoretically and numerically, comparing our results with experimental data. Although the direct effect on swimming velocity of the specific deformation of the swimmer’s body is acknowledged, particularly as it can reduce the drag force Gray (1936); Li et al. (2021), we would like to emphasize that we are not attempting to provide a detailed and precise analysis of a particular mode of locomotion, but rather a more universal description in terms of the forces applied to the fluid. While our swimmer model is minimal, the motion of the surrounding fluid is accurately captured using the full numerical resolution of the 3D Navier-Stokes equation, and our approach encompasses the different swimming regimes of a wide variety of aquatic species. Our model can also remarkably reproduce the characteristic wake vortices observed behind fish due to the flapping of their tails Drucker and Lauder (2002).

Our approach furthermore exhibits universal scaling laws which link the swimming Reynolds number $Re$ to a new dimensionless group, the thrust number defined below. We identify three different regimes: the Stokes regime ($Re<1$), a laminar regime ($1<Re<10^{3}$) and a turbulent regime ($Re>10^{3}-10^{4}$), and calculate the theoretical exponents of the scaling laws in the three regimes using simple scaling analyses (independently of the space dimension). These compare very well with the numerical simulations produced by our generic swimmer model. We also validate our results with experimental data presented in Gazzola et al. (2014) for the laminar and turbulent regimes, and further extend this validation with data collected on micro-swimmers for the Stokes regime.

Figure 1: The swimmer is an ellipsoidal rigid body of length $L$ and width $L/a_{r}$, with $a_{r}$ the aspect ratio. a) Time-dependent dipoles of forces $\{-\bm{{F}},\bm{{F}}\}$ and torques $\{-\bm{{T}},\bm{{T}}\}$ are applied at the center-of-mass of the body $\bm{{X}}_{\mathcal{B}}$ and the “phantom tail” $\bm{{X}}_{t}$, with $\bm{F}(t)=(\pi/2)\bm{F}_{0}\lvert\cos(\omega t)\rvert$ and $\bm{T}(t)=\bm{T}_{0}\cos(\omega t)$. The position of $\bm{{X}}_{t}$ can be controlled to steer the swimmer, but is kept such that $\|\bm{{X}}_{\mathcal{B}}-\bm{{X}}_{t}\|=L/2+L/a_{r}$. b) Static model obtained by averaging the time dependent model.

Figure 2: $Re$ as a function of $Th$ from 3D numerical simulations (crosses), obtained by solving the full Navier-Stokes equations with our swimmer model. We clearly obtain three regimes: $Re\sim Th^{1.0}$ for $Re\lesssim 20$; $Re\sim Th^{0.66}$ for $20\lesssim Re\lesssim 1000$ and $Re\sim Th^{0.51}$ for $Re\gtrsim 1000$. The three lines correspond to fitting curves and give the numerical scaling exponents. Note that the crossovers between the different regimes depend on the geometry of the swimmer; this curve corresponds to $a_{r}=4$.

Swimmer model.

The model uses a time-dependent force dipole combined with a torque dipole (Fig. 1(a)), both attached to a rigid body $\mathcal{B}$ of ellipsoidal shape. An autonomous swimming body creates its own motion, therefore the total sum of forces and torques must cancel out, due to the third law of Newton of the fluid-body system. Feynman et al. (2006). The model developed by Filella et al. (2018) presents some conceptual similarities. It represents each fish as a point-like active particle bearing a dipole in a potential 2D inviscid fluid, which allows consideration of the hydrodynamic interactions between fish in the far-field limit. However, since our aim is to explore a wide spectrum of Reynolds numbers (from viscous to inertial regimes), we solved the full incompressible Navier-Stokes equations in 3D and 2D (see SM ven (a)).

The swimmer exerts on the fluid a force dipole $\{\bm{-F},\bm{F}\}$ like that generally used for a micro-swimmer at a low $Re$ Hernandez et al. (2005); Mehandia and Nott (2008); Jibuti et al. (2014) (see Fig. 1). We used a pusher-like model that reproduces the force distribution of a fish at high Reynolds numbers. This approach can easily be extended to a more detailed model by using more complex force distributions (e.g. Mehandia and Nott (2008)). As shown in Fig. 1(a), the force dipole is composed of one force applied in the fluid at the rear of the body $\bm{{X}}_{t}$, mimicking a swimming organ, and an opposite force exerted inside the body at the center-of-mass $\bm{{X}}_{\mathcal{B}}$. The force is time-dependent with pulsation $\omega$ and pusher-like: $\bm{F}(t)=(\pi/2)\bm{F_{0}}\lvert\cos\,\omega t\rvert$ with $\mathbf{F}_{0}\cdot(\bm{{X}}_{\mathcal{B}}-\bm{{X}}_{t})>0$. The absolute value in the expression of $\bm{F}(t)$ enforces the pusher nature of the swimmer. We also consider a torque dipole: $\bm{T}(t)=\bm{T_{0}}\cos\,\omega t$, collocated with the force dipole. The torque at the back represents the stroke of the swimming organ and causes the vortex street Drucker and Lauder (2002) in the fish’s wake at high $Re$. An opposite torque is applied in the body (Fig. 1(a)), and represents the counter-reaction of the rest of the body. For practical reasons in the scaling analyses below, we also introduce force and torque densities $f_{0}\equiv F_{0}/L^{3}$ and $\tau_{0}\equiv T_{0}/L^{3}$ respectively. Averaging these dipoles over one period of time ($2\pi/\omega$) results in a simple static force dipole $\{\bm{-F_{0}},\bm{F_{0}}\}$ (Fig. 1(b)), while the average torque cancels out. This time averaging over one beating period is very similar to models of pushers and pullers beating at low $Re$ Hernandez et al. (2005); Mehandia and Nott (2008); Jibuti et al. (2014).

Scaling laws. The hydrodynamic nature of our model allows for simple scaling arguments, inspired by Gazzola et al. (2014) for inertial flows but translated to a more generic framework and extended to the non-inertial Stokes regime. In the following, we present $3D$ arguments, but they remain valid in $2D$ (see SM ven (a)). We also consider that all the lengths scale as $L$.

The body of the swimmer is submitted to different dominant drag forces depending on the Reynolds number. To describe this effect across all swimming regimes, we introduce the thrust number $Th$ as the ratio between the applied force density $f_{0}$ multiplied by inertial forces $\rho\lvert D\mathbf{v}/Dt\rvert\sim\rho v^{2}L^{-1}$, and the square of viscous forces $\lvert\eta\Delta\mathbf{v}\rvert^{2}\sim(\eta vL^{-2})^{2}$, which gives

$$Th\equiv\frac{\rho f_{0}L^{3}}{\eta^{2}}.$$ (1)

The thrust number appears naturally at all scales of Reynolds numbers, as shown below. It contains the force term at the origin of the motion, which is characterized by the Reynolds number. It therefore provides a convenient method for evaluating velocity as a function of force. Let us consider classic scaling arguments:

• At high Reynolds numbers, the boundary layer around the body is turbulent and pressure drag dominates. The corresponding force scales as $f_{0}L^{3}\sim\rho v^{2}L^{2}$. Since $v\sim\eta Re/(\rho L)$, we obtain $Re\sim(\rho f_{0}L^{3}/\eta^{2})^{1/2}=Th^{1/2}$.

• For small but finite Reynolds numbers, the regime is laminar and the viscous force in the boundary layer dominates: $f_{0}L^{3}\sim(\eta v/\delta)L^{2}$ where $\delta$ is the thickness of the boundary layer, which obeys the Blasius law Landau and Lifshitz (2013) $\delta\sim L\,Re^{-1/2}$. This finally leads to $Re\sim Th^{2/3}$.

• At low Reynolds numbers, the Stokes drag force dominates and the force applied on the body compensates the drag: $f_{0}L^{3}\sim\eta vL$. This gives $Re\sim\rho f_{0}L^{3}/\eta^{2}=Th$.

From the Stokes to the turbulent regime, we observe that the exponent $\alpha$ of the scaling $Re\sim Th^{\alpha}$ is always below one and decreases. It suggests a diminishing swimming performance as the Reynolds number of the swimmer increases. To confirm these three successive regimes, we present below numerical simulations with our swimmer model, exploring a large range of values for the thrust and Reynolds numbers.

Numerical simulations.

We performed direct numerical simulations of the incompressible Navier-Stokes equations in the presence of our model swimmer, a rigid ellipsoid body with the force and torque dipoles attached. The corresponding fluid momentum balance equation writes:

$$\rho\frac{D\bm{v}}{Dt}-\nabla\cdot\left(2\eta\bm{{E}}(\bm{{v}})\right)+\nabla p=\bm{{\mathcal{F}}}+\bm{{\mathcal{T}}}$$ (2)

where $\bm{{v}}$ and $p$ are respectively the velocity and pressure fields, $D/Dt$ denotes the material derivative $D\bm{{v}}/Dt\equiv\partial\bm{{v}}/\partial t+(\bm{v\cdot\nabla})\bm{v}$, $\bm{{E}}(\bm{{v}})\equiv(\bm{\nabla v}+\bm{\nabla v}^{t})/2$ is the strain-rate tensor, $\rho$ and $\eta$ denote the density and viscosity fields, and $\bm{{\mathcal{F}}}(t)$ and $\bm{{\mathcal{T}}}(t)$ respectively represent the – time-varying – force and torque dipoles attached to the swimmer, as introduced above. The rigid body $\mathcal{B}$ of the swimmer is accounted for with a fictitious domain penalty method inspired by (Janela et al., 2005); in practice, this can be implemented simply with a spatially-variable viscosity (Tanaka and Araki, 2000): $\eta=\eta_{f}+(\eta_{b}-\eta_{f})\,H_{\mathcal{B}}$ where $H_{\mathcal{B}}$ is the indicator – or Heaviside – function of $\mathcal{B}$. In practice, a viscosity ratio $\eta_{b}/\eta_{f}=10^{3}-10^{6}$ is applied between the fluid and the swimmer body to ensure that the rigid motion constraint $\bm{{E}}(\bm{{v}})=0$ is satisfied within $\mathcal{B}$. Density $\rho$ is defined as constant inside and outside $\mathcal{B}$, thus making the swimmer neutrally buoyant. Note that this fictitious domain penalty method allows the swimmer to be treated directly as part of the Navier-Stokes equations through the viscosity field, thereby avoiding the need to deal with moving boundary conditions and potential remeshing issues at the body interface in the discrete setting. The incrompressible Navier-Stokes equations are solved numerically using an implicit $\mathcal{P}2-\mathcal{P}1$ finite element method (Metivet et al., 2018) implemented in the parallel FEEL++ library Prud’homme et al. (2012). The position and orientation of the swimmer are updated at each time-step using a first order Euler scheme with the translational and rotational velocities computed from the fluid velocity field in $\mathcal{B}$. A comprehensive derivation of the numerical model and technical details are provided in SM ven (a).

This numerical framework is used to explore a large range of Reynolds numbers ($10^{-5}<Re<10^{4}$) and Thrust numbers ($10^{-3}<Th<10^{6})$, in order to evaluate the dependence of $Re$ as a function of $Th$ while varying each of the model’s different parameters ($L$, $\eta_{f}$, $f_{0}$, $\omega$ and $\tau_{0}$) separately (see SM ven (a)). As shown in Fig. 2, the numerical simulations are in perfect agreement with the scaling laws presented above, displaying the Stokes regime in the range $10^{-5}\lesssim Re\lesssim 10^{2}$, the laminar regime for $10^{2}\lesssim Re\lesssim 10^{3}$ and the turbulent regime for $10^{3}\lesssim Re$,

Figure 3: 2D Vorticity field for different $Re$. The simulations were performed with $\omega=2\pi$, $L=32$ and $\|\bm{{F_{0}}}\|=\|\bm{{T_{0}}}\|$. (a): $\lVert\bm{T}_{0}\rVert=4000$ and $Re=960$. (b): $\lVert\bm{T}_{0}\rVert=600$ and $Re=300$. (c): $\lVert\bm{T}_{0}\rVert=100$ and $Re=50$. (d): $\lVert\bm{T}_{0}\rVert=0.25$ and $Re=0.1$.

with fitting exponents that match the predictions. Note that the $Re$ ranges of each regime depend on the aspect ratio of the swimmer $a_{r}$, which is kept constant. We also found that $\omega$ and $\tau_{0}$ do not play any role in the $Re(Th)$ dependency, which confirms that all the important parameters are embedded in the $Th$ number.

Wake and vortices. Although torque plays no role in the scaling of $Re$ as a function of $Th$, it is essential to reproduce the wake at the rear of the swimmer in the inertial regime. The torque dipole can be used to account for flagellum, body undulation or tail beating to generate a reverse von Karman vortex street, as observed in the wake of a fish at high $Re$ number Kern and Koumoutsakos (2006); Ko et al. (2023).

Figures 3(a,b) illustrate the typical wakes obtained with $(T_{0},F_{0})=(4000,4000),(600,600)$ in 2D simulations.

Figure 4: Reynolds number $Re$ as a function of the swimming number $Sw$. Crosses are new experimental data available in the database ven (b). Dots are experimental data collected by Gazzola et al. (2014). The three lines correspond to fitted curves and give the numerical scaling exponents of the three different regimes: $Re\sim Sw^{0.98}$ for $Re<10$; $Re\sim Sw^{1.33}$ for $10<Re<10^{4}$; and $Re\sim Sw^{1.0}$ for $10^{4}<Re$.

The vortices created at the same frequency are spaced further apart as the Reynolds number (i.e. velocity) increases from $Re=100$ to $Re=960$. The length over which they dissipate becomes shorter towards the viscous regimes (Fig. 3(c)), vanishing completely at low $Re$ (Fig. 3(d)). Indeed, no vortices are present behind micro-organisms Drescher et al. (2010).

Comparison with experimental data. To compare our scaling results with experimental data, forces must be expressed in terms of observable data (Fig. 4), such as undulation or tail beating frequency. Although force $f_{0}$, torque $\tau_{0}$ and pulsation $\omega$ are independent quantities in the model, physical constraints exist between these quantities in living organisms. We made the reasonable assumption that the size of the swimming organ scales with the size of the body $L$ Gazzola et al. (2014). Force and torque generated by the swimming organ are such that $\tau_{0}\sim f_{0}L$.

In inertial regimes, i.e. excluding the Stokes regime that we address separately, the instantaneous force creates a transient acceleration of the fluid, which scales with $L\omega^{2}$, i.e. $f_{0}\sim\rho L\omega^{2}$. Introducing the observation-based swimming number $Sw=\rho\omega L^{2}/\eta$ Gazzola et al. (2014), the theoretical force drive can be related to experimentally measurable data as $Th\sim Sw^{2}$. The scaling laws previously derived from our model can thus be reformulated in terms of $Sw$: in the laminar regime $Re\sim Th^{2/3}\sim Sw^{4/3}$, while in the turbulent regime $Re\sim Th^{1/2}\sim Sw$, in accordance with the results of Gazzola et al. (2014).

Figure 5: Strouhal number $St=A\omega/v$ as a function of the Reynolds number $Re$. Dots correspond to the data collected in Gazzola et al. (2014), and crosses are new experimental data. The three lines correspond to guides for the eyes with the three different regimes: $\tfrac{St}{2\pi}=1$ for $Re<10$; $St\sim Re^{-0.25}$ for $10<Re<10^{4}$; and $\tfrac{St}{2\pi}\approx 0.3$ for $10^{4}<Re$.

In the Stokes regime, the force exerted by the swimming organ is balanced by viscous drag, so that $f_{0}L^{3}\sim\eta Lv\sim\eta L^{2}\omega$. Reintroducing again the swimming number $Sw$, we obtain $Th\sim Sw$, leading to $Re\sim Sw$. Note that this is a natural consequence of the absence of inertia: each stroke creates a net displacement that scales with $\sim L$, inducing a swimming velocity $v\sim L\omega$.

Figure 4 shows the experimental data collected from micrometer- to meter-size aquatic organisms along with the corresponding fitted scaling laws. In addition to the results from Gazzola et al. (2014), an excellent agreement between experimental data and hydrodynamic scaling laws is also obtained in the Stokes regime.

Note that $v\sim L\omega$ in both the Stokes and turbulent regimes, so that the corresponding Strouhal number $St=Sw/Re$ is constant, as illustrated in Figure 5. The transition between $St/2\pi\approx 0.3$ Gazzola et al. (2014) and $St/2\pi\approx 1$ occurs in the laminar regime where $St\sim Re^{-1/4}$.

Conclusion. By avoiding direct fluid-structure coupling, our generic swimmer model provides an efficient model for hydrodynamic propulsion, while retaining the salient features of swimming organisms across several decades of Reynolds numbers. High numerical stability and efficiency ensure fully tractable 2D and 3D simulations, thereby paving the way to large scale simulations with hundreds of agents. It also proposes a methodology for progressive refinement of the hydrodynamic field, by retaining higher moments of forces and torques, resulting in more complex propulsion models. This approach broadens our understanding of the swimming of aquatic organisms by revealing the universal relationship between the velocity of a swimmer and the force exerted by its swimming organ. The sub-linear dependence demonstrated between $Re$ and $Th$ suggests diminishing swimming performance as the swimmer’s Reynolds number increases. The scaling laws obtained also match the experimental data obtained from thousands of aquatic animals, ranging from large mammalians to micro-organisms. Our results shed new light on the general mechanisms underlying swimming and provide an efficient and robust numerical framework to investigate the collective behavior of swimmers in complex environments.

Acknowledgements. This project received financial support from the French National Research Agency (ANR-21-CE45-0005, FISHSIF project).