# SIMPLE algorithm

### From CFD-Wiki

(Difference between revisions)

m (Fixed the meaning of the SIMPLE acronym) |
|||

Line 1: | Line 1: | ||

- | ==SIMPLE [Semi Implicit Method Pressure Linked | + | ==SIMPLE [Semi-Implicit Method for Pressure-Linked Equations]== |

If a steady-state problem is being solved iteratively, it is not necessary to fully resolve | If a steady-state problem is being solved iteratively, it is not necessary to fully resolve | ||

Line 16: | Line 16: | ||

#Set the boundary conditions. | #Set the boundary conditions. | ||

- | # | + | #Compute the gradients of velocity and pressure. |

#Solve the discretized momentum equation to compute the intermediate velocity field . | #Solve the discretized momentum equation to compute the intermediate velocity field . | ||

#Compute the uncorrected mass fluxes at faces . | #Compute the uncorrected mass fluxes at faces . |

## Revision as of 18:41, 4 October 2006

## SIMPLE [Semi-Implicit Method for Pressure-Linked Equations]

If a steady-state problem is being solved iteratively, it is not necessary to fully resolve the linear pressure-velocity coupling, as the changes between consecutive solutions are no longer small. The SIMPLE algorithm:

- An approximation of the velocity field is obtained by solving the momentum equation. The pressure gradient term is calculated using the pressure distribution from the previous iteration or an initial guess.
- The pressure equation is formulated and solved in order to obtain the new pressure distribution.
- Velocities are corrected and a new set of conservative fluxes is calculated.

## SIMPLE Solver Algorithm

The algorithm may be summarized as follows:

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.

* Return to Numerical Methods *