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

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

RETURN

ENTRY PHI

DO 300 J=2,M2

DO 300 1=2,L2

ALAM(I,J)=RHOCP GAM(I,J)=COND

IF(Y(J).LT.0.25) SC(I,J)=SOURCE 300 CONTINUE

COME HERE TO SPECIFY BOUNDARY CONDITIONS DO 310 J=2,M2 KBCIl(J)=2 KBCLl(J)=2 FLXCIl(J)=QB 310 CONTINUE

DO 320 1=2,L2 KBCMl(I)=2 FLXCMl(I)=HE*TINF FLXPMl(I)=-HE 320 CONTINUE RETURN END

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

8.5.5. Результаты расчетов

RESULTS OF CONDUCT FOR POLAR COORDINATE SYSTEM UNSTEADY CONDUCTION WITH HEAT GENERATION

ITER

TIME

T(4,4) T{10,4) T{7,10) TdClO) QBTM

QTOP

.OOE+00

.OOE+01

.OOE+01

. OOE-hOI

.OOE+01

.OOE+00

.ODE+00

.20E-

. 13E-H01

. 12Е-Ю1

.0,0Е-Ю1

.OOE+01

.7IE+02

.88E+02

. 64E-

.29E+01

.29E+01

.OOE-HOl

.OOE+01

.02E+02

.71E+02

.37E-

.52E+01

. 51Е-Ю1

. 99Е-Ю1

.99E+01

.21E+02

.55E+02

.44E-

. 80E-H01

.77E+D1

. 98E+01

.98E+0I

.OOE+02

.40E+02

. 93E-

. I3E-b01

.09E+01

.96E+01

.96E+01

.18E+02

.26E+02

. 19E-

. 50Е-Ю1

.43E+01

.94Е-Ю1

,94E+01

.60E+02

.I3E+02

.55E-

.90ЕЧ-01

. 79Е-Ю1

. 90E1-01

, 90E+01

.15E+02

.OOE+02

.98E-

.30E+01

.16E+01

.85E+01

.85E+01

. 95E+01

.88E+02

.50E-

.71E+01

.53E+01

. 79E-H01

.79E+01

.93E+01

.77E+02

. 12E-

.llE+01

.88E+01

.73Е-Ю1

.72E+0I

.33E+01

.67E+02

.86E-

.50E+01

.22Е-Ю1

. 67E-H01

.66E+01

.30E-01

.57E+02

.75E-

.86E+01

. 53Е-Ю1

.62E-H01

.61E+01

.97E+01

.49E+02

.82E-

.20E+01

,81E+0I

.59E+01

,58E+01

.74E+01

.42E+02

. lOE-

.51Е-Ю1

. 06E+01

,58E+01

,56E+01

,29E+01

.37E+02

.64E-

.79ЕЧ-01

.29E-H0I

, 60Е-Ю1

,57E+01

, 63E+01

. 34E+02

.05E-

.00E-f02

.48Е-Ю1

. 66E+01

.61E+01

-7 .

.77E+01

.34E+02

.27E-

.03E+02

.64E+01

.73E+01

.66E+01

.74E+01

.36E+02

.54E-

. 04Е-Ю2

.78E-H01

.83E-H01

.74E+01

.56E+01

.40E+02

. 86E-

,06E+02

.90E+0I

.94E+01

,83E+01

.02E+02

.45E+02

.24E-

,07E+02

.OOE+02

, 06Е-Ю1

,92E+01

.08E+02

.51E+02

.70E-

, 09Е-Ю2

.OIE+02

, 18Е-Ю1

,02E+01

, 13E+02

.57E+02

.25E-

, 10ЕЧ-02

.02Е-Ю2

,28E-l-01

,lOE+OI

, 17E+02

. 62E+02

. 91E-

, 10Е-Ю2

,02E+02

,37Е-Ю1

I8E+01

,20E+02

.67E+02 157



71E-01

llE+02

03E+02

44E+01

24E+01

23E+02

71E+02

ITER

TIME

T(4,4) T(10,4) T(7,10) T(10,10) QBTM

QTOP

66E-01

llE+02

03E+02

49E+01

29E+01

24E+02

74E+02

81E-01

12E+02

03E+02

53E+01

32E+01

26E+02

7 6E+02

18E-01

12E+02

03E+02

55E+01

34E+01

26E+02

77E+02

83E-01

12E+02

03E+C2

56E+01

35E+01

27E+02

78E+02

1.18E+00

12E+02

03E+02

57E+01

36E+01

27E+02

78E+02

42E+00

12E+02

03E+02

57E+01

36E+01

27E+92

78E+02

70E+00

12E+02

03E+02

58E+01

36E+01

27E+02

78E+02

04E+00

12E+02

03E+02

58E+01

36E+01

27E+02

78E+02

46E+00

12E+02

03E+02

58E+01

36E+01

27E+02

78E+02

95E+00

12E+02

03E+02

58E+01

36E+01

27E+02

78E+02

54E+00

12E+02

03E+02

58E+01

36E+01

27E+02

78E+02

TH= 0.OOE+OO 1= 8

TH= 9.37E-01 J= 1

Y= 0.OOE+OO J= 8

y= 5.42E-01

4.96E-02 1.63E-01 2.99E-01 4.47E-01 6.03E-01 7.67E-01

9 10 11 12

l.llE+00 1.29E+00 1.43E+00 1.57E+00

2 3 4 5 6 7

92E-01 3.75E-01 4.58E-01 12 13 14

75E-01 9.58E-01 1.OOE+OO

4.17E-02 1.25E-01 2.08E-01 2.

9 10 11

6.25E-01 7.08E-01 7.92E-01 8.

****** TEMPERATURE ******

J 14

OOE+01

41E+01

99E+01

72E+01

56E+01

45E+01

38E+01

28E+01

90E+01

40E+01

09E+01

88E+01

75E+01

67E+01

18E+01

81E+01

23E+01

82E+01

56E+01

38E+01

27E+01

02E+01

66E+01

04E+01

58E+01

26E+01

05E+01

91E+01

81E+01

47E+01

84E+01

34E+01

98E+01

74E+01

58E+01

56E+01

24E+01

63E+01

llE+01

73E+01

46E+01

28E+01

30E+01

OOE+01

42E+01

90E+01

50E+01

22E+01

03E+01

OOE+02

74E+01

20E+01

70E+01

30E+01

02E+01

82E+01

07E+02

05E+02

98E+01

51E+01

13E+01

86E+01

67E+01

14E+02

12E+02

08E+02

03E+02

OOE+02

75E+01

57E+01

22E+02

19E+02

15E+02

12E+02

09E+02

07E+02

06E+02

21E+02

19E+02

16E+02

13E+02

llE+02

lOE+02

09E+02

llE+02

lOE+02

08E+02

07E+02

06E+02

05E+02

05E+02

OOE+01

OOE+02

OOE+02

OOE+02

OOE+02

OOE+02

OOE+02

1= 8

14 3.33E+01 13 3.61E+01 12 4.20E+01 11 4.82E+01 10 5.47E+01 9 6.16E+01 8 6.90E+01 7.69E+01 8.54E+01 9.45E+01 1.05E+02 1.08E+02 1.05E+02 1.OOE+02

3.30E+01 3.58E+01 4.15E+01 4.76E+01 5.40E+01 6.09E+01 6.82E+01 7.60E+01 8.45E+01 9.37E+01 1.04E+02 1.08E+02 1.05E+02 l.OOE+02

29E+01 56E+01 12E+01 72E+01 36E+01 6.04E+01 6.77E+01 7.55E+01 8.40E+01 33E+01 03E+02 07E+02 05E+02 OOE+02

28E+01 55E+01 llE+01 71E+01 34E+01 6.02E+01 6.75E+01 7.53E+01 8.38E+01 9.31E+01 1.03E+02 1.07E+02 1.04E+02 l.OOE+02

.OOE+01 .55E+01 . llE+01 .71E+01 .34E+01 6.02E+01 6.75E+01 7.53E+01 8.38E+01 9.30E+01 1.03E+02 1.07E+02 1.04E+02 5.OOE+01



8.5.6. Обсуждение результатов

Решение этой задачи демонстрирует некоторые интересные тенденции в изменении локальных температур. Из выводимых на печать температур в четырех точках значения Т{4,4) иТ{10,4) соответствуют точкам, расположенным вблизи внутренней поверхности, в то время как Т(7,10) иТ(10,10) - точкам, находящимся во внешней области. В начале процесса высокая температура внутренней поверхности, наличие теплового потока плотностью и выделение тепла приводят к увеличению температуры во внутренней области, и поэтому Т(4,4) иТ(10,4) монотонно возрастают. В то же время внешняя граница омывается более холодной жидкостью с температурой Too- В результате Т{7,10) иТ{10,10) постепенно уменьшаются в течение некоторого времени. Однако в конце концов эти температуры также начинают возрастать из-за все увеличивающегося влияния высокой температуры во внутренней области.

Тепловой поток QTOP на внешней поверхности всегда отрицателен, что свидетельствует о теплоотдаче в окружающую холодную жидкость. На внутренней поверхности поток QBTM сначала положителен, так как тепло поступает в область от границы с высокой температурой Tj. Однако впоследствии внутреннее тепловыделение приводит к тому, что температуры внутри области начинают превышать значения Г,, в результате чего QBTM становится отрицательным.

Заключительное стационарное распределение температуры демонстрирует ожидаемый эффект от заданных различных граничных условий. Самая высокая температура находится на границе с заданной плотностью теплового потока, в то время как самая низкая - около внешней границы, соприкасающейся с холодной внешней средой.

8.5.7. Заключительные замечания

Описание приложения программы CONDUCT к решению нестационарной задачи показывает, что дополнительные усилия, необходимые для этого класса задач, невелики. В то же время можно получать обширную и интересную информацию о поле температуры, зависящем от времени. В этом примере не рассматриваются какие-либо нерегулярности геометрии, однако нет причин, из-за которых вы бы не смогли применить технику, описанную в предыдущих примерах, для анализа нестационарной теплопроводности в областях сложной геометрической формы. Мы использовали граничные условия, которые остаются постоянными с течением времени. Если граничные температуры или тепловые потоки будут постоянно меняться,





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