Eulerian method-Critical Heat Flux(CHF) Problem!
 

Hello everyone,

I want to estimate the critical heat flux of a two-phase flow in a flat mini channel (2mm x 2mm x 200mm). Mesh number is 234000. More details about the model setting are given below.
Model Settings:
Steady State
Eulerian Multiphase Model- Boiling Model - Critical Heat Flux (Primary Phase = Liquid Water; Secondary Phase = Water Vapor)

Realizable k-epsilon model with Enhanced Wall Function
Phase Interactions:Drag -ishii
Scheme: Phase Coupled Simple
Momentum-Volume Fraction,Turbulent Kinetic Energy,Dissipation Rate and Energy:QUICK

Velocity inlet 0.2 m/s apprx. G=200 kg/m2s

Heat Flux: I applied between 40-100 kW/m2

When I make a solution like this, I can not see the phase change. In addition, vf phase 2 and energy p2 values in the residual graph are always zero.http://hizliresim.com/nRbPG1

I have tried dividing the phase intervals like the VOF model by adapt. When I do this, the system gives the error of divergence.

what should I do?