|
[Sponsors] |
Calculating pressure coefficient on surface of an object in ABL simulation |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 | |
New Member
Faiaz Khaled
Join Date: Oct 2017
Location: United States of America
Posts: 17
Rep Power: 8 ![]() |
Hello Everyone,
I would like to thank all the experts who leave comments and suggestions on CFD problems on this forum. I am grateful to them. I am trying to reproduce atmospheric boundary layer flow using openfoam and to measure surface pressure coefficients on a simple building placed in the domain. I am using Komega turbulence model since this model is more useful for simple algorithm. I tried to follow the guidelines regarding yplus for Komega model. I want to calculate the pressure coefficient on the surfaces of the building placed on the ground in the computational domain. I will briefly write out the boundary conditions below: Code:
p inlet { type zeroGradient; } outlet { type fixedValue; value uniform 0; } Code:
Quote:
Now, my specific question is, while calculating pressure coefficient of the building surface, do I need to subtract the p recorded on the side near the wall of the domain? I am little confused on this part, as I am trying to match with experimental data. When I subtract the side p value at reference height from the building p value, there is little gap from experimental value. Otherwise, if I only use the surface recorded p (without subtracting the side p value) then I have better match with experimental data. Just to make it clear, to obtain pressure coefficient from 'p' (obtained from probe or from paraview), should I just multiply it by 2 as my reference velocity is 1 m/s? Cp=(p-p_ref)/(0.5*rho*U^2) Here, U is 1 m/s. rho is already considered under 'p'. Considering the case I described please let me know, the way I am trying to process pressure coefficient is correct or wrong? Please shed some light on how to obtain the p_ref from the domain I described. Thanks Regards, Faiaz |
||
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Kmeti Rao
Join Date: May 2019
Posts: 145
Rep Power: 8 ![]() |
Hi Faiaz,
As per my knowledge in an incompressible simulation, the pressure is defined as gradient. Therefore, the output pressure is the differential pressure and the reference pressure p_ref can be set to zero. Thus you can directly select the output value of surface pressure given by paraview. As you have already described, don't forget to multiply by density. Krao |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
New Member
Faiaz Khaled
Join Date: Oct 2017
Location: United States of America
Posts: 17
Rep Power: 8 ![]() |
Hello Krao,
thank you for your quick reply. Quote:
I have been reading on different threads regarding the topic. But, the question is how to set or control the reference pressure in OpenFOAM simple algorithm? I have 0 p-value on the outlet and also collecting values on building surface in terms of pressure/rho. So, do you think I need to subtract anything from the reading on building surface? Please let me know. If I need to subtract which value should I consider? So far, I placed two probes near the walls of the domain at roof height of the building and subtracting average of those p/rho from the readings on the building surface. I dont know if I am confusing myself too much. Best Regards, Faiaz |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Kmeti Rao
Join Date: May 2019
Posts: 145
Rep Power: 8 ![]() |
Hi Faiaz,
I also use zero at the outlet for my propeller simulation. I am not sure about your problem, but in my case I select the surface pressure on the propeller for further calculations. I usually extract surface pressure directly using paraview and multiply it with density. Regards |
|
![]() |
![]() |
![]() |
![]() |
#5 | |
New Member
Faiaz Khaled
Join Date: Oct 2017
Location: United States of America
Posts: 17
Rep Power: 8 ![]() |
Quote:
Thank you for your reply, Krao. I Really appreciate it. ![]() I was wondering if somebody else has a similar experience on extracting pressure from paraview in openfoam with incompressible case. Just wanted to be sure.. if I have 0 p-value on the outlet then do I need to to subtract any value from p reading on an object placed in the computational domain. Thanks Regards, Faiaz |
||
![]() |
![]() |
![]() |
![]() |
#6 |
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,906
Rep Power: 33 ![]() |
Hi,
Don't do it like that - please use the forceCoeffs function object, which will use the native data and FVM operators, (eg for the viscous stress) within FOAM. Good luck, Hrv
__________________
Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk |
|
![]() |
![]() |
![]() |
![]() |
#7 | |
New Member
Faiaz Khaled
Join Date: Oct 2017
Location: United States of America
Posts: 17
Rep Power: 8 ![]() |
Quote:
![]() I read a little about the forceCoeffs function. Seems like people use it to obtain the drag, lift and moment coefficient etc. I can dive more into it to find how to obtain pressure coefficient with that. If I use this function, what is the impact of having the outlet P boundary condition to be 0. But, my another question, was the approach mentioned in the earlier post right or wrong? I am also using probes to collect p values for specific locations of interest. So, the p-value I am getting from the probes at the end of the converged simulation, is it the final (pressure/rho) value to calculate the pressure coefficient or I need to add/subtract something from that p value? Probably, a more generic question would be if I have a p-value of 0 at the outlet, then can I consider my reference pressure to be zero and process any pressure in the domain without adding or subtracting anything from it? Thanks Please let me know. Regards, Faiaz ![]() |
||
![]() |
![]() |
![]() |
![]() |
#8 | |||
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quick answers:
Quote:
Quote:
Quote:
... although it's still without the hydrostatic pressure, but you could calculate that field too and add it to the expression, to get the total pressure.
__________________
|
||||
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Simulation of a single bubble with a VOF-method | Suzzn | CFX | 21 | January 29, 2018 00:58 |
simpleFoam - pressure (coefficient) of head shape | GJM1991 | OpenFOAM Running, Solving & CFD | 4 | May 12, 2015 17:15 |
How ot plot pressure over an object surface? | vilius | STAR-CCM+ | 2 | March 25, 2015 15:50 |
Calculating the force exerted on an object from fluid pressure | amrbekhit | CFX | 1 | January 30, 2011 16:38 |
Hydrostatic pressure in 2-phase flow modeling (long) | DS & HB | Main CFD Forum | 0 | January 8, 2000 15:00 |