How to include the recoil pressure caused by vapourization in a model
 

Hi All,

I have a small question to ask the gurus here. Right now I would like to model the evaporation process of metal liquid, generated in a capillary. As a result, this evaporation causes high recoil pressure inside the capillary. I would like to include the effect of this pressure in my model. According to my research, some specify this known pressure (function of surface temperature) as the Dirichtlet value on the related evaporation interface, some include it through the interface condition [P]interface = [(surface tension) + (recoil pressure)]*dirac_delta and some include it through the continuity equation div[u] = (mass_flow_rate/density)*dirac_delta. The dirac_delta is the derivative of the smoothed Heaviside function with respect to the distance function (Level set method).

Which method would be the best? Could there be other better solutions?

Thank you in advance
Komkamol