My question is about implementation alone.
- Consider a square domain with regular square, cell centred finite volumes. This is for the multiscale finite volume method (Jenny and Lunati)
I need to solve the Poisson equation on each cell of the "dual mesh", i.e., the mesh constructed by joining centroids of each cell. (the dual mesh will look like a translation of the original mesh; 4 cells of the dual mesh make up one cell of primary mesh)
I will need to use this to construct the "transmissibilities" for each cell of the primary mesh.
My question is:
- What kind of date structure should I use that will help me go from primary to dual mesh easily
- I am confused what I shld store, the coordinates of the cell centres, or the locations of the faces ?
- How do I store the data that will tell me which 4 dual cells make up a given primary cell ?
It seems to me I will need to do a kind of "assembly" for my transmissibility matrix as I visit each dual cell and solve a local Poisson problem.