CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Deriving pressure from velocity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 14, 2014, 09:42
Default Deriving pressure from velocity
  #1
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
Hello everybody,

I am currently trying to implement a CFD code for deriving pressure from velocity in 2D. It's not a problem which has been solved a lot.

Problem : " I have a velocity field on a structured mesh but not necessarily rectangular..and I have to find the pressure field."

First I have tried to solve this problem on a rectangular fully populated grid and then on any grids with the Finite Element method.

I try different methods for solving this problem: an iterative method and the Pressure Poisson equation.
I haven't good results with both.

But I focused on the pressure poisson equation. Do you think this method could give some good results for different kind of flows ? ( not too turbulent )
I compute F=Ux^2 + 2 UxUy + Uy^2 and then solve the linear system AP=F where P is the pressure and A the 2D-laplacian matrix. I implemented Neumann or dirichlet boundaries inside A and F depending on the problem.

No relevant results..So I tried to do the same with a staggered grid with P on the center of the cell and U and V on the edges. But no good results or the A matrix was singular .

So finally I want to know if I am on the right way or if I should choose an other method ? I am quite lost at the time being..

Thank you very much,
I can clarify some points if necessary !

Samy
samycfd is offline   Reply With Quote

Old   August 18, 2014, 16:41
Default Pressure from velocity
  #2
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
Hi again,

Is everyone had already compute the pressure field from velocity field ?
What the more accurate method ?

Thank you very much,

Best,
Samuel
samycfd is offline   Reply With Quote

Old   August 21, 2014, 13:31
Default
  #3
Member
 
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 93
Rep Power: 7
bluebase is on a distinguished road
Hi samycfd,

unfortunately i don't have any helpful advice for this topic. Iin my studies abroad, i met a Professor of Niigata University, Prof. Fujisawa. His laboraty was working on a simular topic.
I recalled the following paper "Measurement of Pressure Field around a NACA0018 Airfoil from PIV Velocity Data", Journal of Visualization.
http://iospress.metapress.com/content/x13646u580l57uh8/
There is a preview page of this short paper. Maybe all relevant information for you is included in that paragraph.

With regards,
Sebastian
bluebase is offline   Reply With Quote

Old   August 22, 2014, 04:44
Default
  #4
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
Thank you for your answer. I know this method but it doesn't work for me.

Have you ever solved the Pressure Poisson equation ?
samycfd is offline   Reply With Quote

Old   August 22, 2014, 11:20
Default
  #5
Member
 
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 93
Rep Power: 7
bluebase is on a distinguished road
No, i haven't.
bluebase is offline   Reply With Quote

Old   August 25, 2014, 12:03
Default
  #6
Member
 
Shenren Xu
Join Date: Jan 2011
Location: London, U.K.
Posts: 63
Rep Power: 6
Shenren_CN is on a distinguished road
Solving the Poisson equation the way you described seems okay. Could you post the velocity field you used and one pressure solution you computed just for the readers to get an idea how good/bad your pressure solution is?

Poisson eqn. might be stiff to solve, so it'll help if you could show the convergence of your Poisson solve.



Quote:
Originally Posted by samycfd View Post
Hello everybody,

I am currently trying to implement a CFD code for deriving pressure from velocity in 2D. It's not a problem which has been solved a lot.

Problem : " I have a velocity field on a structured mesh but not necessarily rectangular..and I have to find the pressure field."

First I have tried to solve this problem on a rectangular fully populated grid and then on any grids with the Finite Element method.

I try different methods for solving this problem: an iterative method and the Pressure Poisson equation.
I haven't good results with both.

But I focused on the pressure poisson equation. Do you think this method could give some good results for different kind of flows ? ( not too turbulent )
I compute F=Ux^2 + 2 UxUy + Uy^2 and then solve the linear system AP=F where P is the pressure and A the 2D-laplacian matrix. I implemented Neumann or dirichlet boundaries inside A and F depending on the problem.

No relevant results..So I tried to do the same with a staggered grid with P on the center of the cell and U and V on the edges. But no good results or the A matrix was singular .

So finally I want to know if I am on the right way or if I should choose an other method ? I am quite lost at the time being..

Thank you very much,
I can clarify some points if necessary !

Samy
Shenren_CN is offline   Reply With Quote

Old   August 26, 2014, 03:45
Default
  #7
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
I compute the Lid-driven cavity case in FLUENT ( pressure and velocity) and then with my algorithm I get the following pressure field( blue )

https://drive.google.com/folderview?...EU&usp=sharing
samycfd is offline   Reply With Quote

Old   August 26, 2014, 08:28
Default
  #8
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
Do you think it would be better solving the PPE with the Finite Volume method ?

Ps: I often have Neumann BC on each edges so I define one reference point on an edge as Dirichlet value to avoid the matrix to be singular ( as discussed in many other threads... )

And the particular point in the pressure field is the point I defined.. If i change the location we will able to see it in the pressure field. So something is apparently wrong !
samycfd is offline   Reply With Quote

Old   August 26, 2014, 09:37
Default update !
  #9
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
I have found an error ! So my pressure field is

, for this velocity field :

samycfd is offline   Reply With Quote

Old   August 28, 2014, 19:51
Default
  #10
New Member
 
Join Date: Aug 2014
Posts: 2
Rep Power: 0
12er is on a distinguished road
May I ask if the flow problem of the thread's author has an incompressibility constraint? Nobody mentioned so far if incompressibility is assumed or not. So I just wanted to know. (I was in the believe that the Poisson Equation can be only derived, if incompressibility is assumed, but I may be wrong such that my question is of no relevance).

Last edited by 12er; September 2, 2014 at 07:40.
12er is offline   Reply With Quote

Old   August 29, 2014, 02:46
Default
  #11
Member
 
Join Date: Aug 2014
Posts: 58
Rep Power: 3
samycfd is on a distinguished road
Quote:
if the flow problem of the thread's author has an incompressibility constraint?
Yes you're right. I assume that the constraint is verified and I checked if it's true. If not I know that results could be wrong.

Do you know if I can reconstruct the velocity field in order it satisfy the constraint? As in the pressure correction algorithm ?(SIMPLE)
samycfd 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
Boundary Condition for pressure and velocity mecman OpenFOAM Pre-Processing 1 May 7, 2014 16:51
Initial pressure and transverse velocity fields to initialize turbulence model nickvinn Main CFD Forum 0 February 29, 2012 11:11
lid-driven cavity in matlab using BiCGStab Don456 Main CFD Forum 1 January 19, 2012 16:00
How to set pressure BC with mass Velocity Magnitud arwang FLUENT 2 March 12, 2007 21:04
what the result is negatif pressure at inlet chong chee nan FLUENT 0 December 29, 2001 06:13


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