微流体力学
微流体力学的理论基础
微流体力学的基础
老师，微流体力学与宏观流体力学有什么区别？只是尺度更小吗？
尺度的变化改变了物理的"支配者"。在宏观尺度，惯性力占主导，但在微尺度（特征长度 $L \sim 1\text{--}100\,\mu\text{m}$）粘性力和表面张力占主导。
具体用数值比较。对于水流速 $U = 1\,\text{mm/s}$、通道宽度 $D = 100\,\mu\text{m}$：
因为Re << 1，惯性力可以忽略，Stokes流（爬流）占主导。
Re = 0.1时完全由粘性主导。不用担心湍流。
正确。还要看毛细数（Capillary number）：
因为Ca << 1，液滴和气泡的形状由表面张力决定。
Stokes方程
在Re << 1的极限情况，Navier-Stokes方程的惯性项消失，得到Stokes方程。
Stokes方程有以下重要性质：
- 线性：解可以叠加
- 时间可逆：反转外力，流体回到原始状态（混合困难的原因）
- 瞬时响应：无惯性，压力变化立即传播到整个系统
时间可逆是说微尺度搅拌无法混合？
用常规方法无法混合。这是微混合器设计困难的原因。只能依靠扩散或利用混沌对流（如曲折流路）。
微流体特有的物理
在微尺度变得重要的物理现象总结如下。
| 现象 | 控制参数 | 宏观影响 | 微观影响 |
|---|---|---|---|
| 表面张力 | Ca, We, Bo | 通常忽略 | 占主导 |
| 电渗流 (EOF) | $\zeta$ 电位、Debye长 | 忽略 | 用作驱动力 |
| 滑流 | Knudsen数 Kn | 无滑 | Kn > 0.01 时滑流 |
| 扩散混合 | Peclet数 Pe | 对流主导 | 扩散主导 |
| 接触角滞后 | 前进/后退接触角 | 不重要 | 控制装置动作 |
Bond数 $\text{Bo} = \rho g L^2 / \sigma$ 在微尺度也变得非常小（$L = 100\,\mu\text{m}$时 $\text{Bo} \sim 10^{-6}$）。重力完全可以忽略，在太空和地球上行为相同。
重力无效这点违反直觉，但这是尺度效应。
微流体工程的诞生——1990年代的μTAS革命和芯片实验室
微流体工程的黎明期始于1990年代初，瑞士ETH的Manz & Widmer(1990)提出"μTAS（微型全分析系统）"概念。利用半导体制造技术（光刻）在玻璃基板上制作数十到数百μm宽的流路，在芯片上完成试剂混合、分离和检测，这就是"芯片实验室"的概念。在这样的流路尺度，Reynolds数小于1，Stokes流占主导，无法使用湍流混合，因此需要T型或螺旋形等被动混合结构。现在COVID-19快速诊断试剂盒和DNA测序仪中的微通道采用了这项技术，CFD在通道设计优化中成为产品开发的核心。
微流体力学的数值计算手法
微流体的CFD手法
微流体的仿真可以直接用普通CFD求解器吗？
在连续体近似成立的范围内（Kn < 0.01）可以用标准的Navier-Stokes型CFD。但微尺度有特殊的数值问题。
二相流（液滴、气泡）的界面追踪
微流体装置中常见液滴生成和T型合流。对比界面追踪的主要方法。
| 方法 | 优点 | 缺点 | 适用例 |
|---|---|---|---|
| VOF（流体体积） | 质量守恒好，成本低 | 界面扩散，接触线问题 | 液滴生成、栓塞流 |
| Level-Set | 界面形状光滑 | 质量不守恒 | 合流、分裂 |
| CLSVOF | 统合VOF+LS优点 | 实现复杂 | 需要高精度的问题 |
| Phase-Field | 自然描述接触线 | 网格要求严格 | 润湿现象、接触角控制 |
| Front-Tracking | 界面显式追踪 | 拓扑变化弱 | 单液滴变形 |
表面张力的数值计算通常用CSF（连续表面力）模型，但微尺度会产生寄生电流（spurious current）问题。这是表面张力离散化误差导致的非物理流动，Ca越小越明显。
表面张力大但网格粗会产生非物理流动，对吧。
对策包括：
- 界面处网格分辨率至少10个单元
- 用Height Function法计算曲率（比CSF精度提高很多）
- Sharp Surface Force法
- 使用Phase-Field模型（热力学描述界面，降低寄生电流）
电渗流（EOF）的计算
微通道常用电场驱动流体，称为电渗流。控制方程为：
$\rho_e$ 是电荷密度，$\varepsilon$ 是介电常数，$\mathbf{E} = -\nabla\phi$ 是电场。Debye长 $\lambda_D \sim 1\text{--}100\,\text{nm}$ 的区域需要很细的网格。
Debye长是nm级，与流路尺度相差3个数量级。
实际上常用Helmholtz-Smoluchowski的滑移速度条件：
作为壁面边界条件，不需要解析EDL（电双层）内部。$\zeta$ 是Zeta电位，$E_t$ 是壁面切向电场。
微流路的"毛细管力"——重力之外的支配力
通道宽度低于100μm时，重力不再是流动的主导力。反而表面张力成为主角，毛细管压力 ΔP=4γcosθ/D（D为通道直径）将液体吸入。实际上，疾病诊断用纸基微流体中，血液可以在不使用泵的情况下通过毛细管力流动。数值处理这种情况需要结合VOF法或格子Boltzmann法来追踪气液界面，手法选择对精度影响很大。
微流体力学的实务应用
微流体装置的CFD解析流程
微流体装置的仿真具体怎么进行？
以典型的T型液滴生成器为例说明。
- 通道宽度100 μm、深度50 μm是典型值
- CAD设计要精确。包括倒角和制造误差（PDMS设备±5 μm）
- 界面所在区域利用AMR（适应性网格细分化）
- 通道宽度方向最少20个单元，界面近处再细2-3倍
- 接触线附近需特别高分辨率
- 连续相（如油）：$\mu_c = 0.01\,\text{Pa}\cdot\text{s}$, $\rho_c = 800\,\text{kg/m}^3$
- 分散相（如水）：$\mu_d = 0.001\,\text{Pa}\cdot\text{s}$, $\rho_d = 1000\,\text{kg/m}^3$
- 界面张力：$\sigma = 0.005\,\text{N/m}$（添加表活剂时）
第4步：边界条件
- 入口：流量指定（注射泵对应）。典型 $Q = 1\text{--}10\,\mu\text{L/min}$
- 出口：压力指定（大气压）
- 壁面：无滑移 + 接触角指定（PDMS水接触角 ≈ 110°）
接触角是宏观CFD中没有的因素。
接触角对液滴大小和生成频率影响很大。用动态接触角模型（如Cox-Voinov法则）会更准确：
$\theta_D$ 是动态接触角，$\theta_S$ 是静态接触角，$L_H$ 是宏观尺度，$L_m$ 是分子尺度。
解析结果的评价指标
微流体装置性能的主要评价指标：
| 指标 | 定义 | 用途 |
|---|---|---|
| 液滴直径 | $d/w$（用通道宽度无次元化） | 液滴生成装置设计 |
| 生成频率 | $f \cdot w / U$（无次元） | 吞吐量评估 |
| 混合度 | $1 - \sigma_c / \sigma_{c,0}$（浓度标准差） | 微混合器性能 |
| 压力损失 | $\Delta p / (\mu U L / w^2)$（无次元） | 泵需求 |
| 分离效率 | 目标粒子的回收率 | 分离装置 |
液滴直径用通道宽度无次元化是微流体的约定。
正是。T型液滴生成器中，$d/w$ 是Ca数的函数。在Squeezing范围（$\text{Ca} < 0.01$）中，$d/w \propto Q_d / Q_c$（依赖流量比），在Dripping范围（$\text{Ca} > 0.1$）中，$d/w \propto \text{Ca}^{-1/3}$，标度关系会改变。
微流体装置的"失败常见问题"——流路堵塞的真实原因
微流体实验中报告的"流动停止"多数是气泡堵塞。1μm宽的通道中1μm气泡会导致毛细管压力超过泵压，流动完全停止。解决方法是进行通道内壁亲水化处理（接触角θ低于5°）或在出口设置气泡陷阱。CFD要事先预测这种情况，需要包含动态接触角模型的多相流解析。实验遇到困难前，先检查表面润湿性的数值参数，这是实践铁则。
微流体力学的软件比较
微流体对应工具的比较
用什么软件做微流体仿真？
微流体涉及多种物理，工具选择很重要。按用途整理。
COMSOL Multiphysics（Microfluidics模块）
微流体领域最广泛使用的工具之一。
- 优势：EOF（电渗流）、DEP（介电泳）、AC EOF、扩散反应耦合容易
- 界面追踪：Level-Set、Phase-Field、Moving Mesh三种方法
- 接触角：静态和动态接触角设置直观
- EDL解析：完整Poisson-Nernst-Planck方程可解
- 弱点：大规模3D计算容易内存不足
Ansys Fluent
Fluent通用性强，也能用于微流体。
- VOF模型：液滴生成仿真有大量应用实绩
- UDF：电场驱动、接触角可自定义
- AMR：界面附近自动网格细分
- 弱点：没有专用EOF模块，需UDF实现
OpenFOAM
开源的自由度在这个领域发挥作用。
- interFoam：基于VOF的二相流求解器（液滴生成有大量论文应用）
- interPhaseChangeFoam：包含相变（蒸发、凝聚）的微流体
- 定制求解器：EOF、DLP、声场等多物理可用C++自由实现
- 社区：GitHub上有多个微流体定制求解器开源代码
专用工具
| 工具 | 开发商 | 特点 |
|---|---|---|
| FLOW-3D | Flow Science | VOF法精度高。表面张力处理有口碑 |
| Gerris/Basilisk | S. Popinet | AMR+VOF。学术界评价非常高。免费 |
| COMSOL | COMSOL AB | 多物理耦合容易 |
| Palabos | FlowKit | 格子Boltzmann法。复杂形状微流体适合 |
Basilisk没听说过。
Basilisk是Stephane Popinet开发的自适应直交网格（Quadtree/Octree）求解器，表面张力问题精度世界一流。Height Function法曲率计算和AMR结合，寄生电流极小。学术论文中频繁引用。
选定指南
按应用推荐工具：
| 应用 | 推荐工具 | 理由 |
|---|---|---|
| EOF/DEP | COMSOL | 电化学耦合容易 |
| 液滴生成 | Basilisk、Fluent、interFoam | 表面张力精度 |
| 扩散混合 | COMSOL、Fluent | 多成分输运 |
| Lab-on-a-chip全面 | COMSOL | GUI易操作，学习成本低 |
| 学术研究（高精度） | Basilisk | AMR+高精度表面张力 |
| 大规模3D计算 | OpenFOAM、Fluent | 并行计算可扩展性 |
COMSOL适合入门，精度追求用Basilisk或OpenFOAM。
微流体商用工具的"陷阱"
微流体芯片设计中直接用汎用CFD工具会遭遇问题。通道宽度低于10μm时，Knudsen数超过0.001，连续体假设可疑。ANSYS和COMSOL默认解Navier-Stokes连续体，不设定滑流修正，细管中流量会比实验多出2-3倍。专用微流体工具或薄膜流选项的有无是汎用和专用工具的大差异。
微流体力学的先端研究
格子玻尔兹曼法（LBM）的微流体计算
听说格子Boltzmann法对微流体最适合？
LBM与微流体相性非常好。有3个原因。
1. 多相流的自然描述：Shan-Chen模型或free-energy模型自动追踪界面。无需VOF那样的界面重构
2. 滑流扩展：Knudsen数大的区域也能自然导入滑移边界条件
3. 复杂形状适配：反弹边界条件简单表现任意形状的壁面
基本BGK（Bhatnagar-Gross-Krook）模型的格子Boltzmann方程：
$f_i$ 是粒子分布函数，$\tau$ 是松弛时间（粘度相关：$\nu = c_s^2(\tau - 0.5)\Delta t$）。
粒子方法与微观物理的亲和力高。
Digital Microfluidics（EWOD）
利用Electrowetting-on-Dielectric（EWOD）的数字微流体是近年快速发展的领域。
通过加电压改变局部接触角，自如操纵液滴。Lippmann-Young方程：
$\theta_V$ 是电压 $V$ 下的接触角，$\theta_0$ 是初始接触角，$d$ 是绝缘膜厚度。
CFD建模主要用Phase-Field模型耦合电场。COMSOL的AC/DC+Phase-Field模块可计算。
用电压驱动液滴，可用于PCR或快速诊断。
Organ-on-a-Chip（器官芯片）
微流体最瞩目的应用是Organ-on-a-Chip。
在微通道中培养血管内皮细胞，重现器官微环境。CFD的作用：
- 剪切应力控制：血管内皮需要 $\tau_w = 1\text{--}10\,\text{dyne/cm}^2$ 的生理剪切
- 氧气和营养运输：评价 $\text{Pe} = UL/D_m$。$D_m \sim 10^{-9}\,\text{m}^2/\text{s}$（氧在水中的扩散）
- 药物浓度梯度设计：用Pe数判断扩散还是对流主导
微流体重现生体物理。
研究动向
微流体力学最前沿：
- 纳流体力学：通道直径低于10 nm时连续体近似破裂。需要MD仿真或DPD（耗散粒子动力学）
- 声波流（Acoustofluidics）：超声定在波进行粒子和细胞的大小分选。BAW/SAW装置CFD建模
- 喷墨高精度化：控制皮升级液滴，抑制卫星液滴。CFD优化
- 3D打印（双光子聚合）：微流体装置直接成形与仿真联动
- 惯性微流体：Re $\sim$ 10-100的中间范围，利用惯性升力的粒子分离
微流体是医疗、生物与CFD的交叉前沿。
器官类器官芯片与微流体——"体外重现体内"的挑战
微流体的前沿之一是"Organ-on-a-chip"（器官芯片）。肝、肠、肺细胞在10-100μm宽的微通道中培养，在接近生体的流动环境（Re≈0.01-0.1）中试验药物反应。与普通细胞培养不同，能精密控制剪切应力、扩散和温度分布，在某些情况下动物实验的精度被超越。流体分析上这是"超低Re、复杂生体物质"的极难问题，格子Boltzmann法和分子动力学混合模型备受关注。
微流体力学的故障排除
微流体计算的典型问题
微流体CFD特别要注意什么问题？
与宏观CFD不同的问题很多。按类型分类。
1. 寄生电流（Spurious Current）导致液滴崩裂
症状：静止液滴发生变形分裂。界面附近出现非物理涡。
原因：表面张力的CSF模型曲率计算的离散化误差。Ca << 1时明显。
定量判据：寄生电流速度超过物理流速10%就不可信。
对策：
- 增加界面处网格分辨率（最少10个单元/液滴直径）
- 切换到Height Function法计算曲率
- 使用Phase-Field模型（Cahn-Hilliard方程）
- OpenFOAM的
interFoam调整interfaceCompression
静止液滴作为测试检查寄生电流大小再进行实际计算。
2. 接触线被钉住，液滴无法移动
症状：壁面液滴虽有驱动力但不移动
原因：无滑移条件与界面追踪矛盾。接触线处速度为零但界面移动的悖论（Huh-Scriven悖论）引起数值问题。
对策：
- Navier滑移条件：壁面引入微小滑移长度 $\beta$。$u_{\text{slip}} = \beta \partial u/\partial n$
- 物理上 $\beta = 1\text{--}10\,\text{nm}$，但网格无法解析时用数值上较大的值
- Phase-Field模型能自然回避接触线特异性
3. 时间步长极端微小
症状：Courant数限制导致 $\Delta t \sim 10^{-10}\,\text{s}$，计算停滞
原因：表面张力毛细波传播速度限制。Brackbill极限：
以 $\Delta x = 1\,\mu\text{m}$、水的物性，$\Delta t < 1.2 \times 10^{-8}\,\text{s}$。
对策：
- 网格不要过度细分
- 隐式表面张力处理（部分求解器支持）
- 用AMR只在界面处细化，远处粗化
表面张力限制时间步是微尺度特有的制约。
4. 二相流质量守恒差
症状：计算过程中液滴体积减少或增加
原因：Level-Set法追踪界面为 $\phi = 0$ 等值面，但重初始化步骤中质量丧失。
对策：
- 使用VOF法（质量守恒好）
- CLSVOF（Level-Set和VOF结合）
- 加入质量修正步骤
- Phase-Field模型（Cahn-Hilliard：质量守恒型）
调试步骤
微流体计算调试按下列顺序。
1. 单相流测试：先用Poiseuille流验证圧力損失 $\Delta p = 12\mu Q L/(w h^3)$（矩形断面）是否匹配
2. 静止液滴测试：球形液滴静置，Laplace压力 $\Delta p = 2\sigma/R$ 正确重现？寄生电流足够小？
3. 接触角测试：壁面液滴平衡接触角与设置值一致？
4. 动态测试：T型合流液滴生成与实验定量比较
分阶段逐一验证也是微流体基本。
特别是Laplace压力测试最低限度要做。这个测不过，二相流结果完全不可信。
微流体解析"数值发散"——元凶九成是网格
微流体解析不收敛，首先怀疑网格长宽比。通道截面10μm×100μm时，如用1:10的长宽比网格，压力梯度数值误差易膨胀。理想是各向最少5个单元，长宽比3以下。另一个麻烦源是入口边界条件的"急剧速度变化"。不假定充分发展流而设置均匀速度时，入口附近会振荡且不收敛。在入口前设置发展长度 L≈0.06·Re·D（微流体通常数十μm）即改善。
相关专题
更多细节
错误