Modelling and controllability of the motion of a slender, flexible
micro-swimmer

arXiv:1803.08795v1 [cs.SY] 21 Mar 2018

Sudin Kadam1 and Ravi N. Banavar1
Abstract— The mechanism of swimming at very low Reynolds
number conditions is a topic of interest to biologists and
engineering community. We develop a novel kinematic model of
a slender flexible swimmer which locomotes in a low Reynolds
number regime. In contrast to existing techniques that model
such systems as a connected set of straight, rigid links, the
novelty of our technique stems from the fact that we model the
swimmer with two components - one is a straight, rigid body
(the head) and the other is a flexible member (the tail). Using
Cox theory we model the gradient of the forces as a function
of the instantaneous shape of the swimmer and its velocity.
By virtue of the low inertia conditions, an expression for the
translational and rotational velocity of the head is obtained
for the planar motion in the form of a Lie algebra of the
Special Euclidean group. We explain the principal fiber bundle
structure of the configuration space of the swimmer and use
that to show a weak controllability result for a type of slender
flexible swimmer where the shape space is the space of all
continuous curves of a given length. A set of simulation results
is presented showing the variation of the swimmer head velocity
for a bump function moving along the swimmer length.

Most of our intuition about locomotion stems from inertiadominant systems used by larger animals. As opposed to
such systems, microorganisms resort to a creeping motion. This motion occurs in a fluid medium with very
low Reynolds number, which is the ratio of the inertial to
viscous forces acting on the swimmer’s body. The essential
mechanism of the motion in this regime is based on the fact
that the inertial effects are almost negligible as compared to
the viscous forces due to Reynolds number of the order of
10−4 . To get a relative sense of the numbers, the Reynolds
number for a man swimming in water is of the order of 104
[4], [5]. Organisms such as motile bacteria and Escherichia
Coli, commonly found in the human intestines, use one or
more flagella, which are long helical filaments and act as
propeller. Figure 1 shows a few such microbes with flagella.

I. INTRODUCTION
Swimming at micro scales is a topic of growing interest.
There has been a lot of research in exploring new and
efficient ways to generate propulsion at these scales, see [1],
[2]. A better understanding of the mechanism of swimming
can lead to many applications in several fields such as
micro-machining, nano technology and medicine for targeted
drug delivery. Micro-robotics is one of the recently evolving
fields in miniature robotics, especially mobile robots with
characteristic dimensions to the scale of a micron. Figure 1
shows an example of a microswimmer developed at Monash
University. Apart from applications in robotics, since cell
or microbial locomotion is an essential part of biological
systems, understanding the means of locomotion in this
regime is of a great academic interest to biologists.

Fig. 1: A flagellum propelled microbot [3]
1 Sudin Kadam and Ravi N. Banavar are with Systems and Control Engineering Department, Indian Institute of Technology Bombay, Mumbai, India,
400072 sudin@sc.iitb.ac.in, banavar@iitb.ac.in

(a) Sperm cells [6]

(b) Spirillum bacteria [7]

There has been a lot of research on the mechanism of
motion of micro-organisms and micro-robotic swimmers.
The Purcell’s swimmer is a famous example of the simplest possible low Reynolds number swimmer [8], [9]. This
swimmer has three straight, slender and rigid links connected through 2 rotary joints through which the actuation
is affected. Through non-reciprocal motion of the outer
links, the Purcell’s swimmer can achieve any point-to-point
motion on its configuration space. Although there has been
a lot of research on the Purcell’s swimmer and a few
other low Reynolds number swimmers which have finitely
many control inputs from mechanics and control theoretic
perspective [5], [8], [10], [11], the dynamic behaviour of
these swimmers differ significantly from most of the microorganisms found in the nature whose limbs are not rigid
links. Most of the micro-organisms have flexible limbs,
whose actuation results in the net motion of the swimmers.
Clearly, there is a need to understand the mechanism of
motion of such microswimmers with flexible link, which is
far more challenging than such analysis with rigid links. The
approach to model these systems is through Stokes equations.
One of the earlier works formulates this problem of

flexible swimmers’ modelling in terms of a gauge field
on the space of the shapes [12]. Existence and uniqueness
of the solution to the model obtain using exterior Stokes
problem when a deformable body is placed in a threedimensional fluid is presented in [13]. The recent works
have extensively studied control theoretic problems such as
controllability and optimal control of such self propelling
deformable bodies [14], [15]. All of these works target to
solve the problem when the deformable body is a smooth
embedding of sphere in the ambient fluid. Our work relates to
a similar control theoretic modelling and analysis of flexible
micro-swimmers, but which are slender. As discussed earlier,
the motivation for this is that a lot of microorganisms found
in the nature are slender. The approach makes use of Cox
theory for modelling, which makes use of the slenderness of
the swimmer to simplify the computation while solving the
exterior Stokes problem.
A. Contribution
To the best of our knowledge, this is the first attempt at
getting a mathematical model of a slender, flexible swimmer’s motion at low Reynolds number conditions. The model
obtained gives an expression for the velocity of the head of
the swimmer as a function of the shape and velocity of the
rest of the body. For many micro-robotics applications, since
the payload of the swimmer is located in the frontal part of
the body, having control over the velocity and position of
the swimmer’s head is of importance. The model presented
in this paper is in a form suitable to do an analysis of the
swimmer head motion. Furthermore, we present a geometric
interpretation of the configuration space of the swimmer. A
weak controllability analysis is presented using the strong
controllability of the planar Purcell’s swimmer which gives
a result on the ability to do point to point reconfiguration of
the swimmer head. Moreover, a set of simulation results for
a class of shape curves and their velocities gives an insight
into the relationship of the parameters of these motion curves
with the swimmer motion.
B. Organization of the paper
Section 2 gives a brief explanation of the Stokes equations
and Cox theory to get the expression of the force gradient
acting on a slender swimmer. This is followed by using
this force expression and the condition of strong dominance
of the viscous forces over inertial forces to obtain the
kinematic model of the swimmer. The geometry of the
configuration space of this swimmer along with a result
on weak controllability is presented in section 3. Section 4
analyses the behaviour of the kinematic model for a class of
spatio-temporal curves of the shape of the swimmer through
a set of simulation results. Section 5 concludes the paper with
a discussion on the possible avenues for the future work.

Fig. 3: Slender flexible swimmer

shows an arbitrary configuration of such a slender body
in a 3-dimensional space along with the inertial reference
frame defined by axes {iR , jR , kR }. We do arc-length
parametrization of the body along its length measured from
the head. So a point on the body is characterised by the
dimensionless parameter s ∈ [0, 1]. The cross sectional
radius at a given point of the curve is taken as bλ(s), where
λ(s) is a dimensionless quantity function of s. Let u∗ (s)
be the velocity of the point at arc-length s with respect
to the reference frame. u(r) be the velocity of fluid at a
general point at r with respect to the reference frame. The
velocity field u and pressure field p in the fluid satisfy the
dimensionless Stokes equations
µ∇2 u − ∇p = 0,

∇◦u=0

with the boundary conditions,
u = u∗ (s) on the body at point s and,
u → 0 as r → ∞
∂
∂
∂
iR + ∂y
jR + ∂z
kR
Where ∇ is the operator defined as ∂x
2
2
2
∂
∂
∂
2
and ∇ = ∂x2 + ∂y2 + ∂z2 is the Laplacian operator.
The complete velocity and pressure field is the solution of
the Stokes equations subject to the boundary conditions. To
write the expression of the force density at each point of
the swimmer we set a reference Cartesian coordinate frame
at each point of the swimmer as shown in figure 4. In
these frames, the axis i is along the tangent direction, j is
orthogonal to i and in the plane formed by u∗ (s) − u and
the i axis, and k completes the right handed frame.

II. K INEMATIC MODEL OF THE FLEXIBLE SWIMMER
Consider a long slender body of a circular cross-section,
the length of the body be l and the cross-sectional radius b along with a reference coordinate frame. Figure 3

(1)

Fig. 4: Slender flexible swimmer

A. Cox theory and fluid force gradient
According to Cox theory [17], the density of the fluid
forces acting at a point at arc-length s are obtained in this
local coordinate frame in terms of the swimmer’s local shape,
its slenderness ratio h = bl and the difference between the
local body velocity and fluid velocity in the following form.
F (s) = 2π(E(h, s)i − 4C(h, s)j − 4D(h, s)k)

(2)

Fig. 5: Forces and moments on the swimmer

where,
1
1
1
(uj − u∗j (s)) +
{ (uj −
(log(h))−1
(log(h))−2 2
1
1
u∗j (s))( + log(2) − log(λ) + log())} + Jj + ...
2
2
1
Jk
D(h, s) =
2(log(h))−2
1
1
1
(ui − u∗i (s)) +
{ (ui −
E(h, s) = −
(log(h))−1
(log(h))−2 2
1
1
u∗i (s))( + log(2) − log(λ) + log())} + Ji + ...
2
2
and J is a vector given by
Z
Z 1
δij
(Ri − R̂i )(Rj − R̂j )
1 s−
+
]{
+
×
Ji = [
2 0
|R − R̂|3
s+ |R − R̂|
1
{δjk − t̂j t̂k }{uk (R̂) − u∗k (ŝ)}}dŝ
(3)
2
where δij is the Kronecker delta defined as
(
0 if i 6= j,
δij =
(4)
1 if i = j.
C(h, s) =

and R is the position vector of the center-line of the infinitesimal segment of length 2 and R̂ is the position r at the point
on the center line with s = ŝ, both in the inertial frame. t̂ is
the unit tangent vector to the curve formed by the swimmer at
arc length ŝ in the inertial frame. The force gradient equation
(2) can be compactly written as the following vector equation
in terms of the tangent vector t̂ at the arclength s to the curve
formed by the swimmer’s body [17]
u − u∗ (s) J + (u − u∗ (s))log(2) λ )
+
]◦
log(k)
(log(h))2
(u − u∗ )
1
[t̂t̂ − 2I] +
◦ [3t̂t̂ − 2I] + O{
}] (5)
2
(log(h))
(log(h))3
where, I is the identity operator defined such that for any
vector in v ∈ R3 , v ◦ I = v. For  → 0, the value of J can
be written as
F (s) = 2π[[

J = −{u(R) − u∗ (s)}log() + O(1)

(6)

We substitute this value of J in the equation 5. We also
assume that the ambient fluid is stationary, u = 0, and that
body has a circular cross-section of a constant radius, i.e.
λ(s) = 1. A a consequence, we get the following form of
force gradient which is independent of 
−u∗ (s) −u∗ (s)log(2)
F (s) = 2π[[
+
] ◦ [t̂t̂ − 2I]−
log(h)
(log(h))2
−u∗ (s)
1
◦ [3t̂t̂ − 2I] + O{
}]
2
2(log(h))
(log(h))3

(7)

B. Model of the swimmer with head and tail
We now obtain the expression for the total force and
moment acting on the swimmer body with respect to the
reference frame at the head of the swimmer (s = 0). We
refer to this frame as the head frame in the rest of the paper.
We define Q(s) to be an appropriate transformation which
transforms the forces and moments in the body frame at the
point at arc-length s to those in the head frame, and r(s) is
the position vector of the point at arc-length s with respect to
the head frame. Thus the net force and moment in the head
frame is obtained by integrating the force gradient from the
equation 7 over the entire body length as follows Z 1
F = 2π
Q(s)F (s)ds
(8)
0
Z 1
M = 2π
r(s) × [Q(s)F (s)]ds
(9)
0

By virtue of the low Reynolds number conditions, the net
forces and moments acting on the body should be zero in
the body coordinate frame at the tip of the head (s = 0).
1
= c and neglect higher
Along with this, we substitute log(h)
order terms in equation 7 to get the following net force and
moment expressions in the head frame Z 1

2π
Q(s) −cu∗ (s) − log(2)c2 u∗ (s) ◦ [t̂t̂ − 2I]−
0

u∗ (s)c2 ◦ [3t̂(s)t̂(s) − 2I] ds = 0
(10)
Z 1

2π
Q(s) r(s) × [−cu∗ (s) − log(2)c2 u∗ (s) ◦ [t̂(s)t̂(s)−
0

2I] − u∗ (s)c2 ◦ [3t̂(s)t̂(s) − 2I]] ds = 0
(11)
We now assume that the swimmer’s head is a small rigid
straight link of length δ and of the same cross sectional
radius as that of the rest of the body. We split the forces
and moments acting on the body as those acting on the head
and the tail to get the force and moment balance equations
as follows Fhead + Ftail = 0

(12)

Mhead + Mtail = 0

(13)

The forces and moments are represented in the frame
attached to the tip of the swimmer head. Figure 5 shows the
schematic. For the further analysis we consider the motion of
the swimmer restricted in the plane formed by the x, y axes
of the head frame. We consider the tip of the head to have

translational velocity as v = [v0x , v0y , 0]T with respect to
the reference frame and angular speed of this segment as ω0
about this tip with respect to the reference frame at the head.
We note that the transformation Q(0) for the head segment
is just the identity transformation. Thus we get the following
force equation for the motion of the head of the swimmer -

where,
t0 = [t0x , t0y ]T is the tangent vector at the head s = 0
a11 = −cδ − c2 δlog(2)(t0x t0x + 1) − c2 δ(3t0x t0x − 2),
a12 = 0,
δ2
a13 = (−t0y − c2 log(2)t0y t0x t0x − c2 log(2)+
2
3c2 t0y t0x t0x − c2 t0y ),

δ

Z


−c(v + st̂(s) × ω) − log(2)c2 (v + st̂(s) × ω)◦

a21 = 0,

0
2

a22 = −cδ − c2 δlog(2)(t0y t0y + 1) − c2 δ(3t0y t0y − 2),
δ2
a23 = (−t0x − c2 log(2)t0x t0y t0y + log(2)c2 t0x −
2
3c2 t0x t0y t0y − 2c2 t0x )



[t̂(s)t̂(s) − 2I] − c (v + st̂(s) × ω) ◦ [3t̂(s)t̂(s) − 2I] ds
Z 1

Q(s) −cu∗ (s) − log(2)c2 u∗ (s) ◦ [t̂(s)t̂(s) − 2I]−
=−
δ

c2 u∗ (s) ◦ [3t̂(s)t̂(s) − 2I] ds
(14)
and the moment equation is obtained as follows-

Similarly, the left hand side of the moment equation 15 gets
simplified as follows Z

Z

δ

r(s) × [−c(v + st̂ × ω) − log(2)c2 (v + st̂ × ω)◦

0

[t̂(s)t̂(s) − 2I] − c2 (v + st̂ × ω) ◦ [3t̂(s)t̂(s) − 2I]]ds
Z 1

=−
r(s) × Q(s)[−cu∗ (s) − log(2)c2 u∗ (s)◦
δ

[t̂(s)t̂(s) − 2I] − c2 u∗ (s) ◦ [3t̂(s)t̂(s) − 2I]] ds (15)
The objective is to first determine the translational and
rotational velocities of the swimmer head v and ω for a given
profile of u∗ (s), s ∈ (δ, 1], and then to design u∗ (s, t) to
achieve a desired velocity trajectory v(t) of the head of the
swimmer. Using these equations we will now obtain the head
velocity of the swimmer as a function of the shape position
and shape velocity u∗ (s). We consider the planar motion.
Thus at each point the z component of the forces acting on
the body is zero and only the moment about the z-axis is
non-zero. We denote by v0 = [v0,x , v0,y ]T the translational
velocity of the tip of the swimmer head (s = 0) and ω0 as
the angular velocity of the swimmer head about the its tip.
Let t̂0 is the tangent vector to the head. Thus the left hand
side of the force equation 14 is computed as follows
Z

δ


−c(v0 + st̂0 × ω0 ) − log(2)c2 (v0 + st̂0 × ω0 )◦
0

t̂0 t̂0 − 2I) − (v0 + st̂0 × ω) ◦ (3t̂0 t̂0 − 2I) ds


 
 
δ 2 t0y
δ2
v0x
v
2
= −cδ
+ ω0
− log(2)c {δ 0x − ω0 ×
v0y
v0y
2 −t0x
2


  
 
t0y
t
t0x
v
} ◦ { 0x
− 2I} − c2 (δ 0x +
−t0x
t0y t0y
v0y


  
2
δ
t0y
t
t0x
ω0
) ◦ {3 0x
− 2I}
−t0x
t0y t0y
2
 

 v
a
a12 a13  0x 
v0y
= 11
(16)
a21 a22 a23
ω0

δ

st̂0 × [−c(v0 + st̂0 × ω0 ) − log(2)c2 (v0 + st̂0 × ω0 )◦

0

[t̂0 t̂0 − 2I] − c2 (v + st̂0 × ω) ◦ [3t̂0 t̂0 − 2I]]ds

   
  
Z δ
t0x
0
−cv0x
t0x
=
s t0y  × −cv0y  + s t0y  ×  0  − log(2)c2
0
0
0
ω0
0
 

        
v0x
t0x
t0x
t0x
0
× v0y  + s t0y  ×  0  ◦ t0y  t0y  − 2I
ω0
0
0
0
0
  
   
v0x
t0x
0
− c2 v0y  + s t0y  ×  0 
0
0
ω0
    

t0x
t0x
◦ 3 t0y  t0y  − 2I ds
0
0

T  
Z δ
cst0y − (2log(2) − 2)st0y
v0x

 v0y  ds
−sct0x + (2log(2) − 2)st0x
=
0
−s2 t20x − 2s2 t20y − (2log(2) − 2)st0x
ω0
 

 v0x
= a31 a32 a33 v0y 
(17)
ω0
where in the last 2 equations just the third row of equations is
written since the the first 2 rows corresponding to the x and
y components would be zero after taking the cross products.
The terms aij in these equations are
δ2
δ2
− (2log(2) − 2) ty ,
2
2
δ2
δ2
a32 = − ctx + (2log(2) − 2) tx ,
2
2
δ3 2
δ3 2
δ2
a33 = − tx − 2 ty − (2log(2) − 2) tx
3
3
2
Thus, combining equations 16 and 17, we get 3 equations in
v0x , v0y , ω0 . These 3 equations are the forces and moments
acting on the swimmer head of length δ, which are the
left hand sides of equations 14 and 15. Substituting these
a31 = c

in the left hand sides of equations 14 and 15, respectively,
and combining the 2 equations, we get the following single
equation
 
v0x
A(k, t̂, δ) v0y  = Ftail (u∗ , t, k)
(18)
ω0
where, A(k, t̂, δ) is a 3×3 matrix such that A(i, j) = aij and
the total forces and moments acting on the tail is Ftail which
is the function of the shape velocity u∗ : (δ, 1] ∈ R 7→ R2
and the instantaneous shape of the body which is the space of
all continuous curves in R2 . Substituting ξ = [v0x , v0y , ω0 ]T
as the resultant velocity of the swimmer head, we get the
following kinematic equation
ξ = −A−1 (k, t̂, δ)Ftail (u∗ , t, k)

(19)

This is the expression for the head velocity of the slender,
flexible microswimmer in the coordinate frame attached to
the tip of the head. In the following section we shall have a
look in detail the form of this equation and the topological
structure of the configuration space of this type of swimmers.
III. G EOMETRY OF THE CONFIGURATION SPACE
The geometry of the configuration space requires attention
for elegant and insightful solutions while studying the problem of locomotion using shape change. For such systems,
the configuration space is usually written as the product of
two manifolds. One is the base space or the shape space
M which describes the configuration of the internal shape
variables of the mechanism, and the other is a Lie group G
which represents the macro-position of the locomoting body
and is usually SE(3) or one of its submanifolds. The total
configuration space of the robot Q then naturally appears as
a product M × G. Such systems follow the topology of a
trivial principal fiber bundle, see [18]. Figure 6 shows an
explanatory figure of a fiber bundle. With such a separation
of the configuration space, locomotion is readily seen as
the means by which changes in shape affect the macro
position. We refer to [19], [20] for a detailed explanation
on the topology of locomoting systems. In particular, [21]
discusses how the flexible locomotion systems such as the
low Reynolds swimmers’ configuration space admit this
topological structure. We now show that the configuration
space of a slender flexible swimmer is a trivial principal
fiber bundle, defined as follows [20].
Definition : For Q a configuration manifold and G a Lie
group, a trivial principal fiber bundle with base M and
structure group G is a manifold Q = M × G with a free
left action of G on Q given by left translation in the group
variable: φh (x, g) = (x, hg) for x ∈ M and g ∈ G.
We model the slender flexible swimmer of length L as an
embedding ψ of interval I = [0, L] into R3
ψ : I 7→ R3

(20)

The configuration space Q for a swimming microorganism is
the space Q of embeddings of the cell membrane, represented

Fig. 6: Fiber Bundle

by I, into R3 . The structure group G in our case is the
Special Euclidean group SE(3). The shape space for the
organism is space of unlocated curves of finite length in R3 ,
which is the quotient of Q by the group SE(3) of rigid
motions acting on Q by composition. Since all the points
q ∈ Q are represented by (x, g) with x ∈ M and g ∈ G, Q
has global product structure of the form M × G. Moreover,
SE(3) acts via the left action as a matrix multiplication, and
has a single identity element, which is a 4×4 identity matrix.
Hence left action of the group, defined by Φh : (x, g) ∈
Q −→ (x, hg) is free, for x ∈ M and h, g ∈ G. Thus, the
configuration space of the slender flexible swimmer satisfies
a trivial principal fiber bundle structure.
A. Controllability analysis
In this section we prove a result on the controllability of
the flexible swimmer whose shape space M is the space of all
continuous curves in a plane. We refer to this swimmer as the
C 0 -planar flexible swimmer. The proof goes by construction
through the result of strong controllability of the planar
Purcell’s swimmer already presented in the literature. As
discussed in the introduction, the planar Purcell’s swimmer is
a the swimmer which has 3 straight, rigid and slender links
with 2 rotary joints, see figure 7. The angular velocity at
these 2 joints is the control input to the system. Clearly, the

Fig. 7: Purcell’s swimmer
Purcell’s swimmer is an example of a C 0 − flexible swimmer,
whose shape variable takes values in the space of all the
continuous curves having 3 interconnected line segments.
This shape space can also be written as SO(2) × SO(2),
which is of dimension 2. The system dynamics for the
Purcell’s swimmer can be written using the Cox theory itself
in a driftless control affine form as follows [8].
  

ẋ
I
=
u
(21)
ξ
−A(x)

Here A(x) is the local connection form defined at each x ∈
M . We recall that for a shape space M , its tangent space at a
point x ∈ M is denoted by Tx M , and the shape velocity ẋ =
(α̇1 , α̇2 ) ∈ Tx M . Thus A(x) : T M 7→ g, is the local form of
connection mapping the points in the tangent bundle T M to
the shape space M to the Lie algebra g of the structure group
G. I is a 2 × 2 identity matrix and u = [α̇1 , α̇2 ]T ∈ Tx M is
the control input at each point (α1 , α2 ) ∈ SO(2) × SO(2).
The connection form and the other notions mentioned here
have roots in geometric mechanics, see [22], [23] for details.
Since our configuration space is naturally split into a shape
and a structure group, we write a point in the configuration
space as q = (x, g) ∈ M × G = Q. We recall that for a time
parametrized shape x(t) ∈ M , the horizontal lift x∗ (t) ∈ Q
is a curve which projects to x(t) under the projection map
defining the principal fiber bundle and the components of its
tangent vectors ẋ∗ (t) ∈ Tq Q satisfy the equation (21). We
define following 2 controllability notions [20].
• A locomotion system is said to be strongly controllable
if, for any initial q0 = (x0 , g0 ) and final qf = (xf , gf ),
there exists a time T > 0 and a curve passing through
q0 satisfying x∗ (0) = q0 and x∗ (T ) = qf .
• A locomotion system is said to be weakly controllable if,
for any initial position g0 ∈ G, and final position gf ∈
G, and initial shape x0 ∈ M , there exists a time T > 0
and a curve in the base space x(t) satisfying x(0) =
x0 such that the horizontal lift of x(t) passing through
(x0 , g0 ) satisfies x∗ (0) = q0 and x∗ (T ) = (x(T ), gf ).
The principal fiber bundle structure in controllability analysis
gives rise to such strong and weak controllability notions
which define finer ideas of controllability for locomotion
systems. These notions are of practical relevance since many
times just reaching the desired group component without
strict requirement on shape of the system is sufficient. We
define following vector spaces in terms of the local connection form A(x), its curvature DA(x), its Lie derivative
Lz DA(x) and their successive Lie brackets [20]

Using the analytical computations, the literature has already
proven that for the planar Purcell’s swimmer the rank of
h2 ⊕ h3 is 3 on the entire space Q, which is the dimension
of the Lie algebra of its structure group SE(2). Hence the
swimmer satisfies the strong controllability conditions at all
the points [11], [24].
Now if we restrict the shape space of the flexible swimmer
to all C 0 curves, then all the shapes that the rigid 3-link
Purcell’s swimmer can take are contained in the shape space
M of C 0 −planar flexible swimmer. Moreover, the Purcell’s
swimmer mechanism when one of its outer link is considered
as the head is similar to the C 0 -planar flexible swimmer.
Also, weak controllability notion concerns with this head’s
motion on SE(2). The global strong controllability condition
for the planar Purcell’s swimmer means that the swimmer
with 3 straight rigid links can be manoeuvred from any
point q0 to q1 on its configuration space Q, where q0 , q1 ∈
Q = SO(2) × SO(2) × SE(2). Hence, the Lie algebra g of
the structure group G of the C 0 −planar flexible swimmer
is spanned by the Lie algebra of the control vector fields,
satisfying equation 23. This proves that the C 0 −planar
flexible swimmer is weakly controllable.
We highlight that the strong controllability of the flexible swimmer is neither guaranteed nor ruled out by this
approach of proof by construction since the shape space of
the C 0 −planar flexible swimmer by definition consists of all
C 0 curves, of which the C 0 curves corresponding to nonstraight link shapes are excluded in the proof. Nonetheless,
the weak controllability result is of practical relevance since
it means such a swimmer can be manoeuvred to achieve any
arbitrary group displacement of its head on SE(2).

Whereas the system is locally strongly controllable if and
only if
g = h2 ⊕ h3 ⊕ · · ·
(23)

V. CONCLUSIONS AND FUTURE WORK
We have presented a kinematic model of a slender, flexible
swimmer at low Reynolds number conditions. Cox theory

IV. S IMULATION RESULTS
In this section we present the response of the head motion
of the planar flexible swimmer subject to a bump function
ψ, which defines displacement of the swimmer body in
y direction of the with respect to the head frame. The
bump function travels from the swimmer head to tail over
simulation period of 15 seconds.
h1 = span{A(x)(X) : X ∈ Tx M },


(
−c2
x ∈ (−1 + c3 t, 1 + c3 t)
c1 exp 1−(x−c
2
h2 = span{DA(x)(X, Y ) : X, Y ∈ Tx M },
3 t)
Ψ(x, t) =
0
otherwise
h3 = span{LZ DA(x)(X, Y ) − [A(x)(Z), DA(x)(X, Y )],
(24)
[DA(x)(X, Y ), DA(x)(W, Z)] : W, X, Y, Z ∈ Tx M }
The values of the parameters of the bump function as c1 =
..
1
106 , c2 = 15, c3 = 15
. The resulting variation in the shape
.
space is shown in figures 8a to 8f, where the head is at
hk = span{LX ξ − [A(x)(X), ξ], [η, ξ] : X ∈ Tx M,
the left tip. This function defines the control input on the
ξ ∈ hk−1 , η ∈ h2 ⊕ · · · ⊕ hk−1 }
tangent bundle of the base space T M of the swimmer, which
Then a system defined on a trivial principal bundle Q is includes the shape Ψ and shape velocity Ψ̇ of the swimmer.
locally weakly controllable near q ∈ Q if and only if the The swimmer head length δ is taken to be 0.05 times the
space of the Lie algebra g of the structure group G is spanned swimmer length. The motion of the head is obtained using
the kinematic equation 19. The velocities and position of the
by the vector fields h1 , h2 , · · · as follows
head obtained are shown in figures 9a to 9d, where quantities
g = h1 ⊕ h2 ⊕ · · ·
(22) corresponding to the translational motion are dimensionless.

shape space. Also, given the weak controllability result in
this paper, motion planning strategies for point to point
reconfiguration is also a very relevant problem to explore
and has applications in robotics and in understanding the
microbial motion.
R EFERENCES
(a) Time = 0 sec

(b) Time = 3 sec

(c) Time = 6 sec

(d) Time = 9 sec

(e) Time = 12 sec

(f) Time = 15 sec

Fig. 8: Swimmer simulation against the moving bump function

(a) Translational velocity

(b) Translational position

(c) Rotational velocity

(d) Rotational position

Fig. 9: Swimmer head motion response to the bump function

is used to model the viscous forces as a function of the
instantaneous shape and shape velocity of the swimmer.
These forces along with the low inertia conditions are used
to obtain the expression of the velocity of the swimmer
head. We highlighted the peculiar topology of the principal
fiber bundle that this swimmer admits. The notion of weak
controllability is presented and it is shown that the planar,
slender flexible micro-swimmer is weakly controllable for
shape space consisting of the C 0 curves. An interesting
avenue for the future work is to explore if the system is
strong and weak controllable for C k -smooth curves in the

[1] L. E. Becker, S. A. Koehler, and H. A. Stone, “On self-propulsion of
micro-machines at low reynolds number: Purcell’s three-link swimmer,” Journal of fluid mechanics, vol. 490, pp. 15–35, 2003.
[2] R. Dreyfus, J. Baudry, M. L. Roper, M. Fermigier, H. A. Stone, and
J. Bibette, “Microscopic artificial swimmers,” Nature, vol. 437, no.
7060, pp. 862–865, 2005.
[3] Technovelgy, “Microrobot to swim like bacteria with flagellar
propeller: Science fiction in the news,” 2016, accessed 13August-2016. [Online]. Available: ”http://www.technovelgy.com/ct/
Science-Fiction-News.asp?NewsNum=845”
[4] A. Najafi and R. Golestanian, “Simple swimmer at low reynolds
number: Three linked spheres,” Physical Review E, vol. 69, no. 6,
p. 062901, 2004.
[5] N. Cohen and J. H. Boyle, “Swimming at low reynolds number: a
beginners guide to undulatory locomotion,” Contemporary Physics,
vol. 51, no. 2, pp. 103–123, 2010.
[6] S.
Lawler,
“Sex
cells:
sperm
that
crawl
instead
of
swim,”
http://theconversation.com/
sex-cells-sperm-that-crawl-instead-of-swim-7219, 2012, [Accessed
1-August-2016].
[7] Website, “”look for diagnosis”,” $https://lookfordiagnosis.com/mesh
info.php?term=spirillumlang=1$, 2016, [Accessed August-2016].
[8] R. L. Hatton and H. Choset, “Geometric swimming at low and high
reynolds numbers,” IEEE Transactions on Robotics, vol. 29, no. 3, pp.
615–624, 2013.
[9] E. M. Purcell, “Life at low reynolds number,” American journal of
physics, vol. 45, no. 1, pp. 3–11, 1977.
[10] H. C. Berg, E. coli in Motion. Springer Science & Business Media,
2008.
[11] S. Kadam and R. N. Banavar, “Geometric controllability of the purcells
swimmer and its symmetrized cousin,” IFAC-PapersOnLine, vol. 49,
no. 18, pp. 988–993, 2016.
[12] A. Shapere and F. Wilczek, “Self-propulsion at low reynolds number,”
Physical Review Letters, vol. 58, no. 20, p. 2051, 1987.
[13] A. L. Silvestre, “On the slow motion of a self-propelled rigid body
in a viscous incompressible fluid,” Journal of mathematical analysis
and applications, vol. 274, no. 1, pp. 203–227, 2002.
[14] G. P. Galdi, “On the steady self-propelled motion of a body in a
viscous incompressible fluid,” Archive for rational mechanics and
analysis, vol. 148, no. 1, pp. 53–88, 1999.
[15] J. Martı́n, T. Takahashi, and M. Tucsnak, “A control theoretic approach
to the swimming of microscopic organisms,” Quarterly of applied
mathematics, vol. 65, no. 3, pp. 405–424, 2007.
[16] A. N. Pressley, Elementary differential geometry. Springer Science
& Business Media, 2010.
[17] R. G. Cox, “The motion of long slender bodies in a viscous fluid part
1. general theory,” Journal of Fluid mechanics, vol. 44, no. 04, pp.
791–810, 1970.
[18] S. Kobayashi and K. Nomizu, Foundations of differential geometry.
New York, 1963, vol. 1, no. 2.
[19] A. M. Bloch, P. Krishnaprasad, J. E. Marsden, and R. M. Murray,
“Nonholonomic mechanical systems with symmetry,” Archive for
Rational Mechanics and Analysis, vol. 136, no. 1, pp. 21–99, 1996.
[20] S. D. Kelly and R. M. Murray, “Geometric phases and robotic
locomotion,” Journal of Robotic Systems, vol. 12, no. 6, pp. 417–431,
1995.
[21] R. Montgomery, A tour of subriemannian geometries, their geodesics
and applications. American Mathematical Soc., 2006, no. 91.
[22] A. M. Bloch, “Nonholonomic mechanics,” in Nonholonomic Mechanics and Control. Springer, 2003, pp. 207–276.
[23] D. D. Holm, T. Schmah, C. Stoica, and D. C. Ellis, Geometric
mechanics and symmetry: from finite to infinite dimensions. Oxford
University Press London, 2009.
[24] L. Giraldi, P. Martinon, and M. Zoppello, “Controllability and optimal
strokes for n-link microswimmer,” in Decision and Control (CDC),
2013 IEEE 52nd Annual Conference on. IEEE, 2013, pp. 3870–
3875.

