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

TwoPhaseEulerFOAM application

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 27, 2008, 16:05
Default Why should it be better with 1
  #21
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Why should it be better with 1e-15 instead than 1e-4 from a numerical point of view? This limitation is useful when alpha goes to zero, because you divide by it, so 1e-4 should have a better "stabilizing effect".

Moreover, you can set the inferior limit of alpha, in the case of a particle flow to something under 1e-4 (1e-6 is what is used in the kinetic theory framework, and it's conservative) without problems from a physical point of view, because the flow is so dilute that you can't consider the two-fluid equations valid in those conditions, and, as a consequence, you should not use them anyway.

Regards,
Alberto
__________________
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 28, 2008, 02:15
Default What I meant is that max(alpha
  #22
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 80
Rep Power: 7
juho is on a distinguished road
What I meant is that max(alpha,anumber) only affects alpha values smaller than anumber (clips the scale at the bottom) while (alpha+anumber) affects all the alpha values. (moves the whole scale upwards)

To me, the clipping seems like a more elegant solution. Assuming the max() function doesn't introduce new problems - like the dimension mentioned by Rasmus.

Thank you all very much for your responses on the whole issue!
juho is offline   Reply With Quote

Old   May 28, 2008, 10:23
Default Hello, yes, sorry. I was no
  #23
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Hello,

yes, sorry. I was not discussing about clipping. I agree with you on that.

I was just pointing out that you can happily clip alpha at a numerically safe value around (like 10^-6) instead than to a very small value, that in a division would return a huge coefficient.

Regards,
Alberto
__________________
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 28, 2008, 18:42
Default Hi Rasmus, I have implement
  #24
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
Hi Rasmus,

I have implemented your changes in our 1.5 release candidate and all seems well.

Thanks

Henry
henry is offline   Reply With Quote

Old   May 29, 2008, 08:45
Default Thanks Henry!
  #25
Senior Member
 
Rasmus Hemph
Join Date: Mar 2009
Location: Sweden
Posts: 108
Rep Power: 7
hemph is on a distinguished road
Thanks Henry!
hemph is offline   Reply With Quote

Old   August 29, 2008, 20:05
Default I think that in "liftDragCoeff
  #26
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 7
danielle is on a distinguished road
I think that in "liftDragCoeffs.H" file for blended case
K should be calculated as:
K = (alpha*Ka + beta*Kb);
not:
K = (beta*Ka + alpha*Kb);
danielle is offline   Reply With Quote

Old   September 9, 2008, 09:34
Default are you agree with me !!!
  #27
Member
 
Danielle PRL
Join Date: Mar 2009
Posts: 42
Rep Power: 7
danielle is on a distinguished road
are you agree with me !!!
danielle is offline   Reply With Quote

Old   September 9, 2008, 10:12
Default assume phase a to be the discr
  #28
Super Moderator
 
niklas's Avatar
 
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 692
Rep Power: 18
niklas will become famous soon enough
assume phase a to be the discrete one and b to be the continous.

standard approach would be to calculate the drag as K = Ka, right?

alpha = 0 -> pure gas.

so...
as alpha -> 0 we should have K -> Ka

but with your approach K -> Kb

so no, i dont agree.
niklas is offline   Reply With Quote

Old   May 28, 2011, 20:53
Default
  #29
New Member
 
Prashant Gupta
Join Date: Mar 2011
Location: Edinburgh
Posts: 29
Rep Power: 5
Prash is on a distinguished road
Hey Guys,

I think I am on the right thread. But I am trying to write a new drag law for my calculations.

I found this :

http://openfoamwiki.net/index.php/Ad...lerFoam_solver

I really cannot understand why -


The Syamlal and O'Brien drag model can be implemented in OpenFOAM by dividing the β' < math > bytheproductofthetwovolumefractions < math > αβ to obtain K.
N.B. Every drag formulation has to be divided by the product αβ before beeing implemented in twoPhaseEulerFoam because this product has been extracted for numerical reasons.



Why do we devide by the product of voulme fractions and where in the source code do we multiply them back ???





Please help me with this.





Best Wishes

Prashant
Prash is offline   Reply With Quote

Old   May 29, 2011, 04:54
Default
  #30
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Hi Prashant,

you divide by the product of the phase fractions because of how the momentum equation is treated (phase-intensive form). You might want to check Henrik Rusche's Ph.D. thesis for some more detail on the implementation.

P.S. Notice that all the models listed there are already part of OpenFOAM now.

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.

Last edited by alberto; May 29, 2011 at 04:56. Reason: Added P.S.
alberto is offline   Reply With Quote

Old   June 6, 2011, 06:25
Default
  #31
New Member
 
Prashant Gupta
Join Date: Mar 2011
Location: Edinburgh
Posts: 29
Rep Power: 5
Prash is on a distinguished road
Hey,

Thank you very much for useful info, I had another question on calculations and use of solid volume fractions for different species. I am trying to do a segregation study on fluidized bed , hence I require total volume fraction(alpha) and volume fraction of each of two species of particles ( different sizes).

I am not sure where to look for them. I am trying to access the species volume fraction in order to calculate drag on individual species.

Waiting eagerly for a reply

Best Wishes
Prashant
Prash is offline   Reply With Quote

Old   June 11, 2011, 16:55
Default
  #32
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
The code deals with mono-disperse flows. If you need more than one dispersed phase, you'll have to extend it.

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   November 5, 2011, 16:03
Default more than one discrete phase simulation
  #33
New Member
 
pramod kumar
Join Date: Sep 2011
Posts: 3
Rep Power: 5
pramod308k is on a distinguished road
hi,
I want to simulate binary fluidized bed problem with OpenFoam, In my case more than one discrete phase are present. My question is can I use twoPhaseEulerFoam for more than one discrete phase.???? If yes than please give me some hint.

thank you
Pramod Kumar
pramod308k is offline   Reply With Quote

Old   November 5, 2011, 16:08
Default
  #34
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
Hi Pramod,

twoPhaseEulerFoam can deal with one particle phase.

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   November 6, 2011, 01:54
Default
  #35
New Member
 
pramod kumar
Join Date: Sep 2011
Posts: 3
Rep Power: 5
pramod308k is on a distinguished road
hi Alberto,
can you please tell me how to set a case for more than one disperse phase in openFoam. I'm very new with this software.

thank you
Pramod
pramod308k is offline   Reply With Quote

Old   November 6, 2011, 02:06
Default
  #36
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,888
Rep Power: 25
alberto will become famous soon enoughalberto will become famous soon enough
There is no Euler-Euler multi-fluid solver in OF at the moment.
__________________
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

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
TwoPhaseEulerFoam sara OpenFOAM Running, Solving & CFD 2 November 6, 2008 20:26
Bug in twoPhaseEulerFoam alberto OpenFOAM Bugs 2 May 20, 2008 22:25
TwoPhaseEulerFoam Bug alondono OpenFOAM Bugs 1 February 19, 2008 21:01
TwoPhaseEulerFOAM application hemph OpenFOAM Bugs 0 November 16, 2006 08:27
TwoPhaseEulerFoam newbee OpenFOAM 0 March 27, 2006 08:41


All times are GMT -4. The time now is 16:54.