MahdiNabil / CFD-PC

interThermalPhaseChangeFoam - CFD simulation platform for liquid-vapor thermal phase change flows
http://sites.psu.edu/mtfe/
67 stars 43 forks source link

strange behaviour of interThermalPhaseChangeFoam #3

Open jerzyDudek2 opened 8 years ago

jerzyDudek2 commented 8 years ago

Hello Mahdi Nabil, I created simple 1D test case for evaporation, which is basically the Stefan problem for which analytical solution exists. The test case is: the left wall T = 398.15K, the right wall is adiabatic and for the rest of the surfaces, empty boundary condition is appllied (alpha1 = 1; only in the vicinity of the left wall alpha1 = 0, TSat = 373.15). In this configuration solver works well and the water evaporates.

However, when I changed the boundary conditions, i.e. the left wall is adiabatic, the right wall T = 398.15 K (alpha1 = 1; only in the vicinity of the right wall alpha1 = 0), solver gives bad solution and water does not evaporate!?? Can you explain me this? In my opinion it should be the mirror image. Probably I am doing something wrong. Here is my test case:

ShuStefanLeeEvap.tar.gz

best regards Przemek