Stabilized edge‐based finite element simulation of free‐surface flows

Free-surface flows occur in several problems in hydrodynamics, such as fuel or water sloshing in tanks, waves breaking in ships, offshore platforms, harbours and coastal areas. The computation of such highly nonlinear flows is challenging since free-surfaces commonly present merging, fragmentation and breaking parts, leading to the use of interface-capturing Eulerian approaches. In such methods the surface between two fluids is captured by the use of a marking function which is transported in a flow field. In this work we present a three-dimensional parallel edge-based incompressible SUPG/PSPG finite element method to cope with free-surface problems with volume-of-fluid (VOF) extensions to track the evolving free surface. The pure advection equation for the scalar marking function was solved by a fully implicit parallel edge-based SUPG finite element formulation. We studied variants of this formulation, considering the effects of discontinuity capturing and a particular tangent transformation designed to increase interface sharpness. Global mass conservation is enforced adding or removing mass proportionally to the absolute value of the normal velocity of the interface. We introduce a parallel dynamic deactivation algorithm to solve the marking function equation only in a small region around the interface. The implementation is targeted to distributed memory systems with cache-based processors. The performance and accuracy of the proposed solution method were tested with several validation problems. Copyright © 2007 John Wiley & Sons, Ltd.

[1]  S. Mittal,et al.  Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements , 1992 .

[2]  Alvaro L. G. A. Coutinho,et al.  Inexact Newton-type methods for the solution of steady incompressible viscoplastic flows with the SUPG/PSPG finite element formulation , 2006 .

[3]  Tayfun E. Tezduyar,et al.  Moving‐interface computations with the edge‐tracked interface locator technique (ETILT) , 2005 .

[4]  Jean-François Remacle,et al.  A quadrature-free discontinuous Galerkin method for the level set equation , 2006, J. Comput. Phys..

[5]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

[6]  Rainald Löhner,et al.  Dynamic deactivation for advection-dominated contaminant transport , 2004 .

[7]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[8]  Eugenio Oñate,et al.  The ALE/Lagrangian Particle Finite Element Method: A new approach to computation of free-surface flows and fluid–object interactions , 2007 .

[9]  Damien Violeau,et al.  Numerical modelling of complex turbulent free‐surface flows with the SPH method: an overview , 2007 .

[10]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

[11]  Alvaro L. G. A. Coutinho,et al.  EdgePack: A Parallel Vertex and Node Reordering Package for Optimizing Edge-Based Computations in Unstructured Grids , 2006, VECPAR.

[12]  S. Koshizuka.  A particle method for incompressible viscous flow with fluid fragmentation , 1995 .

[13]  Tayfan E. Tezduyar,et al.  Stabilized Finite Element Formulations for Incompressible Flow Computations , 1991 .

[14]  Tayfun E. Tezduyar,et al.  Enhanced-discretization Selective Stabilization Procedure (EDSSP) , 2006 .

[15]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[16]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[17]  Brian L. Smith,et al.  New finite-element/finite-volume level set formulation for modelling two-phase incompressible flows , 2006, J. Comput. Phys..

[18]  Deborah Greaves,et al.  Simulation of viscous water column collapse using adapting hierarchical grids , 2006 .

[19]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[20]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[21]  S. Mittal,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. II: Computation of free-surface flows, two-liquid flows, and flows with drifting cylinders , 1992 .

[22]  Marek Behr,et al.  Enhanced-Discretization Interface-Capturing Technique (EDICT) for computation of unsteady flows with interfaces , 1998 .

[23]  Alex M. Andrew,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid Mechanics, Computer Vision, and Materials Science (2nd edition) , 2000 .

[24]  Tayfun E. Tezduyar,et al.  Adaptive implicit-explicit finite element algorithms for fluid mechanics problems , 1990 .

[25]  T. Tezduyar.  Computation of moving boundaries and interfaces and stabilization parameters , 2003 .

[26]  T. Tezduyar,et al.  A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests , 1992 .

[27]  Alvaro L. G. A. Coutinho,et al.  Performance comparison of data‐reordering algorithms for sparse matrix–vector multiplication in edge‐based unstructured grid computations , 2006 .

[28]  Alvaro L. G. A. Coutinho,et al.  Parallel edge-based solution of viscoplastic flows with the SUPG/PSPG formulation , 2006 .

[29]  Alvaro L. G. A. Coutinho,et al.  Control strategies for timestep selection in finite element simulation of incompressible flows and coupled reaction–convection–diffusion processes , 2005 .

[30]  J. C. Martin.  An experimental study of the collapse of liquid column on a rigid horizontal plane , 1952 .

[31]  Feng Xiao,et al.  Description of complex and sharp interface with fixed grids in incompressible and compressible fluid , 1995 .

[32]  Arthur Veldman,et al.  A Volume-of-Fluid based simulation method for wave impact problems , 2005 .

[33]  T. Tezduyar,et al.  Computation of mould filling processes with a moving Lagrangian interface technique , 2002 .

[34]  Tayfun E. Tezduyar,et al.  Computation of free-surface flows and fluid–object interactions with the CIP method based on adaptive meshless soroban grids , 2007 .

[35]  Alvaro L. G. A. Coutinho,et al.  Parallel Edge-Based Inexact Newton Solution of Steady Incompressible 3D Navier-Stokes Equations , 2005, Euro-Par.

[36]  S. Osher,et al.  An improved level set method for incompressible two-phase flows , 1998 .

[37]  A. Galeão,et al.  A consistent approximate upwind Petrov—Galerkin method for convection-dominated problems , 1988 .

[38]  Tayfun E. Tezduyar,et al.  Finite element methods for flow problems with moving boundaries and interfaces , 2001 .

[39]  Catherine Donati-Martin,et al.  A note on the change of variable formula with local time-space for bounded variation Levy processes , 2007 .

[40]  Tayfun E. Tezduyar,et al.  Finite element stabilization parameters computed from element matrices and vectors , 2000 .

[41]  Alvaro L. G. A. Coutinho,et al.  Edge‐based adaptive implicit/explicit finite element procedures for three‐dimensional transport problems , 2005 .

[42]  Kornelia Marchien Theresa Kleefsman.  Water impact loading on offshore structures. - A numerical study , 2005 .

[43]  Tayfun E. Tezduyar,et al.  CFD methods for three-dimensional computation of complex flow problems , 1999 .

[44]  Rainald Löhner,et al.  On the simulation of flows with violent free surface motion , 2006 .

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

[46]  J. Smagorinsky,et al.  GENERAL CIRCULATION EXPERIMENTS WITH THE PRIMITIVE EQUATIONS , 1963 .

[47]  R. Lahey,et al.  Computation of incompressible bubble dynamics with a stabilized finite element level set method , 2005 .

[48]  Tayfun E. Tezduyar,et al.  A moving Lagrangian interface technique for flow computations over fixed meshes , 2001 .