4. 量纲分析与无量纲数

在生物反应器的尺度转换中，最危险的做法是“按几何比例放大一切”。反应器体积放大1000倍时，若把搅拌转速也按比例放大，叶端速度、剪切率、混合时间、氧传递速率会各自走向完全不同的方向。量纲分析的价值就在于：它不依赖具体设备细节，而是从物理量的量纲出发，找出控制过程的无量纲群，从而判断哪些效应可以保持相似、哪些必然发生偏离。

Buckingham π定理：从变量清单到无量纲群

Buckingham π定理指出：若一个物理过程涉及 n 个有量纲变量，这些变量涉及 m 个基本量纲（在生物反应器流体力学中通常为质量 M、长度 L、时间 T 三个），则可用 n−m 个独立的无量纲数描述该过程。操作上分四步：列出影响过程的全部变量；写出每个变量的量纲；确定基本量纲数 m；选取 m 个包含全部基本量纲的核心变量，将其余变量逐一与核心变量组合，消去量纲，得到 π 群。

以搅拌槽中的氧传递为例，影响体积传质系数 k_La 的变量可列为：功率输入 P、液体体积 V、液体密度 ρ、黏度 μ、表面张力 σ、气泡直径 d_b、扩散系数 D、氧浓度差 ΔC、重力 g。取 ρ、V、P 为核心变量，可导出若干 π 群，其中最重要的包括：

Re  = ρ N d² / μ          （惯性力 / 黏性力）
Fr  = N² d / g             （惯性力 / 重力）
We  = ρ N² d³ / σ          （惯性力 / 表面张力）
Sh  = k_L d / D            （对流传质 / 分子扩散）
Pe  = N d² / D 或 u L / D  （对流传质 / 分子扩散）
Da  = 最大反应速率 / 最大传质速率

其中 N 为搅拌转速，d 为桨径，u 为特征速度，L 为特征长度。这些 π 群不是数学游戏，而是把“哪些力在主导”“传质与反应谁更快”这类问题压缩成了可比较的数值。

雷诺数 Re：流态与混合的判据

Re = ρNd²/μ 表征惯性力与黏性力之比。在生物反应器中，Re 决定流态：Re < 10 为层流，10～10⁴ 为过渡流，Re > 10⁴ 为充分湍流。对典型细菌发酵液，ρ≈1000 kg/m³、μ≈1～10 mPa·s；对高黏度丝状真菌或植物细胞培养液，μ 可达 100～1000 mPa·s，同样转速下 Re 会下降一到两个数量级，混合时间显著延长。因此放大时若保持 Re 相同，意味着 N 需按 d⁻² 缩放，大罐转速会低到无法悬浮细胞；若保持叶端速度 N·d 相同，则 Re 随 d 增大而增大，剪切又可能损伤细胞。这正是尺度转换中“相似性冲突”的典型来源。

舍伍德数 Sh 与佩克莱数 Pe：传质的两条路径

Sh = k_Ld/D 是传质系数 k_L 的无量纲化形式，物理意义是“对流传质速率与分子扩散速率之比”。Sh 越大，说明界面更新越快、传质越强。Pe = uL/D 则描述同一竞争在整体流动尺度上的表现：Pe 很大时，对流远快于扩散，混合主要靠涡流输运；Pe 很小时，扩散成为主导，混合缓慢且均匀。

在生物反应器中，Sh 与 Pe 常通过经验关联式连接，例如湍流鼓泡塔中 Sh ∝ Re^0.7Sc^1/3，其中 Sc = μ/(ρD) 为施密特数。对氧气在水中的 D≈2×10⁻⁹ m²/s，Sc≈500，因此 Sh 对 Re 的变化相当敏感。放大时若 Re 下降，Sh 随之下降，k_La 可能不足以维持细胞耗氧，表现为溶氧迅速跌至临界值以下。

丹姆科勒数 Da：反应与传质的竞争

Da = 最大反应速率 / 最大传质速率。在生物反应器中，它直接回答一个关键问题：过程是“传质限制”还是“反应限制”。若 Da ≪ 1，反应速率远低于传质速率，溶氧和底物在液相中基本均匀，反应器可按理想混合处理；若 Da ≫ 1，传质成为瓶颈，细胞实际接触到的底物或氧浓度远低于主体浓度，表现为“表观动力学”偏离本征动力学。

以耗氧为例，若细胞比耗氧速率 q_O₂≈5×10⁻⁴ mol/(g·s)、细胞浓度 X≈20 g/L，则体积耗氧速率 OUR≈0.01 mol/(L·s)。若 k_La≈0.05 s⁻¹、饱和溶氧 C*≈0.2 mmol/L，最大传质速率 OTR≈k_La·C*≈1×10⁻⁵ mol/(L·s)，远小于 OUR，Da≫1，说明该条件下必然出现氧限制。此时提高转速或通气量以增大 k_La，比修改菌株动力学更直接有效。

放大中的相似性冲突与取舍

把上述无量纲数放在一起看，尺度转换的本质就清楚了：不可能同时保持所有 π 群不变。工程上通常按过程瓶颈选择“主导相似准则”。例如：

• 溶氧限制型（高密度好氧发酵）：优先保持 k_La 或 OTR 相同，允许 Re 和混合时间偏离。
• 剪切敏感型（动物细胞、丝状真菌）：优先保持叶端速度 N·d 或剪切率相同，接受传质能力下降。
• 混合限制型（高黏度、大体积）：优先保持混合时间或 Pe 相同，往往需要牺牲部分传质效率。
• 反应限制型（Da≪1）：反应动力学主导，放大时重点控制温度、pH 和底物浓度分布，传质相似性要求可放宽。

实际计算时，可先由小试数据拟合出 k_La 与 P/V 或 N 的关联式，再结合目标罐的几何参数反推所需转速与通气量，最后用 Re、Sh、Da 复核是否落入可接受区间。例如小试 5 L 罐在 N=300 rpm、d=0.07 m 下 k_La=0.08 s⁻¹，放大到 500 L 罐 d=0.25 m，若要求 k_La 不低于 0.06 s⁻¹，按 k_La ∝ (P/V)^0.7 估算，P/V 需从约 1.5 kW/m³ 提高到约 2.5 kW/m³，对应转速约 120～150 rpm，此时 Re 仍处于湍流区，但叶端速度已接近细胞耐受上限，需要结合剪切实验数据做最终判断。