Velocity-pressure coupling
From CFD-Wiki
Line 9: | Line 9: | ||
we have at each cell descretised equation in this form, <br> | we have at each cell descretised equation in this form, <br> | ||
:<math> a_p \vec v_P = \sum\limits_{neighbours} {a_l } \vec v_l - \frac{{\nabla p}}{V} </math> ; <br> | :<math> a_p \vec v_P = \sum\limits_{neighbours} {a_l } \vec v_l - \frac{{\nabla p}}{V} </math> ; <br> | ||
- | we have <br> | + | According to [[Rhie-Chow interpolation]], we have <br> |
:<math> \vec v_P = \frac{{\sum\limits_{neighbours} {a_l } \vec v_l }}{{a_p }} - \frac{{\nabla p}}{{a_p V}} </math> <br> | :<math> \vec v_P = \frac{{\sum\limits_{neighbours} {a_l } \vec v_l }}{{a_p }} - \frac{{\nabla p}}{{a_p V}} </math> <br> | ||
Revision as of 00:07, 24 October 2005
If we consider the discretised form of the Navier-Stokes system, the form of the equations shows linear dependence of velocity on pressure and vice-versa. This inter-equation coupling is called velocity pressure coupling. A special treatment is required in order to velocity-pressure coupling. The methods such as:
- SIMPLE
- SIMPLER
- SIMPLEC
- PISO
provide an useful means of doing this for segregated solvers. However it is possible to solve the system of Navier-Stokes equations in coupled manner, taking care of inter equation coupling in a single matrix.
Contents |
Formulation
we have at each cell descretised equation in this form,
- ;
According to Rhie-Chow interpolation, we have
For continuity :
so we get:
this gives us:
defining
from this a pressure correction equation could be formed as:
This is a poisson equation.
Here the gradients could be used from previous iteration.
SIMPLE
See SIMPLE algorithm
SIMPLER
SIMPLEC
PISO
See PISO algorithm
Return to: