[SaufiSangi]

Hello everybody,
I am working on an axisymmetric mesh, built from a 2D mesh and then made axisymmetric by "makeAxialMesh".
For my CFD code I need the calculation of the divergence of a volVectorField (in particular I need to calculate a curvature, therefore a divergence of normal vectors to an interface). This seems to be well calculated everywhere, but not for the cells on the symmetry axis ( the "empty" one). It just gives me a clearly wrong result for these cells.
The only thing that comes to my mind is the different shape of the cells: on the axis they are prisms, while all over the domain are hexahedrals. Can this be a problem?

Thanks guys,

Sangi