2. 液滴生成物理基础

液滴在微通道中能否稳定生成、生成多大、以什么频率生成,本质上由两相流体的惯性、黏性、界面张力三者的相对强弱决定。工程上不会凭感觉调泵速,而是先把这些力折算成几个无量纲数,再用流型图判断当前工况落在哪个区域。本节围绕四个核心量展开:界面张力 σ、毛细数 Ca、雷诺数 Re、韦伯数 We,并给出如何用它们读懂液滴生成流型图。

一、界面张力:液滴生成的“弹性皮”

界面张力 σ(单位 N/m 或 mN/m)是两相界面单位长度上的自由能,物理上表现为界面像一层有弹性的薄膜,倾向于把界面面积缩到最小。液滴之所以能保持球形而不是被连续相冲散,靠的就是 σ。常见体系参考值:水—空气约 72 mN/m,水—矿物油(含表面活性剂)约 5~20 mN/m,水—氟油约 15~30 mN/m。加入表面活性剂会显著降低 σ,同时抑制液滴聚并,这是液滴工程中“先选体系、再调流速”的原因。

界面张力在液滴生成中扮演两个角色:一是抵抗连续相剪切,使分散相在断裂前能形成清晰的颈部;二是决定断裂后液滴的回缩速度。σ 越大,液滴越难被拉断,需要更高的连续相流速或更小的通道尺寸;σ 越小,界面越易变形,液滴更容易生成,但也更容易发生聚并和润湿问题。

二、三个无量纲数:Ca、Re、We

把黏性力、惯性力、界面张力两两相比,就得到三个关键无量纲数。设连续相黏度 μ_c、密度 ρ_c、特征流速 U、特征尺寸(通道水力直径)L、界面张力 σ,则:

Ca = μ_c · U / σ        (黏性力 / 界面张力)
Re = ρ_c · U · L / μ_c   (惯性力 / 黏性力)
We = ρ_c · U² · L / σ    (惯性力 / 界面张力)

三者关系:We = Re · Ca

毛细数 Ca 是液滴生成中最重要的参数。Ca 小(约 10⁻³ 以下)时界面张力占主导,液滴在颈部被“挤断”,生成尺寸较均匀,属于挤压(squeezing)区;Ca 增大到 10⁻² 量级,黏性剪切开始主导,进入滴落(dripping)区;Ca 继续增大到 10⁻¹ 以上,连续相把分散相拉成细长射流,进入喷射(jetting)区,液滴由射流表面不稳定波破碎产生,尺寸分布变宽。

雷诺数 Re 判断流动是层流还是湍流。微通道中 L 通常为 50~500 μm,U 为 mm/s~cm/s 量级,Re 常在 10⁻³~10 之间,几乎总是层流。这意味着液滴生成不依赖湍流混合,而依赖层流界面变形与断裂,惯性效应通常可忽略,但在高速喷射区 Re 会升到 1~100,惯性对射流断裂有贡献。

韦伯数 We 衡量惯性力相对界面张力的大小。We ≪ 1 时界面张力主导,液滴接近球形;We 接近或大于 1 时惯性力使界面明显变形,液滴被拉长甚至破碎成卫星滴。在液滴生成中 We 常作为辅助判据,用来判断高速工况下是否会出现卫星液滴。

工程速记:先算 Ca 定流型,再用 Re 确认层流假设是否成立,最后用 We 检查高速下是否出现卫星滴。三者不是并列关系,而是有主次的判据链。

三、液滴生成流型图

流型图(regime map)通常以连续相毛细数 Ca_c 为横轴、分散相与连续相流量比 Q_d/Q_c(或分散相毛细数 Ca_d)为纵轴,把挤压、滴落、喷射三个区域画在图上。同一芯片、同一流体体系,改变两相流量就能在图上移动工作点。

以 T 型通道为例,通道宽 200 μm、深 100 μm,水力直径约 133 μm,连续相为含 2% Span 80 的矿物油(μ_c ≈ 0.02 Pa·s,σ ≈ 8 mN/m),分散相为去离子水。当 Q_c = 2 μL/min、Q_d = 0.5 μL/min 时,U_c ≈ 0.0017 m/s,Ca_c ≈ 4×10⁻³,落在挤压区,液滴长度约为通道宽的 2~3 倍;把 Q_c 提到 20 μL/min,Ca_c ≈ 4×10⁻²,进入滴落区,液滴长度降到约 1 倍通道宽;继续提到 100 μL/min,Ca_c ≈ 0.2,进入喷射区,可观察到细长射流和尺寸分布变宽的液滴群。以上为教学示例参数,实际值需按具体体系标定。

液滴生成流型图(Ca_c — Q_d/Q_c) 连续相毛细数 Ca_c(对数) 流量比 Q_d / Q_c 10⁻³ 10⁻² 10⁻¹ 10⁰ 0.1 1 10 100 挤压区 Squeezing 界面张力主导 液滴尺寸均匀 滴落区 Dripping 黏性剪切主导 液滴尺寸随 Q_c 减小 喷射区 Jetting 惯性/剪切拉出射流 尺寸分布变宽 A 挤压 B 滴落 C 喷射 Ca_c < 10⁻² 10⁻² ~ 10⁻¹ Ca_c > 10⁻¹

四、从参数到操作:如何用这套判据调工况

实际调机时,建议按以下顺序推进。第一步,确定流体体系并测量或查表得到 μ_c、ρ_c、σ,这是所有无量纲数的基础。第二步,根据目标液滴尺寸估算所需 Ca_c:若希望液滴尺寸接近通道宽度且分布窄,把 Ca_c 控制在 10⁻³~10⁻²;若希望液滴远小于通道宽度,把 Ca_c 提到 10⁻²~10⁻¹。第三步,由 Ca_c = μ_c·U_c/σ 反推所需连续相流速 U_c,再换算成泵流量 Q_c = U_c·A(A 为通道截面积)。第四步,用 Re 和 We 复核:若 Re > 100 或 We > 1,说明惯性效应不可忽略,需重新评估流型或降低流速。

需要特别注意的是,流型图上的分界线并非固定不变,它随通道几何(T 型、十字型、流动聚焦)、润湿性(连续相是否完全润湿通道壁)、表面活性剂浓度而变化。因此,工程上更可靠的做法是:先用上述无量纲数估算初始工作点,再在目标芯片上做一次小范围流量扫描,记录液滴长度、生成频率和尺寸变异系数 CV,把实测点标回流型图,形成该芯片专属的工作窗口。这样得到的流型图才是可直接用于生产的工艺文件。

常见误区:只看流量不看 Ca。同样的 Q_c 和 Q_d,换成低黏度连续相或高界面张力体系,Ca 会显著变化,流型可能从滴落跳到挤压,液滴尺寸随之改变。调参时务必以无量纲数为准,而不是以泵读数为准。