GitLab

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