Внутренняя задача теплообмена
Внутренняя задача теплообмена – это определение температуры по всей толщине защитного слоя. Существует три вида теплопередачи — лучистая теплопередача, конвекция и теплопроводность. Лучистая теплопередача и конвекция не учитываются в силу малости их величин. При этих упрощающих допущениях проведем решение задачи в одномерной модели:
С - теплоемкость
Область определения:
L – толщина материала Зададим начальные условия (распределение поля температур Т):
Граничные условия на внешней границе определим из решения внешней задачи (22):
Внутренние граничные условия:
При этом на температуру поверхности накладываем ограничения: Tmax
Для решения системы дифференциальных уравнений (25) используем метод конечных разностей (нестационарная задача). Этот метод основан на том, что искомая характеристика ищется не в виде непрерывной функции, а в узлах области поиска решения:
Для аппроксимации системы (25) осуществляется по неявной схеме. Геометрическое представление неявной схемы приведено на рисунке 1. Шаги по пространству и времени рассчитываются по формулам (2.5.4).
0 (рис.1)
N – количество узлов по пространству М – количество узлов по времени Нумерация узлов по пространству: Индекс i – по времени j, тогда в соответствующих узлах обозначение температуры и координат будут иметь вид:
Переходим от дифференциальных операторов к алгебраическим:
За счет подстановки формул (27) и (28) в систему (25) переходим к ее алгебраической записи:
Приведем систему уравнений (29 - 31) к виду 3-х диагональной матрицы:
Для этого необходимо расчисть разностные коэффициенты входящие в (32)-(34). Формулы для расчета разностных коэффициентов уравнения (33) могут быть получены путем перегруппировки слагаемых в уравнении (30):
Для расчета разностных коэффициентов в уравнении (33) выразим T3 из уравнения (30) для второго узла.
и подставим в (30). После преобразования получим:
Сопоставим полученное уравнение с (37) и перегруппировав слагаемые в уравнении получим формулы для расчета разностных коэффициентов уравнения (32)
Подобным образом преобразуем уравнение (34) и получаем формулы для расчета разностных коэффициентов:
Для расчета температуры в узлах разностной сетки на текущем временном слое используется метод прогонки. В методе прогонки для расчета температуры применяется формула:
где
Уравнения (40) и (32) с математической точки зрения формально одинаковы. Сопоставив их между собой, найдем выражения для расчета прогоночных коэффициентов
Для расчета прогоночных коэффициентов
В этом уравнении выразим Т1 через формулу (41) с уже рассчитанными прогоночными коэффициентами. Произведя группировку слагаемых при Т2 и Т3 получим формулы для вычисления прогоночных коэффициентов во втором узле: Аналогичные соображения применяются для получения формул во всех последующих узлах. В этом случае формулы для расчета прогоночных коэффициентов можно записать в универсальном виде:
Для расчета значения температуры ТN произведем подстановку в уравнение (34) формулы (41) для N-1 узла:
Рассчитываем поле температур: Ti = Переходим к расчету профиля температуры на следующем временном слое. Таким образом, получим значение Tij, i =1, N, j = 1, M во всех узлах области поиска решения.
Воспользуйтесь поиском по сайту: ![]() ©2015 - 2025 megalektsii.ru Все авторские права принадлежат авторам лекционных материалов. Обратная связь с нами...
|