Главная  Журналы 

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

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