
[Sponsors] 
September 14, 2005, 10:44 
AMG method?

#1 
Guest
Posts: n/a

Hi
I can't find any information about method using for AMG technique. I mean, is it GaussSeidel method or CG or something else? As i understand, AMG is general technique of solving system of linear equation, but it should includes methods mentioned above. If my opinion is wrong, please correct me Thanks in advance 

September 14, 2005, 17:42 
Re: AMG method?

#2 
Guest
Posts: n/a

AMG means algebraic multi grid. It's a method to solve equations systems effectively. You solve your equation system on different meshes (fin > coarse). Since you have different frequency of errors on differ mesh density, you can solve low frequency errors on coarse grids more efficently. Mayby you should have a look at the theory behind it. Look for litterature by Achi Brandt and Hackbusch. It normally needs a "pre conditioner" e.g. CG (used in STARCD). If you can solve your equation system farily good after only a few CG iterations, it would be too expensive to start a AMG process, hence in STARCD AMG is only applied for the pressure equation which normally needs a lot of itersations. And yes, it is different from Gauss Seidel or CG.


September 15, 2005, 13:10 
Re: AMG method?

#3 
Guest
Posts: n/a

Thank you very much for your detailed answer!
Could you, please, answer some more questions 1. You wrote: "You solve your equation system on different meshes (fin > coarse)" But what the method is used for solving the system of linear equation on, for example, coarse grid in STARCD? Is it any of iterative methods or direct? 2. Could you give any link for Achi Brandt and Hackbusch books? Thanks in advance 

September 15, 2005, 15:43 
Re: AMG method?

#4 
Guest
Posts: n/a

1) It's allways an iterative method. Direct methods would need too much memory. The coarser meshes are virtual and are created authomatically by STARCD, hence you need to allocate more memory during the solver process. You can observe this in v3.15a, where the approximate memory usage is give after "starlink x..."; during the run, more memory will be allocated.
2) http://www.wisdom.weizmann.ac.il/~achi/ and (e.g.) Hackbush, W., MultiGrid Methods and Applications, SpringerVerlag, Berlin(1985). Look for "Hackbush multi" at google. 

September 15, 2005, 18:01 
Re: AMG method?

#5 
Guest
Posts: n/a

imagine that you take some neighbouring cells and join them and make one single volume. So this way you can create coarser control volume(u get a coarse mesh). Now you can again merge these coarse volumes and create another coarser mesh. You can do it and create some coarser meshes, say 45 meshes.
A typical multigrid cycle starts at a finest level. The fine level solution is then transferred to next coarser level, this is called restriction. After some relaxation cycles this coarse level the solution is then restricted to next coarser level until the coarsest level is reached. The solution obtained at the coarsest level is than interpolated back to finer level, this step is called prolongation. The solution from this finer level is interpolated to next finer level after some relaxation iterations, called post multigrid sweeps. The solution is prorogated till the finest level is reached. The whole process is repeated until satisfactory convergence is reached. 

September 16, 2005, 11:19 
Re: AMG method?

#6 
Guest
Posts: n/a

Thanks a lot for yours explanations! its help me very much!
But I don't understand, what concrete iterative method is used in AMG algorithm ( within STARCD) ? For example, in FLUENT documentation I read: " FLUENT solves this linear system using a point implicit (GaussSeidel) linear equation solver in conjunction with an algebraic multigrid (AMG) method" 

September 16, 2005, 16:23 
Re: AMG method?

#7 
Guest
Posts: n/a

i can not speak of starCD, but starCCM+ also uses gaussseidel, that it, at each level you solve by gaussseidel solver. Its easy to use with unstructured grids. Fluent also gives an option of using BiCGStab.


September 20, 2005, 17:14 
Re: AMG method?

#8 
Guest
Posts: n/a

STARCD uses CG (conjugate gradient) for preconditioning.


September 20, 2005, 22:08 
Re: AMG method?

#9 
Guest
Posts: n/a

conjugate gradient are for symmetric metrices, that is for pressure correction equation, this preconditioner is valid, what does it do for other scalars, does it use biCGStab, or anyother similar method?????


September 21, 2005, 09:56 
Re: AMG method?

#10 
Guest
Posts: n/a

Multigrid is a accelleration scheme / procedure for a wide range of quit normal matrix solvers, first of all Gauss Seidel Relaxation. This normal matrix solver is called in multigrid a "smoother".
"Numerical Recipes" gives a good and relativly short explanation (and example fortran and c sources) of what is multigrid. To explain it here even roughly is not possible. If you can understand german language I can give you a documantation, which very easy to understand. But you don't realy need this to know. What you should know about multigird is:  only reasonable for steady state problems.  in 2D solver acceleration of possibly 100 times (depends on the mesh size)  in 3D acceleration of about 5 to 10 times in optimal configuration and for very big meshes. In StarCD perhaps only 1.5 times (about 1 mill. cells, with good luck). Again this depends strongly on the mesh size.  has some stability advantage for compressible flows.  nearly linear convergence behaviour, so it is almost independend of the mesh size (dependent on the smoother). Therefore it becomes more and more advantagous for bigger and bigger meshes. In the past it was good to overcome weak CPUperformance, but in the future it will be necessary to get convergence at all on extremly big meshes, because every solver becomes weak for a certain mesh size. E.g. the limitation of Gauss Seidel is about 2 mio cells. The solution progress is so small that the values of the cells react mainly only to the changes of the neihgbour cell. Milovan Peric and Muzaferia are experts in applied Multigrid. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Comparing between the Fractional step method and the SIMPLE method  ghlee  Main CFD Forum  1  April 10, 2012 16:59 
Finite volume method vs finite difference method?  superfool  Main CFD Forum  4  October 21, 2006 14:37 
Pressure Correction method, Finite Volume Method  Seeker01  Main CFD Forum  2  January 13, 2003 03:49 
hesssmith method and fvm method  yangqing  FLUENT  0  March 20, 2002 20:25 
Projection method and Blockoff method  Leo  Main CFD Forum  0  June 14, 2001 07:22 