March 12, 2006, 14:26
|
PSOR cfd code - Convergence problem
|
#1
|
Guest
|
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
|
|
|