On the calculation of normals in free‐surface flow problems

The use of boundary-conforming finite element methods is considered for the solution of surface-tension-dominated free-surface flow problems in three dimensions. This class of method is based upon the use of a moving mesh whose velocity is driven by the motion of the free surface, which is in turn determined via a kinematic boundary condition for the normal velocity. The significance of the method used to compute the normal direction at the finite element node points for a C0 piecewise-polynomial free surface is investigated. In particular, it is demonstrated that the concept of mass-consistent normals on an isoparametric quadratic tetrahedral mesh is flawed. In this case an alternative, purely geometric, normal is shown to lead to a far more robust numerical algorithm. Copyright © 2004 John Wiley & Sons, Ltd.

[1]  Stojan Petelin,et al.  Numerical errors of the volume‐of‐fluid interface tracking algorithm , 2002 .

[2]  Jeffrey J. Derby,et al.  An assessment of a parallel, finite element method for three‐dimensional, moving‐boundary flows driven by capillarity for simulation of viscous sintering , 2001 .

[3]  P. Jimack,et al.  On the stability of viscous free–surface flow supported by a rotating cylinder , 2001, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[4]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[5]  P. A. Sackinger,et al.  A finite element method for free surface flows of incompressible fluids in three dimensions. Part I. Boundary fitted mesh motion , 2000 .

[6]  P. A. Sackinger,et al.  A Finite Element Method for Free-Surface Flows of Incompressible Fluids in Three Dimensions, Part II: Dynamic Wetting Lines , 2000 .

[7]  Peter K. Jimack,et al.  The solution of two-dimensional free-surface problems using automatic mesh generation , 1999 .

[8]  Sandra M. Troian,et al.  The development of transient fingering patterns during the spreading of surfactant coated films , 1999 .

[9]  Yousef Saad,et al.  An arbitrary Lagrangian-Eulerian finite element method for solving three-dimensional free surface flows , 1998 .

[10]  Joachim Schöberl,et al.  NETGEN An advancing front 2D/3D-mesh generator based on abstract rules , 1997 .

[11]  A. Schmidt.  Computation of Three Dimensional Dendrites with Finite Elements , 1996 .

[12]  N. Ashgriz,et al.  A spine-flux method for simulating free surface flows , 1995 .

[13]  Philip H. Gaskell,et al.  Modelling and analysis of meniscus roll coating , 1995, Journal of Fluid Mechanics.

[14]  Thierry Coupez,et al.  Lagrangian finite‐element analysis of time‐dependent viscous free‐surface flow using an automatic remeshing technique: Application to metal casting flow , 1993 .

[15]  R. Mattheij,et al.  A boundary element solution for 2-dimensional viscous sintering , 1992 .

[16]  E. O. Tuck,et al.  Unsteady spreading of thin liquid films with small surface tension , 1991 .

[17]  H. Kuiken.  Viscous sintering: the surface-tension-driven flow of a liquid form under the influence of curvature gradients at its surface , 1990, Journal of Fluid Mechanics.

[18]  M. Kawahara,et al.  LAGRANGIAN FINITE ELEMENT ANALYSIS APPLIED TO VISCOUS FREE SURFACE FLUID FLOW , 1987 .

[19]  Philip M. Gresho,et al.  The implementation of normal and/or tangential boundary conditions in finite element codes for incompressible fluid flow , 1982 .

[20]  William G. Gray,et al.  Finite element simulation of flow in deforming regions , 1980 .

[21]  Daniel R. Lynch,et al.  Continuously deforming finite elements for the solution of parabolic problems, with and without phase change , 1981 .

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