TY - JOUR

T1 - A pressure correction-volume of fluid method for simulation of two-phase flows

AU - Lin, San Yih

AU - Chin, Ya Hsien

AU - Wu, Chun Muh

AU - Lin, Jeng Feng

AU - Chen, Yi Cheng

PY - 2012/1/20

Y1 - 2012/1/20

N2 - A pressure correction method coupled with the volume of fluid (VOF) method is developed to simulate two-phase flows. A volume fraction function is introduced in the VOF method and is governed by an advection equation. A modified monotone upwind scheme for a conservation law (modified MUSCL) is used to solve the solution of the advection equation. To keep the initial sharpness of an interface, a slope modification scheme is introduced. The continuum surface tension (CST) model is used to calculate the surface tension force. Three schemes, central-upwind, Parker-Youngs, and mixed schemes, are introduced to compute the interface normal vector and the gradient of the volume fraction function. Moreover, a height function technique is applied to compute the local curvature of the interface. Several basic test problems are performed to check the order of accuracy of the present numerical schemes for computing the interface normal vector and the gradient of the volume fraction function. Three physical problems, two-dimensional broken dam problem, static drop, and spurious currents, and three-dimensional rising bubble, are performed to demonstrate the efficiency and accuracy of the pressure correction method.

AB - A pressure correction method coupled with the volume of fluid (VOF) method is developed to simulate two-phase flows. A volume fraction function is introduced in the VOF method and is governed by an advection equation. A modified monotone upwind scheme for a conservation law (modified MUSCL) is used to solve the solution of the advection equation. To keep the initial sharpness of an interface, a slope modification scheme is introduced. The continuum surface tension (CST) model is used to calculate the surface tension force. Three schemes, central-upwind, Parker-Youngs, and mixed schemes, are introduced to compute the interface normal vector and the gradient of the volume fraction function. Moreover, a height function technique is applied to compute the local curvature of the interface. Several basic test problems are performed to check the order of accuracy of the present numerical schemes for computing the interface normal vector and the gradient of the volume fraction function. Three physical problems, two-dimensional broken dam problem, static drop, and spurious currents, and three-dimensional rising bubble, are performed to demonstrate the efficiency and accuracy of the pressure correction method.

UR - http://www.scopus.com/inward/record.url?scp=84255178610&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84255178610&partnerID=8YFLogxK

U2 - 10.1002/fld.2500

DO - 10.1002/fld.2500

M3 - Article

AN - SCOPUS:84255178610

VL - 68

SP - 181

EP - 195

JO - International Journal for Numerical Methods in Fluids

JF - International Journal for Numerical Methods in Fluids

SN - 0271-2091

IS - 2

ER -