Abstract

Discrete versions of Poisson's equation with large contrasts in the coefficients result in very ill-conditioned systems. Thus, its iterative solution represents a major challenge, for instance, in porous media and multiphase flow simulations, where considerable permeability and density ratios are usually found. The existing strategies trying to remedy this are highly dependent on whether the coefficient matrix remains constant at each time iteration or not. In this regard, incompressible multiphase flows with high-density ratios are particularly demanding as their resulting Poisson equation varies along with the density field, making the reconstruction of complex preconditioners impractical. This work presents a strategy for solving such versions of the variable Poisson equation.Roughly, we first make it constant through an adequate approximation. Then, we block-diagonalise it through an inexpensive change of basis that takes advantage of mesh reflection symmetries, which are common in multiphase flows. Finally, we solve the resulting set of fully decoupled subsystems with virtually any solver. The numerical experiments conducted on a multiphase flow simulation prove the benefits of such an approach, resulting in up to 6.6x faster convergences.