# Finite Volume -- SIMPLE Algorithm

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

 December 16, 2008, 11:28 Finite Volume -- SIMPLE Algorithm #1 Roger Guest   Posts: n/a Hello, 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, Roger

 December 17, 2008, 23:30 Re: Finite Volume -- SIMPLE Algorithm #2 Bignose Guest   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.

 December 24, 2008, 15:58 Re: Finite Volume -- SIMPLE Algorithm #3 Jed Guest   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 @article{elman2008tcp, 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}, volume={227}, number={1}, pages={1790--1808}, year={2008}, 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.

 January 2, 2009, 19:30 Re: Finite Volume -- SIMPLE Algorithm #4 samia Guest   Posts: n/a i will traying the same thing good luck for us

 February 6, 2010, 19:04 #5 New Member   Join Date: Feb 2010 Posts: 4 Rep Power: 7 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.

 October 25, 2010, 03:16 You can try SIMPLE-TS #6 New Member   Kiril Shterev Join Date: Jan 2010 Location: Sofia, Bulgaria Posts: 2 Blog Entries: 1 Rep Power: 0 Hi, The SIMPLE-like methods have some disadvantages. You can try SIMPLE-TS, which do not need under-relaxation coefficients. http://www.cfd-online.com/Forums/blo...ble-fluid.html The source code written in c++ is also available. Regards

 February 10, 2011, 04:47 #7 New Member   Join Date: Feb 2011 Posts: 1 Rep Power: 0 Hey 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

 February 11, 2011, 09:18 #8 Member   jk Join Date: Jun 2009 Posts: 64 Rep Power: 8 Hello, 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. thanks jyothish

 June 25, 2011, 22:49 solution #9 Member   kabilan.B Join Date: Jul 2010 Location: chennai Posts: 77 Rep Power: 7 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 me.hareesh.aero@gmail.com thank you

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post gerardosrez Main CFD Forum 0 March 16, 2011 15:49 elankov Main CFD Forum 43 December 18, 2010 17:30 Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00 DSS Main CFD Forum 3 January 15, 2007 04:47 Salman Main CFD Forum 9 January 6, 2006 23:59

All times are GMT -4. The time now is 04:18.