<  Back to the Polytechnique Montréal portal

Numerical simulation of three-fluid Rayleigh-Taylor instability using an enhanced Volume-Of-Fluid (VOF) model: New benchmark solutions

Faroogh Garoosi and Tew-Fik Mahdi

Article (2022)

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


The main objective of the present study is to introduce two novel benchmark solutions namely: two-dimensional three-fluid Rayleigh-Taylor Instability problems, aiming to provide an up-to-date data set and a unique fundamental insight into morphology and hydrodynamic behavior of coupled Rayleigh-Taylor-Kelvin-Helmholtz instability phenomenon. To this end, the Volume-Of-Fluid (VOF) model is adopted to probe the complex configurations and kinetic processes of highly nonlinear multi-fluid flow problems with large topological changes and moving interfaces. However, to improve the performance and accuracy of the classical VOF model and preserve monotonicity for the density and viscosity, a novel high-order bounded advection scheme is first proposed in the context of the Total Variation Diminishing and Normalized Variable Diagram (TVD-NVD) constraints and then is utilized for the discretization of the convection terms in the Navier-Stokes and transport equations. To further increase the accuracy of the numerical simulations, the second-order PLIC-ELVIRA is implemented for the reconstruction of the physical discontinuity between phases and the determination of its curvature. Furthermore, to enhance the consistency and stability of the classical VOF model in handling incompressible multi-fluid flows, a novel semi-iterative pressure-velocity coupling algorithm is constructed by the combination of the standard PISO and SIMPLEC algorithms and is then applied to ensure mass conservation in each grid cell. To demonstrate the versatility and robustness of the proposed model in dealing with the multiphase flows involving large interface deformation and breaking phenomena, a series of canonical test cases such as dam-break over a dry bed with and without stationary obstacle, 2D three-fluid rising bubble, two-fluid and three-fluid Rayleigh-Taylor Instability are adopted. In the last stage, an improved VOF model is applied to solve two new three-fluid Rayleigh-Taylor Instability benchmark problems on a staggered grid system. The results of this study can provide a wide panorama on the improvements of standard VOF model and may be utilized as benchmark solutions for validation of various CFD tools or simply to understand more complex related multi-fluid flows.

Uncontrolled Keywords

High-order TVD-NVD advection scheme; Hybrid PISO-SIMPLEC algorithm; Multi-fluid flows; PLIC-ELVIRA technique; Improved VOF model; Rayleigh-Taylor instability

Subjects: 1000 Civil engineering > 1005 Hydraulic engineering
2950 Applied mathematics > 2955 Numerical analysis
2950 Applied mathematics > 2960 Mathematical modelling
Department: Department of Civil, Geological and Mining Engineering
Grant number: RGPIN-2021-03272
PolyPublie URL: https://publications.polymtl.ca/10411/
Journal Title: Computers & Fluids (vol. 245)
Publisher: Elsevier
DOI: 10.1016/j.compfluid.2022.105591
Official URL: https://doi.org/10.1016/j.compfluid.2022.105591
Date Deposited: 01 Aug 2022 11:41
Last Modified: 09 Jun 2023 16:39
Cite in APA 7: Garoosi, F., & Mahdi, T.-F. (2022). Numerical simulation of three-fluid Rayleigh-Taylor instability using an enhanced Volume-Of-Fluid (VOF) model: New benchmark solutions. Computers & Fluids, 245, 105591 (33 pages). https://doi.org/10.1016/j.compfluid.2022.105591


Total downloads

Downloads per month in the last year

Origin of downloads


Repository Staff Only

View Item View Item