Implementations of source terms into enengy equation
 

Hello Foamer,

I try to incorporate the evaporation model into OpenFOAM based on the interFoam solver. The theory is based on ''Hardt, S., et al. (2008). J Comput Phys 227(11): 5871-5895.''. Actually, I checked the source code of ''evapVOFHardt'' solver which implements the evaporation model very well.

In the evapVOFHardt solver (TEqn.H), they calculated the latent heat term by decomposing -1*(mass source term)*(Enthalpy of vaporization) into hESource*T - hESource*Tsat, and calculate them as implicit and explicit terms respectively.

Actually, I know this method works, however, I tried to calculate the latent heat with a direct way, i.e. using the single explicit term ''-1*(mass source term)*(Enthalpy of vaporization) ''. But I got something strange there when I validate the solver with 1D stefan problem, and the problem is that I got the minimum temperature around 369.677 K (for the BCs in my case (Left wall temperature : 383.15 K, and the rest is 373.15K )), and the simulation was really slow although we can see the interface was moving.

Could you please tell me why I got the strange results? Thank you so much.

Hi,
How did you solve this problem? I am having a similar issue.