Simulations of free surface flows with implementation of surface tension and interface sharpening in the two-fluid model

Two-fluid model used for free surface flows with large characteristic scales is improved; the smeared interface is sharpened with conservative level set method and the surface tension force with wetting angle is implemented. Surface tension force is split between two phases with several models. Detailed analysis showed the splitting of surface tension force with volume averaging as the most appropriate. The improved two-fluid model with interface sharpening and implemented surface tension is validated on several test cases. The pressure jump over a droplet interface test case showed that the pressure jump in simulation converges with grid refinement to the analytical one. The parasitic currents in simulation are one order of magnitude larger than in simulation with volume of fluid model. In the oscillating droplet test case the time period of oscillating droplet with initially ellipsoid or square shape is similar to the analytical time period. In the rising bubble test case, the rising bubble position, terminal velocity, and circularity are similar to the one observed in simulations with level set model. The wetting angle is implemented in the two-fluid model with interface sharpening and surface tension force. Model is tested in the simulation of droplet in contact with wall with different wetting angles.

[1]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

[2]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[3]  M. Ishii,et al.  Thermo-Fluid Dynamics of Two-Phase Flow , 2007 .

[4]  Stojan Petelin,et al.  Coupling of the interface tracking and the two-fluid models for the simulation of incompressible two-phase flow , 2001 .

[5]  Stéphane Vincent,et al.  TEST-CASE NO 10: PARASITIC CURRENTS INDUCED BY SURFACE TENSION (PC) , 2004 .

[6]  Eckhard Krepper,et al.  CFD simulation of convective flow boiling of refrigerant in a vertical annulus , 2008 .

[7]  D. Bestion,et al.  The physical closure laws in the CATHARE code , 1990 .

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

[9]  M. Boucker,et al.  A Two-Phase CFD Approach to the PTS Problem Evaluated on COSI Experiment , 2008 .

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

[11]  James M. Hyman,et al.  Numerical methods for tracking interfaces , 1984 .

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

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

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

[15]  Wei Yao,et al.  A Three-Dimensional Two-Fluid Modeling of Stratified Flow with Condensation for Pressurized Thermal Shock Investigations , 2005 .

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

[17]  Akihiko Minato,et al.  Simulation of Multi-dimensional Heterogeneous and Intermittent Two-Phase Flow by Using an Extended Two-Fluid Model , 2003 .

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

[19]  Iztok Tiselj,et al.  Two‐fluid model with interface sharpening , 2011 .

[20]  D. Kuzmin,et al.  Quantitative benchmark computations of two‐dimensional bubble dynamics , 2009 .

[21]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

[22]  Jean-Marie Seynhaeve,et al.  A first assessment of the NEPTUNE_CFD code: Instabilities in a stratified flow comparison between the VOF method and a two-field approach , 2008 .

[23]  Eckart Laurien,et al.  Experimental and numerical investigation of counter-current stratified flows in horizontal channels , 2008 .

[24]  Clayton T. Crowe,et al.  Multiphase Flow Handbook , 2005 .

[25]  Iztok Tiselj,et al.  Accuracy of the Operator Splitting Technique for Two-Phase Flow With Stiff Source Terms , 2002 .

[26]  Thomas Höhne,et al.  Experimental and numerical prediction of horizontal stratified flows , 2007 .

[27]  Washington.  RELAP5/MOD3 code manual. Volume 4, Models and correlations , 1995 .

[28]  Francesco Saverio D'Auria,et al.  Numerical Simulation of Free Surface Flows With Heat and Mass Transfer , 2007 .

[29]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .