Temperature gradient using wall-functions
 

Hi,
I am calculating a channel flow including heat transfer. I'm using the k-epsilon model and pisoFoam solver. Since the flow is incompressible, I added the temperature energy equation to the solver, similar to the following link:
http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam

The mesh I use is a high-Re mesh, and therefore I have to use wall-functions.
To calculate the heat flux at the wall (constant temperature) I need the temperature gradient at the wall.
Now the question: Is there a (relatively) simple approach in calculating the this temperature gradient using the wall function equations? (something similar to the tool wallShearStress)
I think a similar tool within Fluent would be the wall function heat transfer coefficient.

To compare: the wallShearStress tool gives good results concerning the shear stress, while calculating the shear stress using the wallGradU-tool gives far too small values

Thanks a lot for your help
/Andreas

Dear Andreas
Did you find solution for this problem? I also want to calculate heat flux of a wall in OF-ext4.0, I use wallGradT.C in https://github.com/positroncascade/wallGradT and compile successfully. But when I want to use that in controlDict with bellow code, that give error
wallGradT1
{
type wallGradT;
functionObjectLibs ("libutilityFunctionObjects.so");
enabled true;
outputControl outputTime;
patches (top bottom);
}
if you find any solution, please favor me.

just type wallGradT in terminal at the end of simulation.