Difference between revisions of "OpenFOAM guide/The SIMPLE algorithm in OpenFOAM"
(→The SIMPLE algorithm) |
|||
| Line 13: | Line 13: | ||
The solution of this couple of equations is not straightforward because an explicit equation for the pressure is not available. One of the most common approaches is to derive an equation for the pressure by taking the divergence of the momentum equation and by substituting it in the continuity equation. | The solution of this couple of equations is not straightforward because an explicit equation for the pressure is not available. One of the most common approaches is to derive an equation for the pressure by taking the divergence of the momentum equation and by substituting it in the continuity equation. | ||
| − | The momentum equation can be re-written in a | + | The momentum equation can be re-written in a semi discretized form as follows: |
<math> | <math> | ||
a_p \vec{U_p} = H(\vec{U}) - \nabla p \Longleftrightarrow \vec{U_p} = \frac{H(\vec{U})}{a_p} - \frac{\nabla p}{a_p} | a_p \vec{U_p} = H(\vec{U}) - \nabla p \Longleftrightarrow \vec{U_p} = \frac{H(\vec{U})}{a_p} - \frac{\nabla p}{a_p} | ||
</math> | </math> | ||
| + | |||
| + | where | ||
| + | |||
| + | <math> | ||
| + | H(\vec{U}) = - \sum_n a_n \vec{U}_n + \frac{\vec{U}^o}{\Delta t} | ||
| + | </math> | ||
| + | |||
| + | The first term of <math>H(\vec{U})</math> represents the matrix coefficients of the neighbouring cells multiplied by their velocity, while the second part contains the unsteady term and all the sources except the pressure gradient. | ||
| + | |||
| + | The continuity equation is discretized as: | ||
| + | |||
| + | <math> | ||
| + | \nabla \cdot \vec{U} = \sum_f \vec{S} \vec{U}_f | ||
| + | </math> | ||
| + | |||
| + | where <math>\vec{S}</math> is outward-pointing face area vector and <math>U_f</math> the velocity on the face. | ||
== Implementation of the SIMPLE algorithm in OpenFOAM == | == Implementation of the SIMPLE algorithm in OpenFOAM == | ||
== References == | == References == | ||
Revision as of 13:18, 12 March 2007
1 The SIMPLE algorithm
The Navier-Stokes equations for a single-phase flow with a constant density and viscosity are the following:
The solution of this couple of equations is not straightforward because an explicit equation for the pressure is not available. One of the most common approaches is to derive an equation for the pressure by taking the divergence of the momentum equation and by substituting it in the continuity equation.
The momentum equation can be re-written in a semi discretized form as follows:
where
The first term of 
represents the matrix coefficients of the neighbouring cells multiplied by their velocity, while the second part contains the unsteady term and all the sources except the pressure gradient.
The continuity equation is discretized as:
where 
is outward-pointing face area vector and 
the velocity on the face.
