TY - JOUR
T1 - An efficient stabilized multiple auxiliary variables method for the Cahn–Hilliard–Darcy two-phase flow system
AU - Yang, Junxiang
AU - Kim, Junseok
N1 - Funding Information:
J. Yang is supported by China Scholarship Council ( 201908260060 ). The corresponding author (J.S. Kim) was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education ( NRF-2019R1A2C1003053 ). The authors thank the reviewers for the constructive and helpful comments on the revision of this article.
Publisher Copyright:
© 2021 Elsevier Ltd
PY - 2021/6/15
Y1 - 2021/6/15
N2 - Herein, we develop a totally decoupled, linear, and temporally second-order accurate numerical scheme for the Cahn–Hilliard–Darcy system which models the two-phase incompressible fluid flows in porous medium or in Hele–Shaw cells. Our proposed scheme is based on a simple and efficient stabilized multiple scalar auxiliary variables (S-MSAV) approach. Two time-dependent variables are defined to change the original governing equations to be the equivalent forms and then the Crank–Nicolson type approximation and the explicit Adams–Bashforth approximation are used to temporally discretize the equivalent equations. All nonlinear parts and auxiliary variables are explicitly treated, thus we can decouple the phase-field variable and auxiliary variables in time. The velocity and pressure are decoupled by using a second-order accurate pressure correction method. Therefore, our numerical scheme is very simple and efficient. We analytically prove the discrete energy dissipation law and unique solvability of our scheme with the absence of external force. The benchmark tests are performed to show that our method has desired accuracy and energy stability. Moreover, our method works well in simulating buoyancy-driven pinchoff and viscous fingering.
AB - Herein, we develop a totally decoupled, linear, and temporally second-order accurate numerical scheme for the Cahn–Hilliard–Darcy system which models the two-phase incompressible fluid flows in porous medium or in Hele–Shaw cells. Our proposed scheme is based on a simple and efficient stabilized multiple scalar auxiliary variables (S-MSAV) approach. Two time-dependent variables are defined to change the original governing equations to be the equivalent forms and then the Crank–Nicolson type approximation and the explicit Adams–Bashforth approximation are used to temporally discretize the equivalent equations. All nonlinear parts and auxiliary variables are explicitly treated, thus we can decouple the phase-field variable and auxiliary variables in time. The velocity and pressure are decoupled by using a second-order accurate pressure correction method. Therefore, our numerical scheme is very simple and efficient. We analytically prove the discrete energy dissipation law and unique solvability of our scheme with the absence of external force. The benchmark tests are performed to show that our method has desired accuracy and energy stability. Moreover, our method works well in simulating buoyancy-driven pinchoff and viscous fingering.
KW - Cahn–Hilliard–Darcy system
KW - Decoupled scheme
KW - Efficient S-MSAV approach
KW - Second-order accuracy
UR - http://www.scopus.com/inward/record.url?scp=85103977228&partnerID=8YFLogxK
U2 - 10.1016/j.compfluid.2021.104948
DO - 10.1016/j.compfluid.2021.104948
M3 - Article
AN - SCOPUS:85103977228
VL - 223
JO - Computers and Fluids
JF - Computers and Fluids
SN - 0045-7930
M1 - 104948
ER -