<br clear="all">Hello, Paul,<br><br>Thank you for your answer for the same problem last month.  <br><br>But I still could not figure out this issue. (Domain:6.4/2.0/3.2, Resolution: 12*8*12 elements with 8th order, turbulent Re_tau=178.12, Pr=0.71) ,<br>
piece of my *.rea file <br>*****************************************************************************<br>   1.00000              DENSITY<br>   -178.12               Re_b<br>   0.000000E+00    BETAG<br>   0.00000              GTHETA<br>
   0.000000E+00     PGRADX<br>   0.000000E+00     FLOWRATE<br>   1.0                       RHOCP<br>  -126.47                 CONDUCT<br>*****************************************************************************<br><br>
my *usr file<br><br>*****************************************************************************<br> subroutine userq  (ix,iy,iz,ieg)<br>     include 'SIZE'<br>     include 'TOTAL'<br>     include 'NEKUSE'<br>
     ie   = gllel(ieg)<br><br>     qvol = vx(ix,iy,iz,ie)/15.771<br>     source = 0.0<br><br>     return<br>     end<br>*****************************************************************************<br>Here 15.71 is averaged mean velocity over channel section. <br>
<br><br>After a long time running (dimensionless time > 350 )<br><br> I compared the result ,  the velocity profile looks good.<br> <br>But the temperature is much higher than Kasagi's result. <br><br><br><br>I'm wondering whether the parameters above are correct ?<br>
<br><br> I attached the plot of mean velocity and mean temperature.  <br><br>Thank you again<br>   <br><br><br>-- <br>Zhicheng Wang<br><br>