热表面张力驱动的对流是微重力下浮区法晶体生长中熔体最重要的物质与热输运方式。采用单松弛双分布函数格子Boltzmann模型,自主开发了相应的格子Boltzmann方法的串行和MPI并行程序包,并应用该程序包对开口方腔内流体的二维热表面张力对流进行了数值模拟研究。其中串行程序合并碰撞迁移过程和引入临时数组以连续读入分布函数,相比分开碰撞迁移过程,计算性能提高了二倍;在此基础上,采用单向计算区域分区和非阻塞通信模式,实现了MPI版格子Boltzmann并行程序包开发。对比基于传统有限体积法CFD程序计算结果表明,串行和MPI并行版格子Boltzmann程序包计算结果精确可靠;并行程序具有较好的性能。
Thermocapillary flow driven by an unbalanced surface tension plays the most important role for mass and heat transport in floating zone melt crystal growth under microgravity. The authors develop both a serial and a parallel codes with lattice Bhatnagar-Gross-Krook (BGK) model using two distribution functions, and numerically study the thermocapillary flow in a two dimensional square cavity with a single free surface under microgravity condition. The serial code is developed by combining collision and propagation step, using a temporary array to continuously read distribution functions, and its performance is improved two times faster than the code with separating the collision and propagation step. The MPI parallel code is proposed utilizing one dimensional partitioning and non-blocking communication. The accurate and reliable results are achieved with both the serial and parallel codes by comparing with CFD results simulated by the finite volume method; the MPI parallel code has high performance.