A VOF-Based Conservative Interpolation Scheme for Interface Tracking (CISIT) of Two-Fluid Flows

A simple interface-reconstruction scheme (CISIT), based on the volume-of-fluid (VOF) formulation and applicable to unstructured grids with arbitrary topology, is developed to track the interface in free-surface flows. The interface is represented by the contour surface of VOF value f = 0.5, which is constructed simply by interpolation. The advancing of the interface is made by imposing mass conservation in the cells containing the interface. However, problems arise when the interface moves across the cell face and into another cell. The cell originally occupied by the interface may be overfilled (f > 1) or underfilled (f < 1) as the interface advances (moving in the positive direction of the interface) to leave the cell. In contrast, the cell may be overdepleted (f < 0) or underdepleted (f > 0) as the interface retreats (moving in the negative direction of the interface) to move out of the cell. These situations are remedied via adjusting the transported volume of fluid across the cell face by following the conservation law so that the VOF becomes 1 or 0 in these cells, depending on advancing or retreating of the interface. The resulting VOF distribution is uniform, either in 1 or 0, in the region outside the interface cells, and the interface occupies only one cell in its width. To smooth the fluid properties across the interface region in the velocity calculation, an average smoothing technique is adopted. Application to a number of test cases, including model cases with known velocities and a number of real flows, reveal that this method is accurate and robust.

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

[2]  Farzad Mashayek,et al.  A hybrid finite‐element–volume‐of‐fluid method for simulating free surface flows and interfaces , 1995 .

[3]  W. Tao,et al.  A coupled volume-of-fluid and level set ( VOSET ) method for computing incompressible two-phase flows , 2009 .

[4]  Yeng-Yung Tsui,et al.  Use of Characteristic-Based Flux Limiters in a Pressure-Based Unstructured-Grid Algorithm Incorporating High-Resolution Schemes , 2009 .

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

[6]  Yeng-Yung Tsui,et al.  A Pressure-Based Unstructured-Grid Algorithm Using High-Resolution Schemes for All-Speed Flows , 2008 .

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

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

[9]  Nahmkeon Hur,et al.  A COUPLED LEVEL SET AND VOLUME-OF-FLUID METHOD FOR THE BUOYANCY-DRIVEN MOTION OF FLUID PARTICLES , 2002 .

[10]  L YoungsD,et al.  Time-dependent multi-material flow with large fluid distortion. , 1982 .

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

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

[13]  P. Davami,et al.  Numerical simulation of liquid/gas phase flow during mold filling , 2006 .

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

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

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

[17]  Farzad Mashayek,et al.  Advection of axisymmetric interfaces by the volume‐of‐fluid method , 1995 .

[18]  Nasser Ashgriz,et al.  FLAIR: fluz line-segment model for advection and interface reconstruction , 1991 .

[19]  P. Wesseling,et al.  A mass‐conserving Level‐Set method for modelling of multi‐phase flows , 2005 .

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

[21]  Yeng-Yung Tsui,et al.  A Pressure-Correction Method for Incompressible Flows Using Unstructured Meshes , 2006 .

[22]  Gang Wang,et al.  A computational Lagrangian–Eulerian advection remap for free surface flows , 2004 .

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

[24]  P. Woodward,et al.  SLIC (Simple Line Interface Calculation) , 1976 .

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

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

[27]  Suthee Wiri,et al.  On improving mass conservation of level set by reducing spatial discretization errors , 2005 .

[28]  C. K. Thornhill,et al.  Part IV. An experimental study of the collapse of liquid columns on a rigid horizontal plane , 1952, Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences.

[29]  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..

[30]  Marius Paraschivoiu,et al.  Second order accurate volume tracking based on remapping for triangular meshes , 2003 .

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

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

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

[34]  C. Peskin Numerical analysis of blood flow in the heart , 1977 .

[35]  Xiaofeng Yang,et al.  An adaptive coupled level-set/volume-of-fluid interface capturing method for unstructured triangular grids , 2006, J. Comput. Phys..

[36]  J. Monaghan Simulating Free Surface Flows with SPH , 1994 .

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

[38]  P. Gaskell,et al.  Curvature‐compensated convective transport: SMART, A new boundedness‐ preserving transport algorithm , 1988 .