波动薄板在平行平板/通道中推进的文献逐条介绍
编制说明
本文档由文献清单 literature/confined_waving_sheet_refs.md 整理扩充而成，共收录 24 篇文献，按“受限形式 + 流体类型”分为四组：第一组（A）为牛顿流体中单壁/双壁受限行波板的解析研究（1965–2009），是本项目牛顿基线的直接来源；第二组（B）为黏弹性流体中的受限/近壁波动薄板（2007–2017），与本项目直接对标；第三组（C）为实验与其他受限几何（1978–2025），提供外部参照；第四组（D）为基线与方法学文献（1951–2025）。
每一条目按【发表时间】【第一作者】【题目】标注，正文为约 350 字的研究内容与主要结论介绍；公式一律以 LaTeX 源码书写、并以 Word 原生公式对象（OMML）排版，数学字体为 STIX Two Math。全部题名、卷期、页码与 DOI 已于 2026 年 10 月 18 日经出版方页面或权威索引（Cambridge
Core、AIP、IOPscience、Springer、Elsevier、arXiv）逐一核对，与原清单的差异已同步更正（Sauzade 等 2011 页码改为 1567–1573；CMES 2025 一文第一作者更正为 Rehman Ali Shah）。
一、牛顿流体中的受限波动薄板（解析研究）
本组为 1965–2009 年间用 Stokes 方程（含小惯性修正）解析处理“单壁或双壁平行板间行波板”的经典工作，是本项目牛顿受限基线的直接来源。
A1. Reynolds（1965）：近壁行波板的惯性修正
【发表时间】 1965 年
【第一作者】 A.
J. Reynolds
【题目】 The
swimming of minute organisms
【出处】 Journal
of Fluid Mechanics, 23(2): 241–260；DOI: 10.1017/S0022112065001337
研究内容。 该文首次把 Taylor 的无限长行波薄板模型用于“平行于固体平壁游动”的构型，在小振幅假设下求解壁面附近的速度场与推进速度，并把惯性作为修正纳入分析，讨论壁距、波数与振幅对结果的影响。它与 Tuck（1968）一起构成近壁游动问题最早的一批解析工作。
主要结论。 1. 壁面通过镜像作用改变薄板周围的涡量分布，使推进速度偏离无界值，壁距是控制修正量级的关键参数。 2. 计入惯性后，推进速度随雷诺数增大而下降；Reynolds
与 Tuck 因此被公认为“惯性抑制游动”一说的源头。 3. 该文的小振幅近壁解后来被 Shack & Lardner（1974，长波极限）与 Katz（1974，双壁）分别推广，是本项目牛顿受限基线的直接前驱。
A2. Tuck（1968）：近壁解的小惯性修正
【发表时间】 1968 年
【第一作者】 E.
O. Tuck
【题目】 A
note on a swimming problem
【出处】 Journal
of Fluid Mechanics, 31(2): 305–308；DOI: 10.1017/S0022112068000169
研究内容。 该文针对 Reynolds（1965）留下的近壁游动问题给出系统的小惯性修正，以无量纲黏度ζ=ηk2/(ωρ)
主要结论。 1. 惯性修正使推进速度低于 Stokes 预测，且修正量随ζ
减小（雷诺数增大）单调增长。 2. 文中给出推进效率E2
随ζ 的完整曲线；同一模型下挤压板解与势流解分居其上下两侧，可作为界定数值方法误差带的参照。 3. 该文确立的“壁面 + 小惯性”框架后来被 Felderhof（2009、2015）直接继承，用于两平行壁间的自推进与蠕动泵送研究，也是本项目检验求解器惯性残差的标准解析解。
A3. Katz（1974）：单壁与双壁通道的行波板解
【发表时间】 1974 年
【第一作者】 D.
F. Katz
【题目】 On
the propulsion of micro-organisms near solid boundaries
【出处】 Journal
of Fluid Mechanics, 64(1): 33–49；DOI: 10.1017/S0022112074001984
研究内容。 Katz 用无限长小振幅行波薄板同时求解两种构型：薄板平行于单一平壁游动，以及薄板位于两平壁之间的通道中央游动；薄板既可代表单根鞭毛，也可代表密排纤毛尖端的包络面。按波幅与壁距的相对大小给出小振幅与长波两支解析解，并计算推进速度与薄板做功速率。
主要结论。 1. 推进速度被波传播速度设定上界，不可能超过波速。 2. 壁距减小时，推进速度与薄板做功速率同时上升；这正是本项目牛顿基线Ua/c
随h/b
减小（壁距变小）而上升、由 0.064（h/b=12
）升至 0.487（h/b=2
）的解析来源。 3. 通过调节波速或波形，可以在保持功耗不变的情况下提高推进速度，其中降低拍击频率（波速）最为有效，是“能效约束下优化游动”的早期范例。
A4.
Shack & Lardner（1974）：通道内游动的长波解
【发表时间】 1974 年
【第一作者】 W.
J. Shack
【题目】 A
long wavelength solution for a microorganism swimming in a channel
【出处】 Bulletin
of Mathematical Biology, 36(4): 435–444；DOI: 10.1007/BF02464619
研究内容。 该文把通道内游动问题取长波极限：行波波长远大于通道间距（润滑近似），但允许波幅任意，在两刚性壁之间的低雷诺数 Stokes 流中求解二维自推进薄板的速度场，并给出推进速度与能量耗散率的解析表达式。
主要结论。 1. 推进速度与耗散率均可写成波幅、通道间距、波数与波速的显式函数，便于做渐近标度分析。 2. 长波极限下推进速度的主导项来自窄间隙内的剪切平衡，按(b/h)2
标度（b
为波幅、h
为半通道高），波幅只作次要修正；本项目润滑公式对照图中的红色基线即出自这一支。 3. 该解与 Katz（1974）的小振幅解互为补充，共同构成“通道内波动推进”的两支经典渐近，也是本项目窄通道段（h/b=2∼3
）最近的解析参照。
A5.
Smelser, Shack & Lardner（1974）：主动通道中的精子输运
【发表时间】 1974 年
【第一作者】 R.
E. Smelser
【题目】 The
swimming of spermatozoa in an active channel
【出处】 Journal
of Biomechanics, 7(4): 349–355；DOI: 10.1016/0021-9290(74)90030-X
研究内容。 该文把 Odeblad 提出的“宫颈黏液大分子波动与精子动态相互作用”输运机制数学化：几何理想化为二维通道内的自推进无限薄板，板以横波推进；通道壁按 Odeblad 的设想振动，分别取右行波、左行波与驻波三种形式，计算薄板推进速度、能量消耗与通道内流量。
主要结论。 1. 壁面振动改变通道内压强分布与流量，从而改变薄板的自推进速度，壁波与板波之间存在真实的双向耦合。 2. 模型结果与 Odeblad 的定性推断部分相符、部分相悖，说明仅靠“壁面波动助推”不足以解释精子在宫颈内的快速输运。 3. 这是“主动通道”（壁与游动体同时运动）问题的开创性工作，其分析框架被 Shukla 等（1988）与 Philip & Chandra（1995）沿用于不同流变模型。
A6. Felderhof（2009）：两平行壁间的游动与蠕动泵送
【发表时间】 2009 年
【第一作者】 B.
U. Felderhof
【题目】 Swimming
and peristaltic pumping between two plane parallel walls
【出处】 Journal
of Physics: Condensed Matter, 21(20): 204106；DOI: 10.1088/0953-8984/21/20/204106
研究内容。 Felderhof 用迁移率方法研究两平行平板之间、位于通道中央的无限大行波薄板的低雷诺数游动，并指出该问题与具有双波壁的平面通道内蠕动泵送在数学上同构，可纳入同一理论处理。
主要结论。 1. 给出受限条件下薄板游动速度与迁移率的统一表述，可与无界解直接比较，是“双壁等距薄板”构型最贴近的解析参照。 2. 证明“两条平行游动的柔性薄板”之间的流体动力学相互作用，等价于“双波壁平面通道内的蠕动泵送”，从而把自推进与泵送两类问题统一起来。 3. 文中方法后来被推广到含惯性的情形（Felderhof &
Jones 1994；Felderhof 2015），构成从 Stokes 到有限雷诺数壁面游动修正的完整链条。
A7. Blake（1971）：纤毛包络的无限长模型
【发表时间】 1971 年
【第一作者】 J.
R. Blake
【题目】 Infinite models for ciliary propulsion
【出处】 Journal
of Fluid Mechanics, 49(2): 209–222；DOI: 10.1017/S0022112071002027
研究内容。 Blake 提出“可伸展包络面”概念：用瞬时覆盖全部纤毛的表面代替密集分布的纤毛，把纤毛游动化为包络面上的行波问题。文章给出平面与圆柱两个无限长模型的小振幅推进速度解，并与作者此前的有限球形包络模型对比，还给出纤毛尖端的运动轨迹。
主要结论。 1. 两个无限长模型的推进速度都约为有限球形模型的两倍略多，说明“平面性”比“有限性”更重要；真实盘形生物（如 Opalina）的速度应落在两类模型之间。 2. 无限模型的最大推进速度约为波速的四分之一，与 Opalina 实测的约五分之一相符。 3. 该文与 Taylor（1951）共同构成“行波面推进”的理论基础，其包络思想也是本项目无限长薄板模型的物理依据。
二、黏弹性流体中的受限/近壁波动薄板
本组与项目直接对标：同一类无限长行波薄板，同一类 Oldroyd-B 本构，差别在受限形式（单壁/双壁/无壁）、有限长与否与数值方法。
B1. Lauga（2007）：黏弹性流体中的推进
【发表时间】 2007 年
【第一作者】 E.
Lauga
【题目】 Propulsion in a viscoelastic fluid
【出处】 Physics
of Fluids, 19: 083104；DOI: 10.1063/1.2751388
研究内容。 Lauga 用互易定理把 Taylor 的二维波动薄板推广到非线性非牛顿流体，重点计算 Oldroyd-B 黏弹性流体中的小振幅推进速度，推导黏弹性修正因子随 Deborah 数De=λω
的变化，并与牛顿结果并列比较。
主要结论。 1. 在波形给定的条件下，黏弹性总是使推进速度下降，U/U0=(1+βDe2)/(1+De2)
，De→∞
时趋于βU0
，即只剩溶剂贡献的速度。 2. 该结论是后续“黏弹性抑制游动”论断的解析源头，但其前提是“波形固定”。 3. 文章同时点明：一旦允许流体弹性改变游动体的波形（柔性体、有限振幅或多模态步态），结论可以反转；这一条件后来被 Riley & Lauga（2014）、Thomases
& Guy（2017）与 Salazar 等（2016）充分证实，也是解读本项目
De 扫描结果时必须先排除的混淆项。
B2.
Elfring & Lauga（2015）：复杂流体中游动的理论综述
【发表时间】 2015 年
【第一作者】 G.
J. Elfring
【题目】 Theory of locomotion through complex fluids
【出处】 载 S. E. Spagnolie 编 Complex Fluids in Biological
Systems（Biological
and Medical Physics, Biomedical Engineering），Springer，第 283–317 页；DOI: 10.1007/978-1-4939-2065-5_8
研究内容。 这是 Springer 专著的第 8 章，系统综述复杂流体中小尺度游动的理论：从牛顿 Stokes 游动的基本定理（扇贝定理、阻力矩阵）出发，引入各类本构关系，再给出简单几何模型游动体在复杂流体中的解析与数值结果，并直面当时相互矛盾的实验报道。
主要结论。 1. 对双壁之间的小振幅行波板，给出逐模态的黏弹性缩放因子(1+n2βDe2)/(1+n2De2)
：受限（壁）效应与黏弹性效应可因子化，壁面只改变牛顿基准，黏弹性只是再乘一个只依赖模态数n
与De
的因子。 2. 该章把当时“有的加速、有的减速”的矛盾几乎全部归结为“波形是否被弹性改变”，为整个领域提供了分类框架。 3. 本项目h/b×De
扫描的直接判据即来自该章的因子化预言；实测该比值从h/b=2
的 1.002 漂到h/b=12
的 1.026（约 2.5%），违反因子化，是当前最明确的不自洽信号。
B3.
Ives & Morozov（2017）：近壁黏弹性游动的机制
【发表时间】 2017 年
【第一作者】 T.
R. Ives
【题目】 The
mechanism of propulsion of a model microswimmer in a viscoelastic fluid next to
a solid boundary
【出处】 Physics
of Fluids, 29: 121612；DOI: 10.1063/1.4996839
研究内容。 该文首次把“近壁”与“黏弹性”两个效应放进同一模型：Taylor 波动薄板在 Oldroyd-B 类流体中平行于单一固体边界游动。作者重新审视 Elfring &
Lauga（2015）的小振幅结果并给出机制解释，同时建立可处理任意波幅与任意波形的谱方法数值解。
主要结论。 1. 小振幅下两种效应可以解耦解释：黏弹性使薄板减速，壁面使薄板加速，二者近似独立相乘。 2. 有限波幅下解耦失效：对某些参数组合，壁面的存在会使黏弹性反而提高游速，出现“超牛顿”加速。 3. 壁距、波幅与 Deborah 数共同决定加速/减速的符号，单凭无界结论外推受限问题不可靠；这与本项目在h/b=2∼3
观测到的“曲线变平甚至微升”定性一致。
B4.
Salazar, Roma & Ceniceros（2016）：IB-Oldroyd-B 有限长游动体
【发表时间】 2016 年
【第一作者】 D.
Salazar
【题目】 Numerical study of an inextensible, finite swimmer in Stokesian viscoelastic
flow
【出处】 Physics
of Fluids, 28(6): 063101；DOI: 10.1063/1.4953376
研究内容。 该文用浸入边界（IB）方法数值模拟 Stokes 流 Oldroyd-B 流体中的准不可伸长有限长二维薄板：目标波形由弹性罚力施加，允许薄板在流体应力作用下偏离目标波形，并采用非刚性 IB 格式消除强切向力带来的时间步限制，是离本项目方法最近的一篇文献。
主要结论。 1. 波形给定时，平均推进速度随De
增大单调下降并趋近渐近最小值，与空间周期小振幅渐近理论一致；作者同时提醒应力平流项的离散精度不足会显著改变游动结果。 2. 允许薄板足够柔软时单调性被打破，黏弹性可以增强推进，但仍慢于牛顿流体。 3. 最关键的是：在构象方程中人为加入黏弹性应力扩散项后，推进行为被“戏剧性改变”，可导致游速超过牛顿值。这与本项目构象松弛系数ωb
的敏感性实验完全对应——ωb=1
引入的人工扩散能把“弹性减速”翻转成“弹性加速”，说明数值格式细节与物理效应在同一量级上竞争。
B5.
Thomases & Guy（2017）：身体柔性与划程增强
【发表时间】 2017 年
【第一作者】 B.
Thomases
【题目】 The
role of body flexibility in stroke enhancements for finite-length undulatory
swimmers in viscoelastic fluids
【出处】 Journal
of Fluid Mechanics, 825: 109–132；DOI: 10.1017/jfm.2017.383
研究内容。 该文用渐近分析与数值模拟相结合，研究有限长起伏游动体在 Stokes–Oldroyd-B 流体中被动身体柔性如何改变运动学：分别改变身体弹性与流体弹性，预测波形变化，并计算只计入“波形改变”而不计入流体弹性直接作用的划程速度。
主要结论。 1. 对足够柔软的游动体，流体弹性导致的波形变化会产生更大振幅的划程（stroke boost），从而带来加速。 2. 但在所考察的步态中，流体弹性的直接作用总体上是减速的；加速来自划程增强而非弹性助推。 3. 在大振幅、强弹性流中会出现定性不同的机制：弹性应力在游动体附近高度聚集，造成剧烈减速。该文确立的“两条通道”分解是判断数值结果属于哪一类效应的标准工具。
B6.
Riley & Lauga（2014）：黏弹性抽吸与主动游动
【发表时间】 2014 年
【第一作者】 E.
E. Riley
【题目】 Enhanced active swimming in viscoelastic fluids
【出处】 EPL
(Europhysics Letters), 108(3): 34003；DOI: 10.1209/0295-5075/108/34003
研究内容。 针对“实验常见加速、理论总说减速”的矛盾，该文把游动体当作主动系统：以 Taylor 二维薄板为基础，把薄板形状写成外部流体应力、内部驱动力矩与被动弹性阻力三者的平衡，而不是像传统做法那样施加给定波形。
主要结论。 1. 上述平衡导致一个普适转变：身体偏刚性时黏弹性阻碍游动，足够柔性时黏弹性增强游动。 2. 物理图像是“黏弹性抽吸”：非牛顿正应力增大游动振幅，其增益克服黏弹性阻尼，从而加快游动。 3. 该文说明“减速/加速”之争在很大程度上是“给定波形
vs 柔性响应”之争，为 Espinosa-Garcia 等（2013）的实验提供了机制解释，也界定了本项目（波形固定的行波板）结论的适用范围。
B7.
Espinosa-Garcia, Lauga & Zenit（2013）：弹性增强游动的实验证据
【发表时间】 2013 年
【第一作者】 J.
Espinosa-Garcia
【题目】 Fluid
elasticity increases the locomotion of flexible swimmers
【出处】 Physics
of Fluids, 25: 031701；DOI: 10.1063/1.4795166
研究内容。 该文用实验直接检验黏弹性对柔性游动体的影响：游动体由磁驱动旋转的头部与柔性尾组成，分别置于黏度相近的牛顿流体与葡萄糖基 Boger 流体中，用粒子图像测速（PIV）测量流场与游速，并改变磁场频率以扫描 Deborah 数。
主要结论。 1. 流体弹性系统地提高游速，且增强幅度随 Deborah 数增大而增大。 2. PIV 显示游动体头部前后部分的剪切量存在显著差异，作者据此推断黏弹性正应力产生了一个指向游动方向的净弹性力。 3. 这是“弹性增强游动”最常被引用的实验证据之一，与 Lauga（2007）“固定波形必减速”形成对照，再次说明实验体系（有限长、柔性、大振幅）与理想化模型之间的差异正是矛盾来源。
三、实验与其他受限几何
本组把“通道/壁面”从刚性平板推广到波形壁、柔壁与微流控通道，并引入微极流体、电渗等多物理场，是本项目窄通道段与多物理场拓展的外部参照。
C1.
Gagnon 等（2023）：受限 C. elegans 的黏弹性游动实验
【发表时间】 2023 年 11 月（arXiv 预印本）
【第一作者】 D.
A. Gagnon
【题目】 Undulatory swimming in viscoelastic fluids under confinement
【出处】 arXiv:2311.05027 [physics.flu-dyn]；DOI: 10.48550/arXiv.2311.05027
研究内容。 该文用 C. elegans 线虫在微流控通道中游动，同时扫描“通道宽度（受限强度）”与“流体类型（黏性牛顿流体与黏弹性流体）”两个维度，测量游动学参数与周围速度场，并用弯曲弹性丝在黏弹性流体中的理论解释观测，是唯一同时覆盖“通道宽度 × De”的实验研究。
主要结论。 1. 在黏弹性流体中，弱或中等受限可以提高推进速度，体现壁面与弹性的协同。 2. 强受限下增强消失，游速反而低于无受限情形，对通道宽呈非单调依赖。 3. 理论分析表明，受限通过 Weissenberg 数控制线性响应区向非线性（指数）响应区的转变；速度场计算显示强受限下弹性应力在游动体周围增长，是减速的可能原因。该文可作本项目窄通道段（h/b<2.5
）的外部实验参照。
C2.
Shah 等（2025）：波浪通道内行波板的电渗-非牛顿耦合
【发表时间】 2025 年（6 月收稿、9 月录用、10 月刊出）
【第一作者】 R.
A. Shah（通讯作者 Z. Asghar）
【题目】 Non-Newtonian Electroosmotic Flow Effects on a Self-Propelled Undulating Sheet
in a Wavy Channel
【出处】 Computer
Modeling in Engineering & Sciences, 145(1): 753–778；DOI: 10.32604/cmes.2025.069177
研究内容。 该文研究波浪形通道内、非牛顿电解质溶液中的自推进行波板在电渗驱动下的运动：流体按 Carreau–Yasuda 剪切变稀模型本构，电势由 Poisson–Boltzmann 方程导出电渗体力，在 Stokes 流区求解流函数的四阶常微分方程，分析薄板速度与能量耗散。
主要结论。 1. 电渗效应与流体黏弹性共同决定薄板推进速度，两种机制通过改变通道内压强梯度相互耦合，不能简单叠加。 2. 文章系统讨论了电渗参数与流变参数对速度及能量耗散的联合影响，给出可供工程设计使用的参数依赖关系。 3. 该文把“行波板 + 剪切变稀流变 + 波形通道
+ 电渗驱动”结合起来，代表该方向向多物理场耦合拓展的趋势，可作为本项目后续引入电渗或剪切变稀效应的路线参考。
C3.
Shukla, Rao & Parihar（1978）：宫颈黏液的周边层效应
【发表时间】 1978 年
【第一作者】 J.
B. Shukla
【题目】 Swimming of spermatozoa in cervix: effects of dynamical interaction and
peripheral layer viscosity
【出处】 Journal
of Biomechanics, 11: 15–19；DOI: 10.1016/0021-9290(78)90039-8
研究内容。 该文把精子在宫颈内的游动理想化为二维波动薄板在充液通道中的自推进，重点考察宫颈黏液的“外周层”结构：将黏液处理为具有不同黏度的分层流体，分析外周层黏度以及壁面与游动体之间的动力相互作用对推进速度的影响。
主要结论。 1. 外周层黏度是控制推进速度的关键参数：低黏度外周层有利于薄板推进，与排卵期高水化黏液中精子游速更快的生理观察一致。 2. 壁面—游动体动力相互作用使速度明显偏离无界解，即使在中等通道宽度下也不可忽略。 3. 该文开创的“分层流体 + 宫颈几何”路线直接启发了后续一系列宫颈黏液精子输运研究（Sinha 等 1982；Shukla 等 1988；Philip &
Chandra 1995），也是把“通道”与“生物黏液”结合起来的最早一批工作。
C4.
Shukla 等（1988）：壁面蠕动波与纵波的作用
【发表时间】 1988 年
【第一作者】 J.
B. Shukla
【题目】 Effects
of peristaltic and longitudinal wave motion of the channel wall on movement of
micro-organisms: application to spermatozoa transport
【出处】 Journal
of Biomechanics, 21: 947–954；DOI: 10.1016/0021-9290(88)90133-9
研究内容。 该文在 Smelser 等（1974）“主动通道”框架下引入更真实的壁面运动学：通道壁同时作蠕动波与纵波运动，研究壁面运动对通道内微生物推进与输运的作用，并以精子在生殖道内的输运为应用背景，给出速度与流量的参数化结果。
主要结论。 1. 壁面蠕动波改变通道内时均压强梯度与流量，与薄板推进速度呈非平庸耦合。 2. 纵波与横波作用机制不同：纵波直接改变对流输运，横波通过改变有效间隙影响推进。 3. 文章给出不同壁面波参数下的速度与流量曲线，为评估“壁面主动运动能否解释精子在生殖道内的快速输运”提供了定量依据。
C5.
Philip & Chandra（1995）：微极流体中的柔壁自推进
【发表时间】 1995 年
【第一作者】 D.
Philip
【题目】 Self-propulsion of spermatozoa in microcontinua: effect of transverse wave
motion of channel walls
【出处】 Archive
of Applied Mechanics, 66: 90–99；DOI: 10.1007/BF00786692
研究内容。 该文把宫颈黏液当作微极（micropolar）流体——允许流体微元具有独立转动，从而刻画黏液大分子链的结构效应——研究柔壁通道内二维波动薄板的自推进：壁面因肌肉活动作横向波，惯性略去，并假设稳态下壁波与板波同步。
主要结论。 1. 推进速度随流体偏离牛顿行为（微极性增强）而增大，即微结构本身可以“帮助”推进。 2. 壁面蠕动波则降低推进速度；作者据此提出一种生育控制思路：用生化手段在宫颈壁诱导蠕动波即可降低精子推进速度。 3. 该文与 Sinha 等（1982）共同确立“微极流体 + 宫颈管 + 行波板”支线；其“非牛顿性可加速”的结论与 Oldroyd-B 框架下的“黏弹性减速”方向相反，提示流变模型的选择本身就能决定加速/减速的符号。
C6.
Sinha, Singh & Prasad（1982）：微连续介质分析的开创工作
【发表时间】 1982 年
【第一作者】 P.
Sinha
【题目】 A
microcontinuum analysis of the self propulsion of the spermatozoa in the
cervical canal
【出处】 International
Journal of Engineering Science, 20: 1037–1048；DOI: 10.1016/0020-7225(82)90037-4
研究内容。 该文是微连续介质（micropolar）方法用于宫颈管精子自推进的开创性工作：把宫颈黏液建模为微极流体，精子建模为两刚性壁之间的二维行波薄板，在低雷诺数下取“波长远大于通道间距”的长波近似，但允许波幅任意，求推进速度与能量消耗的解析表达式。
主要结论。 1. 推进速度与薄板做功速率均可写成微极参数N
、L
以及波幅、间距、波数、波速的显式函数，并给出系统的参数依赖图。 2. 推进速度与做功速率均随N 增大、L 减小而上升，即微极性越强推进越快。 3. 该文首次把“微结构流变”引入受限波动推进问题；其结论与十年后 Oldroyd-B 框架下的“黏弹性减速”形成鲜明对照，说明对宫颈黏液这类流体，本构选择对预测方向具有决定性影响。
四、基线与方法学文献
本组提供本项目全部无量纲数与基准值的标准来源。
D1. Taylor（1951）：波动薄板模型的起点
【发表时间】 1951 年
【第一作者】 G.
I. Taylor
【题目】 Analysis of the swimming of microscopic organisms
【出处】 Proceedings
of the Royal Society of London A, 209: 447–461；DOI: 10.1098/rspa.1951.0218
研究内容。 Taylor 提出用二维无限长行波薄板模拟鞭毛/纤毛推进：薄板作横向正弦行进波，在无惯性（Stokes）流场中因往复剪切产生净推力。文章以波幅ε=kb
为小参数作扰动展开，计算稳态推进速度，并讨论波形与推进方向的关系。
主要结论。 1. 给出经典首项结果U/c=12(kb)2
：推进速度与流体黏度无关、与波幅平方成正比，且薄板逆着波传播方向前进。 2. 时间可逆（扇贝型）运动无法在 Stokes 流中产生净位移，行波是低雷诺数游动的必要条件。 3. 该文是全部后续工作（含本项目）的零阶基准；需注意12(kb)2
只是ε→0
的首项，有限波幅（本项目kb=0.377
）的真值为0.4997ε2-0.4704ε4≈0.0615
，比首项 0.0711 低约 13%，直接套用首项会系统性高估基准。
D2.
Sauzade, Elfring & Lauga（2011）：Taylor 级数的改进
【发表时间】 2011 年
【第一作者】 M.
Sauzade
【题目】 Taylor’s swimming sheet: Analysis and improvement of the perturbation series
【出处】 Physica
D, 240(20): 1567–1573；DOI:
10.1016/j.physd.2011.06.023
研究内容。 该文把 Taylor 的扰动级数系统化：证明展开可以形式地做到任意高阶，用符号与数值手段求得各阶系数，研究级数的收敛性，并以边界积分法的数值解为基准检验不同阶数截断的精度。
主要结论。 1.
Taylor 原始级数在波幅达到O(1)
时发散，直接截断不可靠。 2. 经适当变换后得到收敛性质大为改善的级数，与边界积分数值解在中等乃至较大波幅下均吻合良好。 3. 文章给出的有限波幅精确级数U/c=0.4997ε2-0.4704ε4
是当前最实用的无界基准：对本项目kb=0.377
给出U/c≈0.0615
，比 Taylor 首项低 13.4%——该差额是有限波幅的物理效应而非数值误差，构成全部新基线的零点。
D3. Li（2025）：惯性区波动薄板的瞬态运动
【发表时间】 2025 年
【第一作者】 G.
Li
【题目】 Transient motion of a swimming sheet in the inertial regime
【出处】 Journal
of Fluid Mechanics, 1011: A2；DOI: 10.1017/jfm.2025.382
研究内容。 该文把波动薄板问题从稳态推广到惯性区的瞬态：同时考虑流体惯性与游动体自身惯性，在 Laplace 域推导纵波与横波联合作用下薄板瞬态速度的解析表达，再通过数值反演得到起动全过程的速度时程，系统分析加速度性能。
主要结论。 1. 瞬态速度由三部分构成：“滑移”速度、流体对流贡献与水动力贡献，三者可在 Laplace 域分离。 2. 以加速率1/T
、板质量M
与雷诺数Re
分别表征非定常、对流与游动体惯性：强惯性下起动时间按TM2Re
（纵波）与TM2
（横波）标度；弱惯性下纵波趋于常数、横波按T
标度。 3. 横波可诱发速度过冲，提升爆发游动性能。对本项目（Re=0.02
）而言，该文确认瞬态惯性效应可忽略，同时也提供了检验求解器瞬态段的标准算例。
D4.
Lauga & Powers（2009）：微生物游动综述
【发表时间】 2009 年
【第一作者】 E.
Lauga
【题目】 The
hydrodynamics of swimming microorganisms
【出处】 Reports
on Progress in Physics, 72: 096601；DOI: 10.1088/0034-4885/72/9/096601
研究内容。 这是《Reports on Progress in Physics》的邀请综述，系统梳理微米及以下尺度游动的生物物理与力学原理：从低雷诺数流动的基本性质（阻力矩阵、流动奇点、净平动判据）讲起，回顾经典游动运动学计算、阻力理论（resistive-force theory）与细长体理论（slender-body theory），再讨论鞭毛驱动机制与前沿方向。
主要结论。 1. 确立了“低雷诺数游动 = 阻力各向异性 + 非互易运动”的物理框架，是理解全部行波板模型的入口。 2. 明确指出复杂流体中的游动在当时尚无统一结论，实验与理论互相矛盾，并把原因归结为模型理想化程度的差异。 3. 该综述是本领域引用最多的入门文献，本项目全部无量纲数定义（Re、De、β
）与模型选择（Taylor 板、Oldroyd-B）均可在此找到标准表述。
五、与本算例的对照要点
把上述 24 篇映射到本项目（双壁通道内无限长行波板、h/b=2∼12
、β=1/2
、kb=0.377、Re=0.02），有五条要点：
1. 牛顿基线来源：Katz（1974，A3）给出单壁/双壁两支解并明确“壁距减小，推进速度与功耗同时上升”，解释了本项目牛顿Ua/c
从h/b=12
的
0.064 升到h/b=2
的 0.487；Shack & Lardner（1974，A4）的长波解给出窄间隙段(b/h)2
标度的红色基线。
2. 最硬的不自洽信号：Elfring & Lauga（2015，B2）预言黏弹性修正与壁距因子化。本项目实测该比值从h/b=2
的 1.002 漂到h/b=12
的 1.026（约 2.5%），违反因子化；这是当前唯一无法用“定义差异”解释的偏差。
3. ωb
问题的文献先例：Salazar 等（2016，B4）证明构象方程的应力扩散项开/关可以直接翻转加速/减速符号，并提醒应力平流项精度不足会显著改变游动结果。本项目实测ωb=1→1.9
使Ua/c 变化
+6.8% 并翻转符号，与该文结论定量对应。
4. 结论适用边界：Lauga（2007，B1）的“固定波形必减速”与 Riley & Lauga（2014，B6）、Thomases
& Guy（2017，B5）的“柔性可加速”并不矛盾，区别在波形是否被弹性改变；本项目波形给定，应参照前者，但需报告ωb 敏感性，以免与格式效应混淆。
5. 外部实验参照：Gagnon 等（2023，C1）是唯一同时扫描“通道宽度 × De”的实验：弱/中受限增强、强受限失去增强。本项目h/b<2.5
段的曲线变平与该文定性一致，可作为窄通道段合理性的一项旁证。
附录：文献速查表
|
编号 |
年份 |
第一作者 |
来源 |
DOI |
|
A1 |
1965 |
A. |
J. |
10.1017/S0022112065001337 |
|
A2 |
1968 |
E. |
J. |
10.1017/S0022112068000169 |
|
A3 |
1974 |
D. |
J. |
10.1017/S0022112074001984 |
|
A4 |
1974 |
W. |
Bull. |
10.1007/BF02464619 |
|
A5 |
1974 |
R. |
J. |
10.1016/0021-9290(74)90030-X |
|
A6 |
2009 |
B. |
J. |
10.1088/0953-8984/21/20/204106 |
|
A7 |
1971 |
J. |
J. |
10.1017/S0022112071002027 |
|
B1 |
2007 |
E. |
Phys. |
10.1063/1.2751388 |
|
B2 |
2015 |
G. |
Springer 专著第 8 章, pp. 283–317 |
10.1007/978-1-4939-2065-5_8 |
|
B3 |
2017 |
T. |
Phys. |
10.1063/1.4996839 |
|
B4 |
2016 |
D. |
Phys. |
10.1063/1.4953376 |
|
B5 |
2017 |
B. |
J. |
10.1017/jfm.2017.383 |
|
B6 |
2014 |
E. |
EPL |
10.1209/0295-5075/108/34003 |
|
B7 |
2013 |
J. |
Phys. |
10.1063/1.4795166 |
|
C1 |
2023 |
D. |
arXiv:2311.05027（预印本） |
10.48550/arXiv.2311.05027 |
|
C2 |
2025 |
R. |
CMES |
10.32604/cmes.2025.069177 |
|
C3 |
1978 |
J. |
J. |
10.1016/0021-9290(78)90039-8 |
|
C4 |
1988 |
J. |
J. |
10.1016/0021-9290(88)90133-9 |
|
C5 |
1995 |
D. |
Arch. |
10.1007/BF00786692 |
|
C6 |
1982 |
P. |
Int. |
10.1016/0020-7225(82)90037-4 |
|
D1 |
1951 |
G. |
Proc. |
10.1098/rspa.1951.0218 |
|
D2 |
2011 |
M. |
Physica |
10.1016/j.physd.2011.06.023 |
|
D3 |
2025 |
G. |
J. |
10.1017/jfm.2025.382 |
|
D4 |
2009 |
E. |
Rep. |
10.1088/0034-4885/72/9/096601 |