[Rahmani_Zakaria]

Hi, I am simulating a 3D transient multiphase flow using the mixture model with 3 phases and the evaporation-condensation model (Lee model). After few iterations the calculation diverges and show this message:
temperature limited to 1.000000e+00 in 80 cells on zone 4 in domain 1
temperature limited to 1.000000e+00 in 313 cells on zone 4 in domain 1
Stabilizing temperature to enhance linear solver robustness.
temperature limited to 5.000000e+03 in 377818 cells on zone 4 in domain 1
Stabilizing y-momentum to enhance linear solver robustness.
Stabilizing z-momentum to enhance linear solver robustness.
Stabilizing temperature to enhance linear solver robustness.
temperature limited to 1.000000e+00 in 756903 cells on zone 4 in domain 1
temperature limited to 1.000000e+00 on 338 faces
temperature limited to 5.000000e+03 on 6 faces
Stabilizing x-momentum to enhance linear solver robustness.
Stabilizing y-momentum to enhance linear solver robustness.
Stabilizing z-momentum to enhance linear solver robustness.
Stabilizing pressure correction to enhance linear solver robustness.
Stabilizing temperature to enhance linear solver robustness.
temperature limited to 1.000000e+00 in 1151916 cells on zone 4 in domain 1
temperature limited to 5.000000e+03 in 18 cells on zone 4 in domain 1
temperature limited to 1.000000e+00 on 133 faces

reversed flow in 346 faces on outflow 7.
Stabilizing x-momentum to enhance linear solver robustness.
Stabilizing y-momentum to enhance linear solver robustness.
Stabilizing z-momentum to enhance linear solver robustness.
Stabilizing temperature to enhance linear solver robustness.
temperature limited to 5.000000e+03 in 18 cells on zone 4 in domain 1
temperature limited to 1.000000e+00 on 39 faces

reversed flow in 262 faces on outflow 7.
Stabilizing x-momentum to enhance linear solver robustness.
Stabilizing y-momentum to enhance linear solver robustness.
Stabilizing z-momentum to enhance linear solver robustness.
Stabilizing temperature to enhance linear solver robustness.
Stabilizing temperature using GMRES to enhance linear solver robustness.

Divergence detected in AMG solver: temperature
Error at host: floating point exception