Preconditioned Krylov subspace methods used in solving two-dimensional transient two-phase flows

This paper investigates the performance of preconditioned Krylov subspace m ethods used in a previously presented two-fluid model developed for the; si mulation of separated and intermittent gas-liquid flows. The two-fluid mode l has momentum and mass balances for each phase. The equations comprising t his model are solved numerically by applying a two-step semi-implicit time integration procedure. A finite difference numerical scheme with a staggere d mesh is used. Previously, the resulting linear algebraic equations were s olved by a Gaussian band solver. In this study, these algebraic equations a re also solved using the generalized minimum residual (GMRES) and the bicon jugate gradient stabilized (Bi-CGSTAB) Krylov subspace iterative methods pr econditioned with incomplete LU factorization using the ILUT(p, tau) algori thm. The decrease in the computational time using the iterative solvers ins tead of the Gaussian band solver is shown to be considerable. Copyright (C) 1999 John Wiley & Sons, Ltd.