# recover solution from dual grid

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

 July 30, 2010, 15:31 recover solution from dual grid #1 New Member   Join Date: Jul 2010 Posts: 1 Rep Power: 0 Sponsored Links Hi all, 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 21:24.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post gueynard a. Main CFD Forum 19 June 27, 2014 21:22 sandip Phoenics 3 April 15, 2009 13:25 HT Main CFD Forum 13 March 4, 2004 13:26 Art Stretton Phoenics 5 April 2, 2002 05:59 Jens Main CFD Forum 2 August 1, 2000 18:01