微流体阻抗泵中非线性流动产生机制——非线性、共振与波动力学的交互作用

《Annals of Biomedical Engineering》：Nonlinear Flow Generation in Microfluidic Impedance Pumps

　　

专业术语翻译摘要
阻抗泵（Impedance pump）提供了一种无阀流动产生机制，其在微流控生物医学系统中具有潜在应用前景，涵盖芯片实验室（lab-on-chip）与器官芯片（organ-on-chip）等器件，此类系统通常借助膜驱动（membrane actuation）操纵流体流动。以往研究已提出阻抗泵的多重潜在机制，例如线性系统中的共振效应。然而，非线性、波动力学（wave dynamics）与共振对净流量（net flow rate, NFR）产生的相对贡献仍未得到充分理解，制约了此类系统的优化设计。

方法：研究人员在微型矩形膜驱动泵模型上开展了一系列综合性流固耦合（fluid-structure interaction, FSI）数值模拟。每次模拟均联立求解三维、时间相关、层流、不可压缩的Navier-Stokes方程与非线性大变形膜方程。该小尺度模型形成了低雷诺数（low-Reynolds number）与高沃默斯利数（high-Womersley number）条件，这是先前探索不足、但与微流控生物医学应用密切相关的参数区间。

结果：尽管以往研究分别考察了上述机制，但均未建立它们之间的明确关联。研究得出三项关键发现：（1）经由稳态流（steady streaming）的非线性效应直接导致净流量（NFR）的产生；（2）共振与波动力学能够放大速度振荡并增强非线性效应，但并不直接产生NFR；（3）最优挤压器（pincher）位置必须与驻波（standing-wave）波腹对齐，该结论此前未被认识到。上述发现调和了文献中看似矛盾的基于波与基于非线性的两种解释。

结论：本研究为无阀泵送中非线性、共振与波动力学之间的相互作用提供了新颖的机理见解。波腹定位准则以及对非线性、共振与波动力学交互关系的阐明，为微流控生物医学器件中阻抗泵的设计提供了定性指导原则。

论文解读文章

**研究背景与目的**

阻抗泵（Impedance pump），又称Liebau泵，是一种无需内部阀门或叶片即可产生流动的技术，通过在柔性管段的不对称位置进行周期性挤压（pinching），即可获得净流量（net flow rate, NFR）或平均压头。该机制在微流控生物医学系统中具有重要应用前景，包括芯片实验室（lab-on-chip）与器官芯片（organ-on-chip）平台，此类器件通常利用膜驱动（membrane actuation）方式控制流体。尽管已有大量研究探讨阻抗泵的物理机制，但关于非线性、波动力学（wave dynamics）与共振在净流量产生中的相对贡献仍存在显著争议。部分研究将NFR的产生归因于波动力学与共振现象，认为波在阻抗不匹配边界的反射与干涉是主要机制；另一些研究则从理论上证明，任何线性系统均无法产生净流量，非线性效应必不可少。两种解释看似矛盾，缺乏统一的理论框架，严重制约了阻抗泵的优化设计。为此，研究人员开展了系统的三维流固耦合（fluid-structure interaction, FSI）数值模拟研究，旨在阐明非线性、共振与波动力学三者之间的相互作用机制，为微流控生物医学器件中阻抗泵的设计提供定性指导。

**研究内容与总体结论**

研究人员对微型矩形膜驱动泵模型进行了大规模三维FSI数值模拟。该模型总长50 mm、宽10 mm、高1.5 mm，顶部对称布置长度为40 mm、厚度为0.25 mm的乳胶膜。研究采用三维、时间相关、层流、不可压缩的Navier-Stokes方程耦合非线性大变形膜方程，在低雷诺数（Reynolds number, Re）与高沃默斯利数（Womersley number, Wo）条件下，系统考察了挤压频率、挤压位置、挤压强度、挤压器数量及相位差等参数对NFR的影响。研究发现：非线性经由稳态流（steady streaming）直接导致NFR产生；共振与波动力学通过放大速度振荡增强非线性效应，但并不直接产生NFR；最优挤压器位置须与驻波（standing-wave）波腹对齐。这一成果调和了文献中基于波与基于非线性的矛盾解释。

**关键技术方法**

研究人员使用ANSYS求解器包，构建三维、时间相关、层流、不可压缩Navier-Stokes方程与瞬态非线性大变形膜方程的双向流固耦合数值模型。流体域约84,000个六面体单元，结构域5,760个单元，并经网格无关性与时间步独立性验证。NFR通过对10个完整挤压周期的体积流量进行时间平均获得；固有频率通过自由振动测试结合FFT分解识别；通过改变流体黏度调节Re以考察非线性效应。

**研究结果**

*流动场特征*：通过对比f=3.7 Hz（低NFR，0.5 ml/min）与f=53 Hz（高NFR，2.45 ml/min）两个代表性工况，研究人员发现高NFR工况下膜上形成近似驻波模式，速度幅值与膜位移幅值显著增大；低NFR工况下则为传播波。两种工况的平均速度均比均方根（RMS）速度低一个数量级，表明流动高度振荡，泵送效率较低。入口与出口流量在驻波工况下近似同相，在传播波工况下近似反相。

*NFR对挤压频率的依赖*：NFR随频率呈非单调变化，峰值出现在系统固有频率（f_n≈4.1 Hz）的整数倍附近（f/f_n=2.93、4.63、6.10、7.80、9.02、13.17、16.83），表明共振模式下泵送效果最佳。膜位移x-t图显示，大NFR工况通常对应驻波模式，小NFR工况对应传播波模式。例外的是f=19 Hz工况，该频率远离共振频率却产生较大NFR，且膜上呈传播波，证明共振并非产生大NFR的必要条件。

*雷诺数效应*：将流体黏度降低至标准值的1/10（Re×10）后，NFR显著增大且部分频率下方向发生反转；在f=53 Hz工况下，NFR随Re单调增加，在Re/100时NFR几乎消失，在Re×100时NFR约为标准工况的3倍。该结果证实非线性对流项对NFR幅值具有关键影响。

*双挤压器配置*：在f=53 Hz共振工况下加入第二个挤压器（x_p/L=0.12），当两挤压器间相位角约为180°±30°时，驻波模式得以维持，NFR增至2.65 ml/min，较同幅值（p₀=50 Pa）单挤压器的0.32 ml/min实现了每挤压器NFR约四倍提升。在f=3.7 Hz传播波工况下，双挤压器无增益，表明多挤压器的优势取决于其能否强化整体膜波响应，而非数量本身。

*功率与能量分析*：瞬时功率计算显示，大NFR工况下挤压压力波形与膜局部速度同步良好，能量可高效传递至流体；小NFR工况（如f=3.7 Hz）则存在明显的相位失配。f=56 Hz工况虽保持驻波形态，但膜波运动相对53 Hz工况发生时间偏移，导致同步质量下降与振荡幅值减小。NFR与挤压器单位时间做功量在广泛参数范围内呈近似单调关系，为非线性整流机制提供了直接的定量支撑。

**讨论与结论**

讨论部分指出，非线性（经由稳态流）是阻抗泵产生净流量的根本机制（除黏性阻力外），共振与波动力学的作用在于放大速度振荡幅值，为非线性整流提供更强输入。挤压器应放置在驻波波腹附近（如f=53 Hz时约在x_p/L=0.08、0.25、0.42处）以实现最大功传递，该准则此前未被认识。研究提出定性设计准则：（1）选择能产生大幅膜振荡与速度振荡的激励频率，通常接近共振条件；（2）致动器应位于膜位移驻波波腹区域；（3）激励波形应与局部膜运动同步以最大化正功传递；（4）多致动器应通过相位设置强化整体波模式；（5）运行工况应保证足够强的非线性对流效应以将振荡流整流为平均流。

研究结论可归纳为：非线性整流是产生净流量的根本机制，该结论得到理论证据（线性系统无法产流）与数值证据（NFR与挤压器做功量的定量相关性）的双重支持；共振与波动力学通过放大振荡流间接增强非线性整流效果。本研究澄清了文献中三个关键争议：非共振频率下NFR产生的原理（非线性本身即可产流，共振仅起增强作用）；基于波与基于非线性的解释实为互补而非竞争；最优挤压器位置须与驻波波腹对齐这一新的设计准则。