“…The use of such user defined variables has been previously used by the author [36,37] in studying evaporation of pure fluids. As these variables are solvable in the whole flow domain, in the first study [36], saturated vapour mass fraction was assigned in the liquid phase using the concept of ''internal" boundary condition [38].…”