Главная Журналы FRE - произведение JRe; HP - обогреваемый участок периметра; PI - число п; QW - плотность теплового потока на стенке RE - число Рейнольдса; RHOCP - объемная теплоемкость рс,; RIN - радиус внутренней трубы Л,„; ROUT - радиус внешней трубы К,; Т (I, J) - температура Т; ТВ - среднемассовая температура Г,; TEMP - временная переменная для (Г-Г)/(J, - Г); TSUM - jwTdA; TW - температура стенки Т; W (I, J) - продольная скорость w; WBAR - средняя скорость в сечении w; WP - смоченный периметр; WSUM - jw dA . 10.3.4. Листинг подпрограммы ADAPT с сссссссссссссссссссссссссссссссссссссссссссссссссссссссс SUBROUTINE ADAPT с----- с--EXAMPLE 9 - ANNULAR SECTOR DUCT С----- $ INCLUDE: COMMON DIMENSION W(NI,NJ),T(NI,NJ) EQUIVALENCE (F(1,1,1),W(1,1)), (F (1,1,2),T(1,1)) ENTRY GRID HEADER=ANNULAR SECTOR DUCT PRINTF=PRINT9 PL0TF=PL0T9 CALL DATA4 (ROUT,1.,RIN,0.2 , ALPHA, 60. , PI, 3 . 14159) XL=0.5*ALPHA*PI/180. CALL DATA2(R(l),RIN, YL,ROUT-RIN) CALL INTA3(NCVLX,5,NCVLY,10,MODE,3) CALL EZGRID RETURN ENTRY BEGIN TITLE(1)= W/WBAR TITLE(2)= (T-TW)/(TB-TW) CALL INTA6(KSOLVE(1),1,KPRINT(1),1,KPRINT(2) ,1, 1 KPLOT(l),1,KPL0T(2),1,LAST,10) CALL DATA?(AMU,1.,COND,1.,CP,1.,DEN,1.,DPDZ, -1.,TW, 1 0.,DTBDZ,1.) RHOCP=DEN*CP DO 100 J=1,M1 DO 100 1=1,Ll W(l,J)=0. T(I,J)=TW 100 CONTINUE RETURN Q ir - -ic-ic-ic-k-k -k - -k-k-k-k-k-k - -k - -k - -k - -k-k-k - -k-k - -k - *-.-k - -k - -k--k - -k - ENTRY OUTPUT IF(ITER.EQ.3) THEN KSOLVE(1)=0 KSOLVE(2)=1 ENDIF ASUM=0. WSUM=0. TSUM=0. DO 200 J=2,M2 DO 200 1=2,L2 AR=XCV(I)*YCVR(J) ASUM=ASUM+AR WSUM=WSUM+W(I,J)*AR TSUM=TSUM+W(I, J)*T(I,J)*AR 200 CONTINUE WBAR=WSUM/ASUM TB=TSUM/(WSUM+SMALL) WP=XL*(ROUT+RIN)+YL HP=XL*RIN+YL DH=4.*ASUM/WP RE=DH*WBAR*DEN/AMU FRE=-2.*DPDZ*DH/(DEN*WBAR**2+SMALL)*RE QW=DTBDZ*RHOCP*WSUM/HP ANU=QW*YL/(COND*(TW-TB)+SMALL) DO 210 IUNIT=IU1,IU2 IF(ITER.EQ.0) WRITE(lUNIT,220) 220 FORMATdX, 4TER,2X, W(6,8) , 4X, W(4,ll) 1 3X, •T(6,8) ,4X, T (4,11) ,4X, FRE,8X, NU WRITE(lUNIT,230) ITER,W(6,8),W(4,11),T(6,8), 1 T(4, 11) , FRE,ANU 230 FORMAT(2X,12,1P7E10.2) 210 CONTINUE IF(ITER.EQ.LAST) THEN COME HERE TO CALCULATE LOCAL NU DO 240 1 = 2,L2 ANUI = FLUXJ1 (1,2) YL/ (COND* (TW-TB) +SMALL) DO 250 IUNIT=IU1,IU2 IF(I.EQ.2) WRITE(lUNIT,260) 260 FORMAT( , I,8X, TH (I) ,6X, 1 -LOCAL NU (INNER WALL)) WRITE(lUNIT,265) I,X(I),ANUI 2 65 FORMAT(IX,12,5X,1PE9.2,lOX,1PE9.2) 250 CONTINUE 240 CONTINUE DO 270 J=M2,2,-1 ANULFT=FLUXI1(J,2)*YL/(COND*(TW-TB)+SMALL) DO 275 IUNIT=IU1,IU2 IF(J.EQ.M2) WRITE(lUNIT,280) 280 FORMAT( , J,8X, Y(J) , 7X, 1 LOCAL NU (SIDE WALL)) WRITE(lUNIT,285) J,Y(J),ANULFT 2 85 FORMAT(IX,12,5X, 1PE9.2 , lOX, 1PE9 . 2 ) 27 5 CONTINUE 270 CONTINUE DO 290 J=1,M1 DO 290 1=1,LI W(I,J)=W(I,J)/WBAR T(I, J) = (T(I,J)-TW)/(TB-TW) 2 90 CONTINUE CALL PRINT CALL PLOT ENDIF RETURN ENTRY PHI IF(NF.EQ.l) THEN DO 300 J=2,M2 DO 300 1=2,L2 GAM(I,J)=AMU SC(I,J)=-DPDZ CONTINUE ENDIF IF(NF.EQ.2) THEN DO 310 J=2,M2 DO 310 1=2,L2 GAM(I,J)=COND TEMP=(T(I,J)-TW+SMALL)/(TB-TW+SMALL) DTDZ=TEMP*DTBDZ SC(I,J)=-RHOCP*DTDZ*W(I, J) CONTINUE ENDIF COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 320 J=2,M2 KBCLl(J)=2 CONTINUE IF(NF.EQ.2) THEN DO 330 1=2,L2 KBCMl(I)=2 CONTINUE ENDIF RETURN END С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 [ 67 ] 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 |