<  Back to the Polytechnique Montréal portal

Numerical simulation of merging of two rising bubbles with different densities and diameters using an enhanced Volume-Of-Fluid (VOF) model

Faroogh Garoosi, Tarek Merabtene and Tew-Fik Mahdi

Article (2022)

[img] Restricted to: Repository staff only until 3 February 2024
Accepted Version
Terms of Use: Creative Commons Attribution Non-commercial No Derivatives
Request a copy
Show abstract
Hide abstract


In this study, the transient evolution of two rising bubbles with different densities is investigated numerically using an enhanced version of the VOF model, aiming to establish an state-of-the-art benchmark solutions and up-to-date data set for CFD validations. The simulations are performed on the staggered grid system where a novel third‐order accurate monotone convection scheme is applied for the discretization of the convection terms in Navier-Stokes and volume fraction equations while the semi-iterative PISOC algorithm (a combined version of the classical PISO and Chorin's model) are used to solve the pressure-velocity coupling. To reduce the false diffusion errors and mitigate smearing of interface thickness in the regions of physical discontinuities, the interface compression technique is also incorporated into the transport equation. To further enhance the accuracy of the numerical solutions, the idea of Piecewise Linear Interface Calculation (PLIC) based on the ELVIRA technique (Efficient Least-square Volume-of-fluid Interface Reconstruction Algorithm) is also utilized for the interface reconstruction and accurate implementation of surface tension force. The validity and accuracy of the enhanced VOF model is further demonstrated against a series of challenging benchmark cases including draining of liquids from the storage tank (tank draining), single rising bubble, three-phase Rayleigh-Taylor Instability and dam-break flows over dry and wet beds. The comparison of the obtained results with previously published data vividly demonstrates the superiority of the proposed method over the standard VOF/Level-Set models in handling multiphase/multi-fluid flow problems with large topological changes. In the last stage, the morphology and hydrodynamic characteristics of merging of two rising bubbles with different densities and diameters are examined and analyzed in details. The results show that, the initial/final deformations and the subsequent steady-state rise of two bubbles are remarkably influenced by the diameters of leading (upper) and trailing (lower) bubbles.

Uncontrolled Keywords

Monotone convection scheme; PISOC algorithm; Multiphase flows; VOF model; Merging of two rising bubbles

Department: Department of Civil, Geological and Mining Engineering
Grant number: RGPIN-2021-03272
PolyPublie URL: https://publications.polymtl.ca/10169/
Journal Title: Ocean Engineering (vol. 247)
Publisher: Elsevier
DOI: 10.1016/j.oceaneng.2022.110711
Official URL: https://doi.org/10.1016/j.oceaneng.2022.110711
Date Deposited: 09 Feb 2022 11:16
Last Modified: 09 Jun 2023 05:54
Cite in APA 7: Garoosi, F., Merabtene, T., & Mahdi, T.-F. (2022). Numerical simulation of merging of two rising bubbles with different densities and diameters using an enhanced Volume-Of-Fluid (VOF) model. Ocean Engineering, 247, 110711 (27 pages). https://doi.org/10.1016/j.oceaneng.2022.110711


Total downloads

Downloads per month in the last year

Origin of downloads


Repository Staff Only

View Item View Item