CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Post-Processing

Full pressure for rasInterFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 20, 2007, 19:16
Default Did anyone create a tool which
  #1
Senior Member
 
Join Date: Mar 2009
Posts: 222
Rep Power: 8
paka is on a distinguished road
Did anyone create a tool which helps obtain the full pressure including both hydrostatic and dynamic pressures?

Current rasInterFoam solver gives only the modified pressures output in 'pd' file, which does not include hydrostatic terms.

Thanks,
Krystian
paka is offline   Reply With Quote

Old   August 21, 2007, 08:54
Default Yup, I wrote one of those a wh
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,748
Rep Power: 20
hjasak will become famous soon enough
Yup, I wrote one of those a while back. It works in a post-processing mode, but if you wish you can build it into the top-level solver as well. Depends what you want to do with the real pressure...

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   August 21, 2007, 14:19
Default Aloha Hrvoje, Would you min
  #3
Senior Member
 
Join Date: Mar 2009
Posts: 222
Rep Power: 8
paka is on a distinguished road
Aloha Hrvoje,

Would you mind posting a tool over here?

First I just need a real pressure data plot, than I would need to calculate surface forces using calcPressureForce tool.

Does your tool give an output file to something like p or pd file?
Based on one of these files the other tool which I use calculate forces.

Thanks,
Krystian
paka is offline   Reply With Quote

Old   August 21, 2007, 17:17
Default Well, actually, this was neith
  #4
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,748
Rep Power: 20
hjasak will become famous soon enough
Well, actually, this was neither easy to implement nor trivial to work out. Is there anything you are prepared to offer in return?

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak is offline   Reply With Quote

Old   August 21, 2007, 17:34
Default Well, keeping in mind I'm stil
  #5
Senior Member
 
Join Date: Mar 2009
Posts: 222
Rep Power: 8
paka is on a distinguished road
Well, keeping in mind I'm still green, I can only promise my current commitment to OpenFOAM and my input, in-between, like tools development (I had some small input to foamToTecplot postprocessing).

What else could I offer you?

Regards,
Krystian
paka is offline   Reply With Quote

Old   January 30, 2009, 03:40
Default Could someone please help me t
  #6
New Member
 
Richard Kenny
Join Date: Mar 2009
Location: Brisbane, Queensland, Australia
Posts: 4
Rep Power: 7
richard_kenny is on a distinguished road
Could someone please help me to understand the "p" and "pd" fields in "rasInterFoam"?

I am a civil engineer and would like to use OpenFOAM for the design of dam spillways. Specifically, I am trying to calculate the static pressure over a spillway crest. I am new to CFD and have just begun experimenting with OpenFOAM. I have no C++ programming experience. I have built and run a model but am having difficulty interpreting the pressure results.

The following four points summarise my understanding of the "p" and "pd" fields (based on my reading of the message board, OpenFOAM manuals and the rasInterFoam.C file).
1) "p" is the total pressure in Pascals [1 -1 -2 0 0 0]
2) p = pd + rho*gh (from rasInterFoam.C file)
3) rho*gh = static pressure in Pascals [1 -1 -2 0 0 0]
4) "pd" = dynamic pressure in Pascals [1 -1 -2 0 0 0] = 0.5*rho*U^2

I have built a test case with static water (i.e. water level at the crest of the spillway but not flowing over the spillway). The static water depth is 38.63 metres. The total and static pressures should both be 378,960.3 Pascals (p = rho*g*h = 1000*9.81*38.63). The solver is giving me a result of 369,834 Pa which is reasonably close but not exactly the same as the analytical solution (within about 2.4%).

The solver velocities (U) are very close to zero (eg. -0.000287, -0.000288, 0). Therefore, the dynamic pressure (pd) should be very small. However, the model results indicate a "pd" of 1,528,580 Pa. I don't understand this result. I suspect that I don't really understand what the "p" and "pd" fields are.

Any advice you can provide would be greatly appreciated.
richard_kenny is offline   Reply With Quote

Old   January 30, 2009, 17:21
Default Hello Kenny, my
  #7
Member
 
Maruthamuthu Venkatraman
Join Date: Mar 2009
Location: Norway
Posts: 80
Rep Power: 7
maruthamuthu_venkatraman is on a distinguished road
Hello Kenny,
my experience with OpenFoam and Free surface flows is not much , however i would like to mentions some comments.

1. Have you succeed with InterFoam solver? usually RAS modeling needs quite fine grids to predict the mean pressure. I would suggest lets just run interFoam and look at your pressure.

2. Also you shall refine the grid in your test case, Mostly discretisation errors are the main sources in CFD.

May be other Foamers may give better suggestions...
maruthamuthu_venkatraman is offline   Reply With Quote

Old   February 2, 2009, 04:02
Default Hi Richard, find attached a
  #8
New Member
 
Daniel Schmode
Join Date: Mar 2009
Posts: 22
Rep Power: 7
daniels is on a distinguished road
Hi Richard,

find attached a extended interFoam, which computes the pressure including hydrostatic effects. It requires 0/p defining boundary conditions for the pressure. I recommend to use fixedValue at top of your domain, and zeroGradient for all remaining boundaries.

Daniel


daniels is offline   Reply With Quote

Old   February 2, 2009, 04:06
Default here is the attachment again
  #9
New Member
 
Daniel Schmode
Join Date: Mar 2009
Posts: 22
Rep Power: 7
daniels is on a distinguished road
here is the attachment again

interFoamGL.tgz
daniels is offline   Reply With Quote

Old   February 3, 2009, 02:19
Default Maruthamuthu & Daniel, Than
  #10
New Member
 
Richard Kenny
Join Date: Mar 2009
Location: Brisbane, Queensland, Australia
Posts: 4
Rep Power: 7
richard_kenny is on a distinguished road
Maruthamuthu & Daniel,

Thanks for your assistance.

I have tried the interFoam solver and a fine mesh but have the same issue with the "pd" field. I have not yet tried the interFoamGL solver but hope to get time shortly.

In the meantime, could you please confirm whether my understanding of the "p" and "pd" fields is correct? In other words, are the four points listed in my post correct?

Thanks in advance,
Richard
richard_kenny is offline   Reply With Quote

Old   February 3, 2009, 04:23
Default Richard, in principle you a
  #11
New Member
 
Daniel Schmode
Join Date: Mar 2009
Posts: 22
Rep Power: 7
daniels is on a distinguished road
Richard,

in principle you are right, but the p = pd + rho*gh is only valid for simple free surface flows. Have a look at this discussion:

http://www.cfd-online.com/cgi-bin/Op...3280#POST23280

the solver posted above follows this proposal. It is similar to the tool Prof. Jasak mentioned further above.

Daniel
daniels is offline   Reply With Quote

Old   February 17, 2009, 10:13
Default Daniel, Thank you for posting
  #12
Senior Member
 
Kent Wardle
Join Date: Mar 2009
Location: Illinois, USA
Posts: 190
Rep Power: 9
kwardle is on a distinguished road
Daniel,
Thank you for posting your solver here. I have read many of the posts mentioning this issue between static pressure and pd, etc. so your solver is helpful. I have a question about boundary conditions for p, pd, and U for an outlet boundary. For example, say you want to solve for the liquid height of an asymmetrically shaped u-tube manometer with the short side 'submerged', like:

p1
---
| |
| |
| |
|x|
|x| p2
|x| ---
|x| |x|
|x| |x|
|x|_|x|
|x-x-x|

where p1 is atmospheric pressure and p2 is 0 *static pressure*. The initial gamma condition (marked with 'x' in crude diagram) shows the left column has a higher liquid height. Since a zero static p is specified for p2, the equilibrium gamma position in the left column should be equal to the height of the right side of the tube. What are the boundary conditions one would use for p, pd, and U to make this happen? For U, it would seem like you cannot use a pressureInletOutletVelocity BC in this case--one would like to have a BC for U that keeps the total *static pressure* equivalent (i.e. U from bernoulli). Do you (or anyone) have any suggestions on this?
THANK YOU.
kwardle is offline   Reply With Quote

Old   February 17, 2009, 12:31
Default hi foamers! I work on the i
  #13
Member
 
hamdi
Join Date: Mar 2009
Posts: 72
Rep Power: 7
hamcer is on a distinguished road
hi foamers!

I work on the injector (nozzel) of hydraulic turbin , my case is axisymetrique with a RANS turbulence model ,
Can I find a similar case in OpenFoam?, I m already investigated the LesInterfoam, and Rasinterfoam, but, I have a problemes with those cases, some advices will be welcome.
hamcer 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
Parallel rasInterFoam openfoam_user OpenFOAM Running, Solving & CFD 4 November 1, 2008 06:14
Static Pressure BC in rasInterFoam kwardle OpenFOAM Running, Solving & CFD 0 September 19, 2008 17:15
Pressure BCs for rasInterFoam tank fillingdraining problem kwardle OpenFOAM Running, Solving & CFD 8 September 17, 2008 15:37
RasInterFoam or lesInterFoam hsieh OpenFOAM Running, Solving & CFD 2 March 31, 2006 14:42
RasInterFoam cavitation maritozzo OpenFOAM Running, Solving & CFD 2 December 6, 2005 15:09


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