EDEM-CFD耦合的时候，出现计算“卡”在计算体积分数“Calculating volume fraction...”这一步，是什么原因？

EDEM-CFD耦合在计算体积分数的时候，采用的是蒙特卡洛取样点法。就是按照蒙特卡洛算法，在颗粒区域将其分解为若干个取样点（sample point），然后按照取样点的坐标和CFD网格节点的坐标进行位置匹配，查找取样点所在的网格，最后统计每个CFD网格中取样点的数量，以求取颗粒体积分数。如果在计算过程中，有颗粒（或某个取样点）离开了流体域，那么在匹配取样点和网格位置时，就会把流体域所有网格全部循环查找一遍。一旦网格数较多，或者流体域以外的取样点很多，就会花费大量的时间在循环网格和取样点上面，导致耦合计算看起来“卡”在了“Calculating volume fraction...”这一步。一般要注意五种可能出现此问题的情况：1、入口。注意将颗粒工厂放在流体域内，防止颗粒在外面生成。2、出口。注意将EDEM求解域适当缩小，保证颗粒在离开流体域前先离开EDEM的计算域而被删除掉。3、壁面。注意颗粒剪切模量和时间步长的设置，防止颗粒“透”壁面而出。4、对称及周期边界。也需要适当缩小EDEM求解域，防止颗粒不会部分漏出。5、Interface边界。大小不一致的interface边界会自动生成虚拟壁面，这个壁面在EDEM读取msh网格的时候不会被识别，需要手动单独添加。