pressure correction method

ravikant · November 6, 2009, 06:36

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

otd · November 6, 2009, 10:01

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.

ztdep · November 7, 2009, 20:45

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.

ravikant · November 9, 2009, 06:53

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

.

jyothishkumar · November 9, 2009, 09:41

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

ztdep · November 9, 2009, 19:08

Quote:

Originally Posted by jyothishkumar

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.

jyothishkumar · November 10, 2009, 06:18

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

ztdep · November 10, 2009, 07:00

Quote:

Originally Posted by jyothishkumar

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 · November 10, 2009, 07:02

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

atlantis · November 25, 2009, 23:52

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.

November 6, 2009, 06:36	pressure correction method	#1
ravikant New Member A Ravi Kant Join Date: Nov 2009 Posts: 4 Rep Power: 16	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

November 10, 2009, 06:18	momentum interpolation	#7
jyothishkumar Member jk Join Date: Jun 2009 Posts: 64 Rep Power: 16	Hi, Momentum interpolation u mean to say rowuu or rowuv for u momentum eqn. (finding the value between the two nodes) and scalar equation u mean to say the pressure correction equation ? cheers jyothish

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
"Pressure Inlet" Boundary Setup	Wijaya	FLUENT	15	May 18, 2016 10:08
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

November 6, 2009, 10:01		#2
otd Member private Join Date: Mar 2009 Posts: 74 Rep Power: 17	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.

November 7, 2009, 20:45		#3
ztdep Senior Member p ding Join Date: Mar 2009 Posts: 427 Rep Power: 19	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.

November 9, 2009, 06:53		#4
ravikant New Member A Ravi Kant Join Date: Nov 2009 Posts: 4 Rep Power: 16	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 .

November 9, 2009, 09:41		#5
jyothishkumar Member jk Join Date: Jun 2009 Posts: 64 Rep Power: 16	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

November 10, 2009, 07:02		#9
ztdep Senior Member p ding Join Date: Mar 2009 Posts: 427 Rep Power: 19	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

November 25, 2009, 23:52		#10
atlantis New Member Quoit Join Date: Nov 2009 Posts: 1 Rep Power: 0	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.