Squirmers with arbitrary shape and slip: modeling, simulation, and optimization

This paper proposes a Helmholtz decomposition-based framework to model arbitrary-shaped microswimmers, deriving their helical trajectories and rigid-body velocities while identifying optimal slip profiles that minimize power loss through a combination of partial and global optimization.

The Gist

Imagine a tiny, invisible world where microscopic creatures swim through thick, sticky fluid (like honey) rather than water. In this world, inertia doesn't exist; if you stop paddling, you stop moving instantly. These creatures are called microswimmers, and many of them, like bacteria or sperm, move by waving tiny hairs (cilia) on their surface.

This paper is about designing the perfect "swimming robot" for this sticky world, but with a twist: instead of being perfect spheres, these robots can be any shape—dumbbells, triangles, or weird blobs. The authors want to figure out two things:

1. How do they move? (Do they go straight, spin, or spiral?)
2. How can they swim most efficiently? (How do they waste the least amount of energy?)

Here is the breakdown of their findings using simple analogies:

1. The "Slip" Secret: How They Move Without Paddling

In our world, a boat moves by pushing water backward. But these tiny swimmers are too small for that. Instead, they use a "slip" mechanism. Imagine the surface of the swimmer is like a conveyor belt. The belt moves the fluid along the surface, and the reaction pushes the swimmer forward.

The authors created a new mathematical "language" (using something called Helmholtz decomposition) to describe exactly how this conveyor belt moves on any weird shape. Think of it like a universal remote control that can program the surface of any blob to wiggle in any pattern you want.

2. The Helix: Why They Don't Just Go Straight

If you have a perfectly round ball with a symmetrical conveyor belt, it swims in a straight line. But if the swimmer is lopsided or the conveyor belt is uneven (chiral), things get interesting.

The paper proves that any swimmer with a steady, unchanging slip pattern will naturally follow a circular helix path.

• The Analogy: Imagine a corkscrew or a spiral staircase. The swimmer isn't just moving forward; it's also spinning.
• The Result: If the spin and the forward motion are perfectly aligned, it spirals. If they are perfectly perpendicular, it might just spin in a circle. If there is no spin, it goes straight (a "flat" helix).

3. The Shape Matters: The "Squash" Effect

The authors looked at a specific shape: a prolate spheroid (like a rugby ball or a hot dog bun). They found that the shape of the swimmer changes how it reacts to its own "wiggles."

• The Analogy: Think of a long, skinny noodle vs. a short, fat dumpling.
• The Finding: If you make the swimmer very long and skinny (high aspect ratio), it becomes harder to spin it around its long axis, but easier to slide it sideways. The shape acts like a filter, deciding which movements are easy and which are hard.

4. The Optimization Game: Finding the Cheapest Path

The big question: How do we program the conveyor belt to use the least amount of energy?

The authors set up a two-step game:

• Step 1 (The Partial Fix): "Okay, we want the robot to go North. What is the most energy-efficient way to wiggle its surface to get there?"

  • Surprise: If the robot is symmetrical (like a dumbbell) and you ask it to go along its symmetry line, the most efficient way is to not spin at all. It just goes straight.
  • The Twist: If you ask it to go in a direction that breaks its symmetry (like a tilted dumbbell trying to go sideways), the most efficient way is actually to spin while it moves. The rotation helps it cut through the sticky fluid more easily.

• Step 2 (The Global Fix): "Forget the direction. What is the absolute best direction to go, and what is the best way to wiggle?"

  • The Discovery: For some weird, lopsided shapes, the absolute most efficient way to swim is to spin and spiral.
  • The Metaphor: Imagine a car. Usually, you drive straight to save gas. But for a very strange car with a flat tire on one side, it might actually be more fuel-efficient to drive in a circle than to try to drive straight. The authors found that for certain "weird" swimmer shapes, spiraling is the fuel-efficient choice.

5. Why This Matters

This isn't just about math; it's about building better nanobots for medicine.

• Imagine a tiny robot designed to swim through your blood vessels to deliver medicine to a tumor.
• If the robot is shaped like a sphere, it's easy to predict.
• But if we make it shaped like a dumbbell or a triangle to fit through narrow gaps, we need to know how to program its surface movements so it doesn't waste energy.
• This paper gives engineers the "instruction manual" to design these robots so they can swim efficiently, whether they need to go straight or spiral, depending on their shape.

In a nutshell: The authors figured out that for tiny, sticky-world swimmers, shape dictates destiny. If you want a robot to swim efficiently, you can't just pick any shape; you have to match the shape to the movement. Sometimes, the most efficient way to move forward is to spin in a circle.

Technical Summary

1. Problem Statement

The paper addresses the modeling and optimization of microswimmers (squirmers) with arbitrary shapes and arbitrary slip velocity profiles in a low Reynolds number (Stokes) fluid. While previous research has largely focused on spherical or axisymmetric squirmers, biological microorganisms and synthetic Janus particles often exhibit non-spherical geometries and chiral (non-axisymmetric) slip patterns. The authors aim to:

1. Establish a mathematical framework to describe slip velocities on arbitrary surfaces of spherical topology.
2. Characterize the resulting trajectories (specifically proving they are helical).
3. Solve the inverse problem: determining the slip velocity profile that minimizes power loss for a given shape, considering both prescribed and optimal directions of net motion.

2. Methodology

Mathematical Framework:

• Governing Equations: The fluid dynamics are governed by the incompressible Stokes equations. The swimmer is modeled as a rigid body with a tangential slip velocity $\mathbf{v}_S$ on its boundary $\Gamma$.
• Helmholtz Decomposition: To handle arbitrary shapes, the authors generalize the slip velocity definition using the Helmholtz decomposition on the surface. The slip velocity is expanded into a series of tangential basis functions derived from spherical harmonics ($Y_n^m$):
  $$\mathbf{v}_S = \nabla_\Gamma \Phi + \nabla_\Gamma \Psi \times \mathbf{n}$$
  where $\Phi$ and $\Psi$ are scalar potentials expanded in spherical harmonic modes with coefficients $a_n^m$ and $b_n^m$. This allows for the representation of both axisymmetric and chiral (non-axisymmetric) modes.
• Rigid Body Velocities: The translational ($\mathbf{U}$) and angular ($\mathbf{\Omega}$) velocities are determined using the Lorentz Reciprocal Theorem. This reduces the problem to solving a $6 \times 6$ linear system involving auxiliary Stokes flow problems (tractions for unit translation and rotation), avoiding the need to solve the full Stokes equations for every slip profile.

Trajectory Analysis:

• The authors prove that for a time-independent slip velocity, the swimmer's trajectory is a circular helix.
• The trajectory is parametrized by the rigid body velocities $\mathbf{U}$ and $\mathbf{\Omega}$.
  • If $\mathbf{\Omega} = 0$ or $\mathbf{\Omega} \parallel \mathbf{U}$, the path is a straight line.
  • If $\mathbf{\Omega} \times \mathbf{U} \neq 0$ and $\mathbf{\Omega} \cdot \mathbf{U} \neq 0$, the path is a circular helix with radius $|\mathbf{\Omega} \times \mathbf{U}|/|\mathbf{\Omega}|^2$ and pitch $2\pi|\mathbf{\Omega} \cdot \mathbf{U}|/|\mathbf{\Omega}|^2$.

Optimization Strategy:
The objective is to minimize the total power loss $P(\mathbf{v}_S)$ subject to a unit speed constraint along the direction of net motion ($\mathbf{W}$). The optimization is decomposed into two nested problems:

1. Partial Optimization: For a prescribed direction of net motion $\mathbf{W}$, find the optimal slip profile. This is a quadratic minimization problem with a closed-form solution.
2. Global Optimization: Find the optimal direction $\mathbf{W}_{opt}$ that minimizes the power loss obtained from the partial optimization. This is a non-linear optimization problem solved numerically.

Analytical Case Study:
The authors derived analytical expressions for the rigid body velocities of a prolate spheroid considering only the first-order ($n=1$) squirming modes, explicitly showing the dependence on the aspect ratio.

3. Key Contributions

1. Generalized Slip Framework: Introduction of a Helmholtz decomposition-based expansion for slip velocities on arbitrary shapes, extending the classical squirmer model beyond spheres and axisymmetric bodies.
2. Proof of Helical Trajectories: A rigorous proof that time-independent slip on arbitrary shapes results in circular helical trajectories, generalizing previous results limited to spherical chiral squirmers.
3. Efficient Computation: Demonstration that rigid body velocities can be computed via a simple $6 \times 6$ linear system using reciprocal theorem auxiliary problems, enabling rapid evaluation of slip profiles.
4. Optimization Insights:
   • Development of a two-stage optimization algorithm (partial then global) to minimize power loss.
   • Identification of the critical role of geometric symmetry in determining whether optimal motion is rotational or translational.

4. Key Results

• Aspect Ratio Effects (Prolate Spheroids):

  • As the aspect ratio increases (becoming more slender), the translational velocity along the symmetry axis ($U_z$) and rotational velocities ($\Omega_x, \Omega_y, \Omega_z$) decrease in magnitude.
  • Conversely, translational velocities perpendicular to the axis ($U_x, U_y$) increase.
  • First-order modes decouple: $a_1^m$ modes drive translation, while $b_1^m$ modes drive rotation.

• Symmetry and Optimal Motion:

  • Axisymmetric Shapes: The globally optimal motion is always a non-rotating straight line ($\mathbf{\Omega}=0$) along the symmetry axis.
  • Non-Axisymmetric Shapes with Symmetry Planes:
    • If the prescribed net motion direction $\mathbf{W}$ lies within a symmetry plane, the optimal motion is non-rotating ($\mathbf{\Omega}=0$).
    • If $\mathbf{W}$ lies outside a symmetry plane, the optimal motion for that specific direction is helical (rotational) to reduce power loss.
  • Global Optimization:
    • For many shapes (e.g., tilted dumbbell), the globally optimal direction $\mathbf{W}_{opt}$ aligns with a symmetry plane, resulting in a non-rotating straight line that is significantly more efficient (e.g., ~41% power reduction) than the best helical motion in a non-symmetric direction.
    • Crucial Finding: There exist shapes (specifically those with asymmetric lobes, e.g., $\rho(\theta, \phi) = 1 + 0.4 \sin^2\theta \sin(2\phi) + 0.2 \sin^3\theta \cos(3\phi)$) where the globally optimal motion is rotational ($\mathbf{\Omega}_{opt} \neq 0$). In these cases, the optimal direction $\mathbf{W}_{opt}$ does not lie in a symmetry plane.

• Power Loss Landscapes: The authors mapped the power loss landscape ($\hat{P}(\mathbf{W})$) for various shape parameters. They found that breaking cyclic symmetry (making lobes asymmetric) can shift the global minimum from a symmetry plane to an off-plane direction, necessitating rotation for maximum efficiency.

5. Significance

• Biological and Synthetic Relevance: The framework provides a theoretical basis for understanding how real microorganisms (which are rarely perfect spheres) and synthetic swimmers (like Janus particles) utilize shape and slip chirality to navigate efficiently.
• Design of Artificial Swimmers: The results offer specific design guidelines. To minimize energy consumption, engineers should either design swimmers with high symmetry to enable straight-line motion or intentionally introduce specific asymmetries if rotational propulsion is required for maneuverability in complex environments.
• Computational Efficiency: By reducing the optimization of infinite-dimensional slip fields to a finite-dimensional problem via the reciprocal theorem and basis function expansion, the method makes the design and analysis of complex microswimmers computationally tractable.
• Foundation for Future Work: The paper sets the stage for studying swimmers in confined geometries (pipes, near walls) and collective suspensions, where shape-induced chirality and helical trajectories will play a critical role in hydrodynamic interactions.