2. 流体力学基础回顾
精密液路系统的控制精度，最终都要落到“流量是否稳定、压力是否可预测、响应是否可重复”这三个问题上。而这三个问题，本质上由流体力学的基本方程和流态决定。本节不追求完整推导，而是把连续性方程、伯努利方程和纳维-斯托克斯方程还原成精密液路设计中真正会用到的那部分，并给出雷诺数计算与流态判别的可执行方法。
一、连续性方程：流量守恒是液路设计的起点
连续性方程的本质是质量守恒。对不可压缩液体（精密液路中常见的水、乙醇、缓冲液、试剂等，在常温常压下密度变化通常小于0.1%），可简化为体积流量守恒：
Q = A₁·v₁ = A₂·v₂ = 常数其中 Q 为体积流量（mL/min 或 μL/s），A 为流道截面积（mm²），v 为平均流速（mm/s）。这个式子看起来简单，但在精密液路中直接决定了两件事：
- 变径处的流速突变：若管路从内径 1.0 mm 缩到 0.25 mm，截面积缩小 16 倍，流速放大 16 倍。流速突变会带来局部压降和剪切变化，对剪切敏感的生物样本（如细胞悬液、蛋白溶液）可能造成损伤。
- 流量分配的可预测性：在多支路并联的液路中，若各支路截面积和长度不同，仅靠连续性方程无法确定分配比例，必须结合阻力关系（泊肃叶定律）才能求解。
教学示例：某精密加样模块主管内径 1.0 mm，流量设定 500 μL/min，则主管平均流速约为 10.6 mm/s；经三通分流后进入内径 0.25 mm 的加样针，针内平均流速约为 170 mm/s。这个数量级差异说明：同一流量下，针尖处的剪切和压降远高于主管，设计时必须单独校核。
二、伯努利方程：压力、速度与高度的能量账
伯努利方程是沿流线的机械能守恒。对不可压缩、定常、无黏流动：
P + ½ρv² + ρgh = 常数精密液路中，高度差 h 通常很小（桌面级设备），ρgh 项常可忽略；真正重要的是 P 与 ½ρv² 的相互转换。这解释了一个常见现象：管路缩径处流速升高，局部静压下降。若该处静压降到液体饱和蒸气压以下，就会发生气蚀，产生气泡，直接破坏流量稳定性。
实际液体有黏性，伯努利方程需加入损失项 h_f（沿程损失）和 h_j（局部损失）：
P₁/ρg + v₁²/2g + z₁ = P₂/ρg + v₂²/2g + z₂ + h_f + h_j在精密液路中，h_f 和 h_j 往往不是“小修正”，而是主导项。例如内径 0.25 mm、长度 100 mm 的细管，在 500 μL/min 下的沿程压降可达数十 kPa，远大于动压项。因此，精密液路设计不能只看伯努利理想式，必须把黏性损失算进去。
三、纳维-斯托克斯方程：什么时候可以简化
纳维-斯托克斯方程（N-S 方程）是黏性流体运动的完整动量方程。对不可压缩牛顿流体：
ρ(∂v/∂t + v·∇v) = -∇P + μ∇²v + ρg精密液路中，我们几乎从不直接求解完整 N-S 方程，而是根据流态做简化：
- 层流充分发展段：惯性项 v·∇v 可忽略，N-S 方程退化为泊肃叶流动，流量与压降成线性关系 Q = πΔP d⁴ / (128μL)。这是精密液路最常用的设计公式。
- 低雷诺数、非定常：若流量按正弦或阶跃变化，需保留 ∂v/∂t 项，形成非定常泊肃叶流动，用于分析响应滞后。
- 局部突变处：弯头、三通、缩径、阀口附近，惯性项不可忽略，通常用局部阻力系数 ζ 来工程化处理，而不是解 N-S。
四、雷诺数与流态判别
雷诺数表征惯性力与黏性力之比：
Re = ρvd / μ其中 ρ 为密度（kg/m³），v 为平均流速（m/s），d 为特征内径（m），μ 为动力黏度（Pa·s）。圆管流动的工程判别：
| 雷诺数范围 | 流态 | 精密液路中的典型表现 |
|---|---|---|
| Re < 2000 | 层流 | 流量-压降线性，响应可预测，控制精度高 |
| 2000 ≤ Re < 4000 | 过渡区 | 流量脉动、噪声增大，控制重复性下降 |
| Re ≥ 4000 | 湍流 | 压降非线性、脉动强，精密计量基本不可用 |
教学示例计算：25℃水，ρ=997 kg/m³，μ=0.89×10⁻³ Pa·s。内径 0.5 mm 管路，流量 1 mL/min。先算流速：A=π(0.25×10⁻³)²≈1.963×10⁻⁷ m²，Q=1 mL/min=1.667×10⁻⁸ m³/s，v=Q/A≈0.085 m/s。则 Re=997×0.085×0.5×10⁻³ / 0.89×10⁻³ ≈ 47.6，远小于 2000，属于典型层流。若把内径换成 2 mm、流量提到 50 mL/min，v≈0.265 m/s，Re≈594，仍为层流；但若内径 4 mm、流量 500 mL/min，v≈0.663 m/s，Re≈2970，进入过渡区，流量控制开始出现明显脉动。
五、流态对控制精度的影响
层流下，压降与流量成严格线性关系，注射泵或压力控制器给出的指令与实际流量之间可用固定增益描述，重复性好。进入过渡区后，同一压差下流量可能出现 5%～15% 的波动，且波动随温度、气泡、管路振动放大。湍流下，压降近似与流量平方成正比，系统增益随工作点变化，闭环控制容易振荡或稳态误差增大。
因此，精密液路设计的一条硬约束是：把工作雷诺数控制在 2000 以下，并留出至少 2 倍余量。若工艺要求高流量，应优先增大内径而不是提高流速；若必须用小内径，则需降低流量或提高黏度（如更换载液），并校核入口发展段和局部损失。
六、工程落地要点
把上述方程和判别方法落到实际设计，建议按以下顺序执行：
- 先定流量与内径：根据加样体积和节拍确定 Q，再根据 Re < 2000 反推允许的最大流速，得到最小内径。
- 校核入口发展段：计算 L_e ≈ 0.05·Re·d，若实际管长不足，压降按发展段修正，不能直接用泊肃叶公式。
- 累加局部损失：每个弯头、三通、阀口、缩径按 ζ·½ρv² 累加，尤其注意小内径处的 ζ 放大效应。
- 留出压力裕量：泵或压力源的最大输出压力应大于总压降的 1.5～2 倍，避免工作在饱和区。
- 验证流态：在最大流量工况下复算 Re，确认仍低于 2000；若接近过渡区，应通过增大内径或降低流量重新设计。
总结来说，连续性方程管“流量去哪了”，伯努利方程管“压力怎么变”，N-S 方程管“什么时候可以简化成泊肃叶”。而雷诺数是决定简化是否成立、控制是否可靠的开关。把 Re 算准、把流态控住，精密液路的控制精度才有物理基础。