
[Sponsors] 
March 24, 2011, 14:21 
Bug in UEqn of twoPhaseEulerFoam?

#1  
Member
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 79
Rep Power: 7 
Hi All:
I've spent some time interpreting the UEqn.C of twoPhaseEulerFoam. The momentum equation seems OK. however, the drag term seems problematic: See the following link for the standard momentum equation. http://openfoamwiki.net/index.php/BubbleFoam In the UEqn.H, the momentum equation is : Quote:
I've tried to implement that (UbUa), I just added the following line in the rhs: + fvm::Sp(beta/rhoa*K, Ub) For the UbEqns, I've done a similar modification Am I right? Anyone can help? 

March 24, 2011, 15:18 
Thank you!

#3 
Member
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 79
Rep Power: 7 
Yes, Thanks for your reply, it really helps!
The effect of particleparticle interaction is also included if we set g0>0, However, note that there is g0 in radial Model (we get g0 and g0Prime), and we also define g0 in ppProperties. Both alphaEqn.H and pEqn.H are using phi_pp, which involves g0, So which g0 are they using? Thank you in advance. Zhen 

March 25, 2011, 02:19 

#4 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,904
Rep Power: 27 
If you use the kinetic theory model, you should set g0 in ppProperties to zero.
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

May 25, 2011, 21:21 

#5 
New Member
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 7 
Hi Alberto
I am reading the code of kineticTheory in twoPhaseEulerFoam. I find there are some annotation in kineticTheory.C such as "construct the granular temperature equation (Eq. 3.20, p. 44)" in line 61. So, can you tell me which book include these equation. Furthermore, where can I find the theory about g0 in ppProperties! 

May 25, 2011, 23:11 

#6  
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,904
Rep Power: 27 
Quote:
Quote:
Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 26, 2011, 06:11 

#7 
Member
Greg Melling
Join Date: Nov 2010
Posts: 41
Rep Power: 7 
Does anyone know how to get the van Wachem Thesis? There's nothing online. Does anyone have a pdf?


August 26, 2011, 06:41 

#8  
Senior Member
Florian Krause
Join Date: Mar 2009
Location: Munich
Posts: 103
Rep Power: 9 
Quote:
http://repository.tudelft.nl/view/ir...83b1416709a6/ enjoy! Florian 

August 26, 2011, 06:45 

#9 
Member
Greg Melling
Join Date: Nov 2010
Posts: 41
Rep Power: 7 
Oh dear! I was googling for ages but couldnt find it...!
Thanks for proving me wrong! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Something wrong in UEqns.H within twoPhaseEulerFoam  cheng1988sjtu  OpenFOAM  2  June 24, 2011 10:48 
Possible bug in UEqn for raslescavitatingFoam  henry  OpenFOAM Bugs  3  November 14, 2008 00:13 
Bug in twoPhaseEulerFoam  alberto  OpenFOAM Bugs  2  May 20, 2008 21:25 
TwoPhaseEulerFoam Bug  alondono  OpenFOAM Bugs  1  February 19, 2008 21:01 
Bug in twoPhaseEulerFoam wallfunctions  alberto  OpenFOAM Bugs  1  February 9, 2007 15:15 