|
[Sponsors] |
local timestep for agglomeration multigrid method |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
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. |
|
![]() |
![]() |
![]() |
#2 |
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 |
|
![]() |
![]() |
![]() |
#3 |
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.
|
|
![]() |
![]() |
![]() |
#4 |
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 |
|
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Multigrid method for Lax-Wendroff scheme | samira | Main CFD Forum | 1 | July 2, 2010 22:17 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 07:36 |
multigrid method | hawk | Main CFD Forum | 2 | November 17, 2005 16: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 |