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

PCG convergency

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

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

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!

Steve
  Reply With Quote

Old   January 8, 2003, 14:10
Default Re: PCG convergency
  #2
xueying
Guest
 
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
  #3
Anthony Wachs
Guest
 
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,

Anthony
  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
Dynamic moving mesh Pei-Ying Hsieh (Hsieh) OpenFOAM Running, Solving & CFD 64 June 7, 2012 11:04
Floating point exception error Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 15:02
GAMG vs. PCG in fvSolution grandgo OpenFOAM Running, Solving & CFD 0 July 21, 2011 05:39
Velocity blows up suddenly after 30,000+ iterations lordvon OpenFOAM 14 October 27, 2010 12:13
Convergency Difficulty!! Min-Hua Wang CFX 8 November 14, 2001 14:31


All times are GMT -4. The time now is 20:55.