TY - JOUR
T1 - Implementation of two multiphase flow methods in modeling wetting of microporous hydrophobic membranes
AU - Yazgan-Birgi, P.
AU - Arafat, Hassan A.
AU - Hassan Ali, Mohamed I.
N1 - Funding Information:
This publication is based upon work supported by the Khalifa University of Science and Technology under Award No. RC2-2018-009. The authors would like to thank Prof. Gustaf Olsson from Lund University, Sweden, for providing his insight to improve the structure of this article.
Publisher Copyright:
© 2019 Elsevier B.V.
PY - 2019/11/15
Y1 - 2019/11/15
N2 - Pore wetting phenomenon plays a critical role in a porous media and is critical in various processes. For instance, liquid entry pressure (LEP) is one of the critical characteristics of hydrophobic membranes used in membrane distillation (MD) processes. In this study, pore-scale models were developed to assess the accuracy of two multiphase flow computational fluid dynamics (CFD) methods, as modeling tools for predicting two-phase flow in microporous MD membranes. Finite element method (FEM)-based phase field (PF) method (which was applied in the COMSOL package) and finite volume method (FVM)-based volume of fluid (VOF) method (which was applied in Star-CCM+) were the selected CFD tools for the implementation. The boundary conditions of the models were first set based on the experimental procedure for measuring the LEP, as given in the literature. Then, the models were used to capture the LEP under the gradually increased water pressure. Critical tuning of CFD parameters of each tool (such as mesh size, mesh type, and interface thickness) was conducted to investigate their influence on the LEP prediction accuracy and the water/air interface representation at the pore entrance. CFD model results were presented and compared with both experimental LEP data and the calculated value using the Young-Laplace equation (YLE). Both CFD tools were capable of capturing the water/air interface. LEP result from the VOF model showed good agreement with the experimental data, but the PF model overestimated the LEP value closer to the theoretical YLE value. For both approaches, the adjustment of the interface thickness was critical. In the VOF method, a realistic interface thickness could be achieved by adjusting both mesh size and time step simultaneously. In contrast, PF simulations were less mesh sensitive. The accuracy of the VOF model was better due to its mass conservation condition at the interface.
AB - Pore wetting phenomenon plays a critical role in a porous media and is critical in various processes. For instance, liquid entry pressure (LEP) is one of the critical characteristics of hydrophobic membranes used in membrane distillation (MD) processes. In this study, pore-scale models were developed to assess the accuracy of two multiphase flow computational fluid dynamics (CFD) methods, as modeling tools for predicting two-phase flow in microporous MD membranes. Finite element method (FEM)-based phase field (PF) method (which was applied in the COMSOL package) and finite volume method (FVM)-based volume of fluid (VOF) method (which was applied in Star-CCM+) were the selected CFD tools for the implementation. The boundary conditions of the models were first set based on the experimental procedure for measuring the LEP, as given in the literature. Then, the models were used to capture the LEP under the gradually increased water pressure. Critical tuning of CFD parameters of each tool (such as mesh size, mesh type, and interface thickness) was conducted to investigate their influence on the LEP prediction accuracy and the water/air interface representation at the pore entrance. CFD model results were presented and compared with both experimental LEP data and the calculated value using the Young-Laplace equation (YLE). Both CFD tools were capable of capturing the water/air interface. LEP result from the VOF model showed good agreement with the experimental data, but the PF model overestimated the LEP value closer to the theoretical YLE value. For both approaches, the adjustment of the interface thickness was critical. In the VOF method, a realistic interface thickness could be achieved by adjusting both mesh size and time step simultaneously. In contrast, PF simulations were less mesh sensitive. The accuracy of the VOF model was better due to its mass conservation condition at the interface.
KW - Computational fluid dynamics
KW - Liquid entry pressure
KW - Phase field
KW - Porous membrane
KW - Volume of fluid
KW - Wetting
UR - http://www.scopus.com/inward/record.url?scp=85069568595&partnerID=8YFLogxK
U2 - 10.1016/j.scitotenv.2019.07.232
DO - 10.1016/j.scitotenv.2019.07.232
M3 - Article
C2 - 31466205
AN - SCOPUS:85069568595
SN - 0048-9697
VL - 691
SP - 1251
EP - 1261
JO - Science of the Total Environment
JF - Science of the Total Environment
ER -