Dear Pau.
I am a new comer in Nek5000. However, I would like to discuss this post.
If you need to obtain the turbulent scalar (or heat) fluxes, i would like to recommed the subroutine userchk in tutorial case of 'turbChannel'.
The subroutine userchk in 'turbChannel' provides us with an excellent example code for Nek5000 to obtain standard flow and heat (or scalar) transfer statistics, such as Reynolds stress, mean velocity and temperature, turbulent heat (or scalar) fluxes in a channel.
That is to say, we can try to write some code sections on the existing basis of subroutine userchk in 'turbChannel' to meet our needs.
It is worthynoting that the original code of subroutine userchk in 'turbChannel' might have a little bug when calculating the streamwise turbulent heat (or scalar) flux for passive scalar transport in a channel . An attempt to get the correct results can see in this post:
A little bug in the tutorial case "turbChannel"
In summary, i would recommend to make the following operations
(2) Run the tutorial case 'turbChannel'
(3) Read the subroutine 'userchk' in 'turbChannel.usr' file
(4) Try to modify the subroutine 'userchk' to get the statistics which we want.
(5) Use a coarse mesh to validate a benchmark case to check the modfication by ourselves is correct.
Besides, maybe the KTH toolbox is also a good choice to get more statistics, such as budgets of turbulent kinetic energy equation.
I hope this post would help you to some extent.
Best,
Student, Xingguang Zhou
Xi'an Jiaotong Univ.