4. 微通道内细胞流动基础

在微流控芯片中,细胞所处的流场与宏观管道有本质区别。通道特征尺寸通常在10~500 μm量级,流速多在0.1~100 mm/s范围,对应的雷诺数极低,流动几乎总是层流。理解层流的速度剖面、剪切力分布、细胞停留时间以及细胞在流场中的受力,是设计捕获结构、分选通道和评估细胞活性的基础。本节围绕这五个核心概念展开,并给出可计算的示例参数。

一、层流:微通道中的默认流态

层流指流体分层流动、各层之间不发生宏观混合的流态。在圆管中,当雷诺数 Re < 2300 时通常为层流;在矩形微通道中,临界雷诺数更高,实际芯片内 Re 常低于 10,因此层流是默认状态。层流的核心特征是:

  • 速度剖面稳定,不随时间随机脉动;
  • 流体沿流线运动,横向混合仅靠扩散;
  • 通道中心速度最高,靠近壁面速度趋近于零。

对于宽度 w、高度 h 的矩形通道,若 w ≫ h,可近似为两平行平板间的流动,速度剖面为抛物线:

u(y) = 6·Q/(w·h³)·y·(h − y)   (y 从下壁面起算,0 ≤ y ≤ h)
中心最大速度 u_max = 1.5·u_avg
平均速度 u_avg = Q/(w·h)

其中 Q 为体积流量。例如 w = 200 μm、h = 50 μm、Q = 5 μL/min = 8.33×10⁻¹¹ m³/s,则 u_avg = 8.33×10⁻¹¹/(200×10⁻⁶×50×10⁻⁶) = 8.33×10⁻³ m/s ≈ 8.3 mm/s,中心最大速度约 12.5 mm/s。这一速度量级决定了细胞通过观察区的时间。

二、雷诺数:判断流态的无量纲数

雷诺数定义为惯性力与黏性力之比:

Re = ρ·u·D_h/μ
D_h = 2·w·h/(w + h)   (水力直径)

以水基缓冲液为例,ρ ≈ 1000 kg/m³,μ ≈ 1×10⁻³ Pa·s。对上述通道,D_h = 2×200×50/(250) = 80 μm,Re = 1000×8.33×10⁻³×80×10⁻⁶/1×10⁻³ ≈ 0.67。可见 Re 远小于 1,惯性效应可忽略,流动由黏性主导。这意味着:细胞不会因惯性而“冲出”流线,其运动基本跟随流线;混合只能靠扩散,若需快速混合必须设计被动结构或主动扰动。

三、剪切力:细胞承受的机械刺激

壁面剪切率 γ̇ 与剪切应力 τ 是评估细胞损伤的关键参数。对矩形通道,壁面剪切率近似为:

γ̇_wall ≈ 6·Q/(w·h²)
τ_wall = μ·γ̇_wall

代入示例:γ̇_wall = 6×8.33×10⁻¹¹/(200×10⁻⁶×(50×10⁻⁶)²) = 6×8.33×10⁻¹¹/(5×10⁻¹³) ≈ 1000 s⁻¹,τ_wall = 1×10⁻³×1000 = 1 Pa。对于多数哺乳动物细胞,1 Pa 的短时暴露通常可耐受,但若停留时间延长或通道更窄,剪切力会迅速上升。例如将高度减至 20 μm、流量不变,γ̇_wall 将增至约 6250 s⁻¹,τ_wall ≈ 6.25 Pa,可能引起细胞膜损伤或活化。因此设计捕获阵列时,需在捕获效率与剪切损伤之间取平衡。

四、停留时间:细胞在观察/捕获区的暴露时长

停留时间 t_res 指细胞从进入某一功能段到离开该段的时间。对长度为 L 的直通道:

t_res = L/u_avg

若 L = 10 mm,u_avg = 8.3 mm/s,则 t_res ≈ 1.2 s。若用于荧光检测,需保证 t_res 大于相机曝光与信号累积所需时间;若用于捕获,则需保证细胞在捕获结构附近有足够时间被识别并停驻。停留时间与剪切力共同决定细胞受到的累积机械负荷,可用“剪切剂量”τ·t_res 粗略衡量。示例中 τ·t_res ≈ 1 Pa×1.2 s = 1.2 Pa·s,可作为工艺窗口的参考指标。

五、细胞受力分析:流场中细胞受到哪些力

在低雷诺数层流中,细胞主要受到以下几类力:

  • 黏性拖曳力 F_d:与相对速度、细胞尺寸和流体黏度有关,驱动细胞随流线运动;
  • 剪切梯度升力 F_l:由速度梯度引起,使细胞向速度较低或较高区域迁移,取决于细胞与流体的密度差和变形性;
  • 壁面排斥力 F_w:靠近壁面时流体动力学排斥,使细胞倾向于离开壁面一定距离;
  • 惯性升力 F_i:在有限雷诺数下出现,可将细胞推向特定平衡位置,是惯性分选的基础;
  • 外加力:如介电泳、磁力、声辐射力等,用于主动分选。

在 Re ≪ 1 时,惯性升力可忽略,细胞主要受拖曳力和壁面排斥力控制,运动轨迹基本沿流线。若需要利用惯性聚焦,则需将 Re 提高到约 1~100,并设计足够长的直通道使细胞迁移到平衡位置。例如在 50 μm 高、Re ≈ 5 的通道中,细胞可能在数毫米内聚焦到通道中心附近。

设计提示:先用 Re 判断流态,再用 τ_wall 和 t_res 估算细胞机械负荷,最后根据目标(捕获、分选、检测)决定通道尺寸、流量和功能段长度。三者相互耦合,需迭代调整。
微通道内细胞流动与受力分析示意图 矩形微通道(w=200 μm, h=50 μm) 抛物线速度剖面 u(y) u_max ≈ 1.5 u_avg 细胞 壁面剪切应力 τ_wall = μ·γ̇_wall ≈ 1 Pa(示例) 细胞受力分析 F_d 黏性拖曳力(随流) F_l 剪切梯度升力(横向迁移) F_w 壁面排斥力(离开壁面) F_i 惯性升力(Re>1时显著) 关键参数关系 Re = ρ·u·D_h / μ D_h = 2wh/(w+h) γ̇_wall ≈ 6Q/(w·h²) τ_wall = μ·γ̇_wall t_res = L / u_avg 剪切剂量 ≈ τ_wall · t_res 示例:Re≈0.67,τ≈1 Pa,t_res≈1.2 s

六、综合示例:从参数到判断

假设设计一个用于细胞捕获的直通道段,目标是在不损伤细胞的前提下让细胞有足够时间被识别。已知细胞直径约 15 μm,通道 w = 200 μm、h = 50 μm,流量 Q = 5 μL/min。计算得 u_avg ≈ 8.3 mm/s,Re ≈ 0.67,τ_wall ≈ 1 Pa。若捕获区长度 L = 5 mm,则 t_res ≈ 0.6 s。此时剪切剂量约 0.6 Pa·s,对多数细胞可接受。若将流量提高到 20 μL/min,u_avg 增至约 33 mm/s,τ_wall 增至约 4 Pa,t_res 降至 0.15 s,剪切剂量约 0.6 Pa·s 但瞬时剪切力升高,可能触发细胞活化。因此,实际设计中应优先控制壁面剪切应力,再通过延长通道或降低流量来保证停留时间。

对于分选应用,若采用惯性聚焦,则需将 Re 提高到 1~10 范围,并保证通道长度足够使细胞迁移到平衡位置。此时需重新计算 τ_wall 和 t_res,确保细胞在聚焦过程中不被过度剪切。通常的做法是:先根据目标分选精度确定 Re 和通道长度,再校核剪切力是否在细胞耐受范围内,最后通过实验验证细胞活性和回收率。

注意:上述数值均为基于理想矩形通道和牛顿流体的估算示例,实际芯片中入口效应、通道变形、细胞浓度和细胞变形性都会影响结果。设计时应以实验标定为准,不可直接将估算值作为最终工艺参数。