3. 生物反应器多尺度建模基础
生物反应器的数字孪生之所以难做，核心在于它同时跨越了至少四个数量级：胞内代谢反应发生在微米级、毫秒到分钟尺度；细胞群体生长在毫米到厘米级、小时尺度；气泡—液膜传质在毫米级、秒尺度；而反应器整体混合与传热则在米级、分钟尺度。多尺度建模的任务，就是把这四个层次用“状态变量 + 通量耦合”串成一条可计算、可校核的链路，而不是把四个独立模型简单拼在一起。
一、四个尺度各自“管什么”
先明确每个尺度负责的状态量与守恒对象，这是后续耦合的前提。
- 细胞代谢层（微米级）：以胞内代谢网络为对象，状态量是各代谢物浓度与通量。工程上常用简化形式——黑箱产率模型或通量平衡分析（FBA）。黑箱模型给出比生长速率 μ、比底物消耗速率 qS、比产物生成速率 qP 之间的关系，例如 qS = μ/YX/S + mS，其中 YX/S 为细胞得率，mS 为维持系数。
- 反应动力学层（群体级）：把胞内速率放大到反应器内浓度变化，典型形式是 Monod 动力学 μ = μmax·S/(KS+S)，并叠加底物抑制、产物抑制、温度修正（如 Arrhenius 项）。
- 传质传热层（界面级）：描述气—液—细胞之间的氧、二氧化碳、热量传递。核心参数是体积传质系数 kLa 与传热系数，决定反应器能否“供得上、排得掉”。
- 流体力学层（设备级）：通过搅拌桨、通气、挡板产生的流场，决定混合时间、剪切、气泡分散与局部浓度梯度，是前三个尺度在空间上的“搬运者”。
二、耦合关系：谁给谁提供边界条件
多尺度耦合不是单向的，而是“自下而上定速率、自上而下定环境”的双向闭环。下图给出本节的核心耦合结构：流体力学提供局部传质与混合，传质决定溶解氧与底物供给，动力学把供给转化为细胞生长与产物积累，代谢层再反馈出耗氧与产热需求，反过来影响流场与传质负荷。
三、关键耦合参数与量级示例
下面给出一组教学示例参数（属于合理假设，用于说明量级关系，非实测值）。假设一个 5 L 实验室搅拌罐，工作体积 3 L，温度 37 ℃，初始葡萄糖 20 g/L，接种干重 0.5 g/L。
| 参数 | 符号 | 示例值 | 所属尺度 |
|---|---|---|---|
| 最大比生长速率 | μmax | 0.45 h⁻¹ | 动力学 |
| Monod 常数 | KS | 0.15 g/L | 动力学 |
| 细胞得率 | YX/S | 0.50 g/g | 代谢 |
| 维持系数 | mS | 0.03 g/(g·h) | 代谢 |
| 体积传质系数 | kLa | 80 h⁻¹ | 传质 |
| 饱和溶氧 | C* | 7.5 mg/L | 传质 |
| 临界溶氧 | Ccrit | 1.5 mg/L | 传质/动力学 |
用这些参数可以做一个快速校核：当细胞浓度 X = 5 g/L、比耗氧速率 qO2 = 0.20 g O₂/(g·h) 时，需氧量 OUR = X·qO2 = 1.0 g O₂/(L·h)。而最大供氧能力 OTRmax = kLa·C* = 80 × 7.5 = 600 mg O₂/(L·h) = 0.6 g O₂/(L·h)。此时 OUR > OTRmax，溶解氧必然下降，动力学中的氧限制项 CL/(KO+CL) 会压低实际 μ。这说明：传质能力不足会直接改写动力学参数的有效值，这正是多尺度耦合的典型表现。
四、耦合求解的工程做法
实际建模时，通常不追求四个尺度同时全耦合求解，而是采用“分层解耦 + 迭代校正”的策略：
- 先算流场：用 CFD 或经验关联式得到 kLa、混合时间 tm、剪切率分布。若用经验式，常见形式为 kLa = C·(Pg/V)α·vsβ，其中 Pg/V 为单位体积通气功率，vs 为表观气速。
- 再算传质：以 kLa 为参数，建立溶解氧动态方程 dCL/dt = kLa(C* − CL) − OUR。
- 然后算动力学：把 CL、S、T 作为输入，计算 μ、qS、qP，更新 X、S、P。
- 最后回代代谢层：由 μ 和维持需求反算 OUR、CER、RQ，检查是否与传质层假设一致；若不一致，修正 kLa 或通气策略，进入下一轮迭代。
五、常见耦合失配与排查方向
多尺度模型出错，往往不是单层公式写错，而是层间接口对不上。典型失配有三类：
- 时间尺度失配：传质时间常数 1/kLa 约为 0.0125 h（45 s），而细胞生长时间常数 1/μ 约为 2.2 h。若用同一时间步长求解，传质动态会被“抹平”，导致溶解氧预测偏高。应使用自适应步长或分层时间推进。
- 空间尺度失配：CFD 网格若粗于气泡尺度，kLa 只能取平均值，无法反映局部缺氧区。此时应在模型中标注“kLa 为体积平均值”，避免误用为局部值。
- 参数来源失配：μmax、KS 若来自摇瓶实验，而 kLa 来自反应器关联式，两者不在同一物理环境下，直接耦合会系统性偏差。应尽量在同一反应器构型下标定动力学参数，或引入环境修正因子。