View source for OpenFOAM guide/The PISO algorithm in OpenFOAM

{{VersionInfo}}{{Version1.4}}{{Version1.5}}{{Version1.6}}

== The PISO algorithm ==

The PISO (Pressure Implicit with Splitting of Operators) is an efficient method to solve the Navier-Stokes equations in unsteady problems. The main differences from the [[The SIMPLE algorithm in OpenFOAM | SIMPLE]] algorithm are the following:

* No under-relaxation is applied.
* The momentum corrector step is performed more than once.

The algorithm can be summed up as follows:

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

As already seen for the [[The SIMPLE algorithm in OpenFOAM | SIMPLE]] algorithm, the steps 4 and 5 can be repeated for a prescribed number of time to correct for non-orthogonality.

== Implementation of the PISO algorithm in OpenFOAM ==

The PISO algorithm is implemented in OpenFOAM as follows (Details can be found in the icoFoam standard solver provided with OpenFOAM):

* Define the equation for U

<cpp>
 fvVectorMatrix UEqn
 (
   	fvm::ddt(U)
      + fvm::div(phi, U)
      - fvm::laplacian(nu, U)
 );
</cpp>

* Solve the momentum predictor

<cpp>
 solve (UEqn == -fvc::grad(p));
</cpp>

* Calculate the <math>a_p</math> coefficient and calculate U

<cpp>
 volScalarField rUA = 1.0/UEqn().A();
 U = rUA*UEqn().H();
</cpp>

* Calculate the flux

<cpp>
 phi = (fvc::interpolate(U) & mesh.Sf()) 
       + fvc::ddtPhiCorr(rUA, U, phi);
 adjustPhi(phi, U, p);
</cpp>

* Define and solve the pressure equation and repeat for the prescribed number of non-orthogonal corrector steps

<cpp> 
 fvScalarMatrix pEqn
 (
    fvm::laplacian(rUA, p) == fvc::div(phi)
 );
 pEqn.setReference(pRefCell, pRefValue);
 pEqn.solve();
</cpp>

* Correct the flux

<cpp>
 if (nonOrth == nNonOrthCorr)
 {
    phi -= pEqn.flux();
 }
</cpp>

* Calculate continuity errors

<cpp>
# include "continuityErrs.H"
</cpp>

* Perform the momentum corrector step

<cpp>
 U -= rUA*fvc::grad(p);
 U.correctBoundaryConditions();
</cpp>

* Repeat from the calculation of <math>a_p</math> for the prescribed number of PISO corrector steps.

== References ==

J. H. Ferziger, M. Peric, Computational Methods for Fluid Dynamics, Springer, 3rd Ed., 2001.

H. Jasak, Error Analysis and Estimation for the Finite Volume Method with Applications to Fluid Flows, Ph.D. Thesis, Imperial College, London, 1996.


[[Category:Numerical_methods]]