|July 30, 2010, 16:31||
recover solution from dual grid
Join Date: Jul 2010
Posts: 1Rep Power: 0
I'm new to solving pdes. I coded up a finite volume scheme for advection on a dual grid, and I'm seeing some weird behavior. I'll recap my scheme (seems pretty standard), but my real question is on recovering the solution.
So I start out with a triangular mesh, which I then break up into a barycentric dual grid. I'll call the function I wish to represent u(x,y,t), and so initially I set the cell value w for the volume V of the dual equal to the average value of the initial condition
Then my simplified update is (f is velocity field, n_j is normal from edge E_ij between cell i and cell j, alpha is upwind parameter)
Anyways, when I plot the result, I see advection, and my little 2d gaussian moves, but it diffuses as well. But I think that's because I should be looking at the function value u(x,y,t), not w(i,t). Perhaps u(x,y,t) isn't diffusing?
So my question is how to recover u(x,y,t) from w(i,t). Using the barycentric dual scheme, the value of each dual cell is the average value over that cell, so maybe if I do something like
but that seems iffy.
tldr; Is there a standard way of recovering the function value from a dual grid?
edit: in all the scenarios I've tried, the max(w) is decaying exponentially. This makes me think the scheme isn't working as it should. Do you see any problems there?
Last edited by jamesproctor; July 30, 2010 at 22:24.
|Thread||Thread Starter||Forum||Replies||Last Post|
|convergence for higher grid points||sandip||Phoenics||3||April 15, 2009 14:25|
|stagged grid solution is high of numer viscosity?||HT||Main CFD Forum||13||March 4, 2004 13:26|
|Combustion Convergence problems||Art Stretton||Phoenics||5||April 2, 2002 05:59|
|Dual grid structure - enlightment ??||Jens||Main CFD Forum||2||August 1, 2000 19:01|
|grid dependancy||gueynard a.||Main CFD Forum||13||November 20, 1998 13:54|