体积分数算法
EDEM-Fluent耦合接口两种版本中计算颗粒体积分数的方法不同。其中,在欧拉多相流框架(2.0)中使用的是蒙特卡洛方法(图g所示)计算颗粒的体积分数。接口程序将在颗粒内部随机生成N个sample point点,每个point点所代表的颗粒体积为Vp/N,接口程序将根据point点在网格中所处的位置,判定其属于哪个网格。网格中颗粒的体积分数由处在该网格内所有颗粒的point点数决定。
图[1] 体积分数计算公式
基于DPM/DDPM框架的耦合接口(2.2)使用的是DDPM模型计算颗粒体积分数的方法——中心点法(图b所示)。
根据上述对两种模型的简要描述,当需要模拟大颗粒在流场中的行为时,可能会导致无法满足流体网格体积大于颗粒体积10倍以上的条件,DPM模型采用中心点法,无法调整体积分数,可考虑采用欧拉法,更改体积分数计算方法,进而进行模拟。