Physics:Scallop theorem

A scallop swims by slowly opening its two halves, then quickly closing them. This is a successful swimming method because the inertial force of the surrounding water dominates over the viscous force. In a low Reynolds number environment, a scallop would only oscillate in place.

In physics, the scallop theorem states that a swimmer that performs a reciprocal motion cannot achieve net displacement in a low-Reynolds number Newtonian fluid environment, i.e. a fluid that is highly viscous. Such a swimmer deforms its body into a particular shape through a sequence of motions and then reverts to the original shape by going through the sequence in reverse. It does not matter how fast or slow the swimmer executes the sequence. At low Reynolds number, time or inertia does not come into play, and the swimming motion is purely determined by the sequence of shapes that the swimmer assumes.

Edward Mills Purcell stated this theorem in his 1977 paper Life at Low Reynolds Number explaining physical principles of aquatic locomotion. The theorem is named for the motion of a scallop which opens and closes a simple hinge during one period. Such motion is not sufficient to create migration at low Reynolds numbers. The scallop is an example of a body with one degree of freedom to use for motion. Bodies with a single degree of freedom deform in a reciprocal manner and subsequently, bodies with one degree of freedom do not achieve locomotion in a highly viscous environment.

Background

Animation of a 3-sphere swimmer. It has one degree of freedom where the left arm extends and retracts. In low Reynolds number environments, this leads to no net displacement of the whole body as the arm completes a cycle of extension and retraction.

The scallop theorem is a consequence of the subsequent forces applied to the organism as it swims from the surrounding fluid. For an incompressible Newtonian fluid with density [math]\displaystyle{ \rho }[/math] and dynamic viscosity [math]\displaystyle{ \eta }[/math], the flow satisfies the Navier–Stokes equations:

[math]\displaystyle{ \rho\left(\dfrac{\partial}{\partial\mathrm{t}}+\mathbf{u}\cdot\nabla\right)\mathbf{u}=-\nabla p+\eta\nabla^2\mathbf{u},\quad \nabla\cdot\mathbf{u}=0 }[/math],

where [math]\displaystyle{ \mathbf{u} }[/math] denotes the velocity of the fluid. However, at the low Reynolds number limit, the inertial terms of the Navier-Stokes equations on the left-hand side tend to zero. This is made more apparent by nondimensionalizing the Navier–Stokes equations. By defining a characteristic velocity and length, [math]\displaystyle{ u_0 }[/math] and [math]\displaystyle{ L }[/math], we can cast our variables to dimensionless form:

[math]\displaystyle{ \mathbf{\tilde{u}}=\dfrac{\mathbf{u}}{u_0};\quad\mathbf{\tilde{r}}=\dfrac{\mathbf{r}}{L};\quad\tilde{t}=\dfrac{t}{(L/u_0)};\quad\tilde{p}=\dfrac{p}{(\eta u_0/L)} }[/math].

where the dimensionless pressure is appropriately scaled for flow with significant viscous effects. Plugging these quantities into the Navier-Stokes equations gives us:

[math]\displaystyle{ \dfrac{\rho u_{0}^{2}}{L}\left(\dfrac{\partial}{\partial\tilde{t}}+\mathbf{\tilde{u}}\cdot\tilde{\nabla}\right)\mathbf{\tilde{u}}=-\dfrac{\eta u_{0}}{L^{2}}\left( \tilde{\nabla} \tilde{p}+\tilde{\nabla}^2\mathbf{\tilde{u}} \right ),\quad \tilde{\nabla}\cdot\mathbf{\tilde{u}}=0 }[/math]

And by rearranging terms, we arrive at a dimensionless form:

[math]\displaystyle{ \text{Re}\left(\dfrac{\partial}{\partial\tilde{t}}+\mathbf{\tilde{u}}\cdot\tilde{\nabla}\right)\mathbf{\tilde{u}}=-\tilde{\nabla} \tilde{p}+\tilde{\nabla}^2\mathbf{\tilde{u}},\quad \tilde{\nabla}\cdot\mathbf{\tilde{u}}=0 }[/math],

where [math]\displaystyle{ \text{Re} = \rho u_0L/\eta }[/math] is the Reynolds number. In the low Reynolds number limit (as [math]\displaystyle{ \mathrm{Re}\rightarrow 0 }[/math]), the LHS tends to zero and we arrive at a dimensionless form of Stokes equations. Redimensionalizing yields:

[math]\displaystyle{ 0=-\nabla p+\eta\nabla^2\mathbf{u},\quad \nabla\cdot\mathbf{u}=0 }[/math].

Statement

The consequences of having no inertial terms at low Reynolds number are:

• One consequence means that the swimmer experiences virtually no net force or torque.

• A second consequence tells us that the velocity is linearly proportional to the force (same can be said about angular velocity and torque).

• The Stokes equations become are linear and independent of time.

In particular, for a swimmer moving in the low Reynolds number regime, its motion satisfies:

• Independent of time: The same motion may be sped up or slowed down, and it would still satisfy the Stokes equations. More geometrically, this means that the motion of a swimmer in the low Reynolds number regime is purely determined by the shape of its trajectory in configuration space.
• Kinematic reversibility: The same motion may be reversed. Any instantaneous reversal of the forces acting on the body will not change the nature of the fluid flow around it, simply the direction of the flow. These forces are responsible for producing motion. When a body has only one degree of freedom, reversal of forces will cause the body to deform in a reciprocal fashion. For instance, a scallop opening its hinge will simply close it to try to achieve propulsion. Since the reversal of forces does not change the nature of the flow, the body will move in the reverse direction in the exact same manner, leading to no net displacement. This is how we arrive at the consequences of the scallop theorem.

Proof by scaling

This is closer in spirit to the proof sketch given by Purcell. The key result is to show that a swimmer in a Stokes fluid does not depend on time. That is, a one cannot detect if a movie of a swimmer motion is slowed down, sped up, or reversed. The other results then are simple corollaries.

The stress tensor of the fluid is [math]\displaystyle{ \sigma_{i j}=-p \delta_{i j}+\mu (\partial_i u_j + \partial_j u_i) }[/math].

Let [math]\displaystyle{ r }[/math] be a nonzero real constant. Suppose we have a swimming motion, then we can do the following scaling:[math]\displaystyle{ p \mapsto rp; \quad u \mapsto ru;\quad \sigma \mapsto r\sigma }[/math]and obtain another solution to the Stokes equation. That is, if we scale hydrostatic pressure, flow-velocity, and stress tensor all by [math]\displaystyle{ r }[/math], we still obtain a solution to the Stokes equation.

Since the motion is in the low Reynolds number regime, inertial forces are negligible, and the instantaneous total force and torque on the swimmer must both balance to zero. Since the instantaneous total force and torque on the swimmer is computed by integrating the stress tensor [math]\displaystyle{ \sigma }[/math] over its surface, the instantaneous total force and torque increase by [math]\displaystyle{ r }[/math] as well, which are still zero.

Thus, scaling both the swimmer's motion and the motion of the surrounding fluid scales by the same factor, we still obtain a motion that respects the Stokes equation.

Proof by vector calculus

The proof of the scallop theorem can be represented in a mathematically elegant way. To do this, we must first understand the mathematical consequences of the linearity of Stokes equations. To summarize, the linearity of Stokes equations allows us to use the reciprocal theorem to relate the swimming velocity of the swimmer to the velocity field of the fluid around its surface (known as the swimming gait), which changes according to the periodic motion it exhibits. This relation allows us to conclude that locomotion is independent of swimming rate. Subsequently, this leads to the discovery that reversal of periodic motion is identical to the forward motion due to symmetry, allowing us to conclude that there can be no net displacement.

Rate-independence

The reciprocal theorem describes the relationship between two Stokes flows in the same geometry where inertial effects are insignificant compared to viscous effects. Consider a fluid filled region [math]\displaystyle{ V }[/math] bounded by surface [math]\displaystyle{ S }[/math] with a unit normal [math]\displaystyle{ \hat{\mathbf{n}} }[/math]. Suppose we have solutions to Stokes equations in the domain [math]\displaystyle{ V }[/math] possessing the form of the velocity fields [math]\displaystyle{ \mathbf{u} }[/math] and [math]\displaystyle{ \mathbf{u}' }[/math]. The velocity fields harbor corresponding stress fields [math]\displaystyle{ \mathbf{\sigma} }[/math] and [math]\displaystyle{ \mathbf{\sigma}' }[/math] respectively. Then the following equality holds:

[math]\displaystyle{ \iint_S \mathbf{u}\cdot (\boldsymbol{\sigma}' \cdot \hat{\mathbf{n}})~\mathrm{d}S = \iint_S \mathbf{u}' \cdot (\boldsymbol{\sigma} \cdot \hat{\mathbf{n}})~\mathrm{d}S }[/math].

The reciprocal theorem allows us to obtain information about a certain flow by using information from another flow. This is preferable to solving Stokes equations, which is difficult due to not having a known boundary condition. This particularly useful if one wants to understand flow from a complicated problem by studying the flow of a simpler problem in the same geometry.

One can use the reciprocal theorem to relate the swimming velocity, [math]\displaystyle{ \mathbf{U} }[/math], of a swimmer subject to a force [math]\displaystyle{ \mathbf{F} }[/math] to its swimming gait [math]\displaystyle{ \mathbf{u}_S }[/math]:

[math]\displaystyle{ \hat{\mathbf{F}}\cdot\mathbf{U} = -\iint_S \mathbf{u}_S\cdot (\boldsymbol{\hat{\sigma}} \cdot \mathbf{n})~\mathrm{d}S }[/math].

Now that we have established that the relationship between the instantaneous swimming velocity in the direction of the force acting on the body and its swimming gait follow the general form

[math]\displaystyle{ \mathbf{U} = \iint \dot{\mathbf{r}}_S\cdot\mathbf{g}(\mathbf{r}_S)~\mathrm{d}S }[/math],

where [math]\displaystyle{ \mathbf{u}_S\equiv\dot{\mathbf{r}_S}=\mathrm{d}\mathbf{r}_S/\mathrm{d}t }[/math] and [math]\displaystyle{ \mathbf{r}_S }[/math] denote the positions of points on the surface of the swimmer, we can establish that locomotion is independent of rate. Consider a swimmer that deforms in a periodic fashion through a sequence of motions between the times [math]\displaystyle{ t_0 }[/math] and [math]\displaystyle{ t_1 }[/math]. The net displacement of the swimmer is

[math]\displaystyle{ \Delta X=\int_{t_0}^{t_1}\mathbf{U}~\mathrm{d}t }[/math].

Now consider the swimmer deforming in the same manner but at a different rate. We describe this with the mapping

[math]\displaystyle{ t' = f(t),\quad \mathbf{r}_S(t)=\mathbf{r'}_S(t'),\quad \dot{\mathbf{r}}_S(t)=\dfrac{\mathrm{d}\mathbf{r'}_S(t')}{\mathrm{d}t}=\dfrac{\mathrm{d}\mathbf{r'}_S(t')}{\mathrm{d}t'}\cdot\dfrac{\mathrm{d}t'}{\mathrm{d}t}=\dot{\mathbf{r}}_S'(t')\dot{f}(t) }[/math].

Using this mapping, we see that

[math]\displaystyle{ \Delta X'=\int_{t_0}^{t_1}\mathbf{U}'(t')~\mathrm{d}t'=\int_{t_0}^{t_1}\mathbf{U}'(f(t))\dot{f}~\mathrm{d}t= \int_{t_0}^{t_1}\iint \dot{\mathbf{r}}_S'\dot{f}\cdot\mathbf{g}(\mathbf{r'}_S)~\mathrm{d}S\mathrm{d}t=\int_{t_0}^{t_1}\iint \dot{\mathbf{r}_S}\cdot\mathbf{g}(\mathbf{r}_S)~\mathrm{d}S\mathrm{d}t }[/math]

[math]\displaystyle{ =\int_{t_0}^{t_1}\mathbf{U}(t)~\mathrm{d}t\rightarrow\Delta X'=\Delta X }[/math].

This result means that the net distance traveled by the swimmer does not depend on the rate at which it is being deformed, but only on the geometrical sequence of shape. This is the first key result.

Symmetry of forward and backward motion

If a swimmer is moving in a periodic fashion that is time invariant, we know that the average displacement during one period must be zero. To illustrate the proof, let us consider a swimmer deforming during one period that starts and ends at times [math]\displaystyle{ t_0 }[/math] and [math]\displaystyle{ t_1 }[/math]. That means its shape at the start and end are the same, i.e. [math]\displaystyle{ \mathbf{r}_S(t_0)=\mathbf{r}_S(t_1) }[/math]. Next, we consider motion obtained by time-reversal symmetry of the first motion that occurs during the period starting and ending at times [math]\displaystyle{ t_2 }[/math] and [math]\displaystyle{ t_3 }[/math]. using a similar mapping as in the previous section, we define [math]\displaystyle{ t_2 = f(t_1) }[/math] and [math]\displaystyle{ t_3 = f(t_0) }[/math] and define the shape in the reverse motion to be the same as the shape in the forward motion, [math]\displaystyle{ \mathbf{r}_S(t)=\mathbf{r'}_S(t') }[/math]. Now we find the relationship between the net displacements in these two cases:

[math]\displaystyle{ \Delta X'=\int_{t_2}^{t_3}\mathbf{U}'(t')~\mathrm{d}t'=\int_{t_1}^{t_0}\mathbf{U}(t)~\mathrm{d}t=-\int_{t_0}^{t_1}\mathbf{U}(t)~\mathrm{d}t=-\Delta X }[/math].

This is the second key result. Combining with our first key result from the previous section, we see that [math]\displaystyle{ \Delta X'=\Delta X=-\Delta X \rightarrow \Delta X = 0 }[/math]. We see that a swimmer that reverses its motion by reversing its sequence of shape changes leads to the opposite distance traveled. In addition, since the swimmer exhibits reciprocal body deformation, the sequence of motion is the same between [math]\displaystyle{ t_2 }[/math] and [math]\displaystyle{ t_3 }[/math] and [math]\displaystyle{ t_0 }[/math] and [math]\displaystyle{ t_1 }[/math]. Thus, the distance traveled should be the same independently of the direction of time, meaning that reciprocal motion cannot be used for net motion in low Reynolds number environments.

Exceptions

The scallop theorem holds if we assume that a swimmer undergoes reciprocal motion in an infinite quiescent Newtonian fluid in the absence of inertia and external body forces. However, there are instances where the assumptions for the scallop theorem are violated. In one case, successful swimmers in viscous environments must display non-reciprocal body kinematics. In another case, if a swimmer is in a non-Newtonian fluid, locomotion can be achieved as well.

Types of non-reciprocal motion

In his original paper, Purcell proposed a simple example of non-reciprocal body deformation, now commonly known as the Purcell swimmer. This simple swimmer possess two degrees of freedom for motion: a two-hinged body composed of three rigid links rotating out-of-phase with each other. However, any body with more than one degree of freedom of motion can achieve locomotion as well.

In general, microscopic organisms like bacteria have evolved different mechanisms to perform non-reciprocal motion:

• Use of a flagellum, which rotates, pushing the medium backwards — and the cell forwards — in much the same way that a ship's screw moves a ship. This is how some bacteria move; the flagellum is attached at one end to a complex rotating motor held rigidly in the bacterial cell surface.
• Use of a flexible arm: this could be done in many different ways. For example, mammalian sperm have a flagellum which, whip-like, wriggles at the end of the cell, pushing the cell forward. Cilia are quite similar structures to mammalian flagella; they can advance a cell like paramecium by a complex motion not dissimilar to breast stroke.

Geometrically, the rotating flagellum is a one-dimensional swimmer, and it works because its motion is going around a circle-shaped configuration space, and a circle is not a reciprocating motion. The flexible arm is a multi-dimensional swimmer, and it works because its motion is going around a circle in a square-shaped configuration space. Notice that the first kind of motion has nontrivial homotopy, but the second kind has trivial homotopy.

Non-Newtonian fluids

The assumption of a Newtonian fluid is essential since Stokes equations will not remain linear and time-independent in an environment that possesses complex mechanical and rheological properties. It is also common knowledge that many living microorganisms live in complex non-Newtonian fluids, which are common in biologically relevant environments. For instance, crawling cells often migrate in elastic polymeric fluids. Non-Newtonian fluids have several properties that can be manipulated to produce small scale locomotion.

First, one such exploitable property is normal stress differences. These differences will arise from the stretching of the fluid by the flow of the swimmer. Another exploitable property is stress relaxation. Such time evolution of such stresses contain a memory term, though the extent in which this can be utilized is largely unexplored. Last, non-Newtonian fluids possess viscosities that are dependent on the shear rate. In other words, a swimmer would experience a different Reynolds number environment by altering its rate of motion. Many biologically relevant fluids exhibit shear-thinning, meaning viscosity decreases with shear rate. In such an environment, the rate at which a swimmer exhibits reciprocal motion would be significant as it would no longer be time invariant. This is in stark contrast to what we established where the rate in which a swimmer moves is irrelevant for establishing locomotion. Thus, a reciprocal swimmer can be designed in a non-Newtonian fluid. Qiu et al. (2014) were able to design a micro scallop in a non-Newtonian fluid.

See also

• Bacterial motility
• Microswimmer
• Protist locomotion

External links

• E.M. Purcell. Life at Low Reynolds Number, American Journal of Physics vol 45, p. 3-11 (1977)
• Kinematic Reversibility and the Scallop Theorem
• Video of a swimming robot unable to propel in viscous fluid due to the Scallop theorem