CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

interCondensatingEvaporatingFoam vessel wrong T distribution

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By Lorenzo210

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 24, 2020, 13:17
Default interCondensatingEvaporatingFoam vessel wrong T distribution
  #1
Member
 
Petros Ampatzidis
Join Date: Oct 2018
Location: Bath, UK
Posts: 62
Rep Power: 5
petros is on a distinguished road
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 is offline   Reply With Quote

Old   January 24, 2020, 19:32
Default
  #2
Member
 
Petros Ampatzidis
Join Date: Oct 2018
Location: Bath, UK
Posts: 62
Rep Power: 5
petros is on a distinguished road
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
petros is offline   Reply With Quote

Old   September 24, 2020, 12:38
Default Post the Case
  #3
New Member
 
Strumm
Join Date: Oct 2015
Posts: 1
Rep Power: 0
Strumm is on a distinguished road
I am looking at similar problem, if you post the case setup (including your constant file , properties, etc. I will share what I learn.
Strumm is offline   Reply With Quote

Old   September 24, 2020, 13:23
Default
  #4
Member
 
Petros Ampatzidis
Join Date: Oct 2018
Location: Bath, UK
Posts: 62
Rep Power: 5
petros is on a distinguished road
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
petros is offline   Reply With Quote

Old   February 9, 2021, 05:19
Default
  #5
New Member
 
Lorenzo
Join Date: Apr 2020
Location: Italy
Posts: 22
Rep Power: 3
Lorenzo210 is on a distinguished road
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
raj kumar saini likes this.
Lorenzo210 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' muth OpenFOAM Running, Solving & CFD 3 August 27, 2018 05:18
Pressure distribution ok and wrong coefficients (Cl and Cd) Jeanp OpenFOAM Post-Processing 2 March 11, 2017 23:06
meshing of a compound volume in GMSH shawn3531 OpenFOAM 4 March 12, 2015 11:45
What differenceof retention time distribution fuction and residence time distribution zsq FLUENT 0 December 14, 2009 21:29
Rosin Rammler Distribution AndiG FLUENT 1 March 7, 2007 11:56


All times are GMT -4. The time now is 15:06.