1 引言\n\n计算流体力学与离散单元法（CFD-DEM）耦合模拟是研究颗粒-流体两相系统的有力工具，已广泛应用于化工、能源、制药和冶金等工业过程的机理分析与设备优化。LIGGGHTS作为开源DEM软件，因其灵活的扩展性和并行效率而成为耦合模拟的重要选择。耦合模拟产生的海量数据带来了存储、分析和可视化的挑战，促使专业数据处理服务成为该领域不可或缺的技术支撑。本文系统介绍CFD-DEM耦合方法、LIGGGHTS的数据特性及相关数据处理服务，旨在为研究者提供全流程的参考框架。\n\n## 2 CFD-DEM耦合方法\n\n### 2.1 基本原理\n\nCFD-DEM耦合通过双向传递流体与颗粒间的相互作用力实现两相同步求解：流体相满足局部平均的Navier-Stokes方程，在欧拉网格上采用有限体积法离散；颗粒相遵循牛顿运动方程，在拉格朗日框架下由DEM追踪。流体对颗粒的曳力、压力梯度力、虚拟质量力等耦合项构成双向往返的数据链路。\n\n### 2.2 耦合策略（求解层面同步/交替）\n\n按照流体与颗粒求解的同步方式，耦合策略可分为三类：\n\n- 联合求解（紧密耦合）：流体与颗粒方程在每个时间步内同步隐式求解，稳定性好但实现复杂、计算量大；\n- 交替求解（松耦合）：在每一耦合步内分别求解流体与颗粒，通过交换源项迭代达到收敛，适用于颗粒浓度较低的场景；\n- 顺序求解（显式耦合）：为避免时间步悬殊导致的不稳定，采用伪时间推进将颗粒-流体交换量近似而保守地传递。\n\n对于稠密颗粒系统，常采用基于局部平均的CFD-DEM方法；对于涉及四向耦合中颗粒-颗粒及颗粒-壁面碰撞-摩擦效应的描述，多使用软球模型并调用弹性-阻尼接触力模型。\n\n### 2.3 时间尺度比与稳定性\n\n颗粒碰撞时间通常小于流体CFL时间步，二者相差数个量级。为保证数值稳定，常取耦合时间步为min(Coupled\_DT)，并允许颗粒域内��循环加密度权衡计算效率。高效的耦合并行需构建网格-颗粒双向查找（邻域搜索）以提升E‑L (Euler-Lagrange) 映射性能。\n\n### 2.4 耦合工具链与实现软件\n\n当前主流耦合实现包括：\n\n- CFDEM®coupling：首个开源CFD-DEM双向耦合平台，整合OpenFOAM®与LIGGGHTS，支持曳力模型库及IMF稠密颗粒模型等特性，经广泛验证；\n- OpenFOAM + LIGGGHTS原位耦合：针对标准OpenFOAM开发通讯与蒙特卡洛接口，提供了高并行解的支持路径，特别用于气固流化床及颗粒输运特性分析；\n- MESO-Foam与自行开发接口：作为对耦合器件的降阶方案之一。\n\n实际选择时需根据体系颗粒数、流体网格柱状比和精度需求平衡实现成本和硬件。\n\n## 3 LIGGGHTS数据处理需求\n\nDEM/LIGGGHTS的输出在大规模模拟中可达TB(s)级：单个典型案例可取250万含轨迹更新的多重相及变量位置属性，向每一dump原子（一个球团点的输出涉及轨迹高版本下“导出原子截面“>导致额外toc成本以版本编辑和I/O格式要求I/O能力提效20x4来驱动该算法采用具有LAMMP执行结果相同的sttig元素规格标