[seonguk]

Hello, I am working on the pressure increase and the Boil off gas by external heat input using the VOF model in a closed liquid hydrogen tank.
My case is as follows.
The rectangular 2-D axisymmetric model is filled with 50% liquid hydrogen and the heat comes in only on the wall of the liquid hydrogen side.
Multiphase : VOF
Energy : on
Viscous : Realize k-epsilon with enhanced wall treatment
Material
Gas : hydrogen(real-gas-peng-robinson)
Liquid : hydrogen liquid(Boussinesq)
Surface tension : constant
mass transfer : user-input
Cell zone condition :
Operating pressure : 101325pa, operating temperature : 20.268k, operating density : 0
boundary condition : heat wall : 50W/m^2
Method : SIMPLEC, pressure : body force weight, density,momentum,energy : second order
My problem is when the stream function reaches free surface, the residuals continuity goes above 1e-03 and then it causes the problem.
I would really appreciate it if you could tell me how to solve the problem.