strange behaviour of interCondensingEvaporatingFoam for a simple 1D Stefan condensation test case
I created simple 1D test case for condensation, which is basically the Stefan problem for which analytical solution exists. The test case is a rod filled with a water vapour (alpha.liquid = 0) at saturation temperature of 380.26 K. In the first case (StefanCond_left) I set:
- on the left wall T = 370.26 K, the right wall is adiabatic and for the rest of the surfaces, empty boundary condition is appllied . In this configuration one can observe rapid condensation and after time 0.1 s all vapour is condensed (alpha.liquid = 1)
In the second case (StefanCond_right) I changed only the side of the cold wall, namely:
- the left wall is adiabatic, the right wall is at T = 370.26 K and for the rest of the surfaces, empty boundary condition is appllied . In this configuration one can observe very slow condensation and after time 0.1 s only 0.0254598 of the vapour condensed (alpha.liquid = 0.0254598).
Is it a bug? The solution should be the mirror image. Can you explain why the condensation rate decreases so much while it shouldn't chaneg? Here are the two mentioned test cases (35MB):
http://fluid.itcmp.pwr.wroc.pl/~pblasiak/download/StefanCond.tar.gz