CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum

recover solution from dual grid

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 30, 2010, 15:31
Default recover solution from dual grid
New Member
Join Date: Jul 2010
Posts: 1
Rep Power: 0
jamesproctor is on a distinguished road
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
w_i^{0} = \frac{1}{\left|V_i\right|} \int_{V_i} u(x,y,0) dxdy

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)
w_i^{\prime} = -\frac{1}{\left|V_i\right|} \sum\limits_{j\in \mbox{neighbors}_i} \left|E_{ij}\right| \left(n_j \cdot f\right) \left(\alpha_j w_j +(1-\alpha_j)w_i\right)

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
u(\cdot) = \left|V\right| \nabla w_i
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.
jamesproctor is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
grid dependancy gueynard a. Main CFD Forum 19 June 27, 2014 21:22
convergence for higher grid points sandip Phoenics 3 April 15, 2009 13:25
stagged grid solution is high of numer viscosity? HT Main CFD Forum 13 March 4, 2004 12: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 18:01

All times are GMT -4. The time now is 05:37.