[petros]

Hi all,

I have modified the example case condensatingVessel in the interCondensatingEvaporatingFoam folder in order to simulate water evaporation.

Some useful files to start with: geometry, checkMesh

Boundary conditions: T, alpha.liquid, p_rgh, U and p

Also, the necessary fvSchemes and fvSolution

What have I done:
Vessel is half filled with water using setFields. Water is initiated at saturation temperature (298.15K) while the rest of the domain at 313.15K.

The problem:
While I can see that some evaporation is happening at the interface:



the temperature field is blown up after almost 0.12s:



Here is also the pressure field during the same period of time.

Any idea of what I am doing wrong?

Petros

[petros]

Pardon. Images got meshed up.

Here are the alpha.liquid and temperature snapshots during the first 0.5 seconds. The problematic distribution is depicted in the latter.

Petros

[Strumm]

I am looking at similar problem, if you post the case setup (including your constant file , properties, etc. I will share what I learn.

[petros]

I do not have the case at the moment but as I stated in the first post, the base case was the condensatingVessel. Properties in the constant folder remained the same. The links are still working though so you can have a pretty clear idea about the set up.

I dropped the project because I didn't find any solution to the problem.

Best,
Petros

[Lorenzo210]

Hi Petros,
For what I can see inalpha.water image, it seems that your problem concerns with the intrerface compression. Also, I saw that in fvSolution you set cAplha=0, try setting it to 1 (it is a parameter you can find over the description for "alpha." solver).
I recently wrote a reply in [interCondensatingEvaporatingFoam] Film boiling - mixing problem #10 where I explain what I modified in interCondensatingEvaporatingFoam source code in order to have a well defined interface during time process. If cAlpha=1 doesn't change your simulation, you can try the modification I explained in the post above (REMEMBER to copy the original files into a backup folder).

I don't know what's wrong with the temperature, maybe you should define the default value in setFields, or maybe there's something wrong with the thermophhysical properties in the constant folder.

Maybe the next days I'll have the time to simulate your test case in OpenFOAM v2012, which is more stable and accurate in my opinion. Hope to give you some good news.

Lorenzo