Hello DrAmine and Surya,

thank you for this discussion.

I read the thesis and am trying to resimulate that for the first step, with the same options and even mesh.

I wonder if you could provide me with some details of your relevant simulation, to know my mistake.

1) what boundary conditions did you use, are the below correct?

     only the left was wall with temp=5+ T saturated,

     up and down boundaries were symmetric

     and the right side pressureoutlet with the T saturated liquid for the reverse flow.

2) what viscosities did you use for the phases, as in the table and other parts, the viscosities were not mentioned. only this table


which the standard enthalpy of vapor should be 10000 and liquid 0 as the molecular weight is 1.

3) For operating conditions,

where was your reference pressure location and what your Boussinesq parameters were.

4) on the finest grid used in the thesis, my setup was as below, was that correct?

  • VOF with implicit body force option,
  • first phase liquid and vapor as the second phase
  • the built-in Lee model with saturation temp=373.15, from the phase frequency of 0.1 and to phase frequency of 100
  • Simple as the velocity pressure coupling
  • Georeconstruct for interface capturing
  • time step of 1e-5s
  • second-order discretization for momentum and energy
  • node-based green gauss
  • PRESTO for pressure
  • in the material section, the reference temperature for the two phases 298.15
  • Surface tension on (CSF)/off
  • Gravity was on.

In the monitors, I got decreasing liquid mass and increasing vapour mass, but followed by a divergence after a while.

Thank you in advance.