
[Sponsors] 
March 23, 2006, 05:23 
SIMPLE and PISO

#1 
Guest
Posts: n/a

I have a few doubts of implementing SIMPLE and PISO fro transient problems.
1. What do we mean when we say SIMPLE is iterative and PISO is noniterative. Does that indicate that the Momentum and continuity equations are simultaneously satisfied after the "second pressure correction" in PISO? 2, If that is so how does SIMPLE take care of the transient problems. By giving a small delta t and iterating to get converged solutions at that delta t? 3. In PISO why do we use the "correction" mass flux instead of the corrected one (I am not sre whther i am right on that) as the source term of the second pressure corrction equation. Does it not destroy the continuity?? Please answer. Aditya 

March 25, 2006, 11:08 
Re: SIMPLE and PISO

#2 
Guest
Posts: n/a

Hi!
1. That's right, both, momentum and continuity should be fullfilled after the second corrector step if the time step is small enough (Courant Number lower than 1!). You can use more corrector steps if necessary (mostly not more than four). You can derive these equations in the same way like the second correction equation and this answers question 3: if correctly derived you have on the right hand side the divergence of correction massfluxes not of "real massfluxes". Note that in the first corrector step you correct the velocities from the momentum equations in the second corrector step you correct the velocity corrections!!! 2. You already answered! You can also derive a simple relationship between the underrelaxation for momentum equation you need for the iterative process within a time step for SIMPLE and the necessary time for the same calculation with PISO. For this subject try also the following literature: 1) Ferziger and Peric C. Methods for F. Dynamics 2) The papers from the inventor of PISO Issa (Google) 3) The PhD of Mr. Hrvoje Jasak (wrote already very often in this place). Good Luck! 

March 27, 2006, 00:34 
Re: SIMPLE and PISO

#3 
Guest
Posts: n/a

Thank you newposter My doubts 2 and 3 ar now perfectly cleared. However I am still a bit confused on doubt 1. I am referrin Christian Michelsen Research report in which it is clearly stated that PISO is an iterative algorithm. While when I referred Issa's paper, it was explicitly stated that PISO is a noniterative algorithm which makes more sense. So I am bit confused about the nomenclature used by the two authors.
Regards and thanks Aditya 

March 27, 2006, 13:42 
Re: SIMPLE and PISO

#4 
Guest
Posts: n/a

Hi,
I think originally it was developed as a noniterative method. Although many people write that they are using PISO in an iterative variant like SIMPLE or SIMPLEC with a much better convergence rate. I am also using it in the most cases in this variant. The problem in the noniterative variant is the need for a very small timestep! Good Luck! And a question to you: How are you discretizing and calculating right hand side of your PISO (divergence of massfluxcorrection)? Are you using collocated variable arragement? 

March 28, 2006, 00:43 
Re: SIMPLE and PISO

#5 
Guest
Posts: n/a

Yes the Rhie and Chow momentum Interpolation Staggered gris are more or less obsolete I think.
I have a small question about the timestep. I have incompressible equations on my hands. And depending upon whether unsteady or steady my a_p , a_nb and the source term (ref. Prof. Jayathi Murthy's notes) change. So in case of steady equations I nowhere use a timestep. So for using PISO for steady incompressible equations, I dont have an input of time step. Does this sound fundamentally incorrect? Since whenever PISO is mentioned, timestep comes up invariably. By PISO what I do is to correct the first velocity correction at cell through the sum_(a_nb * phi_nb')/a_p which SIMPLE algorithm ignores. Please let me know. Thanks and regards, Aditya 

April 1, 2006, 18:52 
Re: SIMPLE and PISO

#6 
Guest
Posts: n/a

Hi,
That's ok that your APs are changing in the steady formulation of PISO. But then you are running this algorithm like SIMPLE with outer Iterations and not in the original version. It can be faster than SIMPLE or SIMPLEC but must not! Your second question: You are solving a SIMPLE pressure correction equation with the divergence of massfluxes on the right hand side but with neglecting the corrections of neighbourcells. Then you are correcting massfluxes, velocities and pressurefield. After this you have to solve a second pressure correction with the divergence of (a_nb * u_nb')/a_p on the right hand side. You get second pressurecorrections p'' and you have to correct then your pressure field again and your massfluxes and velocities again! After this your pressure field (at the new time step) is: p=p*+p'+p'' Your velocites are now (for example eastface): u=u*dV/A_P*dp'/dxdV/A_P*dp''/dx+(A_nb * u_nb')/A_P And your massfluxes: m=m*A_E*(p'_Ep'_P)A_E*(p''_Ep''_P)+S_e*(A_nb * u_nb')/A_P * means the solution you got after the first momentum solution, S_e means surface on eastside! But look also in literature to derive this equations! Good luck! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Continuity. Stop critertion of PISO and SIMPLE in OpenFOAM solvers.  idrama  OpenFOAM  4  September 21, 2010 03:26 
SIMPLE, SIMPLER or PISO?  Gem Leku  Main CFD Forum  7  October 24, 2007 11:20 
diffrence beetwen SIMPLE and PISO  grzes  FLUENT  0  January 29, 2006 06:57 
PISO vs. SIMPLE  benedikt flurl  Main CFD Forum  2  April 14, 2005 06:54 
advantages/disadvantages of PISO, SIMPISO, SIMPLE  Joel  CDadapco  0  April 16, 2004 06:34 