A refined volume-of-fluid algorithm for capturing sharp fluid interfaces on arbitrary meshes

This paper presents a new volume-of-fluid scheme (M-CICSAM), capable of capturing abrupt interfaces on meshes of arbitrary topology, which is a modification to the compressive interface capturing scheme for arbitrary meshes (CICSAM) proposed in the recent literature. Without resort to any explicit interface reconstruction, M-CICSAM is able to precisely model the complex free surface deformation, such as interface rupture and coalescence. By theoretical analysis, it is shown that the modified CICSAM overcomes three inherent drawbacks of the original CICSAM, concerning the basic differencing schemes, the switching strategy between the compressive downwind and diffusive high-resolution schemes, and the far-upwind reconstruction technique on arbitrary unstructured meshes. To evaluate the performance of the newly proposed scheme, several classic interface capturing methods developed in the past decades are compared with M-CICSAM in four test problems. The numerical results clearly demonstrate that M-CICSAM produces more accurate predictions on arbitrary meshes, especially at high Courant numbers, by reducing the numerical diffusion and preserving the interface shape.

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

[2]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[3]  Robert L. Street,et al.  A computer study of finite-amplitude water waves , 1970 .

[4]  Feng Xiao,et al.  An efficient method for capturing free boundaries in multi‐fluid simulations , 2003 .

[5]  D. Juric,et al.  A front-tracking method for the computations of multiphase flow , 2001 .

[6]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[7]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

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

[9]  Roozbeh Panahi,et al.  Numerical simulation of three‐dimensional interfacial flows , 2007 .

[10]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[11]  D. B. Kothe,et al.  RIPPLE: A NEW MODEL FOR INCOMPRESSIBLE FLOWS WITH FREE SURFACES , 1991 .

[12]  H. Rusche.  Computational fluid dynamics of dispersed two-phase flows at high phase fractions , 2003 .

[13]  Arnaud G. Malan,et al.  Development of a compressive surface capturing formulation for modelling free‐surface flow by using the volume‐of‐fluid approach , 2013 .

[14]  M. Darwish,et al.  Convective Schemes for Capturing Interfaces of Free-Surface Flows on Unstructured Grids , 2006 .

[15]  Chih Hao Chang,et al.  The capturing of free surfaces in incompressible multi-fluid flows , 2000 .

[16]  Roozbeh Panahi,et al.  COMPARISON OF INTERFACE CAPTURING METHODS IN TWO PHASE FLOW , 2005 .

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

[18]  D. K. Walters,et al.  A new interface‐capturing discretization scheme for numerical solution of the volume fraction equation in two‐phase flows , 2009 .

[19]  M. Rudman.  A Volume-Tracking Method for Incompressible Multifluid Flows with Large Density Variations , 1998 .

[20]  M. Sussman,et al.  A Coupled Level Set and Volume-of-Fluid Method for Computing 3D and Axisymmetric Incompressible Two-Phase Flows , 2000 .

[21]  Yeng-Yung Tsui,et al.  Flux-blending schemes for interface capture in two-fluid flows , 2009 .

[22]  Berend van Wachem,et al.  Volume of fluid methods for immiscible-fluid and free-surface flows , 2008 .

[23]  Cornelis Vuik,et al.  A mass conserving level set (MCLS) method for modeling of multi-phase flows , 2003 .

[24]  Chen Yang,et al.  Assessment of different reconstruction techniques for implementing the NVSF schemes on unstructured meshes , 2014 .

[25]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[26]  David R. Emerson,et al.  An Implicit Algorithm for Capturing Sharp Fluid Interfaces in the Volume of Fluid Advection Method , 2006 .

[27]  T. Nakayama,et al.  AN EULERIAN FINITE ELEMENT METHOD FOR TIME‐DEPENDENT FREE SURFACE PROBLEMS IN HYDRODYNAMICS , 1996 .

[28]  Koulis Pericleous,et al.  FREE SURFACE FLOW AND HEAT TRANSFER IN CAVITIES: THE SEA ALGORITHM , 1995 .

[29]  Kensuke Yokoi,et al.  Efficient implementation of THINC scheme: A simple and practical smoothed VOF algorithm , 2007, J. Comput. Phys..

[30]  R. I. Issa,et al.  A Method for Capturing Sharp Fluid Interfaces on Arbitrary Meshes , 1999 .

[31]  Richard H. Pletcher,et al.  The Development of a Free Surface Capturing Approach for Multidimensional Free Surface Flows in Closed Containers , 1997 .

[32]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[33]  A. Gosman,et al.  High resolution NVD differencing scheme for arbitrarily unstructured meshes , 1999 .

[34]  M. Kawahara,et al.  TWO-DIMENSIONAL SLOSHING ANALYSIS BY LAGRANGIAN FINITE ELEMENT METHOD , 1990 .

[35]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .

[36]  Brian L. Smith,et al.  A novel technique for including surface tension in PLIC-VOF methods , 2002 .

[37]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[38]  J. Thuburn.  Multidimensional Flux-Limited Advection Schemes , 1996 .

[39]  Ling Qian,et al.  A free-surface capturing method for two fluid flows with moving bodies , 2006, Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[40]  Woo Il Lee,et al.  A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part I: New free surface‐tracking algorithm and its verification , 2003 .

[41]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

[42]  Jacques Magnaudet,et al.  An interface-capturing method for incompressible two-phase flows. Validation and application to bubble dynamics , 2007 .

[43]  W. B. VanderHeyden,et al.  A general-purpose finite-volume advection scheme for continuous and discontinuous fields on unstructured grids , 2002 .

[44]  Luis Gustavo Nonato,et al.  A front-tracking/front-capturing method for the simulation of 3D multi-fluid flows with free surfaces , 2004 .

[45]  B. P. Leonard,et al.  The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection , 1991 .

[46]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[47]  M. Darwish,et al.  NORMALIZED VARIABLE AND SPACE FORMULATION METHODOLOGY FOR HIGH-RESOLUTION SCHEMES , 1994 .