Equations (12)-(14) were discretized using the control volume formulation developed by Patankar . Calculate both the heat and mass fluxes across the boundaries of each control volume. Algebraic equations were solved using an iterative point-by-point method and relaxation parameters. Residual