4. 发酵罐内流体力学基础
搅拌发酵罐的性能，最终取决于罐内流体“怎么动”。同样一套桨叶、同样的转速，若流型不对，可能出现局部死区、溶氧不均、剪切过强或功率虚耗。理解搅拌槽内的流型、循环流与湍流特性，以及挡板、导流筒如何改变流动边界，是把搅拌设计从“凭经验选桨”推进到“按目标定量分析”的基础。本节围绕这些流动现象展开，并给出关键无量纲数的物理含义与工程用法。
一、搅拌槽内的基本流型：轴向流、径向流与切向流
搅拌桨旋转时，叶轮对流体施加的是“推”和“甩”两种作用，由此形成三种基本流动分量。轴向流沿搅拌轴方向上下运动，典型产生者是45°斜叶桨、螺旋桨；径向流从桨叶沿半径方向甩向罐壁，再沿壁面上下折返，典型产生者是六直叶圆盘涡轮；切向流则绕轴做圆周运动，几乎不产生上下交换，是“只转不混”的主要来源。
实际流型是三者叠加的结果，但主导分量决定了混合效率。轴向流擅长整体循环与固体悬浮，功耗较低；径向流擅长气体分散与局部高剪切，常用于通气发酵；切向流在无挡板时占主导，会形成中心下凹的旋涡（俗称“打旋”），液面中心暴露、混合效果急剧下降。因此判断流型是否合理，第一步就是看切向分量是否被有效抑制。
二、循环流：混合时间的决定因素
循环流描述流体在罐内完成一次“上行—下行”大回路的能力，其强弱用循环流量 Qc 与循环次数表示。循环流量近似与叶轮排量成正比，即 Qc ∝ N·D³（N为转速，D为叶轮直径），而罐内液体体积 V 决定一次完整循环所需时间。工程上常用“每分钟循环次数”Nc = Qc/V 来粗略评估整体混合强度。
循环流不足的直接后果是混合时间 tm 变长，罐内出现浓度梯度：补料区局部糖浓度过高、pH 探头附近与主体不一致、溶氧在死区偏低。对于高径比大的发酵罐，单层桨往往难以覆盖全罐，需要多层桨组合，让下层桨负责气体分散、上层桨负责整体循环。判断循环是否够用，可结合混合时间与发酵周期比较：若 tm 已接近或超过关键反应时间尺度，说明循环流是瓶颈。
三、湍流特性：局部剪切与能量耗散
湍流不是“越猛越好”。叶轮附近是强湍流区，湍动能耗散率 ε 最高，可达到罐平均值的数十倍；罐壁与底部则是低湍流区。湍流的作用有两面：一方面促进气泡破碎、提高传质系数 kLa；另一方面过高的局部剪切会损伤丝状菌、剪切敏感细胞或絮体。
工程上关注两个量：一是平均能量耗散率 ε̄ = P/(ρV)，反映整体输入功率密度；二是最大局部耗散率 εmax，与桨叶尖端线速度 utip = πND 密切相关。设计时通常先由传质或悬浮需求确定 ε̄，再用 utip 校核剪切是否超过菌体耐受上限。这也是为什么“提高转速”不能无限制使用——它同时抬高了传质和剪切。
四、挡板与导流筒：用边界条件塑造流型
挡板的作用是打断切向流。在罐壁安装4块宽度约为罐径 1/10～1/12 的垂直挡板，可把原本绕轴的圆周运动转化为上下翻腾的轴向—径向组合流，显著缩短混合时间、提高功率准数。挡板过多或过宽会过度增加功耗而不改善混合，通常4块标准挡板已足够；对于高黏体系，挡板效果减弱，需要靠桨型本身提供轴向推力。
导流筒则是“主动引导”而非“被动阻挡”。在导流筒反应器中，筒体把流体约束成明确的上升通道与下降环隙，形成定向内循环。其优点是循环路径确定、剪切分布可控、气体停留时间可调；缺点是结构复杂、清洗难度增加。导流筒常用于气升式或强制循环发酵罐，配合低剪切桨实现高循环倍率。
五、关键无量纲数：从流态判断到传质校核
无量纲数把几何、物性与操作条件压缩成可比较的判据，是搅拌放大的核心工具。下面几个数在发酵罐分析中出现频率最高。
| 无量纲数 | 定义式 | 物理含义 | 工程用法 |
|---|---|---|---|
| 雷诺数 Re | ρND²/μ | 惯性力与黏性力之比 | Re<10 层流；10～10⁴ 过渡；>10⁴ 湍流 |
| 弗劳德数 Fr | N²D/g | 惯性力与重力之比 | 判断液面下凹与打旋倾向，无挡板时尤其重要 |
| 功率准数 P₀ | P/(ρN³D⁵) | 无量纲功率输入 | 湍流区近似为常数，用于由转速反算功率 |
| 努塞尔数 Nu | kLd/DAB | 对流传质与分子扩散之比 | 关联传质系数，评估溶氧能力 |
需要说明的是，Nu 在搅拌传质中常以经验关联式出现，例如 kLa 与 (P/V)α·usβ 的幂律关系，其中 us 为表观气速。不同桨型、不同体系指数不同，使用时必须匹配相近几何与物性范围，不能跨体系直接套用。
六、一个教学示例：如何用这些量做初步判断
假设一台示例发酵罐：罐径 T = 1.0 m，液深 H = 1.8 m，装液量 V ≈ 1.4 m³，采用六直叶圆盘涡轮，D = 0.33 m，转速 N = 3 r/s（180 r/min），发酵液密度 ρ = 1000 kg/m³，黏度 μ = 0.002 Pa·s。以下数值为教学假设，用于演示判断流程，并非实测结果。
Re = ρND²/μ = 1000 × 3 × 0.33² / 0.002 ≈ 1.63 × 10⁵ → 充分湍流
u_tip = πND = 3.14 × 3 × 0.33 ≈ 3.11 m/s → 剪切需校核
Fr = N²D/g = 9 × 0.33 / 9.81 ≈ 0.30 → 无挡板时打旋明显
P₀（六直叶涡轮，湍流区）≈ 5.5
P = P₀ρN³D⁵ = 5.5 × 1000 × 27 × 0.0039 ≈ 580 W
ε̄ = P/(ρV) ≈ 580 / 1400 ≈ 0.41 W/kg判断逻辑是：Re 已进入湍流区，说明可以用湍流关联式；Fr 偏大提示必须装挡板，否则液面下凹、混合恶化；utip 约 3.1 m/s 对丝状菌偏高，若菌体剪切敏感，应降低转速或改用大直径低转速桨；ε̄ 约 0.4 W/kg 属于中等功率密度，可结合 kLa 关联式估算是否满足溶氧需求。若计算出的 kLa 低于菌体耗氧速率对应的临界值，则需提高 P/V 或增加通气量，而不是盲目加大桨径。
把上述关系串起来看：挡板与导流筒改变的是流动边界，从而改变流型；流型决定循环流强弱与湍流分布；循环流与湍流又通过 Re、Fr、P₀、Nu 等无量纲数被量化，最终落到混合时间、传质系数与剪切损伤这些可判断的指标上。掌握这条链条，才能在面对具体发酵体系时，有依据地选择桨型、转速与内构件，而不是停留在“经验上差不多”的层面。