CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Using SIMPLE

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 16, 1999, 10:34
Default Using SIMPLE
  #1
david
Guest
 
Posts: n/a
I'm wanting to write a small program to figure a few fluid flow fields for 2-D 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 Line-by-Line method. Then move on to solve p', p (from p' but applying an under-relaxation 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.
  Reply With Quote

Old   July 16, 1999, 20:45
Default Re: Using SIMPLE
  #2
michael fitzmaurice
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.cfd-journal.com, 06/29/1999"An alternative model of the turbulent energy cascade" Michael R. Fitzmaurice.

  Reply With Quote

Old   July 17, 1999, 17:22
Default Re: Using SIMPLE
  #3
John C. Chien
Guest
 
Posts: n/a
(1). I am not going to talk about this 'not so simple SIMPLE METHOD'. (2). As you know that the Navier-Stokes equations are non-linear. To make the thing more complicated, let's take the steady-state equation for incompressible flow for example. You have a first-order continuity equation, momentum equations and energy equation of second-order. (3). Since we don't know where and how to obtain the solution for this set of steady-state 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 steady-state momentum equation with known pressure term is still non-linear. 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 "vorticity-velocity" formulation, or in 2-D case the "vorticity-stream 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 two-equation 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 pressure-based methods also are widely used in commercial codes, it might be a good idea to learn how these pressure-related methods actually work. good luck.
  Reply With Quote

Old   July 20, 1999, 06:54
Default Re: Using SIMPLE
  #4
Nuray Kayakol
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 body-fitted 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.
  Reply With Quote

Old   July 20, 1999, 22:59
Default Re: Using SIMPLE
  #5
Andrew
Guest
 
Posts: n/a
Could you give me some more info about ftp.springer.de? I got lost looking for caffa.for.

Thanks
  Reply With Quote

Old   July 21, 1999, 02:38
Default Re: Using SIMPLE
  #6
Nuray Kayakol
Guest
 
Posts: n/a
please look at /pub/technik/peric/2dgl/sg/*.for
  Reply With Quote

Reply

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


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