2. 灌流培养的流体力学基础
灌流培养的核心,是让培养液在微通道或培养腔中持续流动,在带走代谢废物、补充营养的同时,对细胞施加可控的力学与传质环境。要理解并设计这样的系统,必须先掌握几个基本概念:黏性流、层流、雷诺数、剪切力与压力驱动流。它们共同决定了流速分布、细胞受到的机械刺激以及物质输运的效率。
黏性流与层流:微尺度下的主导行为
在微流控芯片中,通道特征尺寸通常在几十到几百微米,流速多为毫米每秒量级。此时流体的惯性力远小于黏性力,流动由黏性主导,表现为稳定的层流。层流的特点是流线平行、互不掺混,同一截面上的速度呈抛物线分布:通道中心速度最大,靠近壁面速度趋近于零。这一速度梯度正是剪切力的来源。
判断流动状态的关键参数是雷诺数 Re:
Re = ρ · v · D_h / μ
其中 ρ 为流体密度,v 为平均流速,D_h 为水力直径,μ 为动力黏度。对于矩形通道,水力直径 D_h = 2wh/(w+h),w、h 分别为通道宽和高。以水基培养液为例,ρ≈1000 kg/m³,μ≈1×10⁻³ Pa·s。若通道宽 400 μm、高 100 μm,则 D_h≈160 μm;当平均流速为 1 mm/s 时,Re≈0.16,远小于 2000 的临界值,属于典型层流。这意味着在灌流培养中,几乎不会出现湍流,流动可预测、可解析。
剪切力:细胞感受到的机械信号
细胞贴附在通道底面时,感受到的是壁面剪切应力 τ_w。对于充分发展的层流,可用下式估算:
τ_w = 6 μ Q / (w h²)
其中 Q 为体积流量。该式表明,剪切力与流量成正比,与通道高度平方成反比。高度越小,同样的流量会产生更大的剪切。例如,w=400 μm、h=100 μm、Q=24 μL/min(即 4×10⁻¹⁰ m³/s)时,τ_w≈0.036 Pa,即 0.36 dyn/cm²。这一量级接近许多内皮细胞、肝细胞培养中常用的生理剪切范围。若把高度减半到 50 μm,剪切力将增至约 1.44 dyn/cm²。因此,设计灌流参数时,通道几何与流量必须联合考虑。
压力驱动流与沿程压降
微流控灌流通常由注射泵、压力控制器或重力驱动。压力驱动流中,流量与压降的关系由哈根–泊肃叶方程描述。对矩形通道,近似有:
ΔP = 12 μ L Q / (w h³) · f
其中 L 为通道长度,f 为形状修正因子,宽高比越大 f 越接近 1。仍以上述通道为例,L=3 cm、Q=24 μL/min,可估算 ΔP 约为数百帕量级。这个压降看似不大,但在多通道并联、长管路或高黏度培养液中会显著累积,进而影响各支路流量分配的均匀性。实际系统中,应尽量缩短连接管路、增大主管截面,或采用独立压力源分别控制各培养腔。
传质:流速与扩散的竞争
灌流培养不仅要提供合适的剪切,还要保证氧气和营养物质及时到达细胞。传质由对流和扩散共同完成,可用佩克莱数 Pe 判断二者相对强弱:
Pe = v · h / D
其中 D 为溶质扩散系数,小分子营养物在水中约为 10⁻⁹ m²/s 量级。若 v=1 mm/s、h=100 μm,则 Pe≈100,说明沿流动方向对流远强于扩散;但在垂直流动方向的截面上,混合仍主要靠扩散。因此,通道高度不宜过大,否则中心流与壁面附近会形成营养梯度。降低高度、提高流速可减薄扩散边界层,但会同时提高剪切力,需要在传质与力学之间取平衡。
尺度效应:为什么微通道如此不同
宏观管道中,惯性主导,流动容易发展为湍流,混合快但剪切不可控。而在微通道中,比表面积大、黏性效应被放大,层流稳定、剪切可精确计算。这带来两个直接后果:一是细胞受到的力学刺激高度可预测,便于建立定量模型;二是混合主要依赖扩散,若需要快速混合,必须借助蛇形通道、微柱阵列或主动扰动结构。灌流培养正是利用这种“可控的层流”,在细胞周围构建稳定的营养与力学微环境。
参数选择与计算示例
假设要在一个宽 500 μm、高 80 μm、长 2 cm 的矩形培养腔中培养肝细胞,目标壁面剪切力约为 0.5 dyn/cm²(即 0.05 Pa)。由 τ_w = 6μQ/(wh²) 反推流量:
Q = τ_w · w · h² / (6μ)
= 0.05 × 5×10⁻⁴ × (8×10⁻⁵)² / (6×10⁻³)
≈ 2.67×10⁻¹¹ m³/s
≈ 1.6 μL/min
此时平均流速 v = Q/(wh) ≈ 0.67 mm/s,Re 远小于 1,仍为层流。若发现营养供应不足,可适当提高流量,但需重新核算剪切力是否超出细胞耐受范围。实际调试中,建议先以低流量启动,待细胞贴附稳定后再逐步提升至目标值,并观察细胞形态与代谢指标。
掌握这些流体力学基础后,就能在灌流培养中把“流速—剪切—传质”三者联系起来,为后续的泵阀控制、气泡管理和长期稳定运行提供定量依据。