個数密度関数法（PBM）
理论与物理
概述
老师，群体平衡模型（PBM）是什么？
PBM（群体平衡模型）是追踪气泡、液滴、颗粒的尺寸分布在时空上变化的模型。它描述了尺寸分布通过聚并（coalescence）、破碎（breakup）、成核（nucleation）、生长（growth）而发生变化的过程。
在什么场景下使用呢？
在分散相尺寸分布至关重要的各类问题中都会使用，例如鼓泡塔中的气泡径分布、乳化过程中的液滴径分布、结晶中的晶体尺寸分布、气溶胶的粒径变化等。通常与欧拉-欧拉法结合使用。
控制方程
请告诉我PBM的方程。
体积为 $V$ 的颗粒（气泡/液滴）的数量密度 $n(V, \mathbf{x}, t)$ 的输运方程如下。
右边的各项表示以下含义。
- $B_{coal}$: 聚并生成（两个较小颗粒聚并形成尺寸V的颗粒）
- $D_{coal}$: 聚并消亡（尺寸V的颗粒与其他颗粒聚并变得更大）
- $B_{break}$: 破碎生成（大颗粒破碎产生尺寸V的颗粒）
- $D_{break}$: 破碎消亡（尺寸V的颗粒破碎变小）
聚并和破碎的速率是如何建模的？
展示一些典型的闭合模型。
| 过程 | 模型示例 | 驱动机制 |
|---|---|---|
| 聚并频率 | Prince & Blanch (1990) | 湍流碰撞 + 液膜排液 |
| 聚并频率 | Luo (1993) | 湍流能量 |
| 破碎频率 | Luo & Svendsen (1996) | 湍流涡破碎 |
| 破碎频率 | Martínez-Bazán (2010) | 惯性-表面张力平衡 |
Sauter平均直径
如何从尺寸分布求得代表直径？
最常用的是Sauter平均直径 $d_{32}$。
$m_k = \int_0^\infty V^{k/3} n(V) dV$ 是分布的 $k$ 阶矩。$d_{32}$ 对应于体积-表面积比，适用于传质和反应速率的评估。
群体平衡——作为多相流理论“进化论”的工程学
群体平衡方程（PBE）是描述具有尺寸、年龄、组成等“内部坐标”的群体在时空上变化的通用框架。它在数学上与生物群体动态（Lotka-Volterra方程）同构，可以统一处理气泡分布、晶体粒径分布、细胞浓度分布等。化学工程中PBE的应用始于Randolph & Larson（1971年）的结晶论文，50年后的今天，它已进化为与CFD耦合的C-PBE（耦合PBE），成为药品制造、鼓泡塔、液液萃取的设计工具。
各项的物理意义
- 时间项 $\partial(\rho\phi)/\partial t$：想象一下拧开水龙头的瞬间。最初水会不稳定地喷溅出来，过一会儿才会变成稳定的水流，对吧？描述这个“正在变化的过程”的就是时间项。心脏搏动导致血流脉动，发动机阀门每次开闭引起流动变化，这些都是非定常现象。那么定常分析是什么？就是只观察“经过足够时间后流动稳定下来之后”的状态——也就是将此项设为零。由于计算成本大幅降低，先用定常求解是CFD的基本策略。
- 对流项 $\nabla \cdot (\rho \mathbf{u} \phi)$：把落叶丢进河里会怎样？它会随着水流被带到下游，对吧？这就是“对流”——流体运动搬运物质的效果。暖气的热风能到达房间另一端，也是因为空气这个“搬运工”通过对流输送热量。这里有趣的是——这项包含“速度×速度”，因此是非线性的。也就是说，流速变快时此项会急剧增强，变得难以控制。这就是湍流的根本原因。常见的误解：“对流和传导差不多”→ 完全不一样！对流是流动搬运，传导是分子传递。效率有天壤之别。
- 扩散项 $\nabla \cdot (\Gamma \nabla \phi)$：有过在咖啡里倒入牛奶后放置不管的经历吗？即使不搅拌，过一会儿也会自然混合。那就是分子扩散。那么下一个问题——蜂蜜和水，哪个更容易流动？当然是水。因为蜂蜜的粘度（$\mu$）高，所以难流动。粘度越大，扩散项越强，流体的运动就变得“粘稠”。雷诺数小的流动（缓慢、粘稠）中扩散占主导。相反，Re数大的流动中对流占压倒性优势，扩散则成为配角。
- 压力项 $-\nabla p$：按压注射器的活塞，液体就会从针头有力地射出，对吧？为什么呢？因为活塞侧压力高，针头侧压力低——这个压力差成为推动流体的力。水坝放水也是同样原理。天气图上等压线密集的地方会怎样？没错，会刮强风。“有压力差的地方就会产生流动”——这就是纳维-斯托克斯方程压力项的物理意义。这里的误解点：CFD的“压力”多为表压而非绝对压力。切换到可压缩分析时结果突然出错，原因可能就是混淆了绝对压力/表压。
- 源项 $S_\phi$：被加热的空气会上升——为什么呢？因为比周围空气轻（密度低），被浮力推上去了。这个浮力作为源项添加到方程中。此外，燃气灶火焰产生化学反应热、工厂电磁泵对金属熔液施加洛伦兹力……这些都是“从外部向流体注入能量或力”的作用，都用源项表示。忘记源项会怎样？自然对流分析中忘记加入浮力，流体就完全不动——冬天房间里开了暖气，暖空气却不上浮，就会得到这种物理上不可能的结果。
假设条件与适用范围
- 连续介质假设：克努森数 Kn < 0.01（分子平均自由程 ≪ 特征长度）时成立
- 牛顿流体假设：剪切应力与应变速率呈线性关系（非牛顿流体需要粘度模型）
- 不可压缩假设（Ma < 0.3时）：将密度视为常数。马赫数0.3以上需考虑可压缩性效应
- Boussinesq近似（自然对流）：仅在浮力项中考虑密度变化，其他项使用恒定密度
- 不适用情况：稀薄气体（Kn > 0.1）、超音速/高超音速流动（需要捕捉激波）、自由表面流动（需要VOF/Level Set等）
量纲分析与单位制
| 变量 | SI单位 | 注意点·换算备忘 |
|---|---|---|
| 速度 $u$ | m/s | 入口条件从体积流量换算时，注意截面积单位 |
| 压力 $p$ | Pa | 区分表压与绝对压力。可压缩分析使用绝对压力 |
| 密度 $\rho$ | kg/m³ | 空气: 约1.225 kg/m³@20°C、水: 约998 kg/m³@20°C |
| 粘度系数 $\mu$ | Pa·s | 注意与运动粘度系数 $\nu = \mu/\rho$ [m²/s] 混淆 |
| 雷诺数 $Re$ | 无量纲 | $Re = \rho u L / \mu$。层流/湍流转换的判定指标 |
| CFL数 | 无量纲 | $CFL = u \Delta t / \Delta x$。直接关系到时间步长的稳定性 |
数值解法与实现
数值解法详情
请告诉我PBM的数值解法。
PBE是以体积（尺寸）为附加独立变量的积分微分方程，直接求解困难。使用以下离散化方法。
| 方法 | 概要 | 计算成本 | 精度 |
|---|---|---|---|
| MUSIG (多尺寸组) | 将尺寸空间分割为多个区间 | 高（与区间数成正比） | 依赖区间数 |
| QMOM (求积矩法) | 矩法 + 求积公式 | 低（6〜8个矩） | 良好 |
| DQMOM | 直接求积矩法 | 低 | 良好 |
| S-Gamma | 2参数分布假设 | 最低 | 分布形状有约束 |
请告诉我MUSIG和QMOM的区别。
MUSIG将尺寸空间分割为10〜30个区间（尺寸组），将每个区间的数量密度作为输运方程求解。精度高，但需要求解与区间数相同的附加输运方程，计算成本大。
QMOM（McGraw, 1997）不直接求解尺寸分布，而是求解前几个矩（$m_0, m_1, ..., m_5$等）的输运方程。通过Product-Difference Algorithm（PDA）或Wheeler法确定求积公式的节点和权重，计算聚并、破碎的源项。
Fluent中的实现
Ansys Fluent中可通过Eulerian Multiphase + Population Balance Model使用。
| 方法 | Fluent名称 | 区间/矩数 |
|---|---|---|
| 离散方法 | Size Group | 10〜30个区间 |
| QMOM | Quadrature MOM | 6个矩 |
| DQMOM | Direct QMOM | 2〜4个环境 |
| SMM | Standard MOM | 6个矩 |
STAR-CCM+中的实现
STAR-CCM+中标准的是S-Gamma模型（Lo, 2000），用2个参数（平均直径和分散度）追踪分布。计算成本最低，但分布形状受限于伽马分布。MUSIG也可用。
OpenFOAM中的实现
OpenFOAM中，v2006版本后可使用 populationBalanceModel 类。实现了MUSIG（含iMUSIG）和QMOM。在 constant/phaseProperties 内的 populationBalanceCoeffs 中进行设置。
QMOM vs DQMOM——矩量法的权衡
为了将CFD-PBE的计算成本控制在现实可行的范围内，使用了矩量法。QMOM（求积矩法）通过求积点近似分布函数，仅求解矩方程来追踪整个分布的时间演化。但在聚并、破碎复杂的系统中，会出现求积点交叉的“矩反演问题”。DQMOM（直接求积矩法）通过直接求解每个求积点位置和权重的输运方程来避免此问题，但源项的非线性性使收敛变得困难。在鼓泡塔CFD的基准测试中，QMOM和DQMOM预测的气泡径分布存在10〜20%差异的情况很多，模型选择对结果有不可忽视的影响。
迎风格式（Upwind）
一阶迎风：数值扩散大但稳定。二阶迎风：精度提高但有振荡风险。高雷诺数流动中必备。
中心差分（Central Differencing）
二阶精度，但Pe数 > 2时会产生数值振荡。适用于低雷诺数的扩散主导流动。
TVD格式（MUSCL、QUICK等）
通过限制器函数抑制数值振荡同时保持高精度。对捕捉激波或陡峭梯度有效。
有限体积法 vs 有限元法
FVM：自然满足守恒定律。CFD的主流。FEM：对复杂形状、多物理场有利。SPH等无网格法也在发展中。
CFL条件（库朗数）
显式方法：CFL ≤ 1为稳定条件。隐式方法：即使CFL > 1也稳定，但影响精度和迭代次数。LES：推荐CFL ≈ 1。物理意义：一个时间步内信息前进不超过一个网格。
残差监控
连续性方程、动量、能量的各项残差下降3〜4个数量级可判断为收敛。质量守恒的残差尤其重要。
松弛因子
压力：0.2〜0.3、速度：0.5〜0.7为一般初始值。发散时降低松弛因子。收敛后可提高以加速。
非定常计算的内部迭代
在每个时间步内迭代直至收敛到定常解。内部迭代次数：5〜20次为参考值。残差在时间步间波动时需重新审视时间步长。
SIMPLE法的比喻
SIMPLE法是“交替调整”的方法。先假设求解速度（预测步），然后根据该速度修正压力以满足质量守恒（修正步），再用修正后的压力修正速度——重复这种“投接球”过程以接近正确答案。类似于两人调整架子水平的作业：一人调整高度，另一人调整平衡，如此反复。
迎风格式的比喻
迎风格式是“站在河流中重视上游信息”的方法。站在河里的人看下游也无法知道水的来源——反映了上游信息决定下游的物理原理的离散化方法。精度为一阶，但能正确捕捉流动方向，因此稳定性高。
实践指南
实践指南
相关主题
なった
詳しく
報告