悬浮颗粒运动的格子Boltzmann数值模拟

将固体颗粒的牛顿力学和格子Boltzmann方法相结合,研究不规则形状悬浮颗粒在流场中的运动。通过受力分析,精确求得其所受合力、合力矩、合力作用中心等。提出了跟随颗粒运动的动网格计算域技术和模拟悬浮颗粒转动运动的局部数组方法及Euler-Lagrange两套坐标技术。通过对椭圆颗粒运动的数值模拟和对照他人对矩形颗粒的研究,分析了其复杂运动规律,并提供了合理的物理解释。结果表明:运用格子Boltzmann方法和上述特殊技术可以得到与有限元方法相同的模拟精度,且具有计算速度快、对复杂形状边界处理方便灵活、程序简单及特别适合大规模并行计算等优点。