CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

local timestep for agglomeration multigrid method

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 17, 2000, 11:04
Default local timestep for agglomeration multigrid method
  #1
Jian Xia
Guest
 
Posts: n/a
Dear friends,

When I try to use agglomeration multigrid method to solve Euler equations on unstructured grid, I find it difficult for me to determinate local timestep on coarse grid. Could you give me some advice? Thanks.
  Reply With Quote

Old   March 20, 2000, 09:52
Default Re: local timestep for agglomeration multigrid method
  #2
A. Haselbacher
Guest
 
Posts: n/a
Hello,

you can do it in exactly the same way as on the finest grid, i.e., for a given control volume

dt = vol/((u.n+a)ds)

where the symbols have their usual meanings. All the quantities are available on the coarse grids also, the volume is simply the agglomerated volume, and the face vector area ds is also available. Maybe you need to be more specific in your question as to where you encounter problems.

I programmed the agglomeration multigrid method and I used the above definition - there were absolutely no problems with instability in time. For Euler flows, agglomeration multigrid is quite straightforward and works extremely well. I got convergence rates of down to 0.5 for first order schemes and 0.75 for second order schemes (first order on coarse grid levels) with an optimised 4-stage Runge-Kutta scheme. For viscous flows, however, you have to be more clever...

Greetings & good luck

Andreas
  Reply With Quote

Old   March 20, 2000, 11:22
Default Re: local timestep for agglomeration multigrid method
  #3
Jian Xia
Guest
 
Posts: n/a
Dear Haselbacher, Thank you for your help! In fact, the local timestep I am using is just as you say. The residu keep unchanged(about 1.0e-1) after several iterations. Of course, the results are not right. The method I am using is just like Mavriplis(cell-centred symmetric finite-volume spatial discretisation, explicit,multi-stage procedure). In sigle fine grid, the residu 1.0e-6 is achieved after 1000 steps. The CFL can reach 10 in the finest grid, 2 in the second and third grid and 1 in the coarsest grid(for single grid computation). I think maybe you use the same method as I do. Could you give me more information about it? Thanks in advance.
  Reply With Quote

Old   March 20, 2000, 12:20
Default Re: local timestep for agglomeration multigrid method
  #4
A. Haselbacher
Guest
 
Posts: n/a
Hello,

I am happy to give you more information - but I don't really know what else to tell you. The routines which computed the time step, fluxes, residual and enforced boundary conditions and so on were the same on _all_ grid levels. That is, I simply had a wrapper around these routines which drove the multigrid process. The only thing you also need to add is the computation of the forcing function and the restriction and prolongation operators. For inviscid flows, there really is nothing more to it.

The numerical method I used was a finite-volume method based on dual control volumes, with Green-Gauss or least-squares reconstruction and Roe's flux-difference splitting. I don't think that the discretisation should make that much of a difference.

There's a few things you can check:

1. Sum the control-volume areas on all grid levels. Are the sums constant? 2. Accumulate the control-volume face areas? Do they add to zero for all volumes on all grid levels? 3. Try running the solver in single-grid mode on the various coarse levels. This is easily done by restricting the solution down to a given coarse level, ignoring the forcing function, and just running the code to convergence.

I also found it helpful to monitor the residual, correction, and forcing function at a given vertex on the fine grid and the coarse grid. You can see quite easily when something's wrong with your multigrid method if you look at these quantities.

I realise that this is maybe not as specific as you wanted, but I hope it helps all the same.

Greetings & good luck

Andreas
  Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Multigrid method for Lax-Wendroff scheme samira Main CFD Forum 1 July 2, 2010 22:17
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
multigrid method hawk Main CFD Forum 2 November 17, 2005 17:19
Local stepping for steady state Mehdi Main CFD Forum 2 June 17, 2003 09:43
Multigrid - Local Mode Analysis andreas Main CFD Forum 0 March 31, 2003 12:41


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