4. 欧拉-拉格朗日离散相模型（DPM）
在欧拉-拉格朗日框架下，连续相（气体或液体）仍在欧拉网格上求解质量、动量与能量守恒，而离散相（颗粒、液滴、气泡）则被当作拉格朗日质点，在流场中逐个积分运动方程，从而获得每个颗粒的完整轨迹。这一节围绕 DPM 的四个核心问题展开：颗粒怎么被追踪、受哪些力、湍流如何让轨迹发散、以及颗粒如何反过来影响连续相。理解这四点，才能判断一个喷雾或气力输送仿真结果是否可信。
一、颗粒轨迹追踪：从运动方程到积分步长
Fluent 对单个颗粒求解牛顿第二定律，在拉格朗日坐标系下，颗粒速度 up 的变化由单位质量颗粒所受合力决定。最基础的形式是球形颗粒的曳力主导方程：
du_p/dt = F_D (u - u_p) + g(ρ_p - ρ)/ρ_p + F_other
F_D = 18μ/(ρ_p d_p²) · C_D Re_p/24
Re_p = ρ d_p |u_p - u| / μ其中 FD 是曳力系数项，Rep 为颗粒雷诺数，CD 由 Morsi-Alexander 或 Schiller-Naumann 关联式给出。Fother 是可选附加力。轨迹追踪的关键参数是积分时间步长：Fluent 提供自动步长（基于颗粒松弛时间 τp=ρpdp²/(18μ)）和固定步长两种方式。工程经验是步长取 0.1τp～0.5τp，并设置最大追踪步数（如 500～5000）与长度尺度，避免颗粒在回流区无限循环。判断追踪是否收敛，可观察“incomplete particles”数量：若大量颗粒因达到最大步数而终止，说明步长过小或流场存在强旋涡，需要调整。
二、颗粒受力：哪些力必须开，哪些可以忽略
除曳力与重力外，Fluent 允许开启多种附加力，是否开启取决于颗粒尺度与流场特征：
- 虚拟质量力（Virtual Mass）：当 ρp/ρ < 10 时显著，例如气泡在水中；喷雾中水滴与空气密度比约 800，可忽略。
- 压力梯度力（Pressure Gradient）：流场强加速区（如喷嘴出口）需考虑，量级约为曳力的 ρ/ρp 倍。
- Saffman 升力：剪切层中颗粒横向迁移的主因，气力输送弯管处不可忽略。
- Magnus 升力：颗粒自旋时启用，需配合旋转模型。
- 布朗力 / 热泳力：亚微米颗粒（dp<1 μm）才显著。
实操建议：先只开曳力+重力跑通，再逐项加入并对比颗粒轨迹与沉积位置，避免一次性开启全部力导致结果无法归因。
三、湍流扩散：颗粒随机游走与涡相互作用
湍流使颗粒轨迹具有随机性。Fluent 提供两种模型：随机游走模型（Discrete Random Walk, DRW）和颗粒云模型（Cloud Model）。DRW 中，颗粒与湍流涡相互作用，涡的特征时间由 TL=CLk/ε 给出，颗粒穿越涡的时间为 tcross，相互作用时间取二者较小值。每经过一个涡，颗粒速度方向被随机扰动，扰动幅度与 √(2k/3) 成正比。DRW 需要设置“number of tries”，通常 5～20 次，次数越多统计越稳定但计算量线性增长。判断湍流扩散是否合理：对比不同 tries 下颗粒浓度分布，若变化小于 5% 即可认为统计收敛。
四、颗粒与连续相耦合：单向、双向与四向
耦合方式决定颗粒是否影响流场：
| 耦合方式 | 适用条件 | 典型场景 |
|---|---|---|
| 单向耦合 | 体积分数 < 10-6 | 稀疏示踪颗粒 |
| 双向耦合 | 10-6～10-3 | 喷雾、稀相气力输送 |
| 四向耦合 | > 10-3 | 密相输送、流化床 |
双向耦合通过源项将颗粒的动量、质量、能量交换反馈到连续相方程。Fluent 中需设置“Interaction with Continuous Phase”，并指定每多少个连续相迭代步更新一次颗粒源项（通常 1～10）。源项松弛因子不宜过大，否则易发散。判断耦合是否收敛：监控连续相残差与颗粒源项总量，若源项在若干步后趋于稳定，说明耦合已平衡。
五、工业案例：喷雾与气力输送的关键设置
案例一：压力旋流喷嘴喷雾。连续相为空气，离散相为水，dp 用 Rosin-Rammler 分布，最小 10 μm、最大 200 μm、平均 60 μm、分布指数 2.5。注入方式选“pressure-swirl atomizer”，设置喷口直径 0.8 mm、流量 0.05 kg/s、喷雾锥角 60°。开启双向耦合与 DRW（tries=10），曳力用 Schiller-Naumann。关注蒸发模型是否开启（若液滴在高温气流中）。结果判读：观察不同轴向截面的 Sauter 平均直径（SMD）与速度分布，若 SMD 沿轴向快速衰减，说明蒸发或二次破碎设置合理。
案例二：稀相气力输送弯管。颗粒为塑料颗粒，dp=3 mm，密度 950 kg/m³，质量流量 0.2 kg/s，空气速度 20 m/s。弯管处开启 Saffman 升力与双向耦合，壁面边界设为“reflect”并给定恢复系数（法向 0.8、切向 0.9）。重点观察弯管外侧壁面颗粒浓度与磨损率分布。若颗粒大量穿透壁面，检查步长是否过大或恢复系数是否合理。判断输送是否稳定：监测出口颗粒质量流量与入口偏差应小于 2%。
总结本节要点：轨迹追踪的核心是步长与最大步数；受力选择取决于密度比与流场特征；湍流扩散靠 DRW 的 tries 控制统计精度；耦合方式由体积分数决定。把这三者与具体案例参数对应起来，才能让 DPM 仿真既稳定又具有工程参考价值。