I have been waiting for you. At each step the solution get converged. When the normalized residuals of energy equations are below 1e-8 and and for other equation it is below 1e-4. Now I am having 30 iterations in each time step and the solutions get converged at each time step.

I am using phase change model of lee. It loops through each cell and finds cells having both liquid and vapor. As soon as the cell is retrieved mass is transferred. Frequency of evaporation is set now on a value of 100. Please let me know if there is anything needed. Please note that a udf is used for phase change implementation.