[论文解读] Numerical study of 3D vesicles under flow: discovery of new peculiar behaviors

Thierry Biben, Alexander Farutin+1 more|ArXiv.org|Dec 23, 2009

Blood properties and coagulation被引用 4

一句话总结

本研究通过三维数值分析揭示了剪切流中囊泡的两种新行为：随着剪切速率增加，振荡-呼吸（VB）模式区域显著拓宽；且翻滚（TB）与VB振荡周期对剪切速率具有强鲁棒性。结合边界积分模拟与解析建模，作者表明在大形变下，囊泡动力学中的高阶谐波变得至关重要，挑战了在强流条件下经典Keller-Skalak理论的有效性。

ABSTRACT

The study of vesicles under flow, a model system for red blood cells (RBCs), is an essential step in understanding various intricate dynamics exhibited by RBCs in vivo and in vitro. Quantitative 3D analyses of vesicles under flow are presented. The regions of parameters to produce tumbling (TB), tank-treating (TT), vacillating-breathing (VB) {and even Kayaking (K)} modes are determined. New qualitative features are found: (i) a significant widening of the VB mode region in parameter space upon increasing shear rate $\dotγ$ and (ii) a striking robustness of period of TB and VB with $\dotγ$. Analytical support is also provided. These findings shed new light on the dynamics of blood flow.

研究动机与目标

• 系统研究剪切流中囊泡的动力学行为，作为红细胞（RBCs）的模型。
• 识别并表征翻滚（TB）、转动（TT）、振荡-呼吸（VB）和滑行（K）模式的参数区域。
• 揭示囊泡动力学中此前未观测到的行为，特别是VB模式区域随剪切速率变化的展宽，以及TB与VB振荡周期对剪切速率的不变性。
• 评估Keller-Skalak（KS）理论在强流和大形变条件下的局限性。
• 发展并验证一种解析理论，以解释观测到的数值结果，特别是高阶谐波在形状动力学中的作用。

提出的方法

• 采用边界积分（BI）方法数值求解囊泡在简单剪切流中的Stokes方程。
• 引入Helfrich弯曲力以模拟膜弹性，包含弯曲模量κ和面积不可压缩性约束。
• 使用无量纲形变参数$ C_{\bar{k}} $表征囊泡在流场中的形状变化。
• 将模拟结果与Keller-Skalak（KS）理论进行比较，尤其关注囊泡主轴的角速度$ \dot{\theta} $演化。
• 应用摄动解析理论，解释大形变下高阶谐波（如$ \cos(4\theta) $）的出现。
• 进行长时间模拟以评估瞬态动力学及向稳态模式的弛豫过程，特别是在模式转变线附近（如TT-VB、TB-K）。

实验结果

研究问题

• RQ1 随着剪切速率$ \dot{\gamma} $增加，参数空间中振荡-呼吸（VB）模式区域的宽度如何演化？
• RQ2 为何翻滚（TB）与VB模式的振荡周期几乎与剪切速率无关，与以往模拟结果相反？
• RQ3 在强流或大囊泡形变条件下，Keller-Skalak理论的失效程度如何？
• RQ4 高阶谐波（如$ \cos(4\theta) $）在高度形变囊泡的动力学中起什么作用？
• RQ5 在模式转变线附近，囊泡弛豫至稳态动力学需要多长时间？该弛豫过程的时间尺度是什么？

主要发现

• 随着剪切速率$ \dot{\gamma} $增加，参数空间中发生振荡-呼吸（VB）模式的区域显著拓宽，这一行为未被先前的理论或模拟所捕捉。
• 翻滚（TB）与VB模式的振荡周期在剪切速率变化一个数量级的范围内几乎保持恒定，表明其动力学行为对$ \dot{\gamma} $具有强鲁棒性。
• Keller-Skalak理论仅在形变较小或粘度比$ \lambda $较大时能准确描述TB与TT动力学；在高$ \dot{\gamma} $或中等$ \lambda $时出现显著偏差。
• 在大形变（$ C_{\kappa} \geq 8 $）下，高阶谐波（尤其是$ \cos(4\theta) $）在角动量动力学中变得显著，标志着如TB-VB或TB-K分岔等复杂转变的开始。
• 在TT-VB转变线附近，弛豫至稳态动力学的时间尺度约为$ 30 \eta_{in}R^3 / \kappa $，对应约125次振荡，表明存在缓慢的瞬态阶段。
• 本文发展的解析理论成功解释了数值观测结果，特别是高阶项的出现以及振荡周期的鲁棒性。