
[Sponsors] 
December 1, 2005, 14:18 
MULTIGRID 3D

#1 
Guest
Posts: n/a

Hi, every body I've implemented a 3D Multigrid with Neumann and Periodic Boundary condition. I solve a Poisson equation on an irregular mesh of 65*65*65 with finite difference scheme.
I obtain convergence in 12 MultiGrid cycles. Is it a good performance or not ?? PS: I do not use yet the 'Nested iteration'algorithm 

December 1, 2005, 20:28 
Re: MULTIGRID 3D

#2 
Guest
Posts: n/a

That could depend on what is your tolerance of the convergence.


December 2, 2005, 00:39 
Re: MULTIGRID 3D *NM*

#3 
Guest
Posts: n/a



December 2, 2005, 04:45 
Re: MULTIGRID 3D

#4 
Guest
Posts: n/a

This is obviously true. I have impose eps=1.E7. with a Vcycle.


December 3, 2005, 03:03 
Re: MULTIGRID 3D

#5 
Guest
Posts: n/a

it sounds good.


December 3, 2005, 07:51 
Re: MULTIGRID 3D

#6 
Guest
Posts: n/a

Hi Davoche! I am also solving poisson equation with almost similar number of nodes. But I am solving with SOR method and nearly 500 iterations are required for convergence.(1.0E5) Since your solution is getting converged in 12 cycles only, there may be something wrong. Please check the velocity componants used as input in the poisson equation. What about the results?? regards, Jas


December 4, 2005, 03:12 
Re: MULTIGRID 3D

#7 
Guest
Posts: n/a

Jas,
That's the main reason to use multigrid  it accelerates convergence. The improvement over SOR is usually dramatic, so maybe you should consider changing your solver. Rami 

December 4, 2005, 12:30 
Re: MULTIGRID 3D

#8 
Guest
Posts: n/a

For a uniform grid and a smooth solution the answer would be no. A typical number here would be about half that for a properly sorted scheme. Even for a FAS scheme (solving u, v, w & p) for a low Reynolds number driven cavity 6 or 7 iterations would be typical for a sorted scheme.
The important property is less the number of iterations but that it requires the same number on 16^3, 32^3, 64^3, 128^3,... However, when the solution is not smooth and/or your coefficients vary strongly in size (i.e. nonuniform grid) then the performance will drop. So your 12 iterations may be fine but you need to check by running on different grid sizes. Often you can bring back lower numbers of cycles by adopting a better smoother and/or smoothing strategy. Unfortunately, the cost in CPU of such changes is not always worthwhile. 

December 4, 2005, 20:40 
Re: MULTIGRID 3D

#9 
Guest
Posts: n/a

Just a general question to those who are familiar with multigrid. Can you apply a multigrid or multilevel solver to a general linear system of equation, those otherwise could be solved with GaussSeidel etc. I mean lets say its not related to cfd, but just a set of linear equation of the form Ax=Q, can we apply this kind of strategy to accelerate this convergence.


December 5, 2005, 04:58 
Re: MULTIGRID 3D

#10 
Guest
Posts: n/a

It depends on the nature of the equations being solved. Multigrid requires a smoother (a fast simple approximate solver) that strongly reduces local errors so that when you run through the heirarchy of grids all error scales get strongly reduced. The solution also needs to be smooth for a general scheme although there are things one can do to adapt for specific cases.
Wave equations for example do not work out of the box with multigrid. However, transforming the equations from the time domain can produce something that multigrid handles well. Multigrid offers the ideal of scaling linearly with problem size and is often the only viable approach for solving implicit equations with many millions or billions of unknowns. Because of this, most sets of equations have people somewhere working to get multigrid going if it does not work when applied directly. 

December 5, 2005, 05:01 
Re: MULTIGRID 3D

#11 
Guest
Posts: n/a

Hi everybody, > Parth: Sorry but I don't undersatnd your question, can you reformulate it, please.
> Andy: I've tried with different grid size and I have always 12 cycles. But something is very strange. I have even 12 cycles with a UNIFORM grid. Maybe this is due to my Restrication and Prolongation operators. In fact, the transfer from one grid to another is made regularly ( no anisotropic coef ). I have try to implement it but the CPU cost is increased and I have always 12 cycle. No idea about that so I'm trying to use Nested iteration now (Full Multigrid ). 

December 5, 2005, 20:06 
Re: MULTIGRID 3D

#12 
Guest
Posts: n/a

My question was that for a given set of linear equations say AX=q ; and it should not matter how i got them or not neccessarily from CFD formulations, they could just be a set of linear equation. To this set of linear equations can we apply the acceleration due to multigrid kind of scheme.
the reason for asking so was this, i have seen mostly that we apply multigrid for the equations resulting from decretisation of partial differential equations. My question was can we apply it to system of such equations NOT resulting from partial differential equations. Does any such scheme exists. (i am no expert with multigrids so i think people who know better would be able to shed some light on this). 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
LES by Multigrid  mnabi  Main CFD Forum  8  July 11, 2009 20:21 
multigrid  sureshkumar  Main CFD Forum  0  June 9, 2006 01:20 
Multigrid method  Suresh kumar  Main CFD Forum  1  April 24, 2006 20:07 
MULTIGRID at the wall  Jim  Main CFD Forum  0  February 19, 2003 13:22 
Multigrid applied to ke models  Paulo Zandonade  Main CFD Forum  9  May 24, 1999 08:10 