interCondensatingEvaporatingFoam: wrong temperature field inside vapour in phase-transition region near interface
Summary
I initialized a vapor-bubble inside a rectangular liquid domain and set the heat conduction in the vapor to 0 (in order to visualize the problem). The condensation and evaporation coefficients are set to 0. The top wall is hot and the bottom wall cold, gravity is off. A temperature gradient emerges in the liquid as expected.
The wrong behavior occurs in the vapor-bubble near the interface, in the alpha-transition region, where 0 < alpha < 1. The temperature gradient that should only be present in the liquid phase also travels into this transition-region inside the vapor-phase (note that it is not smeared according to alpha in this region, but continued from the liquid) I discussed this problem with colleagues and we couldn't find a physical reason for this behavior. It seams that the temperature-field inside the vapor is only correct if alpha is exactly 0 (or at least smaller than 1e-8).
Steps to reproduce
Run the attached testcase (Allrun), look at the alpha- and temperature-field in paraview and compare their size.
Example case
What is the current bug behaviour?
The temperature-field of the surrounding field is present in the alpha-transition region inside the bubble near the interface.
What is the expected correct behavior?
The temperature should be constant everywhere in the bubble, due to heat conduction kappa = 0.
Relevant logs and/or images
Temperature-field (left) and alpha-field (right)
Temperature field inside the bubble (clipped with alpha = 0.5)
Environment information
- OpenFOAM version : tested in v1912 (fixed version with interface.correct()) and v2006
- Operating system : debian
- Hardware info :
- Compiler : gcc
Possible fixes
Didn't find one yet.