ПОИСК Статьи Чертежи Таблицы Построение подпрограммы из "Численное решение задач теплопроводности и конвективноного теплообмена при течении в каналах " Рассчитаем зависящее от времени распределение температуры в теле, в котором мгновенно изменили граничные условия и увеличили выделение тепла. Решение этой задачи — первый опыт работы в полярной системе координат (0, г). [c.152] Учитывая симметричность геометрических форм и граничных условий, ограничим наши вычисления левой половиной области, как показано на рис. 8.6. [c.153] Зададим значения соответствующих переменных [см. (8.8)— (8.8в)]. Переменной DT присвоено начальное значение At, которое впоследствии будет увеличиваться. Массив температур Т (I, J) заполнен начальным значением Tq. Для стационарных задач эти начальные значения представляют собой лишь начальное приближение, а для нестационарных задач заданные в BEGIN значения Т 1, J) для внутренних точек являются известными начальными температурами. [c.154] Наконец, на внутренней границе зададим температуру Г,. [c.154] В этом примере на каждом временном шаге выводится на печать четыре характерные температуры и суммарные тепловые потоки на внутренней и внешней поверхностях тела. При вычислении этих потоков величины X V I) RV 2 ) и X V I) RV (Ml) представляют собой площади (при единичной глубине ) граней контрольных объемов на внешней и внутренней границе соответственно. [c.154] Определение граничных условий заключается в задании конвективного теплообмена на внешней границе, известного теплового потока на плоской поверхности и условия равенства нулю потока на вертикальной линии симметрии. Все эти величины задаются через соответствующие значения КВС, FLX и FLXP. [c.155] Вернуться к основной статье