颗粒团绕流曳力系数的LBM计算

采用格子Boltzmann方法(LatticeBoltzmannmethod,LBM)中的LBGK(LatticeBhatnagar-Gross-Krook)模型和二阶精度的曲线边界条件处理方法对二维颗粒团绕流现象进行了数值模拟，并同时使用动量交换法计算了两种颗粒团构型中不同颗粒的曳力系数。结果表明：颗粒团曳力系数与颗粒聚团的构型有着密切联系，颗粒聚团的形成将导致颗粒团曳力系数大幅度减小。除颗粒团构型因素外，颗粒间距和流动Reynolds数也是导致颗粒团曳力系数发生改变的主要因素。当颗粒聚团存在时，颗粒团中不同颗粒的受力有较大差异，若忽略颗粒聚团效应，则颗粒团曳力系数的计算必然将产生偏差。