Consistent, essentially conservative and balanced-force Phase-Field method to model incompressible two-phase flows

• Published in: Journal of computational physics 2020-04, Vol.406, p.109192, Article 109192
• Main Authors: Huang, Ziyang, Lin, Guang, Ardekani, Arezoo M.
• Format: Article
• Language: English
• Subjects: Algorithms; Balanced-force scheme; Computational fluid dynamics; Computational physics; Consistency; Consistent scheme; Density; Essentially conservative scheme; Fluid flow; Incompressible flow; Large density ratio; Momentum; Phase-Field; Surface tension; Two phase flow; Viscosity ratio
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2019.109192

•Three consistency conditions for the Cahn-Hilliard Phase-Field model are proposed.•The scheme is consistent & conservative for momentum transport at the discrete level.•The balanced-force algorithm is implemented to the surface force.•The proposed scheme is formally 2nd-order accurate in time and space.•The numerical solution converges to the sharp interface one. In the present work, the Cahn-Hilliard Phase-Field model of incompressible two-phase flows is considered. Conditions needed for consistency of reduction, consistency of mass and momentum transport, and consistency of mass conservation are proposed. The mass flux in the Navier-Stokes equations is defined such that it satisfies the proposed consistency conditions. The analysis in both continuous and discrete levels shows that violation of the consistency conditions result in unphysical solutions and the inconsistent errors are proportional to the density contrast of the fluids. After considering the conservative form of the inertial term, a consistent and conservative scheme for momentum transport is developed. The balanced-force algorithm for the sharp interface model is extended to the surface force derived from the Cahn-Hilliard model. The proposed scheme is formally 2nd-order accurate in both time and space, satisfies the consistency conditions, conserves mass globally and momentum essentially, and is balanced-force, in the discrete level. Its convergence to the sharp interface solution is systematically discussed in cases including large density and viscosity ratios, surface tension, and gravity. Various two-phase flow problems with large density ratios are performed to validate and verify the proposed scheme and excellent agreements with published numerical and/or experimental results are achieved. The proposed scheme is a practical and accurate tool to study two-phase flows, especially for those including large density ratios.