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

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.

[1]  John R. Lister,et al.  Capillary pinch-off in inviscid fluids , 2003 .

[2]  Carlo Cattani,et al.  Boundary integration over linear polyhedra , 1990, Comput. Aided Des..

[3]  Leonard J. Gray,et al.  Numerical simulation of non‐viscous liquid pinch off using a coupled level set boundary integral method , 2007, J. Comput. Phys..

[4]  Peter Hansbo,et al.  Fictitious domain finite element methods using cut elements: I. A stabilized Lagrange multiplier method , 2010 .

[5]  J A Sethian,et al.  Numerical simulations of electrostatically driven jets from nonviscous droplets. , 2014, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  J. Sethian,et al.  A coupled level set-boundary integral method for moving boundary simulations , 2005 .

[7]  Andrew H. Gee,et al.  Regularised marching tetrahedra: improved iso-surface extraction , 1999, Comput. Graph..

[8]  L. Rayleigh.  XX. On the equilibrium of liquid conducting masses charged with electricity , 1882 .

[9]  P. Hansbo,et al.  An unfitted finite element method, based on Nitsche's method, for elliptic interface problems , 2002 .

[10]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods , 1999 .

[11]  O. Vantzos,et al.  Singularities on charged viscous droplets , 2006 .

[12]  Peter Hansbo,et al.  Interior-penalty-stabilized Lagrange multiplier methods for the finite-element solution of elliptic interface problems , 2010 .

[13]  Akos Vertes,et al.  Astable regime in electrosprays. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[15]  William E. Lorensen,et al.  Marching cubes: a high resolution 3D surface construction algorithm , 1996 .

[16]  Marco A. Fontelos,et al.  Evolution of neutral and charged drops in an electric field , 2008 .

[17]  J. Sethian,et al.  Transport and diffusion of material quantities on propagating interfaces via level set methods , 2003 .

[18]  J. Sethian,et al.  A coupled level set projection method applied to ink jet simulation , 2003 .

[19]  J. Strutt.  Scientific Papers: On the Equilibrium of Liquid Conducting Masses charged with Electricity , 2009 .

[20]  J A Sethian,et al.  Simulation of the droplet-to-bubble transition in a two-fluid system. , 2011, Physical review. E, Statistical, nonlinear, and soft matter physics.

[21]  R. Grimm,et al.  Dynamics of field-induced droplet ionization: time-resolved studies of distortion, jetting, and progeny formation from charged and neutral methanol droplets exposed to strong electric fields. , 2005, The journal of physical chemistry. B.

[22]  Mats G. Larson,et al.  $L^2$-error estimates for finite element approximations of boundary fluxes , 2014 .

[23]  P. Hansbo,et al.  Fictitious domain finite element methods using cut elements , 2012 .

[24]  Erik Burman,et al.  A Penalty-Free Nonsymmetric Nitsche-Type Method for the Weak Imposition of Boundary Conditions , 2011, SIAM J. Numer. Anal..

[25]  S Nudurupati,et al.  Electrohydrodynamic removal of particles from drop surfaces. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Geoffrey Ingram Taylor,et al.  Electrically driven jets , 1969, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[27]  Leonard J. Gray,et al.  Axisymmetric boundary integral formulation for a two‐fluid system , 2012 .

[28]  J. M. López-Herrera,et al.  Dynamical behavior of electrified pendant drops , 2013 .

[29]  David L. Chopp,et al.  Some Improvements of the Fast Marching Method , 2001, SIAM J. Sci. Comput..

[30]  Ming Xue,et al.  Effects of microphysical drop size distribution on tornadogenesis in supercell thunderstorms , 2008 .

[31]  Geoffrey Ingram Taylor,et al.  Disintegration of water drops in an electric field , 1964, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[32]  J. Nitsche.  Über ein Variationsprinzip zur Lösung von Dirichlet-Problemen bei Verwendung von Teilräumen, die keinen Randbedingungen unterworfen sind , 1971 .

[33]  Robert T. Collins,et al.  Universal scaling laws for the disintegration of electrified drops , 2013, Proceedings of the National Academy of Sciences.

[34]  P. Hansbo,et al.  A FINITE ELEMENT METHOD ON COMPOSITE GRIDS BASED ON NITSCHE'S METHOD , 2003 .

[35]  H. A. Stone,et al.  Non-coalescence of oppositely charged drops , 2009, Nature.

[36]  Mats G. Larson,et al.  A high order discontinuous Galerkin Nitsche method for elliptic problems with fictitious boundary , 2013, Numerische Mathematik.

[37]  C Guet,et al.  Shape deformations of surface-charged microdroplets. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  Alfonso M Gañán-Calvo,et al.  Revision of capillary cone-jet physics: electrospray and flow focusing. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[39]  T. Russell,et al.  Electrically induced structure formation and pattern transfer , 2000, Nature.

[40]  P. Taborek,et al.  Simulations of coulombic fission of charged inviscid drops. , 2011, Physical review letters.

[41]  Robert E. Benner,et al.  Nonlinear Oscillations and Breakup of Conducting, Inviscid Drops in an Externally Applied Electric Field , 1995 .

[42]  R. Collins,et al.  Electrohydrodynamic tip streaming and emission of charged drops from liquid cones , 2008 .