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

Wiki solver guides

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 25, 2010, 02:57
Default Wiki solver guides
  #1
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi,

a few days ago I updated the following guides on the wiki, for those interested:
If you have comments, suggestions or questions, please use this thread.

This is part of the content that should have gone in the documentation project, and the structure used for the guides is very similar to what was the structure of the documents that should have become part of it. Since it did not work out for the reasons we know, I decided to put on the wiki the part of the material I was working on. Keeping it on my hard drive was not helping anybody.

More will come in some time, since rearranging the format is all but fun (no LaTeX macros can be defined to my knowledge on the wiki. If you know a way, please share ).

Best,
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 16, 2010, 19:42
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi Alberto,

I've seen you're other thread about these same guides, but only just recently saw this thread...
Quote:
Originally Posted by alberto View Post
More will come in some time, since rearranging the format is all but fun (no LaTeX macros can be defined to my knowledge on the wiki. If you know a way, please share ).
It's been a while since I last used LaTeX, so I think there are at least two possible interpretations for "LaTeX macros":
  • wiki macros that interpret LaTeX;
  • "porting" LaTeX macros to wiki, i.e., using the full LaTeX system under wiki.
Uhm... both sound nearly the same. Either way, I did some searching, but haven't gotten around to test either solutions that I've encountered:
  1. Convert LaTeX (and others) to wiki (and others) - Pandoc - by the description in its page, looks to me to be the most promising possibility;
  2. Online and perl script, for basic convertion from LaTeX to wiki - Latex-to-Wiki-converter - this one is pretty basic, but is downloadable and is a single script
  3. Plugin for openfoamwiki - Mediawiki and LaTeX - this would have to be installed in openfoamwiki.net, and might only work for some things, namely equations.
And Thank you for the solver know-how you've shared

In the mean time, I'm still hanging on to the idea of linking up openfoamwiki.net with -extend's SVN in sourceforge... i.e., putting the LaTeX PDF/PS (and even source) versions in the SVN, and making a link to them in the respective solver page in the wiki

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   March 16, 2010, 22:13
Default
  #3
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by wyldckat View Post
Hi Alberto,

I've seen you're other thread about these same guides, but only just recently saw this thread...
Yes, I'm very sorry for the duplicate post. I did not know this thread required moderation to start a new thread, and I simply thought my message would have not been accepted because I posted it in the wrong place, so I opened the other thread. Sorry about that.

Quote:
  1. Plugin for openfoamwiki - Mediawiki and LaTeX - this would have to be installed in openfoamwiki.net, and might only work for some things, namely equations.
Thanks for the tips. The wiki supports LaTeX equations already. The major problem are references, sectioning, command definition and so on, but with soma patience it can be done

Quote:
And Thank you for the solver know-how you've shared
You're welcome.

Quote:
In the mean time, I'm still hanging on to the idea of linking up openfoamwiki.net with -extend's SVN in sourceforge... i.e., putting the LaTeX PDF/PS (and even source) versions in the SVN, and making a link to them in the respective solver page in the wiki
This does not depend on me, and you already know that.

Best,
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 25, 2010, 04:56
Default
  #4
New Member
 
Christian Gross
Join Date: Oct 2009
Location: Germany
Posts: 5
Rep Power: 7
christian_gross is on a distinguished road
Thanks for the guides. They are much appreciated.
christian_gross is offline   Reply With Quote

Old   March 28, 2010, 07:25
Default
  #5
Member
 
Piotr Prusinski
Join Date: Oct 2009
Location: Warsaw, Poland
Posts: 53
Rep Power: 7
piprus is on a distinguished road
First of all, thanks for this contribution. It's a really great job!

One question concerning:

Quote:
Originally Posted by alberto View Post
More will come in some time,
Can I expect such guide for compressibleInterFoam solver in near future?
piprus is offline   Reply With Quote

Old   March 28, 2010, 20:56
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hello,

that solver is not among those I considered to document. Hopefully someone else will contribute the documentation for it, or you could write it if you decide to study how it works.

Best,
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   April 3, 2010, 11:24
Default
  #7
Member
 
Piotr Prusinski
Join Date: Oct 2009
Location: Warsaw, Poland
Posts: 53
Rep Power: 7
piprus is on a distinguished road
Quote:
Originally Posted by alberto View Post
Hello,

that solver is not among those I considered to document. Hopefully someone else will contribute the documentation for it, or you could write it if you decide to study how it works.

Best,

Thanks for your replay, sorry to read it unfortunately, but maybe something will come up due to another thread. Who knows maybe new solver guide
piprus is offline   Reply With Quote

Old   April 3, 2010, 13:29
Default
  #8
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
That would be nice :-)
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   January 27, 2011, 06:29
Default typo, bug or not understood
  #9
Senior Member
 
Wouter van der Meer
Join Date: May 2009
Location: Elahuizen, Netherlands
Posts: 128
Rep Power: 8
wouter is on a distinguished road
Hello Alberto,

Thanks for the many contribution to this forum. I was reading the code in the wiki page on bubblefoam and in Ueqns.H I came across:
UbEqn = ... Cvm*rhob*alpha/rhob .... should this not be:
... Cvm*rhoa*alpha/rhob ... ?
I found in OF-1.7.0 of CAElinux that this statement is also in twoPhaseEulerFoam, same .H file?

Best regards
Wouter
wouter is offline   Reply With Quote

Old   January 27, 2011, 13:49
Default
  #10
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
I think it is correct as it is, the virtual mass force is given by

Cvm * rho_continuous_phase * ...

in the literature.
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   January 27, 2011, 15:52
Default
  #11
Senior Member
 
Wouter van der Meer
Join Date: May 2009
Location: Elahuizen, Netherlands
Posts: 128
Rep Power: 8
wouter is on a distinguished road
Thanks,

I assume I am wrong again if I say that rhob can be removed then altogether?


best regards
Wouter
wouter is offline   Reply With Quote

Old   January 27, 2011, 16:17
Default
  #12
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
You are actually right. I think they left to make the code readable, or because they cut and paste Ua, and replaced ;-)
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   February 16, 2012, 21:44
Default
  #13
Member
 
chenxizhong
Join Date: Feb 2012
Posts: 30
Rep Power: 5
chenxizh is on a distinguished road
hi alberto , I am reading the bubblefoam guides it is very usefull, i thingk there is a mistake in governing equtation of momenton , the Reynolds stress tensor term shuold not have the density rou?
Quote:
Originally Posted by alberto View Post
Hi,

a few days ago I updated the following guides on the wiki, for those interested:
If you have comments, suggestions or questions, please use this thread.

This is part of the content that should have gone in the documentation project, and the structure used for the guides is very similar to what was the structure of the documents that should have become part of it. Since it did not work out for the reasons we know, I decided to put on the wiki the part of the material I was working on. Keeping it on my hard drive was not helping anybody.

More will come in some time, since rearranging the format is all but fun (no LaTeX macros can be defined to my knowledge on the wiki. If you know a way, please share ).

Best,
chenxizh is offline   Reply With Quote

Old   February 17, 2012, 01:05
Default
  #14
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi and thanks for the comment.

I don't think the density should appear, since it is assumed to be constant, and the momentum equation is divided by the phase density (note that the kinematic viscosity is used).

Best,
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   February 17, 2012, 01:21
Default
  #15
Member
 
chenxizhong
Join Date: Feb 2012
Posts: 30
Rep Power: 5
chenxizh is on a distinguished road
hi, I mean that in section 3.2 the third and fourth term in momentum eqution is obvious not dimemsion balance ?
Quote:
Originally Posted by alberto View Post
Hi and thanks for the comment.

I don't think the density should appear, since it is assumed to be constant, and the momentum equation is divided by the phase density (note that the kinematic viscosity is used).

Best,
chenxizh is offline   Reply With Quote

Old   February 17, 2012, 01:33
Default
  #16
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Thanks. I fixed it :-)
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 28, 2012, 11:36
Default
  #17
Member
 
chenxizhong
Join Date: Feb 2012
Posts: 30
Rep Power: 5
chenxizh is on a distinguished road
Hi alberto, recently I run the tutorial bed2(the bubble bed) used twoPhaseEulerFoam, I found that the max volume fraction of the solid will larger than the alphaMax, the alphaMax is 0.62 and alpha actually can reach 0.73? I've turned on the PackingLimter . So what is the problem? can you give me some suggestion? Thank you.
chen
Quote:
Originally Posted by alberto View Post
Thanks. I fixed it :-)
chenxizh is offline   Reply With Quote

Old   March 28, 2012, 14:30
Default
  #18
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Hi, that's an issue on how the particle pressure is managed if you use the kinetic theory. If you take a look at the code, it is treated explicitly.

As discussed elsewhere, I proposed a fix to this problem, but it's not released yet.
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 28, 2012, 20:29
Default
  #19
Member
 
chenxizhong
Join Date: Feb 2012
Posts: 30
Rep Power: 5
chenxizh is on a distinguished road
So do you mean that in this version(2.1.0), the foam cannot avoid this problem ? I used Gidsaspow's radial function, the solid pressure will be large when alpha approach the alphaMax to avoid the problem in physcially. So can you give me a hint about how you fix the problem. Thanks a lot.
Quote:
Originally Posted by alberto View Post
Hi, that's an issue on how the particle pressure is managed if you use the kinetic theory. If you take a look at the code, it is treated explicitly.

As discussed elsewhere, I proposed a fix to this problem, but it's not released yet.
chenxizh is offline   Reply With Quote

Old   March 28, 2012, 21:36
Default
  #20
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by chenxizh View Post
So do you mean that in this version(2.1.0), the foam cannot avoid this problem ? I used Gidsaspow's radial function, the solid pressure will be large when alpha approach the alphaMax to avoid the problem in physcially. So can you give me a hint about how you fix the problem. Thanks a lot.
The code is implemented correctly, and, if you use a small-enough time step, it should indeed work. However its stability is not strong.

The typical approach to enforce the packing limit is to:
  1. Remove the singularity in the particle pressure. One way is to use CarnahanStarling radial distribution function (btw, Gidaspow's radial distribution is not correct, since it is not 1 when alpha is zero. The correct version of that formulation for g0 is the one of Sinclair and Jackson).
  2. Use a continuous function to manage the dense (frictional) limit. In such limit, the kinetic theory is not applicable anyways, because collisions are not binary anymore, and there is sustained contact among particles.
  3. Include the effect of the particle pressure into the equation for alpha. OpenFOAM does this if the "ppMagf" model is used. Note that ppMagf is the "elastic modulus" of the phase, which is nothing but the derivative d(pa)/d(alpha), being pa the particle pressure. The same approach can be extended to manage packing limit with the kinetic theory model.
What kind of problems do you need to solve? I could try to merge some of the developments I have in my code into twoPhaseEulerFoam, and release them. After all, they are all written in my paper
__________________
Alberto Passalacqua

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.
OpenQBMM - An open-source implementation of quadrature-based moment methods
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
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
New solver guides alberto OpenFOAM 7 March 10, 2010 13:13
CFD Wiki - We Need More Help! Jonas Larsson Main CFD Forum 0 September 26, 2005 08:11
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25


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