Divergence on RPI Model
 

Dear all,

I am modeling nucleate boiling on FLuent using the RPI model.

I can detect the vapor phase after some iterations (using steady state) but after that the calculation diverges (normally: turbulent viscosity limited to viscosity ratio 1e+5).

Also, using axisimmetric the vapor phase appears before the saturation temperature (30K before) and the divergence occurs only when Tsat is achieved).

There's a lot of correlations to phase iteration but I am following some papers.

Does anybody have any idea of how to solve this?
I tryed to reduce relaxation factors, change my mesh, 2D and 3D geometries, solution methods (it is more sateble with SIMPLE Presure-Velocity coupling, but everybody says to use Coupled, so I am trying).

Thaks!