4. 搅拌功率计算

搅拌功率是生物反应器放大中最先被追问、也最容易被算错的一个量。它决定了单位体积输入功率、氧传递能力、剪切水平与混合时间，而这些又直接约束细胞或菌体的生长状态。本节不讨论功率与传质之间的换算，而是聚焦于“给定一台反应器和一种流体，如何把搅拌功率算出来”。核心工具是功率准数 N_p，它把叶轮的几何、流体的物性和操作转速统一到一个无量纲数中。

功率准数的基本定义与牛顿流体关联式

搅拌功率 P 的无量纲化形式为功率准数：

N_p = P / (ρ N³ D⁵)

其中 ρ 为流体密度（kg/m³），N 为叶轮转速（r/s，注意不是 rpm），D 为叶轮直径（m）。对牛顿流体，N_p 是雷诺数 Re = ρND²/μ 的函数。在湍流区（Re > 10⁴），N_p 基本恒定，称为湍流功率准数 N_p,t，只取决于叶轮类型与几何比例。这正是生物反应器最常用的工作区间：典型发酵液中 Re 常在 10⁴～10⁶ 量级。

常见叶轮的湍流功率准数（以挡板、标准几何为条件）大致为：六直叶圆盘涡轮（Rushton）约 5.0～5.5；四斜叶桨式约 1.3～1.5；三叶螺旋桨约 0.3～0.5；大直径低剪切轴向流叶轮约 0.3 左右。这些数值不是普适常数，一旦几何偏离标准配置，必须重新取值或修正。

非牛顿流体：表观黏度与 Metzner-Otto 方法

发酵液常表现为假塑性流体，服从幂律模型 τ = Kγ̇ⁿ，其中 K 为稠度系数，n 为流动行为指数（n < 1）。此时黏度不是常数，雷诺数无法直接计算。工程上采用 Metzner-Otto 思路：认为叶轮附近存在一个代表剪切率 γ̇_eff = k_sN，其中 k_s 是与叶轮几何相关的常数（Rushton 涡轮约 11～13，桨式约 10～11）。由此得到表观黏度：

μ_app = K(k_sN)ⁿ⁻¹

再用 μ_app 代替 μ 计算广义雷诺数 Re_g = ρND²/μ_app。若 Re_g 仍落在湍流区，则可沿用湍流功率准数；若落入过渡区，则需用 N_p–Re 曲线插值。注意 k_s 与 N_p 都依赖叶轮类型，不能混用。

通气条件下的功率下降

通气后，叶轮后方空腔被气体占据，液体循环被削弱，功率显著下降。常用修正形式为：

P_g/P₀ = 1 − a·(Q_g/(ND³))^b

其中 P₀ 为不通气功率，Q_g 为通气体积流量，无量纲数 Q_g/(ND³) 称为通气数 N_a。指数 a、b 随叶轮与系统而变，Rushton 涡轮常见 a ≈ 1.0～1.5、b ≈ 0.4～0.6。工程上更常用的是带挡板、有通气时的经验关联式，或直接查 P_g/P₀ 与 N_a、Fr 的关系图。需要强调：通气功率下降幅度可达 30%～60%，放大时若忽略这一点，会高估实际输入功率与 k_La。

叶轮类型与几何配置的影响

叶轮类型决定 N_p,t 的量级与流型：径向流叶轮（Rushton）功率准数高、剪切强、适合分散气体；轴向流叶轮（螺旋桨、斜叶桨）功率准数低、循环量大、剪切温和，适合剪切敏感细胞。几何配置同样关键：D/T（叶轮/罐径比）增大时，同样转速下功率按 D⁵ 增长，但单位体积功率与混合效率的关系并非线性；叶轮离底距离 C、挡板数量与宽度、多层叶轮的间距都会改变实际 N_p。因此查表得到的 N_p,t 只适用于其标定几何，偏离时应以实测或 CFD 校核。

手算案例：0.1 m³ 罐中 Rushton 涡轮的功率

设一台实验室反应器，罐径 T = 0.60 m，装液体积 V_L = 0.10 m³，配单个六直叶圆盘涡轮，D = 0.20 m，转速 N = 300 rpm = 5 r/s，通气量 1.0 vvm。先按牛顿流体、μ = 0.002 Pa·s、ρ = 1000 kg/m³ 计算。

第一步算雷诺数：Re = 1000 × 5 × 0.2² / 0.002 = 1.0 × 10⁵，远大于 10⁴，处于湍流区，可直接取 N_p,t = 5.2。第二步算不通气功率：P₀ = 5.2 × 1000 × 5³ × 0.2⁵ = 5.2 × 1000 × 125 × 0.00032 ≈ 104 W。单位体积功率约 1.04 kW/m³，对一般细菌发酵属于中等偏下水平。

第三步做通气修正。1.0 vvm 对应 Q_g = 0.10 m³/min = 1.67 × 10⁻³ m³/s。通气数 N_a = 1.67×10⁻³/(5 × 0.2³) = 1.67×10⁻³/0.04 ≈ 0.042。取 a = 1.2、b = 0.5（示例参数），则 P_g/P₀ = 1 − 1.2 × 0.042^0.5 ≈ 1 − 1.2 × 0.205 ≈ 0.754，即通气后功率约为 78 W。若换用更激进的关联式（a 更大），下降幅度会更大。这个结果说明：同一台设备，通气前后实际输入功率相差约四分之一，不能直接用 P₀ 去估算 k_La。

若同一体系改为非牛顿流体，K = 0.8 Pa·sⁿ、n = 0.6，则 μ_app = 0.8 × (11.5 × 5)^−0.4 ≈ 0.8 × 57.5^−0.4 ≈ 0.8 × 0.198 ≈ 0.158 Pa·s，Re_g = 1000 × 5 × 0.04/0.158 ≈ 1266。此时已落入过渡区，不能直接取 5.2，需要查 N_p–Re 曲线，过渡区 N_p 通常高于湍流值，功率会明显上升。这正是高黏度假塑性发酵液“看起来转速不高、实际功率却不低”的原因。

把上述过程整理成一张可复用的手算表，字段包括：罐径与叶轮直径、D/T、转速（r/s）、流体模型与参数、Re 或 Re_g、所选 N_p 及来源、P₀、P₀/V_L、N_a、P_g/P₀、P_g 与 P_g/V_L。每次更换叶轮类型或放大规模时，只需替换输入并重走一遍链路，就能快速判断功率是否落在工艺窗口内。