I created 3 planes and 3 surfaces (each case operation is add frozen) and used boolean tool to substract surface 2 from 1,and surface 3 from 1.To apply heat flux I hid surfaces 1 and 2, and applied heat flux at the bottom edge(wall) of surface 3.After 2.389s of simulation the volume fraction of water in the pool is like right side image. However , water in the pool is at saturated state( temp 26.7┬░C and operating pressure 33048Pa) and I am applying heat flux of 120000W/m2 at the bottom of surface 3(copper). From my experimental experience ,the surface should be covered by bubbles/vapor and there should be many bubbles in the system. Do you think water fraction at right side image is still logical? Or, any mistake in process.
Moreover, what might be the process to determine top edge temperature of surface 3? If I know this then I can calculate wall superheat.
During initialization, I selected compute from 'wall_heat flux(the bottom edge of surface_3).
My residual is like below:
iter continuity u-water u-vapor v-water v-vapor energy-p1 energy-p2 k epsilon iac vf-vapor time/iter
47729 1.0868e-06 1.4956e-05 8.6502e-07 1.9258e-05 1.1558e-06 2.0921e-10 3.6352e-19 2.7958e-05 3.6626e-05 6.1474e-07 8.1657e-07 0:00:19 20

Regards Anadi