4. 反应器内流体力学基础

搅拌式生物反应器的性能，最终都要落到流体在罐内的运动状态上。细胞或菌体生活在液体里，它们感受到的氧传递、营养分布、温度均匀性、剪切损伤，全部由流场决定。工程上描述这个流场，最常用的五个量是：剪切速率、雷诺数、功率数、泵送能力和混合时间。它们不是彼此独立的公式，而是一条从“流体怎么动”到“设备要多大、转多快、耗多少功”的推理链。

剪切速率：细胞承受的力学环境

剪切速率 γ̇ 表示流体内部速度梯度，单位是 s⁻¹。在搅拌罐中，最高剪切通常出现在桨叶尖端附近，量级可用叶端速度与桨叶特征尺度估算：γ̇ ≈ πND / δ，其中 N 为转速（r/s），D 为桨径，δ 为剪切层厚度。对剪切敏感细胞（如动物细胞、丝状真菌），工程上更常直接控制叶端速度 u_tip = πND，一般动物细胞培养控制在 0.5～1.5 m/s，植物细胞可放宽到 2～3 m/s，而细菌和酵母通常可承受 3～5 m/s 以上。

需要区分“平均剪切”和“局部峰值剪切”。罐内平均剪切速率可由功率耗散率 ε 估算：γ̇_avg ≈ (ε/ν)^0.5，ν 为运动黏度。但细胞损伤往往由桨叶尖端、挡板边缘、气泡破裂处的瞬时高剪切造成，因此设计时既要看平均 ε，也要看最大 ε。经验上，桨叶区最大局部耗散率可达平均值的 10～100 倍，这是放大时最容易出问题的地方。

雷诺数：判断流态的判据

搅拌雷诺数定义为 Re = ρND²/μ，其中 ρ 为密度，μ 为动力黏度。它衡量惯性力与黏性力之比。对标准六直叶圆盘涡轮，Re < 10 为层流，10～10⁴ 为过渡流，Re > 10⁴ 为充分湍流。注意这里的 N 用 r/s，D 用 m，μ 用 Pa·s。

雷诺数的工程意义在于：它决定了功率数是否进入“恒定区”。在湍流区，功率数 Np 基本不随 Re 变化，此时放大可以按等 P/V 或等叶端速度进行；在层流区，Np 与 Re 成反比，Np·Re 近似为常数，放大逻辑完全不同。因此，计算 Re 是判断“能不能用湍流放大规则”的第一步。

功率数：从转速到功率消耗

功率数定义为 Np = P/(ρN³D⁵)，其中 P 为搅拌轴功率（W）。对几何相似的搅拌系统，Np 是 Re 和弗劳德数 Fr = N²D/g 的函数。在充分湍流且无打旋（有挡板）条件下，Np 近似为常数：六直叶圆盘涡轮 Np ≈ 5.5～6.5，四叶斜桨 Np ≈ 1.3～1.5，三叶螺旋桨 Np ≈ 0.3～0.5。

由此可反推功率：P = Np·ρ·N³·D⁵。这个五次方关系非常关键——桨径增加 10%，功率增加约 61%；转速增加 10%，功率增加约 33%。实际计算时还要除以电机效率 η（通常 0.7～0.9）和机械传动效率，得到装机功率。若通气，还需考虑通气功率下降，常用 Michel-Miller 关联式修正。

泵送能力：流体被“推”出去多少

泵送能力 Q 表示桨叶单位时间排出的液体体积，Q = N_q·N·D³，N_q 为泵送数。六直叶圆盘涡轮 N_q ≈ 0.7～0.9，四叶斜桨 N_q ≈ 0.7～0.8，三叶螺旋桨 N_q ≈ 0.4～0.6。泵送能力决定了大尺度循环时间，也就是流体在罐内走一圈需要多久：t_c = V/Q，V 为装液体积。

泵送能力与剪切是一对矛盾：大直径、低转速桨（如斜叶桨）泵送量大、剪切低，适合剪切敏感体系；小直径、高转速桨（如圆盘涡轮）剪切强、分散气泡好，适合高耗氧发酵。选择桨型，本质上是在“混合效率”和“剪切损伤”之间找平衡。

混合时间：均匀性到底要多久

混合时间 t_m 通常定义为从加入示踪剂到达到最终均匀度 95% 所需的时间。工程上常用经验关联式估算：t_m = 5.9·(T/D)²·(1/N)·(1/N_q)^(1/3)·(ν/ND²)^(1/6)，其中 T 为罐径。更粗略地，湍流条件下 t_m ≈ 5～10 倍循环时间。

混合时间直接影响补料策略、pH 控制、温度均匀性和氧分布。对大型罐（>50 m³），混合时间可达 30～60 s，这意味着局部 pH 或底物浓度可能长时间偏离设定值。因此，放大时不能只看 P/V，还要校核 t_m 是否满足工艺要求。

五个参数如何协同决定设计

实际设计一台反应器，通常按以下顺序推进。先根据细胞类型确定剪切上限，选定叶端速度范围；再根据氧需求确定 P/V 目标（例如大肠杆菌高密度发酵常取 1～5 kW/m³，动物细胞取 0.05～0.2 kW/m³）；然后选择桨型，由 Np 和 N_q 反算转速与桨径；接着用 Re 确认是否处于湍流区，若不在则需调整黏度假设或改用层流放大规则；最后用 t_m 校核混合是否满足补料和 pH 控制要求。

举一个教学示例：装液 2 m³ 的罐，T = 1.2 m，D = 0.4 m，采用六直叶圆盘涡轮，Np = 6.0，N_q = 0.8，ρ = 1000 kg/m³，μ = 0.001 Pa·s。若目标 P/V = 1 kW/m³，则 P = 2000 W。由 P = Np·ρ·N³·D⁵ 得 N³ = 2000/(6.0×1000×0.4⁵) = 2000/(6.0×1000×0.01024) ≈ 32.6，N ≈ 3.2 r/s（约 192 rpm）。此时 Re = 1000×3.2×0.16/0.001 = 5.12×10⁵，充分湍流。Q = 0.8×3.2×0.064 = 0.164 m³/s，循环时间 t_c = 2/0.164 ≈ 12.2 s，t_m 约 60～120 s。这个混合时间对补料分批发酵偏长，若工艺要求 t_m < 30 s，则需要增大桨径或提高转速，但必须重新校核剪切是否超过细胞耐受上限。

工程意义与常见误区

第一个误区是只看 P/V 不看剪切。两台罐 P/V 相同，但桨径和转速不同，叶端速度和局部剪切可能差数倍，对剪切敏感细胞的结果完全不同。第二个误区是忽略黏度变化。发酵过程中菌体浓度升高、代谢产物积累，μ 可能从 1 mPa·s 升到几十甚至几百 mPa·s，Re 随之下降，流态从湍流退到过渡流，Np 不再恒定，原本的放大规则失效。第三个误区是把混合时间当成固定值。t_m 随罐径增大而显著增加，大罐的混合问题往往比小罐严重得多。

因此，在工艺开发阶段就应建立“桨型—转速—P/V—u_tip—t_m”的对应表，并在放大时逐项校核。对剪切敏感体系，优先选择大直径、低转速、高 N_q 的斜叶桨；对高耗氧体系，优先选择高 Np、高剪切分散能力的圆盘涡轮，同时用挡板消除打旋，保证 Np 和 N_q 的关联式适用。只有把五个参数放在同一张设计图上，才能避免“小罐成功、大罐失败”的典型放大问题。