A Parallel Two-Phase Flow Solver on Unstructured Mesh in 3D

摘要

In this paper we introduce a parallel finite element method on 3D unstructured meshes for the two-phase flow problem modelled by a phase-field model consisting of the coupled Cahn-Hilliard and Navier-Stokes equations. A restricted additive Schwarz preconditioned GMRES method is used to solve the systems arising from implicit discretization of the Cahn-Hilliard equation and the velocity equation, and an algebraic multigrid preconditioned CG method is used to solve the pressure Poisson system. Numerical experiments suggest that the overall algorithm scales well on unstructured meshes for problems with up to 150 millions unknowns and on machines with close to 10,000 processors.