CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

SIMPLE algorithm implementation, pressure correction error

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

LinkBack Thread Tools Display Modes
Old   May 12, 2013, 01:25
Question SIMPLE algorithm implementation, pressure correction error
New Member
Join Date: Feb 2013
Posts: 15
Rep Power: 5
northfly is on a distinguished road
I have been implementing this SIMPLE algorithm to a 2D tunnel flow. I use inlet constant velocity boundary condition and outlet constant pressure boundary condition.

I can get very good flow velocity distribution, which have the max velocity in the center and zero on both tunnel wall. The max velocity magnitude also seems reasonable.

But I cannot get good pressure distribution, the pressure inside the tunnel could be high up to 7 times of boundary pressure. And the strange part is the pressure distirbution seems didn't influence the velocity momentum calculation a lot, I can still get good velocity contour.

Is there some one can give me some hint about this problem, which part I should look into? Thanks.
northfly is offline   Reply With Quote

Old   May 12, 2013, 22:56
Default HIs problem is really like me
New Member
Join Date: Feb 2013
Posts: 15
Rep Power: 5
northfly is on a distinguished road
I am also implementing versteeg's book, the pressure correction got similar problem. The velocity distribution profile, pattern are all ok, but the magnitude value, if compared with slit flow theoretical value


the max velocity is always 50% higher than the theoretical number, I believe the problem is at pressure correction and don't know how to solve it.
northfly is offline   Reply With Quote

Old   February 8, 2014, 01:15
Default Simple
New Member
Join Date: Apr 2011
Posts: 10
Rep Power: 7
hadieliasi is on a distinguished road
I think you may use SIMPLER Algorithm that is better than SIMPLE algorithm because in SIMPLER, the pressure equation instead of pressure correction equation is solving directly.
In SIMPLE even velocity distributions is correct but the pressure may not correct and needs to more iterations. but i think 2 way to solve your problem:
1-At first check your code and be sure that is correct. may be the problem is about defining reference point for pressure.
2-Please check b.c and be sure that they are correct. how do you define those b.c?
2-Modify your code to solve simpler algorithm. It is very easy and you should change some parts of your simple code. in patankar's book " heat transfer and fluid flow " these methods are existed.
i have teach and w2d code.If you need, I can send those you.
hadieliasi is offline   Reply With Quote


pressure correction, simple algorithm

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
Native ParaView Reader Bugs tj22 OpenFOAM Paraview & paraFoam 270 January 4, 2016 12:39
Timestep and Pressure Correction Relationship in SIMPLE rks171 Main CFD Forum 23 May 4, 2012 01:04
Problem running perturbUCyl sen.1986 OpenFOAM 14 March 23, 2012 05:12
attach/detach (valve opening/closing) phsieh2005 OpenFOAM Running, Solving & CFD 2 March 21, 2009 06:18
Problem with compile the setParabolicInlet ivanyao OpenFOAM Running, Solving & CFD 6 September 5, 2008 20:50

All times are GMT -4. The time now is 17:39.