Dear all<br>
<br>
we are trying to make calculations of turbulent heat transfer in a<br>
pipe. Everything seems to work fine<br>
with the velocity field, however this is not the case for the<br>
temperature. It seems as if the convection<br>
for the temperature equation is switched off. So, we obtain a<br>
turbulent velocity distribution but a laminar temperature<br>
distribution.<br>
<br>
We have tried several thermal boundary condtions and at the moment we<br>
are just fixing the temperature at the wall<br>
T=0 and adding a heat source qvol=1. In the following we add details<br>
of the case.<br>
<br>
Thanks for your help<br>
<br>
Regards,<br>
Alex<br>
<br>
rea file<br>
 ****** PARAMETERS *****<br>
   2.6099999      NEKTON VERSION<br>
   3 DIMENSIONAL RUN<br>
         103  PARAMETERS FOLLOW<br>
   1.00000     p001 DENSITY<br>
  -750.       p002 VISCOS<br>
   0.00000     p003<br>
   0.00000     p004<br>
   0.00000     p005<br>
   0.00000     p006<br>
   1.00000     p007 RHOCP<br>
  -0.7         p008 CONDUCT<br>
[...]<br>
          13   LOGICAL SWITCHES FOLLOW<br>
 T      IFFLOW<br>
 T      IFHEAT<br>
 T      IFTRAN<br>
 T T F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)<br>
 F F T T T T T T T T T T  IFTMSH (IF mesh for this field is T mesh)<br>
 F      IFAXIS<br>
 F      IFSTRS<br>
 F      IFSPLIT<br>
 F      IFMGRID<br>
 F      IFMODEL<br>
 F      IFKEPS<br>
 F      IFMVBD<br>
 F      IFCHAR<br>
[...]<br>
<br>
usr file<br>
[...]<br>
c-----------------------------<div id=":1il">------------------------------------------<br>
      subroutine userf  (ix,iy,iz,ieg)<br>
      include 'SIZE'<br>
      include 'TOTAL'<br>
      include 'NEKUSE'<br>
<br>
      ffx = 0.0<br>
      ffy = 0.0<br>
      ffz = 2.0<br>
<br>
      return<br>
      end<br>
c-----------------------------------------------------------------------<br>
      subroutine userq  (ix,iy,iz,ieg)<br>
      include 'SIZE'<br>
      include 'TOTAL'<br>
      include 'NEKUSE'<br>
<br>
      qvol   = 1.0<br>
      source = 0.0<br>
<br>
      return<br>
      end<br>
c-----------------------------------------------------------------------<br>
[...]</div>