interCondensatingEvaporatingFoam: interfaceHeatResistance model failing for condensation cases
I initialized a stefans problem case inside a rectangular domain with a small liquid film on the bottom wall and a vapour phase above the film. The film was set to a temperature 10K below the saturation temperature while the vapour was at saturation temperature. The phase change model used was the interfaceHeatResistance model. The top wall is hot and the bottom wall cold, gravity is off. Condensation should appear and the film should grow, while the evaporation massflux should be zero.
The condensation mass flux, which can be viewed in ParaView, seems to produce the correct value but there is also an evaporation mass flux, which is none zero. Also there is a strange behavior within the interface area calculation, which seems to randomly add parts of the walls. If i reversed the temperature setting and looked at an evaporation case the right behavior with growing vapour phase occurred.
Steps to reproduce
Run the attached testcase, look at the mass fluxes and alpha-field in ParaView, especially at the interface.
The example case is attached. Example_case.zip
What is the current bug behaviour?
The solver produces an evaporation mass flux which isn't set to zero for a condensation case (T < Tsat). For information: In the reversed case (evaporation) the condensation mass flux is set to zero. Also the solver uses an wrong interface area with parts of the side walls.
What is the expected correct behavior?
The evaporation mass flux should be set to zero in the case of a condensation case. And the liquid film should grow depending on the analytic solution for the Stefans problem. The interface area, which the solver uses, should only be the interface between the liquid and the vapour phase.