A three-dimensional coupled Nitsche and level set method for electrohydrodynamic potential flows in moving domains

被引:11
作者
Johansson, A. [1 ]
Garzon, M. [2 ]
Sethian, J. A. [3 ]
机构
[1] Simula Res Lab, Ctr Biomed Comp, POB 134, N-1325 Lysaker, Norway
[2] Univ Oviedo, Dept Appl Math, Oviedo, Spain
[3] Univ Calif Berkeley, Lawrence Berkeley Natl Lab, Dept Math, Berkeley, CA 94720 USA
基金
美国国家科学基金会;
关键词
Nitsche finite element method; Level set method; Potential flow; Electrostatic driven flows; FINITE-ELEMENT-METHOD; PINCH-OFF; BOUNDARY; DROPS;
D O I
10.1016/j.jcp.2015.12.026
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this paper we present a new algorithm for computing three-dimensional electrohydrodynamic flow in moving domains which can undergo topological changes. We consider a non-viscous, irrotational, perfect conducting fluid and introduce a way to model the electrically charged flow with an embedded potential approach. To numerically solve the resulting system, we combine a level set method to track both the free boundary and the surface velocity potential with a Nitsche finite element method for solving the Laplace equations. This results in an algorithmic framework that does not require body-conforming meshes, works in three dimensions, and seamlessly tracks topological change. Assembling this coupled system requires care: while convergence and stability properties of Nitsche's methods have been well studied for static problems, they have rarely been considered for moving domains or for obtaining the gradients of the solution on the embedded boundary. We therefore investigate the performance of the symmetric and non-symmetric Nitsche formulations, as well as two different stabilization techniques. The global algorithm and in particular the coupling between the Nitsche solver and the level set method are also analyzed in detail. Finally we present numerical results for several time-dependent problems, each one designed to achieve a specific objective: (a) The oscillation of a perturbed sphere, which is used for convergence studies and the examination of the Nitsche methods; (b) The break-up of a two lobe droplet with axial symmetry, which tests the capability of the algorithm to go past flow singularities such as topological changes and preservation of an axi-symmetric flow, and compares results to previous axi-symmetric calculations; (c) The electrohydrodynamical deformation of a thin film and subsequent jet ejection, which will account for the presence of electrical forces in a non-axi-symmetric geometry. (C) 2016 Elsevier Inc. All rights reserved.
引用
收藏
页码:88 / 111
页数:24
相关论文
共 43 条
[1]   Transport and diffusion of material quantities on propagating interfaces via level set methods [J].
Adalsteinsson, D ;
Sethian, JA .
JOURNAL OF COMPUTATIONAL PHYSICS, 2003, 185 (01) :271-288
[2]   NONLINEAR OSCILLATIONS AND BREAKUP OF CONDUCTING, INVISCID DROPS IN AN EXTERNALLY APPLIED ELECTRIC-FIELD [J].
BASARAN, OA ;
PATZEK, TW ;
BENNER, RE ;
SCRIVEN, LE .
INDUSTRIAL & ENGINEERING CHEMISTRY RESEARCH, 1995, 34 (10) :3454-3465
[3]   Singularities on charged viscous droplets [J].
Betelú, SI ;
Fontelos, MA ;
Kindelán, U ;
Vantzos, O .
PHYSICS OF FLUIDS, 2006, 18 (05)
[4]   A PENALTY-FREE NONSYMMETRIC NITSCHE-TYPE METHOD FOR THE WEAK IMPOSITION OF BOUNDARY CONDITIONS [J].
Burman, Erik .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 2012, 50 (04) :1959-1981
[5]   Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method [J].
Burman, Erik ;
Hansbo, Peter .
APPLIED NUMERICAL MATHEMATICS, 2012, 62 (04) :328-341
[6]   Ghost penalty [J].
Burman, Erik .
COMPTES RENDUS MATHEMATIQUE, 2010, 348 (21-22) :1217-1220
[7]   Fictitious domain finite element methods using cut elements: I. A stabilized Lagrange multiplier method [J].
Burman, Erik ;
Hansbo, Peter .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2010, 199 (41-44) :2680-2686
[8]   Interior-penalty-stabilized Lagrange multiplier methods for the finite-element solution of elliptic interface problems [J].
Burman, Erik ;
Hansbo, Peter .
IMA JOURNAL OF NUMERICAL ANALYSIS, 2010, 30 (03) :870-885
[9]   Simulations of Coulombic Fission of Charged Inviscid Drops [J].
Burton, J. C. ;
Taborek, P. .
PHYSICAL REVIEW LETTERS, 2011, 106 (14)
[10]   BOUNDARY INTEGRATION OVER LINEAR POLYHEDRA [J].
CATTANI, C ;
PAOLUZZI, A .
COMPUTER-AIDED DESIGN, 1990, 22 (02) :130-135