The extended discontinuous Galerkin method adapted for moving contact line problems via the generalized Navier boundary condition

被引:4
作者
Smuda, Martin [1 ,2 ]
Kummer, Florian [1 ,2 ]
机构
[1] Tech Univ Darmstadt, Chair Fluid Dynam, Otto Berndt Str 2, D-64287 Darmstadt, Hessen, Germany
[2] Tech Univ Darmstadt, Grad Sch Computat Engn, Darmstadt, Germany
关键词
extended; unfitted discontinuous Galerkin method; generalized Navier boundary condition; sharp interface formulation; transient two-phase flow; FINITE-ELEMENT-METHOD; FIELD; SINGULARITIES; FLOWS; MODEL; ANGLE;
D O I
10.1002/fld.5016
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this work, an extended discontinuous Galerkin (extended DG/XDG also called unfitted DG) solver for two-dimensional flow problems exhibiting moving contact lines is presented. The generalized Navier boundary condition is employed within the XDG discretization for the handling of the moving contact lines. The spatial discretization is based on a symmetric interior penalty method and the numerical treatment of the surface tension force is done via the Laplace-Beltrami formulation. The XDG method adapts the approximation space conformal to the position of the interface and allows a sub-cell accurate representation within the sharp interface formulation. The interface is described as the zero set of a signed-distance level-set function and discretized by a standard DG method. No adaption of the level-set evolution algorithm is needed for the extension to moving contact line problems. The developed solver is validated against typical two-dimensional contact line driven flow phenomena including droplet simulations on a wall and the two-phase Couette flow.
引用
收藏
页码:2921 / 2945
页数:25
相关论文
共 43 条
[1]   A FAST LEVEL SET METHOD FOR PROPAGATING INTERFACES [J].
ADALSTEINSSON, D ;
SETHIAN, JA .
JOURNAL OF COMPUTATIONAL PHYSICS, 1995, 118 (02) :269-277
[2]   The fast construction of extension velocities in level set methods [J].
Adalsteinsson, D ;
Sethian, JA .
JOURNAL OF COMPUTATIONAL PHYSICS, 1999, 148 (01) :2-22
[3]   Hybrid scheduling for the parallel solution of linear systems [J].
Amestoy, PR ;
Guermouche, A ;
L'Excellent, JY ;
Pralet, S .
PARALLEL COMPUTING, 2006, 32 (02) :136-156
[4]   A fully asynchronous multifrontal solver using distributed dynamic scheduling [J].
Amestoy, PR ;
Duff, IS ;
L'Excellent, JY ;
Koster, J .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 2001, 23 (01) :15-41
[5]   AN INTERIOR PENALTY FINITE-ELEMENT METHOD WITH DISCONTINUOUS ELEMENTS [J].
ARNOLD, DN .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1982, 19 (04) :742-760
[6]   A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems [J].
Bao, Kai ;
Shi, Yi ;
Sun, Shuyu ;
Wang, Xiao-Ping .
JOURNAL OF COMPUTATIONAL PHYSICS, 2012, 231 (24) :8083-8099
[7]   An unfitted finite element method using discontinuous Galerkin [J].
Bastian, Peter ;
Engwer, Christian .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2009, 79 (12) :1557-1576
[8]   A minimization-based finite element formulation for interface-preserving level set reinitialization [J].
Basting, Christopher ;
Kuzmin, Dmitri .
COMPUTING, 2013, 95 (01) :S13-S25
[9]   The velocity field near moving contact lines [J].
Chen, Q ;
Rame, E ;
Garoff, S .
JOURNAL OF FLUID MECHANICS, 1997, 337 :49-66
[10]   Numerical simulation of static and sliding drop with contact angle hysteresis [J].
Dupont, Jean-Baptiste ;
Legendre, Dominique .
JOURNAL OF COMPUTATIONAL PHYSICS, 2010, 229 (07) :2453-2478