西安交通大学：《数值传热学》课程教学资源（课件讲稿）第十一章 通用程序在二维椭圆型流动与换热问题中的应用示例 11.6 Complicated fully developed fluid flow and heat transfer in square duct 11.7 Impinging flow on a rotating disc 11.8 Turbulent flow and heat transfer in duct with a central jet

ENTRYOUTPUTIF(ITER==O) THENPRINT401WRITE(8.401)401FORMAT(1X'ITER,6X'SMAX,8X,SSUM,7X'V(6,4),& 6X'T(2.6).6X,FRE)ELSEPRINT 403,ITER, SMAX, SSUM, V(6,4), T(2,6),FREWRITE(8.403)ITER,SMAX,SSUM,V(6,4),T(2,6),FRE403FORMAT(1X,16,1P5E12.3)ENDIFIF(ITER/=LAST)RETURNDO410J=1.M1DO 411 I=1.L1IDimensionlessoutputforwF(I,J,5)-F(I,J,5)/WBAR411ENDDO410ENDDOCALLPRINTRETURN16/98

16/98 ENTRY OUTPUT IF(ITER==0) THEN PRINT401 WRITE(8,401) 401 FORMAT(1X,' ITER',6X,'SMAX',8X,'SSUM',7X,'V(6,4)', & 6X,'T(2,6)',6X,'F.RE') ELSE PRINT 403, ITER, SMAX, SSUM, V(6,4), T(2,6), FRE WRITE(8,403) ITER,SMAX,SSUM,V(6,4),T(2,6),FRE 403 FORMAT(1X,I6,1P5E12.3) ENDIF IF(ITER/=LAST) RETURN DO 410 J=1,M1 DO 411 I=1,L1 F(I,J,5)=F(I,J,5)/WBAR 411 ENDDO 410 ENDDO CALL PRINT RETURN !Dimensionless output for w

ENTRYGAMSORDO500J=1.M1DO501I=1L1!Ffor velocityGAM(IJ)-AMUIF(NF=4)THEN！TfortempGAM(IJ)=CONDGAM(I,1)=0.!Adiabaticfor southand northboundariesGAM(I,M1)=0.ENDIF501ENDDO500ENDDODO510J=2.M2DO511I=2.L2(i, j+1)IF(NF=-2)THENIF(J/=2)THEN！JST=3forv-eqTM=(T(I.D+T(IJ-1))*O.5Source term of y-eqCON(LJ)=TM*GBR(i,j)(i-1,J)(i+1,j)ENDIFGBR = gPref βTENDIFTMIF(NF=-5)CON(LJ)=-DPDZ511ENDDOV(i,j)! Source term of w-eq510ENDDORETURNEND(i,j-1)17/98

17/98 ENTRY GAMSOR DO 500 J= 1,M 1 DO 501 I= 1,L 1 GAM(I,J)=AMU IF(NF== 4 ) THEN GAM(I,J)=COND GAM(I, 1)= 0 . GAM(I,M 1)= 0 . ENDIF 501ENDDO 500ENDDO DO 510 J= 2,M 2 DO 511 I= 2,L 2 IF (NF== 2 ) THEN IF(J/= 2 ) THEN TM=(T(I,J)+T(I,J - 1))* 0 . 5 CON(I,J)=TM*GBR ENDIF ENDIF IF (NF== 5 ) CON(I,J)= -DPDZ 511 ENDDO 510 ENDDO RETURN END GBR ref = g T   ！Source term of v -eq. ！Source term of w -eq. ！Adiabatic for south and north boundaries ！  for temp. TM ！JST = 3 for v -eq ！  for velocity

11-6-4ResultsanalysisCOMPUTATIONINCARTESIANCOORDINATES*水**水***水水水水*水水水**水水***水****水***水水***水水水*水水水SSUMV(6,4)T(2,6)F.REITER SMAX00.000E+000.000E+000.000E+000.000E+000.000E+000.000E+000.000E+000.000E+001.000E-010.000E+00121.273E+01 -1.907E-061.016E+012.848E-010.000E+0036.308E+001.073E-061.926E+013.445E-010.000E+0042.978E+007.153E-072.076E+013.826E-010.000E+0051.237E+00-5.960E-072.284E+013.854E-010.000E+0063.889E-016.454E-01-4.768E-072.304E+010.000E+0072.911E-017.153E-072.342E+013.894E-010.000E+0081.338E-01-3.278E-072.346E+013.900E-010.000E+0096.046E-02-5.364E-072.352E+013.900E-010.000E+00102.868E-02-5.364E-072.352E+013.900E-010.000E+00111.286E-02-4.321E-072.353E+013.900E-010.000E+0018/98

18/98 COMPUTATION IN CARTESIAN COORDINATES ********************************************** ITER SMAX SSUM V(6,4) T(2,6) F.RE 0 0.000E+00 0.000E+00 0.000E+00 0.000E+00 0.000E+00 1 0.000E+00 0.000E+00 0.000E+00 1.000E-01 0.000E+00 2 1.273E+01 -1.907E-06 1.016E+01 2.848E-01 0.000E+00 3 6.308E+00 1.073E-06 1.926E+01 3.445E-01 0.000E+00 4 2.978E+00 7.153E-07 2.076E+01 3.826E-01 0.000E+00 5 1.237E+00 -5.960E-07 2.284E+01 3.854E-01 0.000E+00 6 6.454E-01 -4.768E-07 2.304E+01 3.889E-01 0.000E+00 7 2.911E-01 7.153E-07 2.342E+01 3.894E-01 0.000E+00 8 1.338E-01 -3.278E-07 2.346E+01 3.900E-01 0.000E+00 9 6.046E-02 -5.364E-07 2.352E+01 3.900E-01 0.000E+00 10 2.868E-02 -5.364E-07 2.352E+01 3.900E-01 0.000E+00 11 1.286E-02 -4.321E-07 2.353E+01 3.900E-01 0.000E+00 11-6-4 Results analysis

SSUMV(6,4)T(2,6)F.RE!ITERSMAX123.901E-016.224E-032.850E-072.353E+010.000E+00133.349E-03-3.660E-072.353E+013.901E-010.000E+003.901E-01141.544E-031.974E-072.353E+010.000E+00：158.407E-04-2.626E-072.353E+013.901E-010.000E+00163.686E-04-1.118E-082.353E+013.901E-010.000E+00171.961E-041.043E-072.353E+013.901E-010.000E+00187.963E-052.775E-073.901E-010.000E+002.353E+01194.327E-053.166E-082.353E+013.901E-010.000E+00202.098E-05-1.825E-072.353E+013.901E-016.000E+01!Solvingw-eq.212.098E-05-1.825E-072.353E+013.901E-015.323E+015.238E+01222.098E-05-1.825E-072.353E+013.901E-01232.098E-05-1.825E-07 2.353E+01 3.901E-015.236E+0124 2.098E-05-1.825E-072.353E+01 3.901E-015.236E+01252.098E-05-1.825E-072.353E+013.901E-015.236E+0119/98

19/98 12 6.224E-03 2.850E-07 2.353E+01 3.901E-01 0.000E+00 13 3.349E-03 -3.660E-07 2.353E+01 3.901E-01 0.000E+00 14 1.544E-03 1.974E-07 2.353E+01 3.901E-01 0.000E+00 15 8.407E-04 -2.626E-07 2.353E+01 3.901E-01 0.000E+00 16 3.686E-04 -1.118E-08 2.353E+01 3.901E-01 0.000E+00 17 1.961E-04 1.043E-07 2.353E+01 3.901E-01 0.000E+00 18 7.963E-05 2.775E-07 2.353E+01 3.901E-01 0.000E+00 19 4.327E-05 3.166E-08 2.353E+01 3.901E-01 0.000E+00 20 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 6.000E+01 21 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 5.323E+01 22 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 5.238E+01 23 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 5.236E+01 24 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 5.236E+01 25 2.098E-05 -1.825E-07 2.353E+01 3.901E-01 5.236E+01 ! Solving w-eq. ! ITER SMAX SSUM V(6,4) T(2,6) F.RE

****************************************************.VEL U.523467I=J70.00E+000.00E+000.00E+000.00E+000.00E+000.00E+0060.00E+00-1.52E+01-1.78E+01-1.77E+01-131E+010.00E+0050.00E+00-8.36E+00-1.40E+01-1.40E+01-9.70E+00.0.00E+0040.00E+007.76E-01.8.31E-02-8.31E-02-7.76E-010.00E+0030.00E+009.70E+00140E+01140E+018.36E+000.00E+0020.00E+001.31E+011.77E+011.78E+01.1.52E+010.00E+0010.00E+000.00E+000.00E+000.00E+000.00E+000.00E+00Naturalconvectionincross section20/98

20/98 ************************** .VEL U. ************************** I = 2 3 4 5 6 7 J 7 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 6 0.00E+00-1.52E+01-1.78E+01-1.77E+01 -1.31E+01 0.00E+00 5 0.00E+00-8.36E+00-1.40E+01-1.40E+01 -9.70E+00 0.00E+00 4 0.00E+00 7.76E-01 8.31E-02 -8.31E-0 2 -7.76E-01 0.00E+00 3 0.00E+00 9.70E+00 1.40E+01 1.40E+01 8.36E+00 0.00E+00 2 0.00E+00 1.31E+01 1.77E+01 1.78E+01 1.52E+01 0.00E+00 1 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 0.00E+00 Natural convection in cross section