 To ensure non-divergence of the flow field, an elliptic equation is solved.The elliptic equation is derived from the continuity equation wherein velocity components are expressed in terms of formula_7.Since the elliptic equation is derived from the discrete form of the continuity equation and the discrete form of the pressure gradient, conservativity is guaranteed (Flassak and Moussiopoulos, 1988).