CFD专栏丨为什么需要CFD+DEM耦合方法分析颗粒两相流？

颗粒相和流体相可以在不同的长度尺度上模拟，通常分为2类：宏观方法（连续介质）和微观方法（离散介质）

DEM离散单元法，就是指每个单元都是离散的，有独立特性的，也就我们常见的颗粒状物料。

离散元法的核心思想就是在拉格朗日坐标体系下，针对每个颗粒进行检索，计算由于接触产生的力，再运用牛顿第二定律进行计算颗粒的加速度/速度和位移的变化，进而得到整个系统的状态。

►Continuity equation:

►Momentum equation:

►Energy equation:

►Newton’s second law of motion (translation):

►Euler’s second law of motion (rotation):

►Energy equation (particle):

1. 颗粒的尺寸可以远远大于CFD的网格尺寸。

2. 支持球形/非球形的颗粒，例如不规则形状的石头在水中的沉降过程。

3. 可模拟能量耦合，例如热的颗粒被空气冷却

4. CFD动网格模拟刚体运动，例如颗粒通过振动的筛子或旋转的风扇叶片。

5. 可模拟传质现象，例如在颗粒表面镀液膜(coating)，或湿颗粒的干燥。

6. CFD和DEM的模型均支持GPU加速计算，对大模型有显著加速效果。

通常可以认为稀疏相的粒子不影响流场，用单向耦合就足够了，计算代价小。

如颗粒是密相，颗粒总体积/流体体积远大于千分之一，则考虑双向耦合。当然这也不是绝对的，主要是看颗粒的存在有没有显著的影响流场。

• 颗粒的存在几乎不影响流场。AcuSolve先完成流场计算后，将稳态结果导出给EDEM即可。
  

• EDEM接着计算颗粒在稳定流场下的运动轨迹。

  
  

CFD稳态流场

颗粒运动的动画

• 颗粒的存在几乎不影响流场。但是 CFD的模型参数可以随时间变化，例如流量，压力，物性参数等等。

• AcuSolve和EDEM需同时求解，AcuSolve只负责数据发送，EDEM只负责数据接收。

CFD流场动画

颗粒运动的动画

• 颗粒的存在显著影响流场。AcuSolve和EDEM需同时求解，在CFD每个迭代步结束的时刻交换流体的速度、压力、温度以及颗粒的位置信息。

• EDEM的模型时间步长通常在1e-5秒的量级。因此在acuSolve每个时间步长中（通常是0.001~0.01秒量级），EDEM已经进行了多轮的迭代。

CFD瞬态流场动画

颗粒运动的动画

下面我们看几个AcuSolve和EDEM耦合仿真的例子。

联合收割机采用气流和筛子相互配合的清选方式。谷粒混合物在筛子中振动，轻的杂质靠气流被吹走，重的杂质由尾筛排除，谷粒通过筛孔进入粮仓。清选室的气流分布对整个筛选效率有重要影响。

联合收割机原理图（图片来自网络）

联合收割机产品图（图片来自网络）

旋风分离器，是用于颗粒分离的一种设备。工作原理为靠气流切向引入造成的旋转运动，使具有较大惯性离心力的固体颗粒或液滴甩向外壁面分开。

• 入口风速13.4 m/s，3000 个颗粒，双向耦合

CFD瞬态流场动画

颗粒运动的动画

流化床，是一种利用流体通过颗粒状固体层而使固体颗粒处于悬浮运动状态，并进行气固相反应过程或液固相反应过程的反应器。

底部风速2.5 m/s，24750 个颗粒，双向瞬态耦合

制粒机，搅拌桨的旋转速度确保运动产生涡流，物料充分混合。同时，粒-刀处交汇区域对翻腾的物料团充分打碎成颗粒。搅拌桨与制粒刀的旋转速度能使物料产生三维运动，颗粒产生碰撞、磨擦、剪切，使其磨擦均匀、细致，最后形定稳定球状颗粒。

入口速度2 m/s，旋转速度600 rpm，5900个颗粒，双向瞬态耦合

气力输送是利用气流的能量，在密闭管道内沿气流方向输送颗粒状物料，是流态化技术的一种具体应用。主要特点是输送量大，输送距离长，输送速度较高；能在一处装料，然后在多处卸料。

• 入口速度15 m/s，一百万粒子，双向耦合

颗粒的气力输运

医用干粉吸入器，可用于模拟干粉药物输送到呼吸道的过程，并识别颗粒积聚的区域。模拟颗粒积聚和吸入器效率。

• 空气体积流量(50 L/min) ，5000 个颗粒

叶片泵模型，泵送带有颗粒的液体，双向耦合。

生产线喷涂和干燥过程，在密闭腔体中向固体表面喷射低温液滴，环境是高温空气 

VolumeAdded参数查看固体表面有多少水分。当开始喷涂时，颗粒表面水分（ VolumeAdded ）增加。一旦停止喷射，涂层开始干燥，水分又逐渐减小。

真空吸尘器，模拟长条形状的杂物被负压吸入清扫的过程。