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

pressure correction method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 6, 2009, 07:36
Angry pressure correction method
  #1
New Member
 
A Ravi Kant
Join Date: Nov 2009
Posts: 4
Rep Power: 7
ravikant is on a distinguished road
Hi all,
I am writing a 2D matlab code based on pressure correction method for lid driven cavity flow problem. When I frame the pressure poisson matrix {grad^2 P'=0} with neumann pressure boundary condition {dp'/dn=0} at all the 4 walls, what I observed is that the determinant [A] of AX=B becomes zero. I would be greatful if anyone can tell me how do I get out of this problem.

Thank you,
Ravi
ravikant is offline   Reply With Quote

Old   November 6, 2009, 11:01
Default
  #2
otd
Member
 
private
Join Date: Mar 2009
Posts: 74
Rep Power: 8
otd is on a distinguished road
The pressure is determined only to a constant with Neumann boundary conditions.

You're free to pin the pressure at one node to solve the matrix problem.

If you solve the pressure equation by relaxation you will get a solution. But only the pressure gradients will be unique - the overall level will depend on your initial guess for the pressure field.
otd is offline   Reply With Quote

Old   November 7, 2009, 21:45
Default
  #3
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Yes , you are right. We can not obatain a solution if we adopte the direct method to solve the Ax=B, we should use the iteration method.
in the incompressible flow, we only care about the pressure difference in the flow field.
ztdep is offline   Reply With Quote

Old   November 9, 2009, 07:53
Default
  #4
New Member
 
A Ravi Kant
Join Date: Nov 2009
Posts: 4
Rep Power: 7
ravikant is on a distinguished road
Finally I am able to solve the lid driven problem by fixing one of the node's pressure correction value=0. Thanks for your suggestions .
ravikant is offline   Reply With Quote

Old   November 9, 2009, 10:41
Default
  #5
Member
 
jk
Join Date: Jun 2009
Posts: 64
Rep Power: 8
jyothishkumar is on a distinguished road
Dear friends,

Thanks, I have written my staggered grid incompressible code using SIMPLE algorithm. Now that i want write it using a collocated grid. Is it direct or should i do some major changes ? Please help me in this regard.

cheers

jyo
jyothishkumar is offline   Reply With Quote

Old   November 9, 2009, 20:08
Default
  #6
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Quote:
Originally Posted by jyothishkumar View Post
Dear friends,

Thanks, I have written my staggered grid incompressible code using SIMPLE algorithm. Now that i want write it using a collocated grid. Is it direct or should i do some major changes ? Please help me in this regard.

cheers

jyo

no major changes are needed. momentum interpolation should be introduced. the coefficient of the U , V and W should be changed at that of the scalar equations.
ztdep is offline   Reply With Quote

Old   November 10, 2009, 07:18
Default momentum interpolation
  #7
Member
 
jk
Join Date: Jun 2009
Posts: 64
Rep Power: 8
jyothishkumar is on a distinguished road
Hi,

Momentum interpolation u mean to say row*u*u or row*u*v for u momentum eqn. (finding the value between the two nodes) and scalar equation u mean to say the pressure correction equation ?

cheers

jyothish
jyothishkumar is offline   Reply With Quote

Old   November 10, 2009, 08:00
Default
  #8
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
Quote:
Originally Posted by jyothishkumar View Post
Hi,

Momentum interpolation u mean to say row*u*u or row*u*v for u momentum eqn. (finding the value between the two nodes) and scalar equation u mean to say the pressure correction equation ?

cheers

jyothish
MIM is used to caculate the interface flux and source term of the pressure correction equations.
please refer to some CFD book.
Ferziger J.H., Peric M. Computational methods for fluid dynamics
ztdep is offline   Reply With Quote

Old   November 10, 2009, 08:02
Default
  #9
Senior Member
 
ztdep's Avatar
 
p ding
Join Date: Mar 2009
Posts: 287
Rep Power: 9
ztdep is on a distinguished road
Send a message via Yahoo to ztdep Send a message via Skype™ to ztdep
MIM is used to caculate the interface flux and source term of the pressure correction equations.
please refer to some CFD book.
Ferziger J.H., Peric M. Computational methods for fluid dynamics
ztdep is offline   Reply With Quote

Old   November 26, 2009, 00:52
Default
  #10
New Member
 
Quoit
Join Date: Nov 2009
Posts: 1
Rep Power: 0
atlantis is on a distinguished road
Hi friends, I, too, am working on this very same problem, and would really appreciate it if someone would be willing to provide their code to me.

Feel free to post it here or send it to me in a PM.

Thanks a lot! I really want to improve what I have by looking at others.
atlantis 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
"Pressure Inlet" Boundary Setup Wijaya FLUENT 14 August 26, 2015 06:02
Pressure Correction Method for Cylindrical Coordinates Eslam Reda Main CFD Forum 4 July 4, 2009 02:15
Streamfunction-vorticity method - Pressure find Perumal Main CFD Forum 1 September 23, 2008 14:43
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02
BC for Pressure Correction on Colocated Grid Matt Umbel Main CFD Forum 1 October 3, 2000 19:44


All times are GMT -4. The time now is 03:57.