[matthias]

Hi,

I'm trying to get the cell centered values surrounding an arbitrary given cell in a certain volField, i.e. volVectorField U. First I have used a search algorithm that finds for a given cell all neighbour cells and returns the cell values but all in all, it's not working in parallel due to missing diagonal dependencies in the mesh.

Afterwards, I got a hint that the extendedStencil stuff could possibly be used for this. So I'm trying to make it work but ...hmm ... no success till now (especially not in parallel).

 CPCCellToCellStencil wideStencil(mesh());

 forAll(mesh().C(), lcellI)
 {
       const labelList& cell27Stencil=wideStencil[lcellI];
 
       forAll(cell27Stencil, cellI)
       {
          // now it becomes tricky ...

         // that's not working since the global addressing of the volField is wrong
         Pout << " center value: " << unFilteredField().[cell27Stencil[cellI]] << endl;
        
       }
  }

Furthermore I have no idea how to get correctly the values from the boundary fields (and processor patches) using the stencil stuff.

So I would be grateful for suggestions how to solve this messy problem! Maybe somebody has already done this or has some experience with the stencil stuff.


Best regards

Matthias

[kmooney]

Hi Matthias,

While I don't have much experience in this matter, there is an extended stencil implementation of leastSquares here:

/OpenFOAM/OpenFOAM-1.6-ext/src/finiteVolume/lnInclude/extendedLeastSquaresVectors.C


You might be able to gain some insight if you haven't already checked it out.

Good luck,

Kyle