Equations (12)-(14) were discretized using thecontrol volume formulation developed by Patankar[4O]. The power law scheme was used to calculate boththe heat and mass fluxes across the boundaries of eachcontrol volume. The simultaneous algebraic equationswere solved using an iterative point-by-point methodand relaxation parameters. The process was repeateduntil the changes in $.,,, T+ and C, between twosuccessive iterations satisfied the criterion