# Driving a Microswimmer with Wall-Induced Flow

Clément Moreau, Kenta Ishimoto

► To cite this version:

Clément Moreau, Kenta Ishimoto. Driving a Microswimmer with Wall-Induced Flow. Micromachines, 2021, 12 (9), pp.1025. 10.3390/mi12091025. hal-05050432

HAL Id: hal-05050432

https://hal.science/hal-05050432v1

Submitted on 29 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Driving a Microswimmer with Wall-Induced Flow

Clément Moreau* and Kenta Ishimoto

Research Institute for Mathematical Sciences, Kyoto University, Kyoto 606-8502, Japan

* Correspondence: cmoreau@kurims.kyoto-u.ac.jp

**Abstract:** Active walls such as cilia and bacteria carpets generate background flows that can influence the trajectories of microswimmers moving nearby. Recent advances in artificial magnetic cilia carpets offer the potentiality to use a similar wall-generated background flow to steer bio-hybrid microrobots. In this paper, we provide some ground theoretical and numerical work assessing the viability of this novel means of swimmer guidance by setting up a simple model of a spherical swimmer in an oscillatory flow and analysing it from the control theory viewpoint. We show a property of local controllability around the reference free trajectories and investigate the bang-bang structure of the control for time-optimal trajectories, with an estimation of the minimal time for suitable objectives. By direct simulation, we have demonstrated that the wall actuation can improve the wall-following transport by nearly 50%, which can be interpreted by synchronous flow structure. Although an open-loop control with a periodic bang-bang actuation loses some robustness and effectiveness, a feedback control is found to improve its robustness and effective transport, even with hydrodynamic wall-swimmer interactions. The results shed light on the potentialities of flow control and open the way to future experiments on swimmer guidance.

**Keywords:** microswimming; active wall; cilia carpet; controllability; optimal control

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Moreau, C.; Ishimoto, K. *Driving a Microswimmer with Wall-Induced Flow. Micromachines* **2021**, *12*, 1025. https://doi.org/10.3390/mi12091025

Academic Editors: Hiroaki Ito and Naoki Takeishi

Received: 4 August 2021
Accepted: 23 August 2021
Published: 27 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. Introduction

The motion of biological and artificial bodies in a fluid at microscopic scale, or microswimmers, has been a subject of growing interest in the last decades, with increasing understanding of the mechanisms at stake in order to efficiently swim at this scale and the rise of promising applications in the biomedical field, such as targeted drug delivery or non-invasive surgery performed by swimming microrobots. To serve these purposes, the need for propulsion and guidance of microdevices and biological particles is becoming increasingly prevalent.

Microswimmers seldom move in a perfectly quiescent fluid, and their motion is largely influenced by surrounding background flows. The general interactions of swimmers with background flows have been extensively considered, both experimentally and theoretically [1,2]. In particular, microswimmers often move in the presence of active biological walls generating an oscillatory flow, thanks to cilia beating in synchrony on walls [3,4] or on the body of microorganisms [5] or the collective motion of a bacterial carpet [6]. Active walls also play significant roles in reproductive processes such as peristaltic and ciliary pumping in the oviduct and are considered to promote sperm and ovum transport [7–12] and induce sperm rheotactic responses [13–15].

Furthermore, inspired by these examples of wall actuation in nature, the last few years have seen the emergence of artificial cilia carpet operated by magnetic fields [16–19], which open the way to practical wall-controlled guidance of both biological and artificial swimmers in the near future. The general idea of using an external flow as a means of guidance of a microswimmer has, to the best of our knowledge, only briefly been explored [20,21], in particular using a variable shear flow as a means of control [22]. It, however, constitutes a promising avenue to deal with situations where the use of more standard propelling
strategies, such as chemical propulsion or direct magnetic guidance, which require active engines or magnetic components embedded into the swimmer, is not possible or suitable.

Therefore, in this paper, we will address the possibility of guiding a generic microswimmer by modulating the amplitude of an oscillatory flow generated by a neighbouring wall, formulating a simple model and analysing its properties with the help of mathematical control theory and numerical simulations. In addition to laying a formal mathematical foundation that fosters understanding of microswimming at a theoretical level, the control-theory standpoint can provide both guidance for designing microdevices and an optimal way of controlling them. Former theoretical studies of microswimmer control produced powerful results of controllability (the ability for a swimmer to reach any target from an initial position) and optimal control (finding the best possible way of reaching the target with respect to a given objective function); to cite only a few, we mention the works about the general deformable body [23–27]; robots made of spheres [28] or rigid segments connected by elastic joints, whether controlled by their deformation [29–31] or by an external magnetic field [32–35]; and more recently, an ellipsoidal swimmer controlled by a shear flow [22]. In line with these studies, the present analysis will thus seek to establish the viability of wall-generated flow control as a means of guidance, while motivating and informing future considerations of complex fluid and swimmer control problems and providing guidelines for future experimental setups.

Hence, in the following, we will first describe our model of a spherical swimmer moving in an oscillatory background flow, as proposed in [36]. Seeing it as a control system, we will move on to examining the theoretical local controllability of the swimmer around its reference trajectories, therefore warranting the soundness of searching for optimal trajectories with classical optimal control tools. We address several practical objectives and provide explicit approximate expressions for the structure of the optimal control and for the minimal time when the swimmer is guided along a direction parallel to the wall. Finally, in the last section, we move on to practical considerations, numerically exploring the advantages of feedback control and studying the influence of the interactions between the wall and the swimmer when assuming it behaves as a squirmer.

## 2. Model

We consider a spherical swimmer of radius $a$, facing a direction $r$ with coordinates in the lab frame ($e_x, e_y, e_z$) and swimming at a constant speed $U^{\text{free}} = vr$, with $v \ge 0$. We assume that the swimmer swims in the reference plane ($e_x, e_z$), so we simply have $U^{\text{free}} = (-v \sin \theta, 0, v \cos \theta)$ and no angular velocity, namely $\Omega^{\text{free}} = \dot{\theta}e_y = 0$, with $\theta$ the angle between $e_z$ and $r$—see Figure 1.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Spherical swimmer in wall-induced oscillatory flow. The swimmer moves in the xz-plane and has propulsive velocity $v$ in the direction indicated by $r$, which makes an angle $\theta$ to the z axis of the laboratory frame.</figcaption></figure>
The swimmer is assumed to be swimming in the presence of a wall on z = 0 that generates an oscillatory flow field of the form $u(t) \sin(\omega t - kx)e_{x}$, typically modelling a metachronal wave created by a carpet of active cilia, with $u(t)$ being a time-varying amplitude that we will prescribe in order to influence the swimmer's motion. Following [36], the reciprocal theorem [37] allows us to nicely calculate the effect on this flow on the swimmer's motion:

$$U_x^{\text{wall}} = u(t)e^{-kz}(1 - kz) \sin(kx - \omega t), \quad (1)$$

$$U_z^{\text{wall}} = -u(t)e^{-kz}kz \cos(kx - \omega t), \quad (2)$$

$$\Omega_y^{\text{wall}} = \frac{1}{2} u(t) e^{-kz} k \sin(kx - \omega t). \quad (3)$$

This wall-generated flow is represented in Figure 1, with flow cells centred at height kz and of width $\pi/k$ rotating in alternate directions. Thus, the dynamics of the swimmer is governed by the following system of equations:

$$\left\{
\begin{aligned}
\dot{x} &= U_x^{\text{free}} + U_x^{\text{wall}} \\
\dot{z} &= U_z^{\text{free}} + U_z^{\text{wall}} \\
\dot{\theta} &= \Omega_y^{\text{wall}}.
\end{aligned}
\right.
\quad (4)$$

In this system, the swimmer-wall hydrodynamic interactions may be neglected, since the wall-induced flow is significantly larger in the practical situation, which will be further discussed in the later section.

Defining the state $X = [x, z, \theta]^T$, we will write System (4) under the generic control system form

$$\dot{X} = f(X, u(t)), \quad (5)$$

where f encodes the dynamics in (4). The time-varying amplitude $u(t)$, which we will refer to as the *control*, allows us to act on the swimmer's state $X$, i.e., its position and orientation. The purpose of this study is to determine the properties of this system from the control theory point of view, i.e., the capacity of $u$ to drive the swimmer to a required state, and how to achieve it in practice. We will start by addressing the *controllability* of system (4); in other words, the existence of controls that allow for steering it to a target state. With this theoretical analysis set up, we will then move on to finding the optimal controls for relevant objectives.

# 3. Results

## 3.1. Controllability Analysis

In this section, we examine the controllability properties of the flow-driven swimmer in a local sense around reference trajectories $X_{\text{ref}}$. The idea is to answer the following question, which constitutes an important theoretical feature for a controlled system: can we target any position and orientation within a neighbourhood of the swimmer path $X_{\text{ref}}$ using a control of “small” amplitude that stays close to a reference value $u_{\text{ref}}$? Here, we will focus on the particular and arguably more practical situation $u_{\text{ref}} = 0$, therefore seeking to control the position and orientation of the free swimmer by generating a flow of “small” amplitude with the neighbouring wall.

Consider an initial state $X_0 = [x_0, z_0, \theta_0]^T$ for the swimmer, and let $X_{\text{ref}} = [x_{\text{ref}}, z_{\text{ref}}, \theta_{\text{ref}}]^T$ be the solution of the system starting at $X_0$ with no flow, on the time interval $[0, T]$, which is simply given by $X_{\text{ref}}(t) = [x_0 - vt \sin \theta_0, z_0 + vt \cos \theta_0, \theta_0]^T$: with no external flow, the swimmer swims along a straight line.

Local controllability around $X_{\text{ref}}$ is schematised on Figure 2 and formally defined as follows [38]: System (4) is said to be locally controllable around the trajectory ($X_{\text{ref}}$) if, for all $\varepsilon > 0$, there exists neighbourhoods $V$ and $W$ of $X_0$ and $X_{\text{ref}}(T)$ such that, for all $X_i$ in $V$, and $X_f$ in $W$, there exists a bounded control $u$ defined on the time interval $[0, T]$, satisfying $|u(t)| \le \varepsilon$ for all $t \in [0, T]$, and such that the solution of System (4) starting at
$X_i$ satisfies $X(T) = X_f$. Figure 2 provides a graphical illustration of this concept, with an example of controlled trajectory (in a dotted line) linking an initial state in $V$ and a target state in $W$ around the reference trajectory (in red).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Schematic representation of local controllability around a reference trajectory $X_{\text{ref}}$, in which the swimmer moves along the dashed trajectory between the neighbourhoods $V$ and $W$ whilst the control $u(t)$ is kept small. Informally, the system is locally controllable around the trajectory $X_{\text{ref}}$ if such a dashed trajectory joining $X_i$ and $X_f$ exists for all states $X_i \in V$ and $X_f \in W$.</figcaption></figure>

To study the local controllability of System (4), we consider the linearised control system of System (4) around $X_{\text{ref}}$, with reference control $u_{\text{ref}} = 0$, that reads

$$ \dot{X} = A(t)X + B(t)u(t), \qquad (6) $$

with $A = \frac{\partial f(X, u)}{\partial X}(X_{\text{ref}}, 0)$ and $B = \frac{\partial f(X, u)}{\partial u}(X_{\text{ref}}, 0)$. It is well-known (see for instance [38] (Theorem 3.6)) that the nonlinear system is locally controllable around $X_{\text{ref}}$ if the linearised system (4) is (globally) controllable. To examine if this is the case, we can apply a simple algebraic “Kalman-type” condition for time-varying linear control systems [39,40]: let $B_1(t) = \dot{B} - AB$ and $B_2(t) = \dot{B}_1 - AB_1$; then, System (6) is controllable in time $T$ if and only if there exists $t \in [0, T]$ such that the matrix $(B(t)|B_1(t)|B_2(t))$ is invertible.

The determinant $\Delta$ of $(B(t)|B_1(t)|B_2(t))$ can be obtained through straightforward calculations:

$$ \Delta(t) = \frac{1}{2} e^{3k(z_0 + tv \sin \theta_0)} k^3 v \omega_c \sin \theta_0 \times [\omega_c(z_0 + tv \sin \theta_0) \cos(kx_0 - \omega_c t) + v \sin \theta_0 \sin(kx_0 - \omega_c t)], \quad (7) $$

with $\omega_c = \omega - kv \cos \theta_0$. We can see that the determinant vanishes if any of the parameters among $k, v, \omega_c$, or $\sin \theta_0$ is equal to 0. In those cases, the linearised system is not controllable. Upon further investigation, we can see that two of these cases are actually trivially not controllable, even globally: if $k = 0$, the swimmer's orientation $\theta$ remains always constant regardless of the control, indicating that a simple sliding motion of a wall is unable to control the swimmer, and the wave-like spatial inhomogeneity is necessary for the wall actuation, and if $v = 0$ (passive particle), the particle can only follow a single flow line. Thus, the wave-like wall actuation is only able to control an active particle. For the two other cases $\omega_c = 0$ (horizontal swimming speed equal to the velocity of the flow wave front) and $\sin \theta_0 = 0$ (vertical initial orientation), the controllability analysis alone does not allow for concluding with certainty that the nonlinear system (4) is not controllable. However, it suggests that these cases are unfavourable for driving the swimmer in an arbitrary direction and call for extra care when trying to steer its trajectory in practice.

If all of the above parameters are nonzero, then the determinant $\Delta(t)$ does not vanish (or, at worst, vanishes on isolated points), and therefore system (4) is locally controllable around $X_{\text{ref}}$.
Overall, the controllability analysis shows that, except in a few unfavourable cases described above, the flow generated by the active wall is able to locally influence the trajectory of the swimmer, both in positional and rotational dynamics. We have therefore theoretically validated the usability of this means of control and will carry on towards constructive determination of optimally controlled trajectories.

## 3.2. Optimal Control

For practical guidance of a swimmer using the active wall flow, we will study optimal control problems minimising the time required to reach a set target. Formally, given an initial state $X_0 = [x_0, z_0, \theta_0]^T$ and a target $X_1$, this problem can be stated as follows:

$$ \left\{ \begin{array}{ll} \min & T \\ \text{s.t.} & \dot{\mathbf{X}} = f(\mathbf{X}, u) \text{ on } [0, T], \\ & u \in [u_{\min}, u_{\max}], \\ & \mathbf{X}(0) = \mathbf{X}_0, \\ & \mathbf{X}(T) = \mathbf{X}_1. \end{array} \right. \qquad (8) $$

We will focus on three situations ("Problems") of interest that together allow for fully steering the trajectory of the swimmer:

1. A displacement parallel to the wall, from $x_0$ to $x_1$ at fixed distance $z_0$: $X_1 = [x_1, z_0, *]^T$;

2. A change of orientation, from $\theta_0$ to $\theta_1$ at fixed distance $z_0$: $X_1 = [*z_0, \theta_1]^T$;

3. A change of distance to the wall, from $z_0$ to $z_1$: $X_1 = [*z_1, *]^T$.

A star (*) denotes the variables for which the final value is not fixed. Problem 1 is a natural goal that one would wish to address in this situation, consisting of enhancing, with the external flow, the transportation of the swimmer along the horizontal direction. Problems 2 and 3 complete this main objective by adding the capacity of steering the swimmer close to or away from the wall, addressing the practical issue of trapping a swimmer near the wall or letting it escape. Note that Problem 3 is less constrained than Problems 1 and 2, as only one variable is fixed at a final time; for that reason, it will be partly contained in Problem 2, as we will see later. As we observe in the following, the capacity to control the swimmer's orientation via Problem 2 is also useful for guidance along the wall (Problem 1) for a long distance by occasionally helping the swimmer to face back in its goal direction.

Using classical optimal control theory, we address the time-optimal problem (8) by writing the associated Hamiltonian $H$:

$$ \begin{aligned} H = & \mathbf{p} \cdot \dot{\mathbf{X}} + p^0, \\ = & u(t)e^{-kz} \left[ \left( p_x(1-kz) + \frac{p_\theta}{2} \right) \sin(kx - \omega t) - p_ykz \cos(kx - \omega t) \right] \\ & + v \left( p_y \cos \theta - p_x \sin \theta \right) + p^0, \end{aligned} \qquad (9) $$

where $p^0 \le 0$, and $\mathbf{p} = [p_x \ p_y \ p_\theta]^T$ is non-trivial and a solution to the adjoint system $\dot{\mathbf{p}} = \frac{\partial H}{\partial \mathbf{X}}$. By virtue of the Pontryagin maximum principle, $u$ must maximise $H$ along the trajectory, from which we deduce that it switches between values $u_{\min}$ and $u_{\max}$, depending on the sign of the term between brackets in (9), which we will denote by $h(\mathbf{p}, \mathbf{X})$. Hence the optimal control is simply expressed as a so-called *bang-bang* function:

$$ u(t) = \begin{cases} u_{\min} & \text{if } h(\mathbf{p}, \mathbf{X}) < 0, \\ u_{\max} & \text{if } h(\mathbf{p}, \mathbf{X}) > 0, \\ 0 & \text{if } h(\mathbf{p}, \mathbf{X}) = 0. \end{cases} \qquad (10) $$

Note that if $h(\mathbf{p}, \mathbf{X})$ vanishes on an interval, one can quickly deduce that $\mathbf{p} = 0$ on this interval, which contradicts the non-triviality of $\mathbf{p}$; thus, the third case of (10) is almost void
here, and the control always switches between $u_{\min}$ and $u_{\max}$. This is valid for any optimal problem including the aforementioned objectives 1, 2, and 3; however, different values and constraints on $X_0$ and $X_1$ yield different dynamics for the Hamiltonian, requiring ad hoc analysis for each problem in order to determine the optimal control.

Let us now set $u_{\min} = -u_{\max}$ and focus on Problem 1: moving the swimmer from $x_0$ to $x_1 > x_0$ in minimal time, while keeping it at a distance $z_0$. To further the analysis, we assume that the wave velocity $\omega/k$ is sufficiently larger than the swimmer speed $v$ and, from Equation (1), the characteristic parallel flow speed at distance $z_0$, defined by

$$V_0 = u_{\max}e^{-kz_0}(1 - kz_0), \quad (11)$$

and that $\theta$ is close to $-\pi/2$, so that the swimmer is moving towards the target $x_1$. These assumptions imply that the dynamics of the adjoint vector $p$, as well as those of $z$ and $\theta$, change slowly with respect to the timescale $2\pi/\omega$, so the sign changes of the bracketed term in (9) are approximately given by those of $\sin(kx - \omega t)$.

In other words, to maximise the Hamiltonian on a short timescale, the flow amplitude $u$ must switch between $u_{\min}$ and $u_{\max}$ every time $\sin(kx - \omega t)$ changes sign, which occurs when the swimmer crosses the border between two of the flow cells generated by the wall. The time $T_{\text{switch}}$ elapsed between two switches can then be estimated as the ratio between the wall wavelength $\frac{\pi}{k}$ and the difference between the wall flow velocity and swimmer effective horizontal velocity, yielding

$$T_{\text{switch}} = \frac{\pi}{\omega - k(v \sin \theta + v_{\text{flow}})}. \quad (12)$$

In (12), the flow velocity term $v_{\text{flow}}$ accounts for the average velocity gained by the swimmer while crossing the flow cell.

In order to estimate $v_{\text{flow}}$, we consider the leading-order approximation of $z$ and $\theta$ at the small timescale $T_{\text{switch}}$, hence simply assuming that the swimmer remains at a constant distance from the wall $z_0$ and orientation $\theta_0$. We then consider a time interval of length $T_{\text{switch}}$ in which the control is set to the constant value $u_{\max}$, which consequently allows us to explicitly solve the dynamics for $x$:

$$x(t) = \frac{1}{k} \left[ \omega t - 2 \arctan \left( \frac{\omega_0}{\omega} \tan \left( \frac{\omega_0(t - T_0)}{2} \right) - \frac{V_0 k}{\omega} \right) \right] + vt \sin \theta_0, \quad (13)$$

with $\omega_0 = \sqrt{\omega^2 - V_0^2 k^2}$, and $V_0$ is defined in (11). At the lowest order, the displacement contributed by the flow between two switches is therefore given by $\Delta x = \frac{2V_0}{\omega}$, and the average velocity $v_{\text{flow}}$ that the swimmer gains from the wall flow satisfies the equation

$$v_{\text{flow}} = \frac{\Delta x}{T_{\text{switch}}} = \frac{2V_0}{\pi} \left( 1 - \frac{k(v \sin \theta_0 + v_{\text{flow}})}{\omega} \right), \quad (14)$$

and hence we obtain, at the lowest order,

$$v_{\text{flow}} = \frac{2V_0}{\pi} \left[ 1 - \frac{k}{\omega} \left( v \sin \theta_0 + \frac{2V_0}{\pi} \right) \right], \quad T_{\text{switch}} = \frac{\pi}{\omega - \left( v \sin \theta_0 + \frac{2V_0}{\pi} \right)}. \quad (15)$$

Finally, the minimal time $T_{\min}$ for the swimmer to reach $x_1$ from $x_0$ can be estimated as

$$T_{\min} = \frac{x_1 - x_0}{v + v_{\text{flow}}}. \quad (16)$$
The expressions obtained above show a linear relationship between the switch time $T_{\text{switch}}^{-1}$ and the frequency $\omega$. Moreover, for small $V_0$, $T_{\text{min}}$ is minimal when $V_0$ is maximal, which is achieved at $z_0 = 2/k$, suggesting that this distance to the wall is optimal for parallel displacement purposes. There, setting $z_0 = 2/k$, we can give an estimation of the optimal horizontal velocity boost $b$ given by the flow

$$b = \frac{v_{\text{flow}}}{v} = \frac{2e^{-2}}{\pi} \cdot \frac{u_{\text{max}}}{v}, \qquad (17)$$

which amounts to approximately 8.5% if $u_{\text{max}} = v$, and up to a boost of 43% if $u_{\text{max}} = 5v$, as in Figure 3a.

To validate this analysis, we numerically computed the optimal controls for a range of parameters, using the CasADi framework [41] and the in-built ordinary differential equation solvers in MATLAB® (The Mathworks, Inc., Natick, MA, USA) [42] for the solutions of the control system (4). For this numerical assessment, we set the swimmer speed $v$ and the wave number $k$ to be of unit magnitude, which is equivalent to nondimensionalise the system (4) with respect to $v$ and the wall wave velocity $\omega/k$; then, in Figure 3, the value of $\omega$, $t$, and $V_0$ respectively represent, in non-dimensional variables, the ratios $\omega/kv$, $tkv$, and $V_0/v$. The range of $\omega$ is taken to match the regime in which we carried out our analysis ($\omega/k$ sufficiently larger than $v$, with a ratio of 4 seeming to be already satisfying when looking at Figure 3). Figure 3 displays the results and shows good agreement with the theory. Figure 3a,b shows the optimal trajectory and associated control function $u(t)$ for different values of $V_0$ and $\omega$, showing a characteristic aspect due to the successive switches in $u$, and confirming our prediction that the optimal time $T_{\text{min}}$ depends very little on the frequency $\omega$. The relationships between the switching time $T_{\text{switch}}$ and $\omega$ and between the optimal time $T_{\text{min}}$ and the characteristic speed $V_0$ are shown in Figure 3c, with remarkable agreement between the numerical values and the curves predicted by Equations (15) and (16).

<figure><img src="image_7.png" /><figcaption>Figure 3. Visualisation of the time-optimal trajectories for a displacement parallel to the wall with units of $k = v = 1$. Panels (a,b) display various trajectories obtained by solving the optimal control problem, with the swimmer going from $x_0 = 0$ to $x_1 = 2\pi/k$. On each pair of graphs, the left panel shows the path of the swimmer within the flow field, and the right panel shows the value of the control along time, with periodic switches. The trajectories for three different frequencies $\omega$ ($\omega = 4, 8, 12$) are displayed along with two different values of the characteristic speed $V_0$ (0.62 and 0.34, respectively, on panels (a,b)), nondimensionalised with respect to the swimmer speed $v$). The top graph of panel (c) shows the relationship between the control switch time $T_{\text{switch}}$ and the frequency $\omega$. The blue and red dots correspond to numerically observed values, while the lines plot the expression obtained in Equation (15). The bottom plot shows the relationship between $V_0$ and the optimal time $T_{\text{min}}$. The dots show the numerical values obtained with various values of $V_0$ (varying $z_0$ and $u_{\text{max}}$), and the line is the prediction given by Equation (16), showing remarkable agreement.</figcaption></figure>
We now move on to our second control problem: changing the orientation of the swimmer from $\theta_0$ to $\theta_1$ with identical initial and final distance to the wall $z_0$. A few numerically computed trajectories realising this objective for various couples $(\theta_0, \theta_1)$ are displayed on Figure 4. While switches in the control occur with roughly the same period $T_{switch}$ as for the parallel displacement problem, in line with the short timescale analysis above, the observation of these trajectories in that case shows a more complex structure emerging for this control problem, with three distinctive phases in the swimmer's trajectory. First, the control makes the swimmer approach the wall and then rotates it once it is close enough to the wall, before driving it away from the wall again to reach the target distance $z_0$. During the first and last phases, the swimmer's orientation almost remains constant, seemingly indicating that an effective change of orientation cannot be obtained while the swimmer is far away from the wall. This can be understood by looking at the rotational component of the flow, $\Omega_y$ (see Equation (4)), along the optimal trajectories, plotted on the right-hand side of Figure 4. Approaching the wall allows $\Omega_y$ to reach large absolute values during the second phase and to generate net rotation by oscillating with non-zero average. This "approach-rotate-escape" strategy is consistently observed regardless of the initial and target orientation as well as the distance from the wall $z_0$. Hence, controlling the swimmer's orientation requires efficiently controlling its distance to the wall, which is the goal of Problem 3. We can therefore conclude that the solution to Problem 3 is contained in Problem 2, with the approach and escape phases in Problem 2 corresponding to the optimal motion of the swimmer towards or away from the wall. This analysis is confirmed by numerical observations: trajectories solving Problem 3 closely resemble the ones obtained for the approach and escape phases in the trajectories solving Problem 2.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Time-optimal trajectories for changing the swimmer's orientation with units of $k = v = 1$ as in Figure 3. In each row is displayed an optimal trajectory from an initial orientation $\theta_0$ to a target orientation $\theta_1$, while remaining at a distance $z_0$ from the wall. The path followed in the $xy$-plane is shown on the left panel, while the evolution of the orientation $\theta$ with respect to time is shown on the middle panel, and the local rotational speed of the flow $\Omega_y$ with respect to time is shown on the right panel. The values of initial and target orientation $(\theta_0, \theta_1)$ for rows (a-c) are respectively given by $(-\pi/4, -3\pi/4)$, $(-\pi, 0)$, and $(-\pi/2, \pi/2)$ and represented by red arrows on the trajectory plots on the left. For each case, we can distinguish three phases: approach to the wall with little variation of $\theta$, rotation close to the wall, and going away from the wall again.</figcaption></figure>
This exploration of controllability and optimal control properties of the wall-driven swimmer shows that the wall-generated flow can be expected to efficiently steer the swimmer's trajectory, both improving the swimmer's speed when swimming parallel to the wall and allowing it to change its orientation. However, this analysis assumes idealised control and neglects the interaction of the swimmer with the wall. In the next section, we address these issues by tackling the robustness of the control and studying the influence of added far-field wall interactions.

### 3.3. Robust Driving Strategies

The optimal control analysis conducted in the previous section showed that a bang-bang periodic control with period $T_{switch}$ allows for increasing the velocity of the swimmer while driving it parallel to the wall. Thinking of the practical application of this result, one may be tempted to use this control policy in an *open-loop* fashion, i.e., without changing the policy depending on the initial or current state of the swimmer, so that $u$ is given by

$$u(t) = u_{\max}\operatorname{sgn}\left(\sin\left(\pi\frac{t}{T_{\text{switch}}}\right)\right). \qquad (18)$$

This strategy has the obvious advantage of not requiring any measurement of the swimmer's position or orientation but will typically be sensitive to measurement errors and initial conditions for that same reason. In the case of the wall-controlled swimmer, we saw in the optimal control study that the periodic control needs to be synchronised with the swimmer position, so that it switches when the swimmer reaches the border between two flow cells. If the synchronisation fails because of a small error on setting $T_{switch}$, the swimmer's motion can become far from optimal. This is illustrated by the numerical simulation displayed in Figure 5a, where the open-loop strategy is applied to the swimmer for a short time and quickly shows desynchronisation, preventing it from gaining extra velocity from the flow.

To recover better synchronisation between the swimmer and the flow, we move on from this naive open-loop strategy to a feedback (or *closed-loop*) control policy, where the flow amplitude $u$ is allowed to depend on the state $X$. In line with the observation made in the optimal control study that the control should switch sign along with $\sin(kx - \omega t)$, we then define

$$u(t, X) = -u_{\max}\operatorname{sgn}(\sin(kx - \omega t)), \qquad (19)$$

and apply this policy to drive the swimmer parallel to the wall. The resulting trajectory is displayed in Figure 5b and resembles the optimal one obtained in the previous section with good synchronisation between the flow and the swimmer and a horizontal speed gain close to the optimal one. However, one can also observe in this figure that this simple feedback strategy defined in (19) tends to rotate the swimmer towards the wall over several periods and reduce its distance to the wall, with the risk of crashing into the wall after some time. Such a phenomenon highlights the necessity of combining both strategies of parallel displacement and change of orientation to compensate for the trajectory deviations emerging even when using a feedback control.

The analysis carried on so far demonstrates the caution needed to use the optimal controls in practice and provides solid grounds for the flow guidance of microswimmers, using a simple model that neglects the influence of the swimmer's motion on the surrounding fluid in the presence of the neighbouring wall. As briefly explained in Section 2, this is qualitatively justified by considering that the scale of these wall-swimmer interactions is typically significantly smaller than the free-swimming velocity of the swimmer, which itself is dominated by the background flow in our study; hence, the wall interactions can be neglected as a first approximation. To further validate our results, we will conclude this section by discussing the effect on the controlled trajectories when taking this interaction into account.
The wall-swimmer fluid interaction depends on the swimming strategy of the swimmer itself, which we will model from now on as a standard “squirmer” [43,44] of radius a with one characteristic parameter $B_2$.

Following [36], at the leading order (far-field approximation), the additional effect of the boundary on the squirmer’s motion is given by

$$U_x^b = -\frac{3a^2B_2}{40z^2} \sin 2\theta, \quad U_z^b = \frac{9a^2B_2}{16z^2} \cos 2\theta, \quad \Omega_y^b = -\frac{3a^2B_2}{40z^3} \sin 2\theta, \quad (20)$$

with a being the radius of the swimmer. By linearity, including these effects (20) in Equation (4) is simply done by adding them together, yielding the new system

$$\left\{
\begin{aligned}
\dot{x} &= U_x^{\text{free}} + U_x^{\text{wall}} + U_x^b, \\
\dot{z} &= U_z^{\text{free}} + U_z^{\text{wall}} + U_z^b, \\
\dot{\theta} &= \Omega_y^{\text{wall}} + \Omega_y^b.
\end{aligned}
\right.
\quad (21)$$

The swimmer is called a puller, a pusher, or a neutral swimmer if the characteristic parameter $B_2$ respectively satisfies $B_2 > 0$, $B_2 < 0$, or $B_2 = 0$. Note that the wall effects (20) vanish at the leading order in the case of a neutral swimmer, and we can reasonably neglect the higher order terms in the context of this study.

When $B_2 \neq 0$, it is well known [45] by cell-wall hydrodynamic interactions that pullers face towards neighbouring walls while pushers tend to swim parallel to them. We subsequently investigated, in Figure 5c-e, the influence of this additional term in the feedback control trajectories. As expected, in the pusher case ($B_2 = -20$), displayed in Figure 5c, the repulsive effect bends the trajectory away from the wall, although relatively weakly; a larger (and possibly unrealistic) absolute value of $B_2$ ($B_2 = -100$) yields a trajectory looking almost parallel to the wall in Figure 5d. Hence, this added repulsive effect seemingly provides increased stability for the control. On the contrary, in the case of a puller ($B_2 = 20$), the trajectory bends faster towards the wall, reducing the average horizontal speed of the swimmer and seemingly making the task of guiding it along the wall more challenging.

It is worth highlighting that the values of the parameter $B_2$ yielding a visible effect on the swimmer’s trajectories in Figure 5 are notably large with respect to the other parameters of the system, suggesting that the flow controls prescribed above are relatively robust with respect to the wall effect on the squirmer’s motion. This observation also corroborates, as predicted above, the choice of neglecting the boundary terms (20) for most of our theoretical study.

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Effect of feedback control and wall interactions on the controlled trajectories. All panels show the trajectory of the swimmer starting at $(0, 3, -\pi/2)$ until $T = 5$ with various parameter control policies (always with $u_{\text{max}} = 5$). The grey circle and red arrow respectively represent the swimmer and its orientation at the end of its motion. Panel (a) displays the trajectory obtained with the open-loop control policy defined in Equation (18), aimed at reproducing the optimal control determined in the previous section. The trajectory shows a bad synchronisation with the background flow, with little improvement of the swimmer’s average horizontal speed compared with the free swimming case. A feedback control policy is used in panel (b), defined as in Equation (19). The feedback helps recover synchronisation, but the swimmer’s orientation slowly changes, with risk of crashing into the wall. Panels (c-e) show the effect of the wall interaction term $B_2$ under the same feedback control policy. When $B_2 < 0$ (puller), as on panels (c,d), repulsion from the wall reduces the risk of crashing, all the more so the larger $|B_2|$ is; on the contrary, if $B_2 > 0$ (pusher), the swimmer is attracted to the wall.</figcaption></figure>
## 4. Discussion

In this article, we have theoretically and numerically addressed the control and guidance of a microswimmer using the flow generated by a wall. We have described a model of a simple spherical swimmer moving within a plane in the presence of a wall, which itself generates an oscillatory flow with a variable amplitude that we control in order to influence the positional and rotational dynamics of the swimmer and steer its trajectory. Using classical tools of control theory, we have examined the local controllability of the system around reference trajectories, thus partially validating its viability as a practical means of control. Further, we studied the minimal time problem for several suitable objectives, including moving along the wall, and obtained predictive results about the structure of the optimal controls and the estimated speed boost in this situation, which are remarkably confirmed by numerical simulations.

We then observed that the open-loop control strategy applied numerically has rather low robustness to initial conditions and measurement errors, due to the high level of accuracy required in the synchronisation between the swimmer position and the control switches. To compensate for the disadvantages of open-loop control, we briefly explored the possibility of using feedback control, with improved results approaching the idealised optimal solutions.

Finally, we discussed the validity of our model by taking into account the influence of the wall at far-field in the squirmer case. With our numerical observations confirming that the influence of the wall on the swimmer's trajectory remains modest, even for high values of the squirming parameters, we conclude that they can be safely neglected for practical purposes, though the repulsive effect of the wall on a puller may slightly help controlling it when moving parallel to the wall.

Though we have restricted this study to planar motion, our model dynamics are straightforwardly adaptable to three-dimensional motion in the presence of a wall along the (x,y)-plane by linearly superposing the two flows $u_1(t) \sin(\omega_1 t - k_1 x)e_x$ along the x direction and $u_2(t) \sin(\omega_2 t - k_2 y)e_y$ along the y direction, as well as the associated solutions for the swimmer, with a suitable set of 3D coordinates. Then, the motion in a given direction parallel to the wall should be obtained by tuning the controls $u_1$ and $u_2$ to generate flow in this direction only, hence reducing to the 2D case. With the goal of this study being to lay some simple ground guidelines for the novel idea of the guidance of swimmers by wall-induced flow, we leave further exploration of more complex, 3D-specific cases to future investigation.

In summary, in this paper, we have assessed the controllability and optimal control properties of a spherical swimmer controlled by the modulation of an oscillatory flow generated by a wall and described the key features of the optimal control for motion along the wall, notably its periodicity. We have shown that the control of a swimmer with such a setup is fully realisable and provided theoretical and practical guidelines, in particular regarding feedback control and the influence of wall interactions, aimed at calibrating future experimental work, with the hope of making efficient use of the emergent artificial cilia carpets operated by a magnetic field [16–19]. Our results on how oscillatory flows can be used to boost a microswimmer's speed could also help understand in more depth the mechanisms at stake in biological contexts featuring microorganisms swimming near active walls.

**Author Contributions:** Conceptualization, C.M. and K.I.; formal analysis, C.M.; investigation, C.M. and K.I.; writing—original draft preparation, C.M.; writing—review and editing, C.M. and K.I. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was partially supported by JSPS-KAKENHI for Young Researchers (18K13456) and JST, PRESTO Grant Number JPMJPR1921 (to K.I.). The APC was funded by K.I.

**Data Availability Statement:** The data presented in this study are openly available in the following Github repository: https://github.com/clementmoreau/wall-controlled-swimmer.
**Acknowledgments:** C.M. is a JSPS International Research Fellow (PE20021).

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Denissenko, P.; Kantsler, V.; Smith, D.J.; Kirkman-Brown, J. Human spermatozoa migration in microchannels reveals boundary- following navigation. *Proc. Natl. Acad. Sci. USA* **2012**, *109*, 8007–8010. [CrossRef] [PubMed]

2. Fu, H.C.; Powers, T.R.; Stocker, R. Bacterial rheotaxis. *Proc. Natl. Acad. Sci. USA* **2012**, *109*, 4780–4785. [CrossRef]

3. Kantsler, V.; Dunkel, J.; Polin, M.; Goldstein, R.E. Ciliary contact interactions dominate surface scattering of swimming eukaryotes. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 1187–1192. [CrossRef] [PubMed]

4. Man, Y.; Ling, F.; Kanso, E. Cilia oscillations. *Philos. Trans. R. Soc. B* **2020**, 375, 20190157. [CrossRef]

5. Gilpin, W.; Prakash, V.N.; Prakash, M. Vortex arrays and ciliary tangles underlie the feeding-swimming trade-off in starfish larvae. *Nat. Phys.* **2017**, *13*, 380–386. [CrossRef]

6. Darnton, N.; Turner, L.; Breuer, K.; Berg, H.C. Moving fluid with bacterial carpets. *Biophys. J.* **2004**, *86*, 1863–1870. [CrossRef]

7. Winet, H.; Bernstein, G.S.; Head, J. Observations on the response of human spermatozoa to gravity, boundaries and fluid shear. *Reproduction* **1984**, *70*, 511–523. [CrossRef]

8. Smith, D.; Gaffney, E.; Blake, J.; Kirkman-Brown, J. Human sperm accumulation near surfaces: A simulation study. *J. Fluid Mech.* **2009**, *621*, 289–320. [CrossRef]

9. Berke, A.P.; Turner, L.; Berg, H.C.; Lauga, E. Hydrodynamic attraction of swimming microorganisms by surfaces. *Phys. Rev. Lett.* **2008**, *101*, 038102. [CrossRef]

10. Aranda, V.; Cortez, R.; Fauci, L. A model of Stokesian peristalsis and vesicle transport in a three-dimensional closed cavity. *J. Biomech.* **2015**, *48*, 1631–1638. [CrossRef]

11. Ishimoto, K.; Gaffney, E.A. Mechanical tuning of mammalian sperm behaviour by hyperactivation, rheology and substrate adhesion: A numerical exploration. *J. R. Soc. Interface* **2016**, *13*, 20160633. [CrossRef] [PubMed]

12. Nosrati, R.; Driouchi, A.; Yip, C.M.; Sinton, D. Two-dimensional slither swimming of sperm within a micrometre of a surface. *Nat. Commun.* **2015**, *6*, 8703. [CrossRef] [PubMed]

13. Miki, K.; Clapham, D.E. Rheotaxis Guides Mammalian Sperm. *Curr. Biol.* **2013**, *23*, 443–452. [CrossRef] [PubMed]

14. Ishimoto, K.; Gaffney, E.A. Fluid flow and sperm guidance: A simulation study of hydrodynamic sperm rheotaxis. *J. R. Soc. Interface* **2015**, *12*, 20150172. [CrossRef]

15. Ishimoto, K. Guidance of microswimmers by wall and flow: Thigmotaxis and rheotaxis of unsteady squirmers in two and three dimensions. *Phys. Rev. E* **2017**, *96*, 043103. [CrossRef]

16. Meng, F.; Matsunaga, D.; Yeomans, J.M.; Golestanian, R. Magnetically-actuated artificial cilium: A simple theoretical model. *Soft Matter* **2019**, *15*, 3864–3871. [CrossRef]

17. Massana-Cid, H.; Meng, F.; Matsunaga, D.; Golestanian, R.; Tierno, P. Tunable self-healing of magnetically propelling colloidal carpets. *Nat. Commun.* **2019**, *10*, 2444. [CrossRef]

18. Matsunaga, D.; Hamilton, J.K.; Meng, F.; Bukin, N.; Martin, E.L.; Ogrin, F.Y.; Yeomans, J.M.; Golestanian, R. Controlling collective rotational patterns of magnetic rotors. *Nat. Commun.* **2019**, *10*, 4696. [CrossRef]

19. Gu, H.; Boehler, Q.; Cui, H.; Secchi, E.; Savorana, G.; De Marco, C.; Gervasoni, S.; Peyron, Q.; Huang, T.Y.; Pane, S.; et al. Magnetic cilia carpets with programmable metachronal waves. *Nat. Commun.* **2020**, *11*, 2637. [CrossRef]

20. Walker, B.J.; Ishimoto, K.; Wheeler, R.J.; Gaffney, E.A. Response of monoflagellate pullers to a shearing flow: A simulation study of microswimmer guidance. *Phys. Rev. E* **2018**, *98*, 063111. [CrossRef] [PubMed]

21. Khalil, I.S.M.; Klingner, A.; Hamed, Y.; Hassan, Y.S.; Misra, S. Controlled Noncontact Manipulation of Nonmagnetic Untethered Microbeads Orbiting Two-Tailed Soft Microrobot. *IEEE Trans. Robot.* **2020**, *6*, 1–13. [CrossRef]

22. Moreau, C.; Ishimoto, K.; Gaffney, E.A.; Walker, B.J. Control and controllability of microswimmers by a shearing flow. *R. Soc. Open Sci.* **2021**, *8*, 211141. [CrossRef] [PubMed]

23. Shapere, A.; Wilczek, F. Efficiencies of self-propulsion at low Reynolds number. *J. Fluid Mech.* **1989**, *198*, 587–599. [CrossRef]

24. Martín, J.; Takahashi, T.; Tucsnak, M. A control theoretic approach to the swimming of microscopic organisms. *Q. Appl. Math.* **2007**, *65*, 405–424. [CrossRef]

25. Lohéac, J.; Scheid, J.F.; Tucsnak, M. Controllability and time optimal control for low Reynolds numbers swimmers. *Acta Appl. Math.* **2013**, *123*, 175–200. [CrossRef]

26. Lohéac, J.; Munnier, A. Controllability of 3D low Reynolds number swimmers. *ESAIM Control Optim. Calc. Var.* **2014**, *20*, 236–268. [CrossRef]

27. Chambrion, T.; Giraldi, L.; Munnier, A. Optimal strokes for driftless swimmers: A general geometric approach. *ESAIM Control Optim. Calc. Var.* **2019**, *25*, 6. [CrossRef]

28. Alouges, F.; DeSimone, A.; Heltai, L.; Lefebvre-Lepot, A.; Merlet, B. Optimally swimming stokesian robots. *Discret. Contin. Dyn. Syst. Ser. B* **2013**, *18*, 1189–1215. [CrossRef]

29. Tam, D.; Hosoi, A.E. Optimal stroke patterns for Purcell’s three-link swimmer. *Phys. Rev. Lett.* **2007**, *98*, 068105. [CrossRef]

30. Alouges, F.; DeSimone, A.; Giraldi, L.; Zoppello, M. Self-propulsion of slender micro-swimmers by curvature control: N-link swimmers. *Int. J. Non-Linear Mech.* **2013**, *56*, 132–141. [CrossRef]
31. Giraldi, L.; Martinon, P.; Zoppello, M. Controllability and optimal strokes for N-link microswimmer. In Proceedings of the 52nd IEEE Conference on Decision and Control, Firenze, Italy, 10–13 December 2013; pp. 3870–3875.

32. Giraldi, L.; Pomet, J.B. Local controllability of the two-link magneto-elastic micro-swimmer. *IEEE Trans. Autom. Control* **2017**, *62*, 2512–2518. [CrossRef]

33. Giraldi, L.; Lissy, P.; Moreau, C.; Pomet, J.B. Addendum to “Local Controllability of the Two-Link Magneto-Elastic Microswimmer”. *IEEE Trans. Autom. Control* **2018**, *63*, 2303–2305. [CrossRef]

34. Moreau, C. Local controllability of a magnetized Purcell’s swimmer. *IEEE Control Syst. Lett.* **2019**, *3*, 637–642. [CrossRef]

35. El Alaoui-Faris, Y.; Pomet, J.B.; Régnier, S.; Giraldi, L. Optimal actuation of flagellar magnetic microswimmers. *Phys. Rev. E* **2020**, *101*, 042604. [CrossRef]

36. Papavassiliou, D.; Alexander, G.P. The many-body reciprocal theorem and swimmer hydrodynamics. *EPL Europhys. Lett.* **2015**, *110*, 44001. [CrossRef]

37. Stone, H.A.; Samuel, A.D. Propulsion of microorganisms by surface distortions. *Phys. Rev. Lett.* **1996**, *77*, 4102. [CrossRef]

38. Coron, J.M. *Control and Nonlinearity*; American Mathematical Society: Providence, RI, USA, 2007; Volume 136.

39. Silverman, L.; Meadows, H. Controllability and time-variable unilateral networks. *IEEE Trans. Circuit Theory* **1965**, *12*, 308–314. [CrossRef]

40. Chang, A. An algebraic characterization of controllability. *IEEE Trans. Autom. Control* **1965**, *10*, 112–113. [CrossRef]

41. Andersson, J.A.E.; Gillis, J.; Horn, G.; Rawlings, J.B.; Diehl, M. CasADi—A software framework for nonlinear optimization and optimal control. *Math. Program. Comput.* **2019**, *11*, 1–36. [CrossRef]

42. Shampine, L.F.; Reichelt, M.W. The MATLAB ODE Suite. *SIAM J. Sci. Comput.* **1997**, *18*, 1–22. [CrossRef]

43. Lighthill, M. On the squirming motion of nearly spherical deformable bodies through liquids at very small Reynolds numbers. *Commun. Pure Appl. Math.* **1952**, *5*, 109–118. [CrossRef]

44. Blake, J.R. A spherical envelope approach to ciliary propulsion. *J. Fluid Mech.* **1971**, *46*, 199–208. [CrossRef]

45. Lauga, E. *The Fluid Dynamics of Cell Motility*; Cambridge University Press: Cambridge, UK, 2020; Volume 62.