CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Bugs

Bug in UEqn of twoPhaseEulerFoam?

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

Like Tree1Likes
  • 1 Post By florian_krause

Reply
 
LinkBack Thread Tools Display Modes
Old   March 24, 2011, 15:21
Default Bug in UEqn of twoPhaseEulerFoam?
  #1
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 77
Rep Power: 6
cheng1988sjtu is on a distinguished road
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:
00022 UaEqn =
00023 (
00024 (scalar(1) + Cvm*rhob*beta/rhoa)*
00025 (
00026 fvm::ddt(Ua)
00027 + fvm::div(phia, Ua, "div(phia,Ua)")
00028 - fvm::Sp(fvc::div(phia), Ua)
00029 )
00030
00031 - fvm::laplacian(nuEffa, Ua)
00032 + fvc::div(Rca)
00033
00034 + fvm::div(phiRa, Ua, "div(phia,Ua)")
00035 - fvm::Sp(fvc::div(phiRa), Ua)
00036 + (fvc::grad(alpha)/(fvc::average(alpha) + scalar(0.001)) & Rca)
00037 ==
00038 // g // Buoyancy term transfered to p-equation
00039 - fvm::Sp(beta/rhoa*K, Ua) //Drag term
00040 //+ beta/rhoa*K*Ub // Explicit drag transfered to p-equation
00041 - beta/rhoa*(liftCoeff - Cvm*rhob*DDtUb)
00042 );
00043
00044 UaEqn.relax();
00045 }
The drag term should involve the difference between Ua and Ub, not just Ua.

I've tried to implement that (Ub-Ua), 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?
cheng1988sjtu is offline   Reply With Quote

Old   March 24, 2011, 15:31
Default
  #2
Senior Member
 
Laurence R. McGlashan
Join Date: Mar 2009
Posts: 370
Rep Power: 13
l_r_mcglashan will become famous soon enough
The explicit part of the drag that contains Ub has been transferred to pEqn.H. Have a look there and see if it makes sense to you.
__________________
Laurence R. McGlashan :: Website
l_r_mcglashan is offline   Reply With Quote

Old   March 24, 2011, 16:18
Default Thank you!
  #3
Member
 
Charlie
Join Date: Dec 2010
Location: 415 Kinross Dr. Newark, DE 19711
Posts: 77
Rep Power: 6
cheng1988sjtu is on a distinguished road
Yes, Thanks for your reply, it really helps!

The effect of particle-particle 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
cheng1988sjtu is offline   Reply With Quote

Old   March 25, 2011, 03:19
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,890
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
If you use the kinetic theory model, you should set g0 in ppProperties to zero.
__________________
Alberto

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
GeekoCFD 32bit - The 32bit edition of GeekoCFD.
GeekoCFD text mode - A smaller version of GeekoCFD, text-mode only, with only OpenFOAM. Available in a variety of virtual formats.
alberto is offline   Reply With Quote

Old   May 25, 2011, 22:21
Default
  #5
New Member
 
NieYongguang
Join Date: Sep 2010
Posts: 27
Rep Power: 6
nygbook is on a distinguished road
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!
nygbook is offline   Reply With Quote

Old   May 26, 2011, 00:11
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,890
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by nygbook View Post
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.
Berend van Wachem 's Ph.D. thesis.

Quote:
So, can you tell me which book include these equation. Furthermore, where can I find the theory about g0 in ppProperties!
In the literature ppMagf is identified as G(alpha) = d(pa)/d(alpha), being pa the particle pressure, and it is called "normal stress modulus". You find it in Gidaspow's early papers for example.

Best,
__________________
Alberto

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
GeekoCFD 32bit - The 32bit edition of GeekoCFD.
GeekoCFD text mode - A smaller version of GeekoCFD, text-mode only, with only OpenFOAM. Available in a variety of virtual formats.
alberto is offline   Reply With Quote

Old   August 26, 2011, 07:11
Default
  #7
Member
 
Greg Melling
Join Date: Nov 2010
Posts: 41
Rep Power: 6
grjmell is on a distinguished road
Does anyone know how to get the van Wachem Thesis? There's nothing online. Does anyone have a pdf?
grjmell is offline   Reply With Quote

Old   August 26, 2011, 07:41
Default
  #8
Senior Member
 
Florian Krause
Join Date: Mar 2009
Location: Munich
Posts: 103
Rep Power: 7
florian_krause is on a distinguished road
Quote:
Originally Posted by grjmell View Post
Does anyone know how to get the van Wachem Thesis? There's nothing online. Does anyone have a pdf?
Not true! I needed 30sec using google twice (because I knew that he got his PhD at TU Delft) to find it

http://repository.tudelft.nl/view/ir...-83b1416709a6/

enjoy!
Florian
ThomasV likes this.
florian_krause is offline   Reply With Quote

Old   August 26, 2011, 07:45
Default
  #9
Member
 
Greg Melling
Join Date: Nov 2010
Posts: 41
Rep Power: 6
grjmell is on a distinguished road
Oh dear! I was googling for ages but couldnt find it...!
Thanks for proving me wrong!
grjmell 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
Something wrong in UEqns.H within twoPhaseEulerFoam cheng1988sjtu OpenFOAM 2 June 24, 2011 11: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 22: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


All times are GMT -4. The time now is 18:12.