[Philipp Dietrich]

Hi,

I am currently trying to add a heat source to the icoFoam solver.

The heat source will be dependent on the Temperature gradient in a specific direction. Therefore I have to access the temperature gradient of each cell.
I tried to use the "fvc::grad(T)" to define a new volVectorField for the gradient.
However, it seems, that especially in the cells at a boundary this gradient is not calculated right. (nonzero temperature gradient at "empty" boundary; unphysical high/low temperature gradient at a wall with fixed temperature).

Does anyone know how to correctly calculate the temperature gradient for the whole internal Field?
Regards
Philipp