Algorithm
The algorithm is iterative. The basic steps in the solution update are as follows:
- Set the boundary conditions.
- Compute the gradients of velocity and pressure.
- Solve the discretized momentum equation to compute the intermediate velocity field .
- Compute the uncorrected mass fluxes at faces .
- Solve the pressure correction equation to produce cell values of the pressure correction .
- Update the pressure field: where urf is the under-relaxation factor for pressure.
- Update the boundary pressure corrections .
- Correct the face mass fluxes:
- Correct the cell velocities: ; where is the gradient of the pressure corrections, is the vector of central coefficients for the discretized linear system representing the velocity equation and Vol is the cell volume.
- Update density due to pressure changes.
Read more about this topic: SIMPLE Algorithm