Numerical and computational efficiency of solvers for two-phase problems

  • Authors:

  • O. Axelsson;P. Boyanova;M. Kronbichler;M. Neytcheva;X. Wu

  • Affiliations:

  • IG-AS, Studentska 1768, 708 00, Ostrava-Poruba, Czech Republic and King Abdulaziz University, Jeddah, Saudi Arabia;Uppsala University, Box 337, 751 05 Uppsala, Sweden and IICT-BAS, Acad. G. Bonchev Str., bl. 25A, 1113 Sofia, Bulgaria;Uppsala University, Box 337, 751 05 Uppsala, Sweden;Uppsala University, Box 337, 751 05 Uppsala, Sweden;Uppsala University, Box 337, 751 05 Uppsala, Sweden

  • Venue:
  • Computers & Mathematics with Applications
  • Year:
  • 2013

Quantified Score

Hi-index 0.09

Visualization

Abstract

We consider two-phase flow problems, modelled by the Cahn-Hilliard equation. In this work, the nonlinear fourth-order equation is decomposed into a system of two coupled second-order equations for the concentration and the chemical potential. We analyse solution methods based on an approximate two-by-two block factorization of the Jacobian of the nonlinear discrete problem. We propose a preconditioning technique that reduces the problem of solving the non-symmetric discrete Cahn-Hilliard system to a problem of solving systems with symmetric positive definite matrices where off-the-shelf multilevel and multigrid algorithms are directly applicable. The resulting solution methods exhibit optimal convergence and computational complexity properties and are suitable for parallel implementation. We illustrate the efficiency of the proposed methods by various numerical experiments, including parallel results for large scale three dimensional problems.