[论文解读] Energy-optimal strokes for multi-link microswimmers: Purcell's loops and Taylor's waves reconciled

Alouges, Fran\c{c}ois, Antonio DeSimone+3 more|arXiv (Cornell University)|Jan 15, 2018

Micro and Nano Robotics参考文献 43被引用 4

一句话总结

本文推导了在小振幅波动下，N链微游泳器的能量最优摆动模式，表明最优摆动在关节角空间中为平面椭圆，且当N较大时趋近于行进波。对于较大位移，数值解揭示了非平面、复杂的环路结构——通过几何控制理论与最优控制优化，统一了Purcell的环路驱动与Taylor的波浪驱动游泳范式。

ABSTRACT

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.

研究动机与目标

• 统一微游泳中两种主导范式：通过闭合环路实现非互易形变（Purcell）与行进波（Taylor）。
• 在小振幅波动下，识别细长多链微游泳器的能量最优摆动模式。
• 将分析扩展至N=3和N=5游泳器的大振幅、完全非线性摆动。
• 表征高维形貌空间中最优摆动的几何结构，特别是当N≥5时。
• 推导大N极限下位移与能量消耗的标度律。

提出的方法

• 利用阻力系数理论与小振幅近似，推导N链微游泳器的一阶动力学模型。
• 求解约束最优控制问题，以在单周期内实现预定位移的同时最小化能量消耗。
• 将最优摆动识别为关节角变化中二次型的最大特征值解。
• 运用几何控制理论，证明最优摆动位于(N-1)维关节角空间中的二维平面内。
• 使用Bocop优化工具箱，对N=3和N=5的情况数值求解完全非线性最优控制问题。
• 通过固定振幅与相位差的行进波步态，推导大N极限下位移与能量的标度律。

实验结果

研究问题

• RQ1 在小振幅驱动下，多链微游泳器的能量最优摆动在几何结构上呈现何种特征？
• RQ2 随着预定位移的增加，最优摆动如何从平面椭圆过渡为复杂非平面环路？
• RQ3 小振幅结果在小位移下与完全非线性解的收敛程度如何？
• RQ4 在大N极限下，位移与能量消耗的标度律是什么？
• RQ5 对于N≥5，最优摆动几何结构是否可超越平面或波浪形态进行表征？

主要发现

• 在小振幅摆动下，能量最优步态在(N-1)维关节角空间中始终为平面椭圆，与N无关。
• 随着N增大，最优摆动渐近趋近于弯曲行进波，振幅从中心向两端对称衰减。
• 对于大预定位移，三链游泳器的最优摆动为非凸平面环路，与先前结果一致。
• 对于五链游泳器，最优摆动为非平面结构，展现出复杂高维几何特征，暗示最优步态中存在新型结构特性。
• 当预定位移较小时，小振幅与完全非线性解均收敛于同一平面椭圆。
• 标度律表明，大N下位移与ε²l成正比，能量与ε²Nl³/T成正比，且在固定链长下能量消耗随N线性增长。