CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Main CFD Forum (
-   -   Convergence pressure correction equation (

michujo September 24, 2012 05:45

Convergence pressure correction equation
Dear all, I finally got my little CFD code to work (FVM, SIMPLE method, steadyflow, 2D incompressible code for the lid-driven cavity problem). However I'm experiencing some convergence issues.

I'm using upwind discretization for the convective terms and Gauss-Seidel + relaxation (SOR) solution method. I am having a hard time making the pressure correction to converge.

In each iteration of the SIMPLE method loop I update the velocity and pressure field like this:


, where u and v are the velocity components, p is the pressure, ustar, vstar and pstar refer to the guessed velocity and pressure fields and alpha_p is a relaxation coefficient I use to avoid instabilities.

For a 20cell x 20cell grid and Reynolds number 100 I am obliged to use an extremely low value for alpha_p (alpha_p = 0.01).

Did anyone found the same problem? Do you think this is a normal value for the relaxation coefficient? If you computed the same problem, which value of alpha_p did you use?

Thank you in advance for your help.


P.S: The code seems to work well for coarse grids, but it seems impossible to me to compute larger grids, it does not converge.

michujo September 25, 2012 03:17

Sorry, I just wanted to bring this up again.

Any hint from anyone?



All times are GMT -4. The time now is 08:33.