[MABH]

I am trying to implement the Lee model (used for phase changes in Fluent) in the solver interPhaseChangeFoam in Open Foam (version 3.0.1).

The interPhaseChangeFoam solver includes phase change models where the phase change depends on pressure differences (specially for cavitation simuations). I have added the Lee mode (in a new file in the file PhaseChangeTwoPhaseMixtures) in order to simulate a phase change wich depends on temperature (saturation temperature more precisely).

I have already added the Lee model in the PhaseChangeTwoPhaseMixtures file but I don't know how to impose that phase change depending on temperature in the PhaseChangeTwoPhaseMixtures.C solver.

I have already tried to change the pressure term by a temperature term but I have an error, apparently I have to initialize the T field before, how can i do that?