To solve the asymmetric systems of linear equations that arise from approximations of momentum, temperature and species equations Equation 1 (see “Form of the Numerical Algorithm”), a preconditioned generalized conjugate gradient method (Ref. 16) is used.
Incomplete LU factorization is used for preconditioning.