CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

full cell stencil with parallelisation

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

LinkBack Thread Tools Display Modes
Old   April 11, 2012, 14:51
Default full cell stencil with parallelisation
Matthias Walter
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 63
Rep Power: 10
matthias is on a distinguished road

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 is offline   Reply With Quote

Old   April 13, 2012, 11:48
Senior Member
kmooney's Avatar
Kyle Mooney
Join Date: Jul 2009
Location: Amherst, MA USA - San Diego, CA USA
Posts: 321
Rep Power: 11
kmooney is on a distinguished road
Hi Matthias,

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


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

Good luck,

kmooney is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
FvMatrix coefficients shrina OpenFOAM Running, Solving & CFD 10 October 3, 2013 14:38
extended cell stencil tomislav_maric OpenFOAM Programming & Development 4 February 1, 2012 10:38
Cells with t below lower limit Purushothama Siemens 2 May 31, 2010 21:58
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15
Warning 097- AB Siemens 6 November 15, 2004 05:41

All times are GMT -4. The time now is 01:41.