CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

p_rgh in VOF solvers in version 1.7

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 2, 2010, 09:59
Default p_rgh in VOF solvers in version 1.7
  #1
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 16
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
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:
  • First thing is that I don't know how VOF "understands" what is the phase at which the buoyant term should be applied. I think it doesn't, as it solves one momentum equation for all.
  • When the density of one phase is much higher than the density of the other phase, the solution of one momentum equation in the domain, taking into account both densities, will give inaccurate results, especially when those two densities have a very big difference between them.
  • In spraying applications like the ones I simulate, most people don't put their domain one metre or two up, but design their (3D) domains in the positive x-axis, with the domains being symmetrical on the y and z axes and usually gravity is applied on one of them (y or z). When p_rgh is used with standard designs, the fluids inside the domain are (assuming for instance a cylindrical domain symmetrical on the y-z plane), half of the domain has a positive potential energy and the other half has a negative potential energy. So is the buoyant effect correct in this case?
Please give me your ideas and I might come back with other philosophical questions on the application of VOF.

Cheers
thekay is offline   Reply With Quote

Old   July 3, 2010, 02:35
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by thekay View Post
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).
The explanation of preferring p_rgh to p is numerical, not physical, and it is aimed at improving the robustness of the algorithm (force balance involving g and p).
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:
  • First thing is that I don't know how VOF "understands" what is the phase at which the buoyant term should be applied. I think it doesn't, as it solves one momentum equation for all.

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:
  • When the density of one phase is much higher than the density of the other phase, the solution of one momentum equation in the domain, taking into account both densities, will give inaccurate results, especially when those two densities have a very big difference between them.
This is what VOF does. Why should it be inaccurate?

Quote:
  • In spraying applications like the ones I simulate, most people don't put their domain one metre or two up, but design their (3D) domains in the positive x-axis, with the domains being symmetrical on the y and z axes and usually gravity is applied on one of them (y or z). When p_rgh is used with standard designs, the fluids inside the domain are (assuming for instance a cylindrical domain symmetrical on the y-z plane), half of the domain has a positive potential energy and the other half has a negative potential energy. So is the buoyant effect correct in this case?
Answered above :-)

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.
alberto is offline   Reply With Quote

Old   July 4, 2010, 11:54
Default
  #3
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 16
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
Quote:
Originally Posted by alberto View Post
The explanation of preferring p_rgh to p is numerical, not physical, and it is aimed at improving the robustness of the algorithm (force balance involving g and p).
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.
[/LIST]
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.




This is what VOF does. Why should it be inaccurate?



Answered above :-)

Best,
Thanks for the answer Alberto. I'm mostly discussing it in a philosophical way to get what is the idea behind adding the term and making it standard for VOF. I think it is more robust having used it once. But I will explain better, in order to clarify why I ask and what kind of discussion I am trying to raise. I think it's an excellent addition to Eulerian solvers but I can't see why it is a must for VOF.

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
thekay is offline   Reply With Quote

Old   July 4, 2010, 12:52
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by thekay View Post
Thanks for the answer Alberto. I'm mostly discussing it in a philosophical way to get what is the idea behind adding the term and making it standard for VOF. I think it is more robust having used it once. But I will explain better, in order to clarify why I ask and what kind of discussion I am trying to raise. I think it's an excellent addition to Eulerian solvers but I can't see why it is a must for VOF.
Hardly a must, but a good addition for sure. Try to think to some extreme situation, with phase separation. Solving for p_rgh gives you the "hydrostatic balance" straight out of the pressure equation.

Quote:
On my first point, I was comparing to Eulerian multiphase cases (e.g. fluidised beds), hence the question.
Well, in fluidised beds nearly at the minimum fluidization the story is very different because what you have there are four terms: the shared pressure p, the particle pressure and the momentum exchange term between the phases.
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:
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?
Why should you do that? (I'm trying to make you think about it )
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.
alberto is offline   Reply With Quote

Old   July 4, 2010, 13:12
Default
  #5
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 16
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
Thank you Alberto, I am getting the idea better now.
thekay is offline   Reply With Quote

Old   July 4, 2010, 13:14
Default
  #6
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
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
alberto is offline   Reply With Quote

Old   July 4, 2010, 13:50
Default
  #7
Member
 
Costas K
Join Date: Nov 2009
Location: UK
Posts: 51
Rep Power: 16
thekay is on a distinguished road
Send a message via ICQ to thekay Send a message via MSN to thekay
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.
thekay is offline   Reply With Quote

Old   July 4, 2010, 15:49
Default
  #8
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 36
alberto will become famous soon enoughalberto will become famous soon enough
Quote:
Originally Posted by thekay View Post
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.
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.
alberto is offline   Reply With Quote

Old   June 10, 2012, 14:13
Question Hydrostatic Pressure
  #9
Member
 
Mohammad Fereshtehpour
Join Date: Jul 2011
Location: Iran
Posts: 61
Rep Power: 15
MOHAMMAD67 is on a distinguished road
Send a message via Skype™ to MOHAMMAD67
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
Attached Images
File Type: jpg Screenshot.jpg (27.6 KB, 30 views)
MOHAMMAD67 is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


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


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