Hi,<br><br>I'm running a DNS of fully developed heat transfer in a 
channel, the same as H. Kawamura et. al, International Journal of Heat 
and Fluid Flow, 1998.<br><br>I use following subroutine to get rid of additional heat<br><br>c-----------------------------------------------------------------------<br>      subroutine userq  (ix,iy,iz,ieg)<br>      include 'SIZE'<br>

      include 'TOTAL'<br>      include 'NEKUSE'<br>      ie   = gllel(ieg)<br>      qvol = -vx(ix,iy,iz,ie)<br>      source = 0.0<br><br>      return<br>      end<br>c-----------------------------------------------------------------------<br>

<br>Other procedures are similar to turbChannel example.<br><br>But I find the mean temperature is not correct,  there should not be negative  value in my case.<br><br>I attached my *.usr and *.rea and means.dat to this email.<br>

<br>is there any one  else did DNS of heat transfer in a channel ?<br><br><br>Thank you very much<br><br><br>zhicheng<br clear="all"><br><br>