Energy-optimal strokes for multi-link microswimmers: Purcell's loops and Taylor's waves reconciled
Autor: | Yizhar Or, François Alouges, Laetitia Giraldi, Oren Wiezel, Antonio DeSimone |
---|---|
Rok vydání: | 2019 |
Předmět: |
media_common.quotation_subject
FOS: Physical sciences General Physics and Astronomy assymptotic analysis Condensed Matter - Soft Condensed Matter Inertia Ellipse 01 natural sciences Displacement (vector) 010305 fluids & plasmas 0103 physical sciences Physics - Biological Physics optimal strokes 010306 general physics energy efficiency Eigenvalues and eigenvectors media_common Physics Plane (geometry) Computer Science::Information Retrieval Mathematical analysis Fluid Dynamics (physics.flu-dyn) Physics - Fluid Dynamics microswimmer dynamics Optimal control Arbitrarily large Biological Physics (physics.bio-ph) Drag Soft Condensed Matter (cond-mat.soft) |
Zdroj: | New Journal of Physics. 21:043050 |
ISSN: | 1367-2630 |
DOI: | 10.1088/1367-2630/ab1142 |
Popis: | Micron-scale swimmers move in the realm of negligible inertia, dominated by viscous drag forces. In this paper, we formulate the leading-order dynamics of a slender multi-link (N-link) microswimmer assuming small-amplitude undulations about its straight configuration. The energy-optimal stroke to achieve a given prescribed displacement in a given time period is obtained as the largest eigenvalue solution of a constrained optimal control problem. Remarkably, the optimal stroke is an ellipse lying within a two-dimensional plane in the (N – 1)-dimensional space of joint angles, where N can be arbitrarily large. For large N, the optimal stroke is a traveling wave of bending, modulo edge effects. If the number of shape variables is small, we can consider the same problem when the prescribed displacement in one time period is large, and not attainable with small variations of the joint angles. The fully nonlinear optimal control problem is solved numerically for the cases N = 3 (Purcell’s three-link swimmer) and N = 5 showing that, as the prescribed displacement becomes small, the optimal solutions obtained using the small-amplitude assumption are recovered. We also show that, when the prescribed displacements become large, the picture is different. For N = 3 we recover the non-convex planar loops already known from previous studies. For N = 5 we obtain non-planar loops, raising the question of characterizing the geometry of complex high-dimensional loops. |
Databáze: | OpenAIRE |
Externí odkaz: |