about bounded alpha equation for bubbleFoam and twoPhaseEulerFoam

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 April 12, 2011, 09:44 about bounded alpha equation for bubbleFoam and twoPhaseEulerFoam #1 Member   Kai Join Date: May 2010 Location: Stockholm Posts: 40 Rep Power: 7 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 \nabla \cdot \bar{U}=0. 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.

 May 5, 2011, 11:00 #2 Member   Kai Join Date: May 2010 Location: Stockholm Posts: 40 Rep Power: 7 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.

 May 6, 2011, 05:33 #3 Member   Sabin Ceuca Join Date: Mar 2010 Location: Munich Posts: 42 Rep Power: 7 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post alberto OpenFOAM Running, Solving & CFD 26 May 26, 2015 03:05 raagh77 OpenFOAM Running, Solving & CFD 98 September 23, 2013 06:31 su_junwei OpenFOAM Running, Solving & CFD 3 November 18, 2011 03:35 piterw82 OpenFOAM Bugs 0 June 21, 2010 08:45 juho OpenFOAM Running, Solving & CFD 2 May 20, 2008 11:31

All times are GMT -4. The time now is 19:27.

 Contact Us - CFD Online - Top