PISO algorithm

PISO algorithm (Pressure implicit with splitting of operator) proposed by Issa in 1985 without iterations and with large time steps and less computing effort. it is extension of SIMPLE algorithm which use in CFD computational fluid dynamics to solve Navier-stock equation.PISO Is a pressure velocity calculation procedure for Navier-stock equation developed originally for noniterative computation of unsteady compressible flow.but it has been adapted successfully on steady state problems.

PISO involves one predictor step and two corrector steps and In this conservation of mass is design to satisfy within the predictor corrector steps.

Algorithm Steps

The algorithm can be summed up as follows:

1. Set the boundary conditions.
   
2. Solve the discretized momentum equation to compute an intermediate velocity field.
   
3. Compute the mass fluxes at the cells faces.
   
4. Solve the pressure equation.
   
5. Correct the mass fluxes at the cell faces.
   
6. Correct the velocities on the basis of the new pressure field.
   
7. Update the boundary conditions.
   
8. Repeat from 3 for the prescribed number of times.
   
9. Increase the time step and repeat from 1.
   

As already seen for the SIMPLE algorithm, the steps 4 and 5 can be repeated for a prescribed number of time to correct for non-orthogonality.
predictor step
guess the pressure field p* and get velocity field component u* and v* using Discretized momentum equation.here gussed pressure may be correct or not.
corrector step 1
velocity component obtain from predictor step not be satisfy in continuity equation ,so we define correction factor p',v',u' pressure field and velocity field. solve momentum equation by put correct pressure p** field and get corresponding correct velocity component u** and v**.
p'=p**-P*
v'=v**-v*
u'=u**-u*
where ;
p**,v**,u** :correct pressure field and velocity component
p',v',u' :correction in pressure field and correction in velocity component
p*,v*,u* :guessed pressure field and velocity component
we define as above p',u',v'. by put correct term of pressure field p** in discretized momentum equation we get correct velocity component v**&u**.once pressure correction p' known we can find v' and u' correction component of velocity.

corrector step2 in piso another corrector step .
p***=p**+p"  ; p"=p*+p'
u***=u**+u"  ; u"=u*+u'
v***=v**+v"  ; v"=v*+v'
where : P***,v***,u*** are respectively correct pressure field and correct velocity component
and p",v",u" are second correction pressure and velocity field.
set p=p***
v=v***
u=u***

where ; p,v,u are correct pressure and velocity field

Flow chart of PISO algorithm

advantage and disadvantage

1. generally gave more stable results and took less CPU time but not suitable for all process.
2. suitable numerical schemes for solving the pressure-velocity linked equation.
3. for laminar backward facing step piso faster than SIMPLE but it slower for concerning flow through heated fin.
4. if momentum and scalar equation have weak or no coupling than PISO better than SIMPLEC.

Reference

1. An Introduction to Computational Fluid Dynamics The Finite Volume Method, 2/e By Versteeg ISBN 978-0131274983
2. Computational Fluid Dynamics for Engineers By Bengt Andersson Ronnie Andersson Love Håkansson Mikael Mortensen Rahman Sudiyo Berend van Wachem ISBN 978-1-107-01895-2
3. Computational Fluid Dynamics in Fire Engineering: Theory, Modelling and Practice By Guan Heng Yeoh, Kwok Kit Yuen ISBN 978-0750685894
4. http://openfoamwiki.net/index.php/OpenFOAM_guide/The_PISO_algorithm_in_OpenFOAM