Numerical Stability and Accuracy of Contact Angle Schemes in Pseudopotential Lattice Boltzmann Model for Simulating Static Wetting and Dynamic Wetting  

作　　者：Dongmin Wang Gaoshuai Lin 

机构地区：[1]Key Laboratory of Multiphase Flow and Heat Transfer in Shanghai Power Engineering,School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai,200093,China [2]Jiangsu Key Laboratory of Micro and Nano Heat Fluid Flow Technology and Energy Application,Suzhou University of Science and Technology,Suzhou,215009,China

出　　处：《Computer Modeling in Engineering & Sciences》2023年第10期299-318,共20页工程与科学中的计算机建模（英文）

基　　金：sponsored by the National Natural Science Foundation of China under Grant No.52206101;Shanghai Sailing Program under Grant No.20YF1431200;the Experiments for Space Exploration Program and the Qian Xuesen Laboratory,China Academy of Space Technology under Grant No.TKTSPY-2020-01-01.

摘　　要：There are five most widely used contact angle schemes in the pseudopotential lattice Boltzmann(LB)model for simulating the wetting phenomenon:The pseudopotential-based scheme(PB scheme),the improved virtualdensity scheme(IVD scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the fluid layer density above the wall(MPB-C scheme),the modified pseudopotential-based scheme with a ghost fluid layer constructed by using the weighted average density of surrounding fluid nodes(MPB-W scheme)and the geometric formulation scheme(GF scheme).But the numerical stability and accuracy of the schemes for wetting simulation remain unclear in the past.In this paper,the numerical stability and accuracy of these schemes are clarified for the first time,by applying the five widely used contact angle schemes to simulate a two-dimensional(2D)sessile droplet on wall and capillary imbibition in a 2D channel as the examples of static wetting and dynamic wetting simulations respectively.(i)It is shown that the simulated contact angles by the GF scheme are consistent at different density ratios for the same prescribed contact angle,but the simulated contact angles by the PB scheme,IVD scheme,MPB-C scheme and MPB-W scheme change with density ratios for the same fluid-solid interaction strength.The PB scheme is found to be the most unstable scheme for simulating static wetting at increased density ratios.(ii)Although the spurious velocity increases with the increased liquid/vapor density ratio for all the contact angle schemes,the magnitude of the spurious velocity in the PB scheme,IVD scheme and GF scheme are smaller than that in the MPB-C scheme and MPB-W scheme.(iii)The fluid density variation near the wall in the PB scheme is the most significant,and the variation can be diminished in the IVD scheme,MPB-C scheme andMPBWscheme.The variation totally disappeared in the GF scheme.(iv)For the simulation of capillary imbibition,the MPB-C scheme,MPB-Wscheme and GF scheme simulate the dynamics of the liq

关 键 词：Pseudopotential lattice Boltzmann model contact angle scheme static wetting dynamic wetting capillary imbibition 

分 类 号：TN91[电子电信—通信与信息系统]