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

Finite Volume -- SIMPLE Algorithm

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

LinkBack Thread Tools Search this Thread Display Modes
Old   December 16, 2008, 10:28
Default Finite Volume -- SIMPLE Algorithm
Posts: n/a

My name is Roger,

I am a grad student writing a lid-driven cavity code in matlab using the Finite-Volume-Method for a class. I am having trouble getting my SIMPLE algorithm to work. Does any body have any sample code for the Simple Algorithm, especially the initialization process for multiple iterations?

(I am using the Intro to Computational Fluid Dynamics by Versteeg as my guide. In example 6.2 I do not get the same value for the discretized momentum equation of node 1 for the second iteration as does Versteeg.)

Thank you,

  Reply With Quote

Old   December 17, 2008, 22:30
Default Re: Finite Volume -- SIMPLE Algorithm
Posts: n/a
I would suggest going to the book written by the guy who basically invented SIMPLE. Patankar's Numerical Heat Transfer and Fluid Flow. Your school library should have a copy (probably multiple copies), your professor should probably have a copy, too. It really is the seminal work on the SIMPLE procedure.
  Reply With Quote

Old   December 24, 2008, 14:58
Default Re: Finite Volume -- SIMPLE Algorithm
Posts: n/a
FWIW, SIMPLE is a rather poor special case of a more general class: Schur complement derived preconditioners for iterative (usually Krylov) methods. For a good review, see

title={{A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Navier-Stokes equations}},
author={Elman, H.C. and Howle, V.E. and Shadid, J. and Shuttleworth, R. and Tuminaro, R.},
journal={Journal of Computational Physics},
publisher={Academic Press}

There isn't much reason to use the classical methods such as SIMPLE and Uzawa, except to show that they're not very good.
  Reply With Quote

Old   January 2, 2009, 18:30
Default Re: Finite Volume -- SIMPLE Algorithm
Posts: n/a
i will traying the same thing good luck for us
  Reply With Quote

Old   February 6, 2010, 18:04
New Member
Join Date: Feb 2010
Posts: 4
Rep Power: 16
SamR is on a distinguished road
Hi Roger,

I am also trying to implement the SIMPLE algorithm (in Matlab) using the example from Versteeg, and have had the same problem with the second iteration. Towards the end of page 208, it says that they should have included under-relaxation measures for aP, Su and d, so I figured that was why I got different values for the discretised momentum equation on the second iteration.

I know its over a year since you originally posted , but was wondering if you found any useful resources to solve your issues?

Hope you can help.
SamR is offline   Reply With Quote

Old   October 25, 2010, 03:16
Default You can try SIMPLE-TS
New Member
kirilhs's Avatar
Kiril Shterev
Join Date: Jan 2010
Location: Sofia, Bulgaria
Posts: 2
Blog Entries: 1
Rep Power: 0
kirilhs is on a distinguished road
The SIMPLE-like methods have some disadvantages.

You can try SIMPLE-TS, which do not need under-relaxation coefficients.

The source code written in c++ is also available.

kirilhs is offline   Reply With Quote

Old   February 10, 2011, 03:47
New Member
Join Date: Feb 2011
Posts: 1
Rep Power: 0
spejbl is on a distinguished road

I have written a code for a driven cavity problem (Finite-Volume-Method / SIMPLE-method for MATLAB) but's not really working. So it would be really nice if someone have a working MATLAB code for that problem which you can send to me so that I can check out where my mistakes are. I really dont know.

Thanks for helping!

spejbl is offline   Reply With Quote

Old   February 11, 2011, 08:18
Join Date: Jun 2009
Posts: 64
Rep Power: 16
jyothishkumar is on a distinguished road

Have you tried the SIMPLE algorithm for 1d nozzle or diffuser. if not please try to do that and it wont take much time. If you are not getting the solution then let me know.


jyothishkumar is offline   Reply With Quote

Old   June 25, 2011, 22:49
Default solution
Join Date: Jul 2010
Location: chennai
Posts: 79
Rep Power: 15
kabilan is on a distinguished road
hai dear
i am trying to write code for 1d nozzle problem in versteeg book example 6.2. in that the program not working properly and i have more doubt in residual solving, so please send me program to
thank you
kabilan is offline   Reply With Quote

Old   September 25, 2023, 12:04
New Member
Join Date: Sep 2023
Posts: 1
Rep Power: 0
panyx is on a distinguished road
do you finish the problem?
panyx is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Control Volume Finite Element Method gerardosrez Main CFD Forum 0 March 16, 2011 14:49
Finite Difference Vs. Finite Volume elankov Main CFD Forum 43 December 18, 2010 16:30
[blockMesh] Axisymmetrical mesh Rasmus Gjesing (Gjesing) OpenFOAM Meshing & Mesh Conversion 10 April 2, 2007 14:00
Finite Volume WENO Scheme DSS Main CFD Forum 3 January 15, 2007 03:47
Finite Element vs. Finite Volume Codes Salman Main CFD Forum 9 January 6, 2006 22:59

All times are GMT -4. The time now is 05:10.