Using SIMPLE
I'm wanting to write a small program to figure a few fluid flow fields for 2D cleanroom applications. I plan on using Patankar's SIMPLE at first. However, I could use a few tips. From what I understand, after guessing p*, the momentum equations are solved using iterations of the LinebyLine method. Then move on to solve p', p (from p' but applying an underrelaxation factor) u,v, and any other discretization eqns as required. Then iterate again using the p as new p*.
I guess I've got a little confused on which procedures require iteration, and which don't. ANyways, any advice would be greatly appreaciated. 
Re: Using SIMPLE
Before choosing the simple method i recomend you investigate the fractional step method.
I have coded both methods and for my applications (time dependent flows) i find fractional step method to be much fa For references of the fractional step method see www.cfdjournal.com, 06/29/1999"An alternative model of the turbulent energy cascade" Michael R. Fitzmaurice. 
Re: Using SIMPLE
(1). I am not going to talk about this 'not so simple SIMPLE METHOD'. (2). As you know that the NavierStokes equations are nonlinear. To make the thing more complicated, let's take the steadystate equation for incompressible flow for example. You have a firstorder continuity equation, momentum equations and energy equation of secondorder. (3). Since we don't know where and how to obtain the solution for this set of steadystate equations, the iterative method is normally used to obtain the solution from the discretized equations. (4). In the incompressible flow, the pressure field is a function of the velocity field. In other word, you don't have to solve the velocity field and the pressure field at the same time. But, in the primitive variable formulation, velocity and pressure are used to be the primitive variable at the same time. (5). This is the case with the SIMPLE METHOD, which deal with the velocity and the pressure directly at the same time. Therefore, in order to solve the momentum equations with two unknown primitive variables ( velocity vector and the pressure), one has to first guess at the value of the pressure field ( guessing is basically iterative) and solve the momentum equations for the velocity field ( How?) (6). The steadystate momentum equation with known pressure term is still nonlinear. So this part has to be solved iteratively in order to obtain the converged final velocity field consistent with the boundary conditions of the particular problem you are trying to solve. (7). But since the pressure field is the identical twin of the velocity field, the guessed initial pressure values are not the right answer. So, you have to update the pressure field through a separate route. For this reason, a separate pressure route, or pressure equation is derived in such a way that an updated pressure can be computed. This is saying that the pressure field ( the identical twin of the velocity field) also has to be iteratively updated. (8). If you are lucky, if your initial guess is good (I don't know how to define it yet), and if the iterative method parameters are properly selected ( I don't know how to set those either), the posibility of obtaining the answer is relatively high. (9). Since the pressure can be decoupled from the original momentum equations easily, the separate route is to use the "vorticityvelocity" formulation, or in 2D case the "vorticitystream function" formulation. These methods can also be used for variable density, turbulent flows. In this case, the pressure field is computed after the converged velocity field is obtained. I have been using the vorticity related approach for turbulent, variable density flow with low Reynolds number twoequation models for over 25 years. It is relatively easy to use with this approach. Unfortunately, you are not going to find my published papers in this area. (10). Since the pressurebased methods also are widely used in commercial codes, it might be a good idea to learn how these pressurerelated methods actually work. good luck.

Re: Using SIMPLE
You can find open source code from ftp.springer.de For example caffa.for (Finite Volume Method using SIMPLE algorithm on colocated bodyfitted grids)
You may also look at H.K Versteeg and W. Malalasekera, " An introduction to computational fluid dynamics: The finite volume, John Wiley & Sons Inc., New York. This book is very usefol for beginners. 
Re: Using SIMPLE
Could you give me some more info about ftp.springer.de? I got lost looking for caffa.for.
Thanks 
Re: Using SIMPLE
please look at /pub/technik/peric/2dgl/sg/*.for

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