compute heatflux (q) in rhoCentralFoam
 

Hi,

I would like to calculate the heat flux in rhoCentralFoam.If i use the below equation will be correct?

volVectorField q = -k * fvc::grad(T);

I tried it but i get some strange results.
Has anyone tried that before?

Thanks