Hi, I am currently evaluating
 

Hi,
I am currently evaluating OpenFOAM in order to develop high order DNS solver, using the interpolation technic of Carpentier (INRIA) (least square, ENO type) using explicit correction.
I am trying to understand the different classes in OpenFoam and I have some questions:
- do the interpolation scheme in OpenFoam take into account the boundary condition and how ?
- I understand how to get all the neighboring cells with mesh.cellCells(), but I don't know how to get all the faces of a given cell, and specially the boundary cell ?
- from a general point of view, what derived class has to be developed for changing the interpolation scheme in OpenFoam

Any help would be appreciated
Thanks,
Marc