• 论文 •

### 不同密度与粘性的多相流移动接触线问题的自适应有限元方法

1. 北京应用物理与计算数学研究所计算物理重点实验室, 北京 100094
• 收稿日期:2015-08-31 出版日期:2015-12-15 发布日期:2015-12-18
• 基金资助:

国家自然科学基金资助项目(1141033, 1141034, 11471049, 11475033);国家高技术研究发展计划(2012AA01A303);中国工程物理研究院院长基金资助项目(201402037);中国工程物理研究院科学技术发展基金资助项目(2015B0202033,2015B0202035);计算物理重点实验室基金资助项目(9140C690104150C69303).

Shi Yi. AN ADAPTIVE FINITE ELEMENT METHOD FOR MOVING CONTACT LINE PROBLEMS WITH DIFFERENT DENSITIES AND VISCOSITIES[J]. Journal of Numerical Methods and Computer Applications, 2015, 36(4): 297-309.

### AN ADAPTIVE FINITE ELEMENT METHOD FOR MOVING CONTACT LINE PROBLEMS WITH DIFFERENT DENSITIES AND VISCOSITIES

Shi Yi

1. Institute of Applied Physics and Computational Mathematics, Beijing, China
• Received:2015-08-31 Online:2015-12-15 Published:2015-12-18

In this paper, we propose an adaptive finite element method for the moving contact line problems of multi-phase flows with different densities and viscosities. The Cahn-Hilliard-Navier-Stokes system with the generalized Navier boundary condition is employed for our simulation. Splitting method is employed for the system in time:for the Cahn-Hilliard equation, we use the convex-splitting method, and for the Navier-Stokes equation, we use the pressure stabilization method. The method is efficient for problems with large density ratios, and is energy stable under certain conditions. We use the LBB stable finite element method in space. Moreover, we propose the adaptive mesh refinement method based on the a posteriori error estimates for the system. Numerical simulations demonstrates the accuracy and efficiency of our method.

MR(2010)主题分类:

()
 [1] Kai Bao, Yi Shi, Xiao-Ping Wang and Shuyu Sun. A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems[J], J. Comput. Phys., 2012, 231:8083-8099.[2] Bonometti T and Magnaudet J. Transition from spherical cap to toroidal bubbles[J], Phys. Fluids, 2006, 18:052102.[3] Hang Ding, Spelt P D M and Chang Shu. Diffuse interface model for incompressible two-phase flows with large density ratios[J], J. Comput. Phys., 2007, 226:2078-2095.[4] Dussan E B and Davis S H. On the motion of a fluid-fluid interface along a solid surface[J], J. Fluid Mech., 1974, 65:71-95.[5] Min Gao and Xiao-Ping Wang. A gradient stable scheme for a phase field model for the moving contact line problem[J], J. Comput. Phys., 2012, 231:1372-1386.[6] Min Gao and Xiao-Ping Wang. An efficient scheme for a phase field model for moving contact line problem with variable density and viscosity[J], J. Comput. Phys., 2014, 272:704-718.[7] Girault V and Raviart P. Finite element methods for Navier-Stokes Equations:Theory and Algorithms[M], Springer, Berlin-Heidelberg-New York, 1981.[8] Guermond J L and Salgado A. A splitting method for incompressible flows with variable density based on a pressure Poisson equation[J], J. Comput. Phys., 2009, 228:2834-2846.[9] Shuo Guo, Min Gao, Xiaomin Xiong, Yongjian Wang, Xiao-Ping Wang, Ping Sheng and Penger Tong, Direct measurement of friction of a fluctuating contact line[J], Phys. Rev. Lett., 2013, 111:026101.[10] Qiaolin He, Glowinski R and Xiao-Ping Wang. A least-squares/finite element method for the numerical solution of the Navier-Stokes-Cahn-Hilliard system modeling the motion of the contact line[J], J. Comput. Phys., 2011, 230:4991-5009.[11] Jacqmin D. Calculation of two-phase Navier-Stokes flows using phase-field modeling[J], J. Comput. Phys., 1999, 155:96-127.[12] Tiezheng Qian, Xiao-Ping Wang, and Ping Sheng. Molecular scale contact line hydrodynamics of immiscible flows[J], Phys. Rev. E, 2003, 68:016306.[13] Tiezheng Qian, Xiao-Ping Wang, and Ping Sheng. Molecular hydrodynamics of the moving contact line in two-phase immiscible flows[J], Commun. Comput. Phys., 2005, 1:1-52.[14] Tiezheng Qian, Xiao-Ping Wang, and Ping Sheng. A variational approach to moving contact line hydrodynamics[J], J. Fluid Mech., 2006, 564:333-360.[15] Jie Shen and Xiaofeng Yang. A phase-field model and its numerical approximation for two-phase incompressible flows with different densities and viscosities[J], SIAM J. Sci. Comput., 2010, 32:1159-1179.[16] Yi Shi, Kai Bao, Xiao-Ping Wang. 3D adaptive finite element method for a phase field model for moving contact line problems[J], Inverse Probl. Imaging, 2013, 947-959.[17] Yi Shi and Xiao-Ping Wang. Modeling and simulation of dynamics of three-component flows on solid surface[J], Jpn. J. Ind. Appl. Math., 2014, 31:611-631.[18] Xiao-Ping Wang, Tiezheng Qian, and Ping Sheng. Moving contact line on chemically patterned surfaces[J], J. Fluid Mech., 2008, 605:59-78.[19] Pengtao Yue, James J. Feng, Chun Liu, and Jie Shen. A diffuse-interface method for simulating two-phase flows of complex fluids[J], J. Fluid Mech., 2004, 515:293-317.
 [1] 刘辉, 冷伟, 崔涛. 并行自适应有限元计算中的负载平衡研究[J]. 数值计算与计算机应用, 2015, 36(3): 166-184. [2] 刘辉, 崔涛, 冷伟. hp自适应有限元计算中一种新的自适应策略[J]. 数值计算与计算机应用, 2015, 36(2): 100-112. [3] 吴斌, 刘志峰, 王晓宏. 油水两相流非活塞式驱替的数值算法探讨[J]. 数值计算与计算机应用, 2012, 33(4): 274-282. [4] 刘辉, 张林波. 自适应有限元常用标记策略及其在PHG中的并行实现[J]. 数值计算与计算机应用, 2009, 30(4): 315-320.