
[Sponsors] 
May 9, 2016, 12:04 
help with algebraic multigrid ineffectiveness

#1 
New Member
Michael Liu
Join Date: May 2016
Posts: 4
Rep Power: 2 
Hi everyone!
I'm working on an projectionmethod based solver. But when I applied AMG(Algebraic Multigrid) to the 4h layer of grid, which has nx/4 * ny/4 cells(the finest mesh is nx*ny), it didn't work. When the error of 4hgrid is prolongated back to 2hgrid, the error for 2hgrid overshoot instead of decreasing as it should. However, it worked perfectly when I only used 2hgrid. Has anyone come across such problem, that multigrid goes ineffective at certain coarse grid? Any pointer would be very helpful!! thanks! 

May 9, 2016, 13:15 

#2  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 504
Rep Power: 13 
Quote:
Also when correction is prolongated error could increase, this increase in error is quickly removed by smoother and this is what smoother's role is. So check if your smoother is working properly. 

May 11, 2016, 10:52 

#3  
New Member
Michael Liu
Join Date: May 2016
Posts: 4
Rep Power: 2 
Quote:
And yes, when I said the level, I was talking about the matrix size. Sorry I wasn't clear. And my smoother works just fine when it's working on the other level. I still didn't find any solution. And I just wonder if anyone has the same problem. I do suspect there could be a very sneaky bug, since multigrid method is garanteed theoretically. 

May 11, 2016, 11:52 

#4  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 504
Rep Power: 13 
Quote:
Well no, I at least dont have and i have written pretty much all types of popular multigrids myself. One thing I have learned over the years (which is quite simple but I still not able to get it right) is that if something works for others it shall work for us too provided we do everything right. With this in mind you boviously doing something wrong. Now that you claim smoother is working fine, then you have problem with one of these things 1. You coarse level matrix is wrong. 2. Your restriction is wrong 3. Your prolongation is wrong. last but most of the time overlooked 4. Your coarest level system is not fully solved. point number 4 is what most of the people take for granted. If you do not have 1 equation at coarsest level then make sure that you use direct solver to get the solution there. Many cases this small thing create AMG inefficiency. Good luck. Do everything right and it shall work. 

May 12, 2016, 06:25 

#5  
New Member
Michael Liu
Join Date: May 2016
Posts: 4
Rep Power: 2 
Quote:
I found out my prolongation made the error overshoot and I wasn't taking the relaxation number into account. Then I changed my relaxation number( made it smaller), and the multigrid worked. I'm not quite sure what relaxation number leads to the fastest convergence rate, though. 

May 12, 2016, 06:36 

#6  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 504
Rep Power: 13 
Quote:
It means you still have problem. Check your coarse level matrix construction. Find a way to make sure it is okay. One way to do is take a poisson problem and agglomorate as if coarsen system is coarse grid problem. (in 2d coarse 4 equations into one as if descretization is of 2H mesh if H was of finer level). In this case coarse level matrix is exactly same as if you descretiye 2H grid problem. If you do not get this coarse system then you have problem. In case you get this coarse system then more relaxations should only improve your AMG solver. 

May 13, 2016, 23:27 

#7  
New Member
Michael Liu
Join Date: May 2016
Posts: 4
Rep Power: 2 
Quote:
I've run my codes with a smaller relaxation number, and the AMG was more robust but the convergence seemed slowed down. 

May 14, 2016, 06:26 

#8  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 504
Rep Power: 13 
Quote:
yes larger relaxation number. Ideally more relaxation would improve convergence. PS: It is not completely true though. After a certain number your convergence rate would not change because smoother could only remove certain frequencies. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
algebraic multigrid vs geometric multigrid  Anna Tian  Main CFD Forum  5  June 18, 2013 05:03 
Algebraic Multigrid for CFX Transient Analysis  vagantes  CFX  3  April 30, 2013 09:36 
GAMG = geometric algebraic multigrid?  schmittp54  OpenFOAM Programming & Development  1  December 11, 2010 08:25 
Black Box Algebraic Multigrid  Neazen  Main CFD Forum  0  June 4, 2007 07:10 
Seeking Algebraic Multigrid Code  Anna  Main CFD Forum  1  May 11, 2006 20:05 