# PCG convergency

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

 January 8, 2003, 12:42 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

 January 8, 2003, 14:10 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.

 January 16, 2003, 04:47 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post grandgo OpenFOAM Running, Solving & CFD 1 November 13, 2015 09:58 lordvon OpenFOAM Running, Solving & CFD 15 October 19, 2015 13:52 Pei-Ying Hsieh (Hsieh) OpenFOAM Running, Solving & CFD 64 June 7, 2012 10:04 Alan OpenFOAM Running, Solving & CFD 10 April 6, 2012 14:02 Min-Hua Wang CFX 8 November 14, 2001 14:31

All times are GMT -4. The time now is 16:35.