On a free-surface problem with moving contact line: From variational principles to stable numerical approximations

Abstract We analyze a free-surface problem described by time-dependent Navier–Stokes equations. Surface tension, capillary effects and wall friction are taken into account in the evolution of the system, influencing the motion of the contact line – where the free surface hits the wall – and of the dynamics of the contact angle. The differential equations governing the phenomenon are first derived from the variational principle of minimum reduced dissipation, and then discretized by means of the ALE approach. The numerical properties of the resulting scheme are investigated, drawing a parallel with the physical properties holding at the continuous level. Some instability issues are addressed in detail, in the case of an explicit treatment of the geometry, and novel additional terms are introduced in the discrete formulation in order to damp the instabilities. Numerical tests assess the suitability of the approach, the influence of the parameters, and the effectiveness of the new stabilizing terms.

[1]  L. R. Scott,et al.  The Mathematical Theory of Finite Element Methods , 1994 .

[2]  Jean-Frédéric Gerbeau,et al.  Simulations of MHD flows with moving interfaces , 2003 .

[3]  P. Sheng,et al.  Molecular Hydrodynamics of the Moving Contact Line in Two-Phase Immiscible Flows , 2005, cond-mat/0510403.

[4]  Erik Burman,et al.  Stabilization of explicit coupling in fluid-structure interaction involving fluid incompressibility , 2009 .

[5]  Shawn W. Walker,et al.  A mixed formulation of a sharp interface model of stokes flow with moving contact lines , 2014 .

[6]  Nicola Parolini,et al.  Fluid-structure interaction problems in free surface flows: application to boat dynamics , 2008 .

[7]  C. H. Bosanquet.  LV. On the flow of liquids into capillary tubes , 1923 .

[8]  R. Temam,et al.  Navier-Stokes equations: theory and numerical analysis: R. Teman North-Holland, Amsterdam and New York. 1977. 454 pp. US $45.00 , 1978 .

[9]  Vivette Girault,et al.  Finite Element Methods for Navier-Stokes Equations - Theory and Algorithms , 1986, Springer Series in Computational Mathematics.

[10]  Y. Shikhmurzaev.  Moving contact lines in liquid/liquid/solid systems , 1997, Journal of Fluid Mechanics.

[11]  E. Bänsch,et al.  An ALE finite element method for a coupled Stefan problem and Navier–Stokes equations with free capillary surface , 2013 .

[12]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[13]  Gretar Tryggvason,et al.  Direct numerical simulations of gas/liquid multiphase flows , 2011 .

[14]  Gretar Tryggvason,et al.  Direct Numerical Simulations of Gas–Liquid Multiphase Flows: Distributions concentrated on the interface , 2011 .

[16]  G. Batchelor,et al.  An Introduction to Fluid Dynamics , 1968 .

[17]  Abner J. Salgado,et al.  A diffuse interface fractional time-stepping technique for incompressible two-phase flows with moving contact lines , 2013 .

[18]  Gunilla Kreiss,et al.  A conservative level set method for contact line dynamics , 2009, J. Comput. Phys..

[19]  Anders Logg,et al.  The FEniCS Project Version 1.5 , 2015 .

[20]  Shawn W. Walker,et al.  Droplet Footprint Control , 2015, SIAM J. Control. Optim..

[21]  J.-F. Gerbeau,et al.  Generalized Navier boundary condition and geometric conservation law for surface tension , 2008, 0804.1563.

[22]  Stéphane Zaleski,et al.  A mesh-dependent model for applying dynamic contact angles to VOF simulations , 2008, J. Comput. Phys..

[23]  L. Richardson.  The Approximate Arithmetical Solution by Finite Differences of Physical Problems Involving Differential Equations, with an Application to the Stresses in a Masonry Dam , 1911 .

[24]  Shawn W. Walker,et al.  A DIFFUSE INTERFACE MODEL FOR ELECTROWETTING WITH MOVING CONTACT LINES , 2011, 1112.5758.

[25]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[26]  J. Li,et al.  Numerical simulation of moving contact line problems using a volume-of-fluid method , 2001 .

[27]  Yan Xu,et al.  Space-time discontinuous Galerkin method for nonlinear water waves , 2007, J. Comput. Phys..

[28]  Roberto F. Ausas,et al.  Variational formulations for surface tension, capillarity and wetting , 2011 .

[29]  F. Brezzi,et al.  On the Stabilization of Finite Element Approximations of the Stokes Equations , 1984 .

[30]  Alain Miranville,et al.  Cahn–Hilliard–Navier–Stokes systems with moving contact lines , 2016 .

[31]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[32]  Ricardo H. Nochetto,et al.  Geometrically Consistent Mesh Modification , 2010, SIAM J. Numer. Anal..

[33]  M. C. Delfour,et al.  Shapes and Geometries - Metrics, Analysis, Differential Calculus, and Optimization, Second Edition , 2011, Advances in design and control.

[34]  Yasufumi Yamamoto,et al.  Numerical simulations of spontaneous capillary rises with very low capillary numbers using a front-tracking method combined with generalized Navier boundary condition , 2013 .

[35]  Wei-Ming Ni,et al.  Global dynamics of the Lotka–Volterra competition–diffusion system with equal amount of total resources, II , 2016 .

[36]  Sashikumaar Ganesan,et al.  Modelling and simulation of moving contact line problems with wetting effects , 2009 .

[37]  A. Huerta,et al.  Arbitrary Lagrangian–Eulerian Methods , 2004 .

[38]  P. Sheng,et al.  A variational approach to moving contact line hydrodynamics , 2006, Journal of Fluid Mechanics.

[39]  A. Sonnet,et al.  Dissipative Ordered Fluids: Theories for Liquid Crystals , 2012 .

[40]  Anders Logg,et al.  DOLFIN: a C++/Python Finite Element Library , 2012 .

[41]  Fabio Nobile,et al.  A Stability Analysis for the Arbitrary Lagrangian Eulerian Formulation with Finite Elements , 1999 .

[42]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[43]  Sofia Guzzetti,et al.  Hierarchical model reduction for incompressible flows in cylindrical domains , 2014 .

[44]  R. Scardovelli,et al.  A variational approach to the contact angle dynamics of spreading droplets , 2009 .