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

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

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

iter

t(3,3)

t(5,4)

t(10,7)

heat flow (left face

l.ooe+02

ooe+02

ooe+02

o.oooe+00

1.29e+02

38e+02

38e+01

-4.948e+02

1.30e+02

41e+02

70e+01

-5.158e+02

1.30e+02

41e+02

70e+01

-5.126e+02

1.30e+02

41e+02

70e+01

-5.127e+02

1.30e+02

41e+02

70e+01

-5.127e+02

1.30e+02

41e+02

70e+01

-5.127e+02

1.30e+02

41e+02

70e+01

-5.127e+02

y(j)

flux(left)

flux(right)

58e-01

-6.77e+02

-2.45e+02

75e-01

-7.15e+02

-4.33e+02

92e-01

-7.97e+02

-6.83e+02

08e-01

-9.42e+02

-1.20e+03

25e-01

-1.21e+03

-1.34e+03

пе-02

-1.81e+03

-1.51e+03

x(i)

flux(bottom)

flux(top)

ooe-02

8.00e+02

2.82e-13

50e-01

8.ooe+02

-5.76e-13

50e-01

8.ooe+02

2.82e-12

50e-01

8.ooe+02

5.40e-13

50e-01

8.ooe+02

-5.97e-13

50e-01

8.ooe+02

5.26e-13

50e-01

8.ooe+02

-8.53e-14

50e-01

8.ooe+02

-1.50e+03

50e-01

8.ooe+02

-5.83e+02

50e-01

8.ooe+02

-1.60e+02



OVERALL HEAT BALANCE

HEAT OUTFLOW GENERATION DIFFERENCE

885E+02

3.885ЕЧ02

-8.809E

O.OOE+00

.OOE-02

.50E-01

.50E-01

.50E-01

.50E-01

50E-01

6.50E-01

.50E-01

.50E-01

.50E-01

.OOE+OO

0.OOE+OO

.17E-02

.25E-01

.08E-01

.92E-01

.75E-01

58S-01

J~ Y=

5.00E-01

**** TEMPERATURE

l.OOE+02

06E+02

17E+02

25E+02

29E+02

28E+02

21E+02

1.00Ef02

06E+02

18E+02

1.25E+02

29E+02

2BE+02

21E+02

l.OOE+02

07E+02

19E+02

27E+02

31E+02

30E+02

25E+02

1.00E+02

08E+02

21E+02

30E+02

35E+02

35E+02

31E+02

l.OOE+02

09E+02

24E+02

35E+02

41E+02

42E+02

38E+02

l.OOE+02

llE+02

30E+02

42E+02

49E+02

51E+02

47E+02

1. OOE+02

16E+02

39E+02

52E+02

60E+02

61E+02

5BE+C2

l.OOE+02

22E+02

45E+02

59E+02

66E+02

68E+02

64E+02

07E+02

OOE+01

OOE+01

OOE+01

OOE+02

08E+02

48E+01

70E+01

6QE+01

73E+01

14E+02

80E+01

42E+01

12E+01

67E+01

22E+02

06E+02

48E+01

88E+01

92E+01

30E+02

17E+02

87E+01

66E+01

48E+01

39E+02

26E+02

08E+02

50E+01

21E+01

50E+02

37E+02

19E+02

50E+01

07E+01

56E+02

43E+02

25E+02

OlE+02

OOE+02

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