3 phases vof evaporation model
 

Hello everybody. I'm trying to run a cfd simulation of evaporation of water in water/ mercury mixture with air as primary phase , water as 1st secondary phase, and mercury as 2nd secondary phase, but cant see any change as if therse no heat source applied, Could anyone tell me what change to make for the three phase evapotion problem.
Thanks for any help!


DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) /*energy source udf*/
{
real enrg;
Thread *pri_th, *sec_th;
pri_th = THREAD_SUB_THREAD(mix_th,0);
sec_th = THREAD_SUB_THREAD(mix_th,1);
enrg = 0.0;
if (C_T(cell,sec_th)>T_SAT && C_VOF(cell, sec_th) != 0.0 && C_VOF(cell, sec_th) != 1.0 )
{
enrg = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*LAT_HT*(C_ T(cell,sec_th)-T_SAT)/T_SAT;
dS[eqn] = -0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*LAT_HT/T_SAT;
}
if (C_T(cell,pri_th)<T_SAT && C_VOF(cell, sec_th) != 0.0 && C_VOF(cell, sec_th) != 1.0 )
{
enrg = 0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*LAT_HT*(T_ SAT-C_T(cell,pri_th))/T_SAT;
dS[eqn] = -0.1*C_VOF(cell,pri_th)*C_R(cell, pri_th)*LAT_HT/T_SAT;
}
return enrg;
}