Xu Xingchun, Hu Yanwei, Dai Bing, Yang Lei, Han Jiecai, He Yurong, Zhu Jiaqi
National Key Laboratory of Science and Technology on Advanced Composites in Special Environments, Harbin Institute of Technology, Harbin 150080, China.
School of Energy Science & Engineering, Harbin Institute of Technology, Harbin 150001, China.
Phys Rev E. 2021 Sep;104(3-2):035305. doi: 10.1103/PhysRevE.104.035305.
Based on the phase-field theory, a multiple-relaxation-time (MRT) lattice Boltzmann model is proposed for the immiscible multiphase fluids. In this model, the local Allen-Chan equation is chosen as the target equation to capture the phase interface. Unlike previous MRT schemes, an off-diagonal relaxation matrix is adopted in the present model so that the target phase-field equation can be recovered exactly without any artificial terms. To check the necessity of removing those artificial terms, comparative studies were carried out among different MRT schemes with or without correction. Results show that the artificial terms can be neglected at low March number but will cause unphysical diffusion or interface undulation instability for the relatively large March number cases. The present modified model shows superiority in reducing numerical errors by adjusting the free parameters. As the interface transport coupled to the fluid flow, a pressure-evolution lattice Boltzmann equation is adopted for hydrodynamic properties. Several benchmark cases for multiphase flow were conducted to test the validity of the present model, including the static drop test, Rayleigh-Taylor instability, and single rising bubble test. For the rising bubble simulation at high density ratios, bubble dynamics obtained by the present modified MRT lattice Boltzmann model agree well with those obtained by the FEM-based level set and FEM-based phase-field models.
基于相场理论,提出了一种用于不混溶多相流体的多松弛时间(MRT)格子玻尔兹曼模型。在该模型中,选择局部艾伦 - 陈方程作为捕捉相界面的目标方程。与先前的MRT格式不同,本模型采用非对角松弛矩阵,从而可以精确恢复目标相场方程而无需任何人工项。为了检验去除这些人工项的必要性,对有无修正的不同MRT格式进行了对比研究。结果表明,在低马赫数下人工项可以忽略不计,但对于相对较大马赫数的情况会导致非物理扩散或界面波动不稳定性。本改进模型通过调整自由参数在减少数值误差方面表现出优势。由于界面传输与流体流动耦合,采用压力演化格子玻尔兹曼方程来描述流体动力学性质。进行了几个多相流基准案例测试本模型的有效性,包括静态液滴测试、瑞利 - 泰勒不稳定性测试和单个上升气泡测试。对于高密度比下的气泡上升模拟,本改进的MRT格子玻尔兹曼模型得到的气泡动力学与基于有限元的水平集模型和基于有限元的相场模型得到的结果吻合良好。