[Nek5000-users] heat transfer in a channel

[nek5000-users at lists.mcs.anl.gov]

  Hi,

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.

I use following subroutine to get rid of additional heat

c-----------------------------------------------------------------------
      subroutine userq  (ix,iy,iz,ieg)
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'
      ie   = gllel(ieg)
      qvol = -vx(ix,iy,iz,ie)
      source = 0.0

      return
      end
c-----------------------------------------------------------------------

Other procedures are similar to turbChannel example.

But I find the mean temperature is not correct,  there should not be
negative  value in my case.

I attached my *.usr and *.rea and means.dat to this email.

is there any one  else did DNS of heat transfer in a channel ?


Thank you very much


zhicheng
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100917/83091fb4/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: turbHeat.box
Type: application/octet-stream
Size: 510 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100917/83091fb4/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: turbHeat.rea
Type: application/octet-stream
Size: 3512 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100917/83091fb4/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: turbHeat.usr
Type: application/octet-stream
Size: 16611 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100917/83091fb4/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: means.dat
Type: application/x-ns-proxy-autoconfig
Size: 2260 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20100917/83091fb4/attachment.dat>