On the rectification of oscillatory flows by flexible leaflets in a confined geometry

问题陈述
柔性通道及受限几何结构中的流固耦合（FSI）现象是生物传输与微流体设计的基础。虽然单个柔性结构诱导的局部非线性已被充分表征，但预测多个相互作用结构在时变流场下的集体行为仍是一个重大挑战。现有文献主要侧重于稳态机制或单体结构，往往忽略了由多个柔性元件相互作用产生的丰富的动态特性。具体而言，关于叶片阵列如何通过调节振荡流（在低雷诺数环境下打破斯托克斯流型的固有时间反演对称性）以产生净单向传输的机制，目前尚未得到充分理解。

研究方法
作者利用一个极简的二维模型来研究这一问题，该模型包含一个底部锚定有非对称叶片床的通道。系统通过顶壁的垂直振荡驱动，模拟挤压流。本研究采用了一种结合以下技术的全耦合数值方法：

1. 格子玻尔兹曼方法 (LBM)： 用于求解低马赫数极限下的不可压缩纳维-斯托克斯方程。为了处理具有时间依赖性的问题而不受苛刻的时间步长限制，研究采用了由多重网格加速的双时间步进（DTS）方案。
2. 浸没边界法 (IBM)： 用于将柔性叶片几何形状集成到流体域中。叶片被建模为绕支点旋转的刚性板，通过连接的旋转弹簧响应流体动力矩。
3. 解析框架： 开发了一个连续极限模型，用以隔离叶片弹性和密度对系统的影响。该模型利用局部力矩平衡和全局角动量守恒来预测稳态极限下的流向整流，并随后通过准静态近似将其扩展至振荡状态。

核心贡献与结果
本研究通过表征叶片弹性、密度与通道几何形状之间的相互作用，阐明了净流体传输的产生机制。主要发现包括：

• 弹性-粘性数 ($\eta$) 的作用： 系统行为受无量纲弹性-粘性数控制，该数定义为粘性流体动力与叶片恢复弹性力之比。净流量表现出对 $\eta$ 的非单调依赖关系：
  • 刚性叶片 ($\eta \ll 1$)： 形变可以忽略不计，导致极小的流向整流效果。
  • 最佳柔性： 中间的 $\eta$ 值使传输效率最大化。在此状态下，叶片发生足够的形变，在不丧失结构完整性的前提下，在收缩阶段与舒张阶段之间创造出几何不对称性。
  • 过度柔软的叶片 ($\eta \gg 1$)： 叶片发生过度形变（通常超过90度），降低了几何不对称性，从而削弱了净传输。
• 集体相互作用与密度： 高叶片密度能最大化集体相互作用和净传输。随着密度的增加，离散的叶片阵列趋向于连续极限，此时有效边界轮廓变得平滑。净流速随密度增加而上升直至饱和，这一行为可通过唯象标度函数 $Q(\eta, \phi) \sim \phi/(\phi + K(\eta))$ 来捕捉。
• 稳态与动态机制： 在稳态极限下，净流量与雷诺数和振荡频率无关，仅取决于 $\eta$ 和密度。然而，在完全振荡的机制下，系统在壁面运动与叶片响应之间存在相位滞后。这种动态延迟受无量次元时间尺度比 $T_r$ 控制，该比值代表叶片的粘性松弛时间与壁面振荡周期的比值。当 $T_r$ 显著时，产生的瞬态构型会削弱净整流传输，而这一现象是纯准静态模型无法捕捉到的。
• 几何依赖性： 虽然准静态模型能准确预测高密度极限下的流动，但它无法捕捉到模拟中所观察到的净流量对通道长度 ($L_x$) 的依赖性，这凸显了动态效应对短通道的重要性。

意义
本文建立了一个分析柔性阵列在粘性环境中如何进行动态相互作用的基础框架。通过开发经过全耦合 FSI 模拟验证的解析连续模型，作者提供了一种用于理解由柔性阵列驱动的流向整流的预测工具。该研究证明了集体相互作用可以打破斯托克斯流型的对称性以产生单向传输，为理解生物流体传输系统以及设计被动式微流体整流器提供了见解。通过确定最佳弹性-粘性数并表征动态相位滞后，本研究为控制柔性几何结构中的流动提供了具体的工程设计参数。