interCondensatingEvaporatingFoam with InterfaceHeatResistance model does not adequately evaluate the volume change due to evaporation
I am testing interCondensatingEvaporatingFoam of openFoam v2312. The openFoam binary is obtained from software.opensuse.org and installed on the PC of Intel cpu. The OS is openSuse Leap 15.5. The test problem is as follows: An open vessel half filled with water is heated from below. A bubble is installed in the water. The bubble will grow as the water evaporates. This should cause the water surface to rise.testCase.zip
When Constant model is employed as the phaseChangeTwoPhaseModel in constant/phaseChangeProperties, the water surface rises with the growth of bubble, as shown in Fig.1.
Fig.1: Twall=Tsat+2K, phaseChangeTwoPhaseModel=constant, coeffC=100, coeffE=100
However, when Interface Heat Resistance model is employed, the water surface does not rise even as the bubble grow, as shown in Fig.2. (No boiling is occurring near both ends of heated bottom wall, but that's not the issue here. The boiling on heated wall can be easily reproduced numerically by putting thin vapor film on the wall initially.)
Fig.2: Twall=Tsat+2K, phaseChangeTwoPhaseModel=InterfaceHeatResistance, R=1e6, maxAlphaRate=1, spread=3.
These results suggest that the interCondensatingEvaporatingFoam with InterfaceHeatResistance model does not adequately evaluate the volume change due to evaporation.