|
[Sponsors] |
July 2, 2010, 09:59 |
p_rgh in VOF solvers in version 1.7
|
#1 |
Member
|
I would like to ask you, if any of you knows, what is the explanation of the use of p_rgh in multiphase solvers that use VOF. One would say that it provides a correction which represents the buoyancy, but this is as far as I see it, minor in terms of VOF applications that involve two different phases (e.g. air and water). I put a few points up for discussion and I would really like to hear your ideas about them:
Cheers |
|
July 3, 2010, 02:35 |
|
#2 | ||||
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
The set of equations solved in 1.6.x and in 1.7.x are formally identical. What changes is how they are treated to obtain a numerical solution. Quote:
This is really basic stuff! ;-) You are solving for a momentum equation for a mixture, but the density of this mixture is not constant, as a consequence the term rho*g in the momentum equation contains the density gradients. 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 (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
|||||
July 4, 2010, 11:54 |
|
#3 | |
Member
|
Quote:
On my first point, I was comparing to Eulerian multiphase cases (e.g. fluidised beds), hence the question. "Inaccurate" wasn't the right expression. I will phrase my question differently: having a big difference in densities, gives big gradients in a single cell. So, does the p_rgh have any effect on the convergence? But this is a case specific problem. I was wondering how other people understand the convergence effect (if any) of rgh in VOF in large density differences. On my last point, in order to get physically accurate results, does one need to "raise" his domain at the height where the actual application is performed in real life? And finally, having a symmetrical domain at h=0, does it have any effect apart from calculating buoyant terms which almost completely counter each other? Thanks |
||
July 4, 2010, 12:52 |
|
#4 | |||
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Quote:
Quote:
In a bed which tends to the packed condition, you have a competition between the particle pressure, the drag force and the gravity. Solve this balance incorrectly (solving it correctly is harder due to the presence of the particle pressure!), and your code won't be even close to the concepts of accuracy and stability :-) Quote:
Think to how the pressure is specified in your system. :-) 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 (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
||||
July 4, 2010, 13:12 |
|
#5 |
Member
|
Thank you Alberto, I am getting the idea better now.
|
|
July 4, 2010, 13:14 |
|
#6 |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
Maybe share why you do not have to raise your domain. Others might find that useful ;-)
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. Last edited by alberto; July 4, 2010 at 13:15. Reason: Grammar |
|
July 4, 2010, 13:50 |
|
#7 |
Member
|
It really depends on what pressure one wants to present. p_rgh can give negative p_rgh values. Somebody might not want to have it like that. That's how I see it.
|
|
July 4, 2010, 15:49 |
|
#8 |
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36 |
You can plot both (the code computes the actual p too), and this has nothing to do with the independence from the "height".
__________________
Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. |
|
June 10, 2012, 14:13 |
Hydrostatic Pressure
|
#9 |
Member
|
Dear Friends
Hi I want to know how I can obtain hydrostatic pressure. The results in a multiphase flow show that p is less than p-rgh, why? I couldn't understand which parameter I should use in comparison with experimental data. ( In reality we can have the hydrostatic pressure using pizometer) Also, What is the atmospheric pressure in these cases? Thanks in advance |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
OpenFOAM-1.6 install cookbook | MadsR | OpenFOAM Installation | 372 | November 20, 2010 11:57 |
[OpenFOAM] Problem with paraFoam on a linux-64 bit | bunni | ParaView | 4 | April 14, 2010 20:55 |
paraFoam reader for OpenFOAM 1.6 | smart | OpenFOAM Installation | 13 | November 16, 2009 21:41 |
Looking for VOF version of Nast2d | invadoria | Main CFD Forum | 0 | June 8, 2009 08:56 |
lack of support for VOF in version 3.26 | Tom | Siemens | 10 | August 6, 2006 08:20 |