CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

GMRES-solver does not converge correctly

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   February 13, 2013, 18:08
Default GMRES-solver does not converge correctly
  #1
New Member
 
Join Date: Feb 2013
Posts: 1
Rep Power: 0
staentz_b is on a distinguished road
Hi there,

my question is quite specific, but i hope someone can help me out.
I have serious trouble with my GMRES-solver.
This solver has to work on a poisson equation, which arises from a finite-volume 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 Jacobi-solver and a BiCGSTAB-solver 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 GMRES-algorithm.


How can that be?
Is there any solution?
Ever heared of such a problem?

Thanks a lot for any kind of help
staentz_b is offline   Reply With Quote

Old   February 14, 2013, 05:19
Default
  #2
Senior Member
 
andy
Join Date: May 2009
Posts: 129
Rep Power: 6
andy_ is on a distinguished road
> 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.
andy_ is offline   Reply With Quote

Old   February 14, 2013, 06:55
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 1,000
Rep Power: 15
FMDenaro will become famous soon enough
Quote:
Originally Posted by staentz_b View Post
Hi there,

my question is quite specific, but i hope someone can help me out.
I have serious trouble with my GMRES-solver.
This solver has to work on a poisson equation, which arises from a finite-volume 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 Jacobi-solver and a BiCGSTAB-solver 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 GMRES-algorithm.


How can that be?
Is there any solution?
Ever heared of such a problem?

Thanks a lot for any kind of help

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 ...
FMDenaro is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Interfoam blows on parallel run danvica OpenFOAM Running, Solving & CFD 16 December 22, 2012 02:09
Unexplained Error during Solver Runs cfb CFX 6 November 9, 2012 15: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 20:02


All times are GMT -4. The time now is 13:52.