|
| |
|
Главная Журналы wbar - средняя скорость в сечении w; wp - смоченный периметр; wsum - {w dA . 10.2.4. Листинг подпрограммы ADAPT с сссссссссссссссссссссссссссссссссссссссссссссссссссссссс SUBROUTINE ADAPT с---- с--EXAMPLE 8 - CIRCULAR TUBE WITH RADIAL INTERNAL FINS С--------- $ 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=CIRCULAR TUBE WITH RADIAL INTERNAL FINS PRINTF=PRINT8 PLOTF=PLOTS NFINS=6 CALL DATA2(HRATIO,0.3,TRATIO,0.15) RAD=1. DIA=2.*RAD H=RAD*HRATIO C=RAD-H PI=3.14159 ALPHA=PI/FLOAT(NFINS) ALPHAT=ALPHA*TRATIO ALPHAC=ALPHA-ALPHAT CALL INTA3(NZX,2,NCVX(1),3,NCVX(2),7) CALL INTA3(NZY,2,NCVY(1),8,NCVY(2),4) CALL DATA4(XZONE(1),ALPHAT,XZONE(2),ALPHAC,YZONE(1),C, 1 YZ0NE(2),H) CALL ZGRID R(1)=0. M0DE=3 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,6) CALL DATA?(AMU,1.,DEN,1.,CP,1.,COND,1.,DPDZ,-1., 1 DTDZ,1.,TW,0.) RHOCP=DEN*CP CNDFIN=C0ND*3.5 DO 100 J=1,M1 DO 100 1=1,Ll W(I,J)=0. Т(I,J)=TW 100 CONTINUE RETURN 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) IF(X(I).LT.ALPHAT.AND.Y(J).GT.C) AR=0. 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=ALPHAC* RAD+H+ALPHAT * С DH=4.*ASUM/WP RE=DH*WBAR*DEN/AMU FRE=-2.*DPDZ*DH/(DEN*WBAR**2+SMALL)*RE HP=ALPHA*RAD QW=DTDZ*RHOCP*WSUM/HP ANU=QW*DIA/(COND*(TW-TB)+SMALL) DO 210 IUNIT=IU1,IU2 IF(ITER.EQ.O) WRITE(IUNIT,220) 220 FORMAT(IX, ITER",2X, "W ( 8,8) ,4X, W (11,4) , 1 3X, T (8,8) ,4X, •T(11,4) ,4X, FRE,8X, NU) WRITE(lUNIT,2 30) ITER,W(8,8),W(11,4),T(8,8), 1 T(11, 4) , FRE,ANU 230 FORMAT(2X,12,1P6E10.2) 210 CONTINUE IF(ITER.EQ.LAST) THEN DO 240 IUNIT=IU1,IU2 WRITE(lUNIT,250) 250 FORMAT( , I,8X,TH(I),6X, 1 -LOCAL NU (TUBE WALL)) 240 CONTINUE DO 2 60 1=2,L2 ANUL0C=FLUXM1(1,2)*DIA/(COND*(TW-TB) ) DO 270 IUNIT=IU1,IU2 WRITE(IUNIT,280) I,X(I),ANULOC 280 FORMAT(IX,12,5X,1PE9.2, lOX,1PE9.2) 27 0 CONTINUE 260 CONTINUE DO 290 J=1,M1 DO 290 1=1,LI W(I,J)=W(I,J)/WBAR IF(W(I,J).LT.l.E-10) W(I,J)=0.0 T(I, J) = (T(I,J)-TW)/(TB-TW) 290 CONTINUE CALL PRINT COME HERE TO FILL IBLOCK(I,J) BEFORE CALLING PLOT DO 2 95 J=2,M2 DO 295 1=2,L2 IF(X(I).LT.ALPHAT.AND.Y(J).GT.C) IBLOCK(I,J)=1 295 CONTINUE 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 IF(X(I).LT.ALPHAT.AND.Y(J).GT.C) GAM(I,J)=BIG 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 IF(X(I).LT.ALPHAT.AND.Y(J).GT.C) GAM(I,J)=CNDFIN SC(I,J)=-RHOCP*DTDZ*W(I, J) CONTINUE ENDIF COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 320 J=2,M2 KBCIl(J)=2 KBCLl (J)=2 320 CONTINUE DO 330 1=2,L2 KBCJl(I)=2 330 CONTINUE RETURN END С CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 10.2.5. Результаты расчетов RESULTS OF CONDUCT FOR POLAR COORDINATE SYSTEM CIRCULAR TUBE WITH RADIAL INTERNAL FINS
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 |