
[Sponsors] 
February 13, 2013, 19:08 
GMRESsolver does not converge correctly

#1 
New Member
Join Date: Feb 2013
Posts: 1
Rep Power: 0 
Hi there,
my question is quite specific, but i hope someone can help me out. I have serious trouble with my GMRESsolver. This solver has to work on a poisson equation, which arises from a finitevolume discretization. The matrix of the linear system is symmetric and positiv definite. The solver does minimize the residuum, the solution does fullfill the linear system with a certain error, which is extremely small, the residuum has values to 10⁻15 or so... But it is not the correct solution (i tested it with an alaytical solution, a Jacobisolver and a BiCGSTABsolver do work on that problem). Plus the solution changes when changing the starting vector, the solution is always "near" to the starting vector. The algorithm definetly works somehow, i calculated the "real" residuum at each iterationstep and it is aquivalent with the residuum, which is calculated by the GMRESalgorithm. How can that be? Is there any solution? Ever heared of such a problem? Thanks a lot for any kind of help 

February 14, 2013, 06:19 

#2 
Senior Member
andy
Join Date: May 2009
Posts: 129
Rep Power: 9 
> How can that be?
What is the determinant? Conclusion? > Is there any solution? Yes but with a qualification. If you plot your different solutions it should be clear what is happening. > Ever heared of such a problem? This is just one example of the many practical insights one gets working with CFD. It is of course much easier performing exercises like this within an experienced CFD group rather than posting to online chat groups. 

February 14, 2013, 07:55 

#3  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,586
Rep Power: 32 
Quote:
but maybe you are solving for the pressure equation? because your solution will be defined apart a constant... and, if you do not fix the value, such value is fixed somehow by the initial vector ... 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Interfoam blows on parallel run  danvica  OpenFOAM Running, Solving & CFD  16  December 22, 2012 03:09 
Unexplained Error during Solver Runs  cfb  CFX  6  November 9, 2012 16:42 
Strange residuals of the Density Based Solver  Pat84  FLUENT  0  October 22, 2012 15:59 
Quarter Burner mesh with periosic condition  SamCanuck  FLUENT  2  August 31, 2011 11:34 
Working directory via command line  Luiz  CFX  4  March 6, 2011 21:02 