|
| |
|
Главная Журналы Н=1. SHALF=0.5*S TH=0.2S X1 = 0.5TH X2=SHALF-X1 HFIN=0.75*H C=H-HFIN CALL INTA4(NZX,3,NCVX{1),2,NCVX{2),8,NCVX{3),2) CALL DATA3(XZONE(1),XI,XZONE(2),X2-XI,XZONE(3) , XI) CALL INTA4(NZY,3,NCVY(1),3,NCVY(2),5,NCVY(3),3) CALL DATA3(YZONE(1),C,YZONE(2),HFIN-C, YZONE(3) ,C) CALL ZGRID RETURN ENTRY BEGIN TITLE(1)= W/WBAR TITLE(2)=(T-TINF)/(TB-TINF) CALL INTA6(KSOLVE(1),1,KPRINT(1),1,KPRINT(2),1, 1 KPLOT(l),1,KPL0T(2),1,LAST,10) CALL DATA9(AMU,1.,COND,1.,CP,1.,DEN,1.,DPDZ, -1., I TINF,0.,DTBDZ,1.,BITOP,5.,BIBOT,0.2) HTOP=BITOP*COND/H HBOT=BIBOT*COND/H RHOCP=DEN*CP CNDFIN=1.8*C0ND DO 100 J=1,M1 DO 100 1=1,Ll W(I,J)=0. T(I,J)=TINF+1. 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)*YCV(J) IF(X(I).LT.Xl.AND.Y(J).LT.HFIN) AR=0. IF(X(I).GT.X2.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=S+2.*HFIN DH=4.*ASUM/WP RE=DHWBAR* DEN/AMU FRE=-2.*DPDZ*DH/(DEN*WBAR**2+SMALL)*RE HP=S QW=DTBDZ*RHOCP*WSUM/HP ANU=QW*H/(COND*(TINF-TB)+SMALL) QTOP=0. QBOT=0. DO 210 1=2,L2 QT0P=QT0P+FLUXM1(1,2)*XCV{I) QB0T=QB0T+FLUXJ1(1,2)*XCV(I) 210 CONTINUE ANUT=QTOP*H/{SHALF*COND*(TINF-TB)+SMALL) ANUB=QBOT*H/(SHALF*COND*(TINF-TB)+SMALL) DO 220 IUNIT=IU1,IU2 IF(ITER.EQ.O) WRITE(lUNIT,230) 230 F0RMAT(1X, ЧТЕРЧгх, W(8,8) , 5X, T(8,8) , 6X, 1 FRE,8X,NU,7X,NU(TOP)V3X,•NU(BOTTOM)) WRITE (lUNIT,240) ITER,W(8,8) ,T(8,8), 1 FRE,ANU,ANUT,ANUB 240 FORMAT(IX,12,1P6E11.3) 220 CONTINUE IF(ITER.EQ.LAST) THEN DO 250 J=1,M1 DO 250 1=1,LI W(I,J)=W(I,J)/WBAR IF(W(I,J).LT.l.E-lO) W(I,J)=0. T(I,J)=(T(I,J)-TINF)/(TB-TINF) 250 CONTINUE CALL PRINT COME HERE TO FILL IBLOCK(I,J) BEFORE CALLING PLOT DO 2 60 J=2,M2 DO 2 60 1=2,L2 IF(X(I).LT.Xl.AND.Y(J).LT.HFIN) IBLOCK(I,J)=1 IF(X(I).GT.X2.AND.Y(J).GT.C) IBLOCK(I,J)=1 CONTINUE 260 С- 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.Xl.AND.Y(J) IF(X(I).GT.X2.AND.Y(J) SC(I,J)=-DPDZ CONTINUE ENDIF .LT.HFIN) GAM(I,J)=BIG .GT.C) GAM(I,J)=BIG IF(NF.EQ.2) THEN DO 310 J=2,M2 DO 310 1=2, L2 GAM{I,J)=COND IF{X{I).LT.X1.AND.Y{J).LT.HFIN) GAM{I,J)=CNDFIN IF{X(I).GT.X2.AND.Y{J).GT.C) GAM{I,J)=CNDFIN TEMP=(T(I, J)-TINF)/(TB-TINF+SMALL) DTDZ=TEMP*DTBDZ SC(I,J)=-RHOCP*DTDZ*W(I,J) 310 CONTINUE ENDIF COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 320 J=2,M2 KBCIl{J)=2 KBCLl{J)=2 320 CONTINUE IF{NF.EQ.2) THEN DO 330 1=2,L2 KBCMl(I)=2 KBCJl(I)=2 FLXCMl{I)=HTOP*TINF FLXPMl(I)=-HTOP FLXCJl(I)=HBOT*TINF FLXPJl{I)=-HBOT 330 CONTINUE ENDIF RETURN END С cccccccccccccccccccccccccccccccccccccccccccccccccccccccc 10.4.5. Результаты расчетов RESULTS OF CONDUCT FOR CARTESIAN COORDINATE SYSTEM STAGGERED FIN ARRAY
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 |
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||