CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

PCG convergency

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

LinkBack Thread Tools Display Modes
Old   January 8, 2003, 12:42
Default PCG convergency
Steve Song
Posts: n/a

I am developing FEM codes for 3-D incompressible fluid flows.

My current problem is that, for large scale problems (1M nodes), LU method is running extremely slow. It takes two days for one LU decomposition with FEM penalty method. An iterative solver (preconditioned conjugate gradient) was considered to be a solution. However, this case is difficult to converge especially in 3-D if using this iterative solver. I am not sure if the problem is caused by the big number (inverse of penalty parameter) in the matrix.

Any ideas and suggestions? Many Thanks!

  Reply With Quote

Old   January 8, 2003, 14:10
Default Re: PCG convergency
Posts: n/a
I think usually the convergence is related with the problem itself, not the size. You can try a small problem to see if it converge or not. If not, try to load the matrix to MATLAB to check the eigenvalues. I use GMRES for usual 3D flow. It converges. But if I add free surface boundary condition, it doesn't converge.
  Reply With Quote

Old   January 16, 2003, 04:47
Default Re: PCG convergency
Anthony Wachs
Posts: n/a
The convergence speed of the PCG algorithm is directly related to the condition number of the matrix of the linear system.

In the penalty method, the greater epsilon (in div(v)=-(1/epsilon)*p), the greater the condition number and the worse the convergence rate of the method.

At that point there is nearly nothing else trivial to do, as far as I know, than trying to enhance the preconditioning of the matrix.

Otherwise, you can switch to augmented lagrangian methods.

Hope it could help,

  Reply With Quote


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
GAMG vs. PCG in fvSolution grandgo OpenFOAM Running, Solving & CFD 2 March 7, 2017 08:38
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM Running, Solving & CFD 15 October 19, 2015 13:52
Dynamic moving mesh Pei-Ying Hsieh (Hsieh) OpenFOAM Running, Solving & CFD 64 June 7, 2012 10:04
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02
Convergency Difficulty!! Min-Hua Wang CFX 8 November 14, 2001 14:31

All times are GMT -4. The time now is 06:19.