CFD Online Discussion Forums

CFD Online Discussion Forums (
-   Main CFD Forum (
-   -   Nonlinear Multigrid: Coarse Grid Operators (

Greg Lewin December 1, 2006 10:39

Nonlinear Multigrid: Coarse Grid Operators
I am implementing a 3D CFD solver, and I have a question about nonlinear multigrid.

I have implemented a 3rd order advection / 2nd order viscous routine on a single grid. The system is solved using a second order implicit scheme, with Roe splitting, similar to Rogers and Kwak (and others):

R(Q)=0 => L(dQ) = R(Q) + some iteration residual

where Q are the flow variables, R is the discretization of the NS equations (large-stencil) and L is a simplified operator of the non-linear operator R. L is constructed so that the system is block-tridiagonal. Fairly standard fare.

Cut to the multi-grid. Following Press et al., on coarser grids I will end up with equations of the form:


where the underscores represent a coarse grid operator or variable and f is a bunch of stuff left over from restriction/prolongation. I*Q is the restricted solution.

The question is thus: Do I have to calculate the coarse grid operator R_ using the same method that I used for the fine grid operator R? This would require evaluating all of the advection terms _at each grid level_, which seems like it would be costly. I would much rather calculate R_ as an agglomeration of the terms in R, or even L:


This way, L, and all the coarse grid L_'s are calculated once per pseudo-time step. But this seems like I might be breaking the spirit of non-linear multigrid, or drastically hurting convergence.


All times are GMT -4. The time now is 10:45.