Level set based parallel computations of unsteady free surface flows

The simulation of free surface ow is a computationally demanding task. Real life simulations require often the treatment of unstructured grids to capture the complex geometries and the use of parallel supercomputers for the ecient solution. For that aim, we have developed an ecient parallel free surface ow solver based on a full tetrahedral discretization. An implicit stabilized nite element method is used for solv- ing the unsteady incompressible two-phase ow in three-dimensions. A pressure stabilized Petrov Galerkin technique is used to avoid spurious pressure modes, while an upwind nite volume discretization is used to discretize the advective uxes in a stable manner. The interface between the uid phases is captured with the level set method implemented with a quadrature-free discontinuous Galerkin method. The parallel implementation is based on the MPI message-passing standard and is fully portable. We show the ee ctiveness of the method in the simulation of complex 3D ows, such as the ow past a cylindrical and the ow in a partially-l led tank of a car that suddenly brakes.

[1]  Björn Engquist,et al.  A finite element based level-set method for multiphase flow applications , 2000 .

[2]  T. Hughes,et al.  The Galerkin/least-squares method for advective-diffusive equations , 1988 .

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

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

[5]  Mark Sussman,et al.  An Efficient, Interface-Preserving Level Set Redistancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow , 1999, SIAM J. Sci. Comput..

[6]  Jean-François Remacle,et al.  A stabilized finite element method using a discontinuous level set approach for solving two phase incompressible flows , 2006, J. Comput. Phys..

[7]  Yannis Kallinderis,et al.  Finite element method for incompressible viscous flows with adaptive hybrid grids , 1994 .

[8]  Wang Hai-bing,et al.  High-order essentially non-oscillatory schemes for Hamilton-Jacobi equations , 2006 .

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

[10]  E. Puckett,et al.  Second-Order Accurate Volume-of-Fluid Algorithms for Tracking Material Interfaces , 2013 .

[11]  Charbel Farhat,et al.  A Minimum Overlap Restricted Additive Schwarz Preconditioner and Applications in 3D Flow Simulations , 1998 .

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

[13]  R. Ibrahim.  Liquid Sloshing Dynamics: Theory and Applications , 2005 .

[14]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[15]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[16]  Stuart E. Rogers,et al.  Steady and unsteady solutions of the incompressible Navier-Stokes equations , 1991 .

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

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

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

[20]  A. Smolianski.  Finite‐element/level‐set/operator‐splitting (FELSOS) approach for computing two‐fluid unsteady flows with free moving interfaces , 2005 .

[21]  W. Rumold.  Modeling and Simulation of Vehicles Carrying Liquid Cargo , 2001 .

[22]  Tayfun E. Tezduyar,et al.  Stabilized-finite-element/interface-capturing technique for parallel computation of unsteady flows with interfaces , 2000 .

[23]  J. Ferziger.  Interfacial transfer in Tryggvason's method , 2003 .

[24]  T. Barth.  Aspects of Unstructured Grids and Finite-Volume Solvers for the Euler and Navier-Stokes Equations , 1994 .

[25]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[26]  M. Rudman.  INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 671–691 (1997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS , 2022 .

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

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

[29]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[30]  Anna-Karin Tornberg,et al.  Interface tracking methods with application to multiphase flows , 2000 .

[31]  Philippe Geuzaine,et al.  Newton-Krylov Strategy for Compressible Turbulent Flows on Unstructured Meshes , 2001 .

[32]  J. Remacle,et al.  Efficient visualization of high‐order finite elements , 2007 .

[33]  A. Chorin.  A Numerical Method for Solving Incompressible Viscous Flow Problems , 1997 .

[34]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: V. Circumventing the Babuscka-Brezzi condition: A stable Petrov-Galerkin formulation of , 1986 .

[35]  S. Osher,et al.  High-order essentially nonsocillatory schemes for Hamilton-Jacobi equations , 1990 .