01

EMMS-DPM 粗粒化模型 

为了兼顾计算精度与效率，中国科学院过程工程研究所团队创新性地将能量最小多尺度（Energy-Minimization Multi-Scale, EMMS）模型与离散颗粒法（DPM）相结合，提出了粗粒化离散颗粒法（EMMS-DPM），该模型将粗颗粒尺寸限定在团聚物尺寸以下，从底层描述团聚物演化过程；同时根据颗粒动理论修正粗颗粒碰撞参数，保证粗粒化体系与真实颗粒体系的碰撞能量损失一致。核心思路是 "粗粒化"—— 不再追踪单个真实颗粒，而是将运动状态相近的一群颗粒打包成一个 "粗粒化颗粒"（CGP），这种处理方式大幅减少了需要模拟的颗粒数量，从根源上降低了计算成本。

粗粒化并非简单的 "合并同类项"，关键在于如何保证模拟的准确性。EMMS 模型可以精准确定粗粒化颗粒的关键参数：其粒径必须小于气固流场中自然形成的颗粒团特征尺寸，既保证计算效率又不丢失介尺度结构信息；内部空隙率则根据颗粒团的疏密程度动态调整，确保能真实反映颗粒群的堆积状态。同时，基于颗粒流动理学理论，修正了粗粒化颗粒的碰撞恢复系数，保证能量耗散特性与真实颗粒群一致。

在相间作用力计算方面，可以采用 EMMS 曳力模型，充分考虑了流场的非均质性，相比传统曳力模型大幅提升了计算精度。而在计算硬件层面，EMMS-DPM 方法巧妙利用 CPU-GPU 异构计算架构：CPU 负责复杂的气相流场求解，GPU 则承担海量粗粒化颗粒的简单碰撞计算，让硬件性能得到最大化发挥。

压降结果高度吻合：EMMS-DPM 各粗粒化比（k2/k3/k4）模拟得到的床层压降，与传统 DPM 模拟结果、经典 Ergun 欧根方程计算值均高度一致；流化起始后，所有模拟结果的压降变化趋势也保持统一，且略低于理论流化临界压降，该偏差与实际流化的非均匀性规律相符，符合工程实际。

颗粒分布统计特征一致：不同粗粒化比下，相同时刻的颗粒具体分布因初始 / 边界条件敏感性存在差异，但稳态下的统计分布特征保持可比，均能捕捉到固定床流化过程中的均匀膨胀、鼓泡等核心流动模式。

粗粒化比存在合理限值：粗粒化比 k 增大会降低固相描述分辨率，使颗粒分布细节变模糊，当 k=4 时虽床层整体膨胀特性仍接近，但颗粒分布已出现定性扭曲；说明粗粒化比并非越大越好，其限值与系统中粗粒化颗粒的总数量相关，小型测试系统的 k 不宜过大。

床层重心高度变化趋势统一：表观气速低于 0.3m/s 时，EMMS-DPM 与传统 DPM 模拟的床层重心高度几乎恒定，变化趋势完全一致；高流速下的微小偏差，主要由粗粒化带来的壁面效应（大粒径颗粒与壁面间距大，床层空隙率略升）和颗粒数量减少导致的流场均匀性波动引起。

粗粒化未改变核心流化特性：在粗粒化比选取合理的前提下，EMMS-DPM 模拟的固定床核心流化特性（压降、床层膨胀、重心高度变化）均与传统 DPM 基本一致，证明该方法的粗粒化颗粒物性关联、相互作用模型是合理的，为后续将其应用于更复杂的循环流化床模拟奠定了基础。

二维 Horio 提升管：DPM vs EMMS-DPM 

传统 DPM 模拟得到的空隙率分布云图

传统DPM方法：

• 真实颗粒数量：420 万个

• 固相计算占比：93%

• 计算极慢，只能 2.14 秒 / 天

• 无法用于更大尺度、更高颗粒量的模拟

案例 e78k6cp70H 在 EMMS–DPM 模拟下，5–10 秒时间平均的空隙率分布及颗粒团演化过程（底部 0.8 米区域）

基于同一装置开展 EMMS-DPM 模拟，先通过 EMMS 模型确定合理粗粒化比并设置不同粗粒化颗粒内部空隙率、弹簧刚度开展多组对比模拟，结果显示，在粗粒化参数符合 EMMS 模型建议值的前提下，EMMS-DPM 的轴向空隙率分布与实验数据的吻合度优于传统 DPM，能精准捕捉颗粒团形成、破碎、融合的动态演化过程，且弹簧刚度对模拟结果影响极小，仅粗粒化比接近 EMMS 模型上限时密相区和过渡区精度略有下降，同时 EMMS-DPM 通过粗粒化处理大幅减少了追踪颗粒数量、增大了颗粒时间步长，实现了计算效率的大幅跃升致。

在计算速度方面，将算例e73k6cp48中EMMS–DPM和DPM模拟进行了比较，其中EMMS-DPM 跟踪的粒子数量减少了18.72倍，粒子时间步长增加了10倍，计算速度从 2.14s/天提升至 32s/天，效率提升一个数量级以上。

中科院过程工程研究所VPE平台的循环流化床提升管

先依据 EMMS 模型确定该提升管体系的颗粒团特征参数，选取粗粒化比 k=6 开展 3.54m/s、5.62m/s 两种表观气速的 2D 模拟，将归一化轴向压降分布与实验数据对比，结果整体高度吻合，且内部空隙率更贴合 EMMS 模型建议值的模拟组，精度表现更优，还能清晰捕捉颗粒团形成、破碎、融合的动态演化规律。

(a) 表观气速 Ug=3.54 m/s下，EMMS–DPM 模拟得到的 10–15 s 时间平均空隙率分布；(b) 对应工况下的颗粒团演化过程（仅展示底部1.5m区域）

(a) 5～10 s 时间平均的轴向空隙率分布；(b) 密相区、过渡区与稀相区内，5～10 s 时间平均的径向空隙率分布；(c) 不同区域、不同时刻下高度为 5 mm 的切片内瞬时颗粒分布（俯视图）

2D与3D模拟得到的 5～10 s 时间平均轴向空隙率分布，与实验结果对比

核心突破为5.62m/s 气速下的 3D 工业级规模模拟：针对模拟域内约 75.6 亿个真实颗粒，EMMS-DPM 通过 k=10 的粗粒化处理，将模拟颗粒量大幅缩减至 1450 万个粗粒化颗粒，结合轴向负载均衡策略与 120 个 CPU-GPU 异构进程并行计算，实现2 小时完成 1 秒真实时间的流场演化，解决了传统 DPM 在工业级规模下计算量过大、难以实现的核心痛点。

同时，该模型在工业级 3D 模拟中仍能精准还原提升管的多尺度非均质性流场特征：轴向清晰呈现底部密相区、中部过渡区、上部稀相区的典型分布，径向捕捉到 “壁面密相 - 中心稀相” 的核心 - 环核结构，还能精准反映该高密相提升管的工业级特殊流动特征，模拟结果与实验数据总体一致。

EMMS-DPM 模型核心优势：

超强的计算效率与工业级可扩展性：通过粗粒化处理实现海量真实颗粒的大幅缩减（75.6 亿→1450 万），结合 CPU-GPU 异构并行计算，突破了传统 DPM 在工业级规模下的计算瓶颈，实现了工业级 3D 气固流模拟的高效推进。

精度保持性与参数适配性：无论是 2D 多工况模拟还是 3D 工业级模拟，模型计算结果均与实验数据高度吻合，且可通过 EMMS 模型精准匹配粗粒化参数（如内部空隙率），进一步提升模拟精度，在大幅提升效率的同时，未丢失流场核心特征与关键规律。

复杂流场的精准捕捉能力：针对更接近工业实际的复杂气固流体系，模型能精准还原轴向、径向的多尺度非均质性结构，还能捕捉到工业级装置特有的流动特征（如高密相提升管的颗粒整体向上运动、底部颗粒加速），体现了对复杂工程流场的良好适配性与描述能力。

2D 到 3D 的灵活拓展性：模型可无缝从 2D 小工况模拟拓展至 3D 工业级规模模拟，且在拓展过程中无需大幅调整核心算法，仅需根据体系特征适配粗粒化参数与计算策略，展现了良好的工程应用灵活性。

内容来源：Liqiang Lu, Ji Xu, Wei Ge, etal. EMMS-based discrete particle method (EMMS–DPM) for simulation of gas–solid flows. [J] Chemical Engineering Science 120 (2014) 67-87.<button></button>