Главная Журналы GEN=0. DO 233 J=2,M2 DO 233 1=2,L2 IF(X(I).GT.0.7.AND.Y(J).GT.0.25) GO TO 233 GEN=GEN+(1000.-4.E-5*T(I,J)**3)*XCV(I)*YCV(J) 2 33 CONTINUE HTBAL=HTOUT-GEN CONSTRUCT FINAL PRINTOUT DO 240 IUNIT=IU1,IU2 WRITE(lUNIT,250) 250 FORMAT(IX,/, J,8X,Y(J),5X,FLUX(LEFT),4X, 1 FLUX(RIGHT) •) DO 260 J=M2,2,-1 WRITE(lUNIT, 27 0) J,Y(J) , FLUXIl(J,1),FLUXLl(J,1) 270 F0RMAT(1X,I2,5X,1PE9.2,3X,1PE9.2,5X,1PE9.2) 2 60 CONTINUE WRITE(lUNIT,280) 280 FORMAT(IX,/,• I,8X,X(I),5X,FLUX(BOTTOM),2X, 1 FLUX(TOP)) DO 290 1=2,L2 WRITE(lUNIT,27 0) I,X(I),FLUXJl(1,1),FLUXMl(1,1) 2 90 CONTINUE WRITE(IUNIT,2 91) HTOUT,GEN,HTBAL 291 FORMAT(/IX,OVERALL HEAT BALANCE/IX,20(*)/IX, 1 HEAT OUTFLOW,5X,GENERATION,6X,DIFFERENCE/2Х, 2 1PE10.3,6X,1PE10.3,5X,1PE10.3) 24 0 CONTINUE CALL PRINT CALL PLOT KST0P=1 ENDIF RETURN ENTRY PHI DO 300 J=2,M2 DO 300 1=2,L2 IF(X(I).GT.0.7.AND.Y(J).GT.0.25) THEN GAM(I,J)=1.+0.01*T(I, J) ELSE GAM(I,J)=AK1 TP2=4.E-5*T(I,J)**2 SC(I,J)=1000.+2.*TP2*T(I, J) SP(I,J)=-3.*TP2 ENDIF 300 CONTINUE COME HERE TO SPECIFY BOUNDARY CONDITIONS 136 do 310 1=2,l2 квел (i) =2 flxcjl(i)=qw if(x(i).lt.0.7) kbcm1(i)=2 310 continue do 320 j=2,m2 kbcll(j)=2 flxcl1(j)=he-*tinf flxpll(j)=-he 320 continue return end с cccccccccccccccccccccccccccccccccccccccccccccccccccccccc 8.2.5. Результаты расчетов results of conduct for cartesian coordinate system ************************************************** steady conduction with mixed boundary conditions
OVERALL HEAT BALANCE HEAT OUTFLOW GENERATION DIFFERENCE
8.2.6. Обсуждение результатов В этой задаче нелинейность была достаточно слабая и источниковый член был оптимально линеаризован. Поэтому видно, что решение сошлось очень быстро. Распечатка результатов дает нам значения заданных и искомых локальных плотностей потоков на всех границах области. Для адиабатического участка верхней границы полученная плотность потока очень мала, но не равна нулю. Эти малые числа являются результатом погрешностей округления при расчетах. При анализе температурного поля видно, что локальные температуры превышают температуры на границах. Это обусловлено внутренним выделением тепла и заданными граничными плотностями тепловых потоков. На нижней границе температура довольно высокая, что связано со значеннем заданной плотности теплового потока на этой границе. 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 |