“…However, further investigations using the initial form of the model report serious numerical difficulties related to the type of numerical schemes, time step size, and meshes methods used (uniform, non-uniform) in converging to a solution (Zhou, 2005), and also an inadequacy in satisfying the mass conservation (Zhang and Datta, 2004). Thus, an improved model was suggested by replacing the instantaneous phase change of liquid to vapour and by introducing an evaporation rate term (Huang et al, 2007). Other improvements were also made (Zhang et al, 2005(Zhang et al, , 2007Zhang and Datta, 2006;Hadiyanto et al, 2007Hadiyanto et al, , 2008 which could be categorized as equilibrium approach and include:…”