Главная Журналы 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
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 ******
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 |