We consider preconditioned iterative solution methods to solve the algebraic systems of equations arising from finite element discretizations of multiphase flow problems, based on the phase-field model.The aim is to solve coupled physics problems, where both diffusive and convective processes take place simultaneously in time and space. To model the above, a coupled system of partial differential equations has to be solved, consisting of the Cahn-Hilliard equation to describe the diffusive interface and the time-dependent Navier-Stokes equation, to follow the evolution of the convection field in time.We focus on the construction and efficiency of preconditioned iterative solution methods for the linear systems, arising after conforming and non-conforming finite element discretizations of the Cahn-Hilliard equation in space and implicit discretization schemes in time. The non-linearity of the phase-separation process is treated by Newton's method. The resulting matrices admit a two-by-two block structure, utilized by the preconditioning techniques, proposed in the current work. We discuss approximation estimates of the preconditioners and include numerical experiments to illustrate their behaviour.
展开▼
