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>