A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems

被引：62

作者：

Bao, Kai ^{[1
,2
]}

Shi, Yi ^{[2
]}

Sun, Shuyu ^{[1
]}

Wang, Xiao-Ping ^{[2
]}

机构：

[1] King Abdullah Univ Sci & Technol, Div Math & Comp Sci & Engn, Thuwal 239556900, Saudi Arabia

[2] Hong Kong Univ Sci & Technol, Dept Math, Kowloon, Hong Kong, Peoples R China

来源：

JOURNAL OF COMPUTATIONAL PHYSICS | 2012年 / 231卷 / 24期

关键词：

Moving contact line; Generalized Navier boundary condition; Finite element method; Convex splitting; Navier-Stokes equations; Cahn-Hilliard equations; CONSERVATION; INTERFACE; SURFACES; CHANNEL; MOTION; FLOWS; MODEL;

D O I：

10.1016/j.jcp.2012.07.027

中图分类号：

TP39 [计算机的应用];

学科分类号：

081203 ; 0835 ;

摘要：

In this paper, a semi-implicit finite element method is presented for the coupled Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn-Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier-Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples. (C) 2012 Elsevier Inc. All rights reserved.

引用

页码：8083 / 8099

页数：17

共 16 条

[1] A column pre-ordering strategy for the unsymmetric-pattern multifrontal method [J].

Davis, TA .

ACM TRANSACTIONS ON MATHEMATICAL SOFTWARE, 2004, 30 (02) :165-195

[2] MOTION OF A FLUID-FLUID INTERFACE ALONG A SOLID-SURFACE [J].

DUSSAN, EB ;

DAVIS, SH .

JOURNAL OF FLUID MECHANICS, 1974, 65 (AUG12) :71-&

[3] Error analysis of a mixed finite element method for the Cahn-Hilliard equation [J].

Feng, XB ;

Prohl, A .

NUMERISCHE MATHEMATIK, 2004, 99 (01) :47-84

[4] Fully discrete finite element approximations of the Navier-Stokes-Cahn-Hilliard diffuse interface model for two-phase fluid flows [J].

Feng, Xiaobing .

SIAM JOURNAL ON NUMERICAL ANALYSIS, 2006, 44 (03) :1049-1072

[5] A gradient stable scheme for a phase field model for the moving contact line problem [J].

Gao, Min ;

Wang, Xiao-Ping .

JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 231 (04) :1372-1386

[6] Generalized Navier boundary condition and geometric conservation law for surface tension [J].

Gerbeau, J. -F. ;

Lelievre, T. .

COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2009, 198 (5-8) :644-656

[7]

Girault V., 1981, Finite Element Method for Navier-Stokes Equations: Theory and Algorithms

[8] 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].

He, Qiaolin ;

Glowinski, Roland ;

Wang, Xiao-Ping .

JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 230 (12) :4991-5009

[9] Displacement of a two-dimensional immiscible droplet in a channel [J].

Kang, QJ ;

Zhang, DX ;

Chen, SY .

PHYSICS OF FLUIDS, 2002, 14 (09) :3203-3214

[10] Moving contact line over undulating surfaces [J].

Luo, Xiongping ;

Wang, Xiao-Ping ;

Qian, Tiezheng ;

Sheng, Ping .

SOLID STATE COMMUNICATIONS, 2006, 139 (11-12) :623-629

← 1 2 →