PSOR cfd code - Convergence problem
 

hi, I have been trying to write a cfd code(in C++) to implement point successive over-relaxation method(PSOR) in a 2D duct flow with a constant pressure gradient. I have calcuated my relaxation factor many times (which is w = 1.87), however this factor leads to indeterminite results in my code. I have also written a code to implement Point Gauss Seidel scheme, which works correctyly, and I just added the following components to Point Gauss Seidel scheme: u(i,j) = (1 - w)*u(i, j)~ + w * PGauss Seidel solution

~ denotes previous iteration value

this application always diverges before it reaches 100 iterations. however solution converges to correct value when w = 1.0(which turns into Poıint Gauss Seidel)..

if anyone can help me with the problem, I'll be very glad Thanks, Taskin