CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

about bounded alpha equation for bubbleFoam and twoPhaseEulerFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 12, 2011, 09:44
Default about bounded alpha equation for bubbleFoam and twoPhaseEulerFoam
  #1
Member
 
Kai
Join Date: May 2010
Location: Stockholm
Posts: 40
Rep Power: 7
kaifu is on a distinguished road
Hi Foamers,
This time I will deal with subcooled boiling flows modeling which includes evaporation and condensation models. Actually the task is the continuous work of Edy, see Modeling Evaporation at a Heated Wall

I start my solver from bubbleFoam and twoPhaseEulerFoam, which use Weller's strategy to solver UEqn, pEqn, and alphaEqn. In Weller's approach, one of the most important things is that he found a way to bound alpha in [0,1]. Here I need to point out that the reason why alpha is bounded at one is due to the "conservative" feature of volume mixed flux, that is <eq>\nabla \cdot \bar{U}=0</eq>.

However in my case, since there is a phase change term (due to evaporation and condensation) in alpha equation, which means alpha could be still bounded at one but not zero any more by Weller's approach, the program crashed. I noticed that the divergence may come from pEqn.H, which should promise the "conservative" feature of volume mixed flux. However it reports as,

time step continuity errors : sum local = 4.22651e+09, global = 5.82749e+07, cumulative = 5.82749e+07

And then the program crashed...

So anybody could give me suggestions? I have actually considered to start my solver from interPhaseChangeFoam, but I do not need to track the interface...

Thanks a lot.
kaifu is offline   Reply With Quote

Old   May 5, 2011, 11:00
Default
  #2
Member
 
Kai
Join Date: May 2010
Location: Stockholm
Posts: 40
Rep Power: 7
kaifu is on a distinguished road
Some updates: this problem comes from pressure equation. In bubbleFoam since there is no phase change, the volume flux is conservative. When dealing with phase change, we have to use the mass flux instead of volume flux as in bubbleFoam.
kaifu is offline   Reply With Quote

Old   May 6, 2011, 05:33
Default
  #3
Member
 
Sabin Ceuca
Join Date: Mar 2010
Location: Munich
Posts: 42
Rep Power: 7
sabin.ceuca is on a distinguished road
Hi Kaifu,
I’m currently working on direct contact condensation with VoF and if I understood your problem correctly then interPhaseChangeFoam is not the right source for you. That solver considers phase change due to cavitation and has a one fluid approach. The only information gain is that it shows you where you’ll have to add source terms in the alphaEqn.
sabin.ceuca is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Rewriting twoPhaseEulerFoam in conservative form alberto OpenFOAM Running, Solving & CFD 26 May 26, 2015 03:05
TwoPhaseEulerFoam and Boundary conditions raagh77 OpenFOAM Running, Solving & CFD 98 September 23, 2013 06:31
Bounded scheme gives unbounded solution su_junwei OpenFOAM Running, Solving & CFD 3 November 18, 2011 03:35
mpi problem with twoPhaseEulerFoam piterw82 OpenFOAM Bugs 0 June 21, 2010 08:45
TwoPhaseEulerFoam viscosity models juho OpenFOAM Running, Solving & CFD 2 May 20, 2008 11:31


All times are GMT -4. The time now is 14:59.