
[Sponsors] 
July 16, 1999, 10:34 
Using SIMPLE

#1 
Guest
Posts: n/a

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. 

July 16, 1999, 20:45 
Re: Using SIMPLE

#2 
Guest
Posts: n/a

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. 

July 17, 1999, 17:22 
Re: Using SIMPLE

#3 
Guest
Posts: n/a

(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.


July 20, 1999, 06:54 
Re: Using SIMPLE

#4 
Guest
Posts: n/a

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. 

July 20, 1999, 22:59 
Re: Using SIMPLE

#5 
Guest
Posts: n/a

Could you give me some more info about ftp.springer.de? I got lost looking for caffa.for.
Thanks 

July 21, 1999, 02:38 
Re: Using SIMPLE

#6 
Guest
Posts: n/a

please look at /pub/technik/peric/2dgl/sg/*.for


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
SIMPLE algorithm in 3D cylindrical coordinates  zouchu  Main CFD Forum  1  January 20, 2014 18:02 
The correction on pressure equation of SIMPLE algorithm in MRFSimpleFOAM solver  renyun0511  OpenFOAM Running, Solving & CFD  0  November 10, 2010 02:47 
How to add correcting pressure equation in SIMPLE of MRFSimpleFOAM?  renyun0511  OpenFOAM Programming & Development  0  November 4, 2010 02:38 
PISO vs. SIMPLE  benedikt flurl  Main CFD Forum  2  April 14, 2005 06:54 
SIMPLE algorithm  Jonathan Castro  Main CFD Forum  3  December 10, 1999 05:59 