PISO algorithm

<h2 id="algorithm-steps">Algorithm steps</h2>

<p>The algorithm can be summed up as follows:
</p>
<ol><li>Set the boundary conditions.</li>
<li>Solve the discretized momentum equation to compute an intermediate velocity field.</li>
<li>Compute the mass fluxes at the cells faces.</li>
<li>Solve the pressure equation.</li>
<li>Correct the mass fluxes at the cell faces.</li>
<li>Correct the velocities on the basis of the new pressure field.</li>
<li>Update the boundary conditions.</li>
<li>Repeat from 3 for the prescribed number of times.</li>
<li>Increase the time step and repeat from 1.</li></ol>
<p>Steps 4 and 5 can be repeated for a prescribed number of times to correct for non-orthogonality.
</p><p>Predictor step
</p><p>Guess the pressure field 
  
    
      
        
          p
          
            ∗
          
        
      
    
    {\displaystyle p^{*}}
  
 and get velocity field components 
  
    
      
        
          u
          
            ∗
          
        
      
    
    {\displaystyle u^{*}}
  
 and 
  
    
      
        
          v
          
            ∗
          
        
      
    
    {\displaystyle v^{*}}
  
 using discretized momentum equation. The initial guess for the pressure may or may not be correct. 
Corrector step 1Velocity component obtained from predictor step may not satisfy the continuity equation, so we define correction factors p',v',u'  for the pressure field and velocity field. Solve the momentum equation by inserting correct pressure field 
  
    
      
        
          p
          
            ∗
            ∗
          
        
      
    
    {\displaystyle p^{**}}
  
 and get the corresponding correct velocity components 
  
    
      
        
          u
          
            ∗
            ∗
          
        
      
    
    {\displaystyle u^{**}}
  
 and 
  
    
      
        
          v
          
            ∗
            ∗
          
        
      
    
    {\displaystyle v^{**}}
  
.
</p><p>
  
    
      
        
          p
          ′
        
        =
        
          p
          
            ∗
            ∗
          
        
        −
        
          p
          
            ∗
          
        
      
    
    {\displaystyle p'=p^{**}-p^{*}}

</p><p>
  
    
      
        
          v
          ′
        
        =
        
          v
          
            ∗
            ∗
          
        
        −
        
          v
          
            ∗
          
        
      
    
    {\displaystyle v'=v^{**}-v^{*}}

u
          ′
        
        =
        
          u
          
            ∗
            ∗
          
        
        −
        
          u
          
            ∗
          
        
      
    
    {\displaystyle u'=u^{**}-u^{*}}

where ;
  
    
      
        
          p
          
            ∗
            ∗
          
        
        ,
        
          u
          
            ∗
            ∗
          
        
        ,
        
          v
          
            ∗
            ∗
          
        
      
    
    {\displaystyle p^{**},u^{**},v^{**}}
  
     :correct pressure field and velocity component

p
          ′
        
        ,
        
          u
          ′
        
        ,
        
          v
          ′
        
      
    
    {\displaystyle p',u',v'}
  
 :correction in pressure field and correction in velocity components

p
          
            ∗
          
        
        ,
        
          u
          
            ∗
          
        
        ,
        
          v
          
            ∗
          
        
      
    
    {\displaystyle p^{*},u^{*},v^{*}}
  
 :guessed pressure field and velocity component
We define 
  
    
      
        
          p
          ′
        
        ,
        
          u
          ′
        
        ,
        
          v
          ′
        
      
    
    {\displaystyle p',u',v'}
  
 as above.
By putting the correct pressure field 
  
    
      
        
          p
          
            ∗
            ∗
          
        
      
    
    {\displaystyle p^{**}}
  
 into the discretized momentum equation we get the correct velocity components 
  
    
      
        
          v
          
            ∗
            ∗
          
        
      
    
    {\displaystyle v^{**}}
  
 and 
  
    
      
        
          u
          
            ∗
            ∗
          
        
      
    
    {\displaystyle u^{**}}
  
. Once the pressure correction 
  
    
      
        
          p
          ′
        
      
    
    {\displaystyle p'}
  
 is known we can find the correction components for the velocity: 
  
    
      
        
          u
          ′
        
      
    
    {\displaystyle u'}
  
 and 
  
    
      
        
          v
          ′
        
      
    
    {\displaystyle v'}
  
.
</p><p>Corrector step 2
In piso another corrector step can be used.

p
          
            ∗
            ∗
            ∗
          
        
        =
        
          p
          
            ∗
            ∗
          
        
        +
        
          p
          ″
        
      
    
    {\displaystyle p^{***}=p^{**}+p''}
  
    ;    
  
    
      
        
          p
          ″
        
        =
        
          p
          
            ∗
          
        
        +
        
          p
          ′
        
      
    
    {\displaystyle p''=p^{*}+p'}

u
          
            ∗
            ∗
            ∗
          
        
        =
        
          u
          
            ∗
            ∗
          
        
        +
        
          u
          ″
        
      
    
    {\displaystyle u^{***}=u^{**}+u''}
  
     ;     
  
    
      
        
          u
          ″
        
        =
        
          u
          
            ∗
          
        
        +
        
          u
          ′
        
      
    
    {\displaystyle u''=u^{*}+u'}

v
          
            ∗
            ∗
            ∗
          
        
        =
        
          v
          
            ∗
            ∗
          
        
        +
        
          v
          ″
        
      
    
    {\displaystyle v^{***}=v^{**}+v''}
  
     ;    
  
    
      
        
          v
          ″
        
        =
        
          v
          
            ∗
          
        
        +
        
          v
          ′
        
      
    
    {\displaystyle v''=v^{*}+v'}

where : 
  
    
      
        
          p
          
            ∗
            ∗
            ∗
          
        
        ,
        
          v
          
            ∗
            ∗
            ∗
          
        
        ,
        
          u
          
            ∗
            ∗
            ∗
          
        
      
    
    {\displaystyle p^{***},v^{***},u^{***}}
  
 are the correct pressure field and the correct velocity components, respectively
and 
  
    
      
        
          p
          ″
        
        ,
        
          v
          ″
        
        ,
        
          u
          ″
        
      
    
    {\displaystyle p'',v'',u''}
  
 are second corrections to the pressure and velocity field.
Set 
  
    
      
        p
        =
        
          p
          
            ∗
            ∗
            ∗
          
        
        ,
        v
        =
        
          v
          
            ∗
            ∗
            ∗
          
        
        ,
        u
        =
        
          u
          
            ∗
            ∗
            ∗
          
        
      
    
    {\displaystyle p=p^{***},v=v^{***},u=u^{***}}
  
 
where; 
  
    
      
        p
        ,
        v
        ,
        u
      
    
    {\displaystyle p,v,u}
  
 are correct pressure and velocity field
</p>
<h2 id="advantages-and-disadvantages">Advantages and disadvantages</h2>
<ol><li>Generally gives more stable results and takes less CPU time but not suitable for all processes.</li>
<li>Suitable numerical schemes for solving the pressure-velocity linked equation.</li>
<li>For laminar backward facing step PISO is faster than SIMPLE but it is slower concerning flow through heated fin.</li>
<li>If momentum and scalar equation have weak or no coupling then PISO is better than SIMPLEC.</li>
<li>PISO is most time effective method</li></ol>
<h2 id="see-also">See also</h2>
<ul><li><a href="/facts/Fluid_mechanics/NIy4tRfE">Fluid mechanics</a></li>
<li><a href="/facts/Computational_fluid_dynamics/zuWdtNJp">Computational fluid dynamics</a></li>
<li><a href="/facts/Algorithm/fnl5NmRt">Algorithm</a></li>
<li><a href="/facts/SIMPLE_algorithm/Wo36xVHU">SIMPLE algorithm</a></li>
<li><a href="/facts/SIMPLER_algorithm/Wo36xVHU">SIMPLER algorithm</a></li>
<li><a href="/facts/SIMPLEC_algorithm/YqP5Bqjt">SIMPLEC algorithm</a></li></ul>

<ol><li>An Introduction to Computational Fluid Dynamics The Finite Volume Method, 2/e By Versteeg  <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0131274983</li>
<li>Computational Fluid Dynamics for Engineers by Bengt Andersson, Ronnie Andersson, Love Håkansson, Mikael Mortensen, Rahman Sudiyo, Berend van Wachem <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-1-107-01895-2</li>
<li>Computational Fluid Dynamics in Fire Engineering: Theory, Modelling and Practice by Guan Heng Yeoh, Kwok Kit Yuen <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0750685894</li>
<li><a href="http://openfoamwiki.net/index.php/OpenFOAM_guide/The_PISO_algorithm_in_OpenFOAM">http://openfoamwiki.net/index.php/OpenFOAM_guide/The_PISO_algorithm_in_OpenFOAM</a></li>
<li>Computational fluid dynamics by T. J. Chung, University of Alabama in Huntsville <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 0 521 59416 2</li>
<li>Computational method for fluid dynamics by Joel H.Ferziger, Milovan Peric <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 3-540-42074-6</li>
<li>Solution of the implicitly discretized fluid flow equations by operator-splitting, <i>Journal of Computational Physics</i> 62 by R. Issa</li></ol>

PISO algorithm open-in-new

PISO algorithm