I need to add a latent heat flux of evaporation to the hot side boundary, which is a function of the pressure difference between the cold side and the hot side boundary.However, the pressure at the cold side boundary is a function of temperature. How does it apply to the hold side? How can I solve this problem?