Parking 3-sphere swimmer II. The long arm asymptotic regime

The paper carries on our previous investigations on the complementary version of Purcell's rotator: a low-Reynolds-number swimmer composed of three balls of equal radii. In the asymptotic regime of very long arms, the Stokes induced governing dynamics is derived and then experimented in the context of energy minimizing self-propulsion characterized in the first part of the paper.


INTRODUCTION
In his seminal paper [16], Purcell explains how at small Reynolds numbers any organism trying to swim using the reciprocal stroke of a scallop, which moves by opening and closing its valves, is condemned to go back to its original position at the end of one cycle. This observation leads to the question of finding the simplest mechanisms capable of self-propulsion at these scales; by this, we mean the ability to moving by performing a cyclic shape change, a stroke, in the absence of external forces. Several proposals have been put forward and analyzed (see, e.g., [6,7,13,15,16,18] and the review paper [12]).
In this paper, we focus on a very specific microswimmer: the complementary version of Purcell's three-sphere rotator (sPr 3 ) introduced in [13] and fully described in Section 2. This swimmer consists of three non-intersecting balls (B i ) i∈AE3 of Ê 3 centered at b i ∈ Ê 3 and of equal radii a > 0 (for n ∈ AE we set AE n := {1, . . . , n}). The three balls can move along three coplanar axes that mutually meet at a point c ∈ Ê 3 , the center, with fixed angles of 2π/3 one to another; this reflects a situation where the balls are linked together by very thin telescopic arms that can elongate (see Figure 1). The swimmer can freely rotate around c in the horizontal plane containing the arms, although owing to the symmetries of the system, it is forced to stay in this plane.
Full controllability of sPr 3 , as well as for a broader class of model swimmers, i.e., the ability of the swimmer to reach any point in the plane with any orientation, has been proved in [4], while analytical investigations on the optimal control problem have been the object of [5]. Compared to [4,5], we propose here a quantitative analysis and we want to stress that, by contrast to the earlier works on the topic (cf. [1,2,3,6,7,8,9,15,16,17]), here the presence of three control variables and three position variables makes the analysis more involved and rich.
The main aim of this second part is to put the optimality results proved in [5] into a concrete setting, specifically: the Stokes induced governing control system for sPr 3 in the asymptotic regime of very long arms. First, we derive closed-form expressions for the dynamics and use asymptotic analysis to simplify the results. Then, we focus on the analysis of energy-minimizing strokes, and we identify the optimal parameters of the control system in terms of the initial length of the arms and the radius of the three balls. Finally, we present numerical simulations that show the qualitative features of the optimal swimming style.

KINEMATICS AND DYNAMICS OF SPR 3
As the three balls are not allowed to rotate around their axes, the shape of the swimmer can be parametrized by the lengths ζ 1 , ζ 2 , ζ 3 of its three arms, measured from c to the center of each of the balls. Therefore, the possible geometrical configurations of the swimmer can be described by introducing two sets of variables: • The vector of shape variables ζ: between the balls are obtained (cf. (2) and (3)). The lower bound on M is imposed to exclude any overlap of the spheres.
• Position and orientation of sPr 3 in the plane are specified by the coordinates of the center c ∈ Ê 2 × {0}, and by the angle θ that one arm, e.g., the arm connected to B 1 , makes with the fixed direction z 1 . We refer to p = (c, θ) ∈ Ê 2 × Ê as the vector of position variables.
Precisely, without loss of generality, we assume that in its initial configuration, the three arms of the swimmer sit in the plane Ê 2 × {0}. In order to compute the position of the three balls, we take the vertices of the equilateral triangle defined as the convex hull of the unit vectors z 1 , z 2 , z 3 ∈ Ê 3 , with z 1 := (1, 0, 0) T , z 2 := R T (2π/3)z 1 , z 3 := R(2π/3)z 1 , and R(φ) the planar rotation through an angle φ ∈ Ê around the vectorê 3 = (0, 0, 1): Then, the center b i of the i-th ball of the swimmer is at position (cf. Figure 2) where, here and in the sequel, we identify Ê 2 with Ê 2 × {0} and, similarly, we identify the action of R in (1) with its twodimensional analog. Since the balls cannot intersect, the matrix The time evolution of the swimmer can be traced through the state variables (ζ, p) ∈ M × Ê 3 . For i ∈ AE 3 , the instantaneous velocity of the i-th sphere is obtained by differentiating relation (2) with respect to time with z ⊥ i :=R(π/2)z i . The viscous resistance of the arms is deemed negligible and, therefore, we assume that the fluid fills up the whole space outside the balls, that is, the exterior domain Ω:=Ê 3 \ ∪ 3 i=1B i . The geometry of Ω is uniquely determined by the common radius a of the three spheres, and by the matrix b = (b i ) i∈AE3 having as where u and p are, respectively, the velocity field and the pressure of the fluid, and µ is its viscosity. As the structure of the swimmer is deformable but made of rigid balls, the governing equations are subject to no-slip boundary conditions on the balls. Because of the linearity of Stokes equations, the vector u := (u 1 , u 2 , u 3 ) collecting the three velocities in (4) can be expressed in the algebraic form (cf. [11,15]) where H is the Oseen tensor (which depends on the viscosity) and f := (f 1 , f 2 , f 3 ) ∈ Ê 6 is the vector collecting the forces acting on the balls. Symmetry arguments show that in the long arm asymptotic regime the hydrodynamic relation (6) takes the form: where the stokeslet represents a fundamental solution of the Stokes system [10], and ν := 6πµa ∈ Ê + is the drag coefficient linking, at small Reynolds numbers, the force to the velocity of a spherical object of radius a ∈ Ê + immersed in a fluid of viscosity µ.
It will be convenient to rewrite (7) in the form where I := diag(I, I, I) is the 6 × 6 identity matrix, and L the mutual interaction matrix defined by

DYNAMICS OF SPR 3 IN THE LIMIT OF VERY LONG ARMS
Due to the negligible inertia, the total viscous force and torque exerted by the surrounding fluid on the swimmer must vanish. In other words, the dynamics is subject to the balance equations Here, the cross product stands for the determinant form on Ê 2 and the b i 's are given by (2). Clearly, for every i ∈ AE 3 , there and, therefore, the balance equations (11) can be expressed in the concise form where the matrix W is defined by .
We assume that the three arms of the swimmer have the same initial length ξ 0 ∈ Ê + with ξ 0 ≫ a, and we set ζ i := ξ 0 + ξ i with |ξ i | ≪ ξ 0 . We want to show that in the limit of very long arms, and at the leading order, the swimming problem for sPr 3 reduces to a control problem of the forṁ with ξ := (ξ 1 , ξ 2 , ξ 3 ), whose structural symmetries have been fully investigated in the first part of the paper (cf. [5]). First, since the vector of the velocities depends linearly both onξ andṗ, we can recast relations (4) in the form where X , Y are the shape matrices given by In the limit of large arms, the mutual interaction matrix becomes a perturbation of the diagonal part (1/ν)I and equation (9) can be inverted to give (at the leading order) by use of (15). Multiplying both members by W, and after simplifying by ν, we infer that (cf. (12)) with the convenient and not dangerous abuse of notation W(ξ, θ) := W(ζ, θ). This is of the desired form (14) with where, to shorten notation, we left understood the parameters ξ 0 , ξ and θ. Moreover, because of the invariance of Stokes equations under the group of rotations, according to [5] (Prop. 1), we can factorize the control system F in the form F (ξ, θ) = R(θ)F (ξ) with F (ξ) := F (ξ, 0), and thereforė Also (cf. [5] (Prop. 4)), in the limit of small strokes, i.e., in the regime |ξ|/ξ 0 < a/ξ 0 ≪ 1 (see also Section 5), we can expand F to leading order in ξ. This gives Here, a straightforward computation shows that F 0 := F (0) is given by with ϕ(a, ξ 0 ) : Instead, the first order correctors (A k ) k∈AE3 (cf. [5] (Corollary 1)) have a special structure which can be fully characterized in terms of four real parameters. Precisely, there exist α:=α(a, ξ 0 ), β:=β(a, ξ 0 ), γ:=γ(a, ξ 0 ), and λ:=λ(a, ξ 0 ), depending only on the radius a of the balls and on the initial common length ξ 0 of the arms of sPr 3 , such that and With the aid of a symbolic computation software and expanding F in terms of ξ in (19) we can identify the entries and get α(a, ξ 0 ) : and γ(a, ξ 0 ) := 1 ξ 2 We remark that only the skew-symmetric parts (M k ) k∈AE3 of the matrices (A k ) k∈AE3 contribute to a net displacement of the swimmer after one stroke (cf. [5]). For any ξ ∈ Ê 3 they can be expressed by the actions forming an orthogonal basis of Ê 3 .

OPTIMAL SWIMMING
Following the notion of swimming efficiency proposed by Lighthill in [14] (cf. also [9,17]), we adopt the following notion of kinematic optimality: energy minimizing strokes are those minimizing the kinetic energy dissipated during one stroke in order to reach a prescribed net displacement δp ∈ Ê 3 . In mathematical terms, the total energy dissipation due to a smooth stroke ζ : I → M, can be evaluated by considering the instantaneous power dissipated at time t ∈ I, defined by P(u) = f · u. We note thatṗ is linear inξ because of (14), and so are f and u due to (15) and (17). Thus P(u) turns out to be a quadratic form inξ that we write in the following form for a suitable matrix-valued function G that, by the rotational invariance of the problem, does not depend on θ.
At the leading order in the limit of small strokes (cf. [5] ( § 5)) the instantaneous power dissipated at time t ∈ I reads as P(u(t)) = G 0ξ (t) ·ξ(t), with G 0 := G(0), and the total energy dissipation associated with a stroke ζ : I → M is given by (recalling that ζ i := ξ 0 + ξ i ) It can be readily checked that, as derived in [5] ( § 5), the matrix G 0 is symmetric, positive-definite, and has the following special structure with the two parameters h, κ depending only on the ratio a/ξ 0 between the radius of the balls of sPr 3 , and on the common initial length of its arms. Again, a symbolic computation shows that It is convenient to denote by g 1 :=(κ − h) and g 2 :=(κ + 2h) the eigenvalues of G 0 . Note that g 1 is of multiplicity two. Their expanded expressions read as In [5] (Theorem 5.1) we proved that the stroke ζ : I → M that produces a prescribed change of position and orientation δp ∈ Ê 3 of the swimmer at the minimal cost G(ζ) is an ellipse of Ê 3 . This optimal stroke is given by where the vectors u, v ∈ Ê 3 can be fully computed from δp, the coefficients α, γ of the skew-symmetric matrices (M k ) k∈AE3 , and the eigenvalues of G 0 . Namely, as shown in [5] (Theorem 5.1), any minimizer is, in ξ, an ellipse of Ê 3 centered at the origin, and the minimum value of G is equal to |ω| where  More precisely, considering two orthogonal vectors ς 1 , ς 2 ∈ Ê 3 in the plane orthogonal to ω and such that |ς 1 | 2 = |ς 2 | 2 = |ω|, we can compute the vectors u and v in (38) via the relations with U = (τ i /|τ i |) i∈AE3 (cf. (31)) and Λ := diag(g 1 , g 1 , g 2 ).
Summarizing, at the leading order in the range of small strokes and very long arms, the governing dynamics of sPr 3 for energy minimizing strokes is given by (cf. (21)) with (A j ) j∈AE2 given by (24), and u, v ∈ Ê 3 given by (40). In particular, the angular velocity of the swimmer is constant in time and is zero when the prescribed net displacement δp is purely translational (δp 3 = 0). It is easily seen that energy minimizing net displacements along the x-axis direction are achieved via elliptic strokes contained in the plane orthogonal to the vector τ 1 . Similarly pure along-y (resp. along-θ) net displacements are achieved via elliptic strokes contained in the plane orthogonal to τ 2 (resp. to τ 3 ).
The results of numerical simulations of (41)-(42) when the control ζ is the optimal swimming strategy for a prescribed net displacement δp along the x, y and θ directions are shown in Figure 3. Although we put some effort in drawing pictures that give a good feeling of how the swimmer performs, we are aware that the dynamics can be better appreciated by watching a video rather than looking at static frames; in that regard, in the supplementary electronic material, it is possible to find a video demonstrating the motion traced by sPr 3 during optimal swimming.

CONCLUDING REMARKS
Note that, for j = 1, 2, we have that lim ξ0→∞ A j (a, ξ 0 ) = 0 and lim a→0 A j (a, ξ 0 ) = 0. However, since γ(0, ξ 0 ) = 1/ 6 √ 3ξ In other words, the asymptotic limit of very small balls differs from one of very long arms. This is understood by the presence of two fundamental geometric scales: the common radius a of the three balls, and the initial length ξ 0 of its arms. In this respect, the two following asymptotic regimes are different: where we have denoted by |ξ| the "average" stroke intensity.
• In the limit a/ξ 0 ≪ |ξ|/ξ 0 ≪ 1 the swimmer offers great resistance to a net displacement in the (x, y) coordinates, but it is strikingly still able to produce net angular displacements in the θ variable. • The second condition in (44) represents the limit of very long arms and is more interesting for the applications as it allows for both translations and rotations.

ACKNOWLEDGMENTS
This work was partially supported by the Labex LMH (grant ANR-11-LABX-0056-LMH) in the Programme des Investissements d'Avenir. Also, the second author acknowledges support from the Austrian Science Fund (FWF) through the special research program Taming complexity in partial differential systems (Grant SFB F65).