Numerical simulation of a single rising bubble by VOF with surface compression

The capability of the direct volume of fluid method for describing the surface dynamics of a free twodimensional rising bubble is evaluated using quantities of a recently published benchmark. The model equations are implemented in the open source computational fluid dynamics library OpenFOAM®. Here, a main ingredient of the numerical method is the so-called surface compression that corrects the fluxes near the interface between two phases. The application of this method with respect to two test cases of a benchmark is considered in the main part. The test cases differ in physical properties, thus in different surface tension effects. The quantities centre of mass position, circularity and rise velocity are tracked over time and compared with the ones given in the benchmark. For test case one, where surface tension effects are more pre-eminent, deviations from the benchmark results become more obvious. However, the flow features are still within reasonable range. Nevertheless, for test case two, which has higher density and viscosity ratios and above all a lower influence of the surface tension force, good agreement compared with the benchmark reference results is achieved. This paper demonstrates the good capabilities of the direct volume of fluid method with surface compression with regard to the preservation of sharp interfaces, boundedness, mass conservation and low computational time. Some limitation regarding the occurrence of parasitic currents, bad pressure jump prediction and bad grid convergence have been observed. With these restrictions in mind, the method is suitable for the simulation of similar two-phase flow configurations. Copyright © 2012 John Wiley & Sons, Ltd.

[1]  Hakan Serhad Soyhan,et al.  Critical evaluation of CFD codes for interfacial simulation of bubble‐train flow in a narrow channel , 2007 .

[2]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[3]  J. Crank,et al.  A practical method for numerical evaluation of solutions of partial differential equations of the heat-conduction type , 1947, Mathematical Proceedings of the Cambridge Philosophical Society.

[4]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

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

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

[7]  Tracie Barber,et al.  Assessment of Interface Capturing Methods in Computational Fluid Dynamics (CFD) Codes — A Case Study , 2009 .

[8]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[9]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[10]  O. Ubbink.  Numerical prediction of two fluid systems with sharp interfaces , 1997 .

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

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

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

[14]  M. Davidson,et al.  An analysis of parasitic current generation in Volume of Fluid simulations , 2005 .

[15]  S. Muzaferija.  Computation of free-surface flows using the finite-volume method and moving grids , 1997 .

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

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

[18]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

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

[20]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[21]  Stéphane Popinet,et al.  An accurate adaptive solver for surface-tension-driven interfacial flows , 2009, J. Comput. Phys..

[22]  M. Shashkov,et al.  Moment-of-fluid interface reconstruction , 2005 .

[23]  Borut Mavko,et al.  Simulations of free surface flows with implementation of surface tension and interface sharpening in the two-fluid model , 2009 .

[24]  Patricio Bohorquez,et al.  Study and Numerical Simulation of Sediment Transport in Free-Surface Flow , 2008 .

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

[26]  Mikhail J. Shashkov,et al.  Reconstruction of multi-material interfaces from moment data , 2008, J. Comput. Phys..

[27]  Bernard Laval.  Book review / Critique de livre : Numerical computation of internal and external flows: the fundamentals of computational fluid dynamics. 2nd ed. , 2008 .

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

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