- •Введение
- •Теория теплопроводности
- •Механизмы переноса тепла
- •Методы изучения физических явлений
- •Температурное поле
- •Тепловой поток. Закон Фурье
- •Коэффициент теплопроводности
- •Дифференциальное уравнение теплопроводности
- •Уравнение теплопроводности
- •Лекция № 2 Условия однозначности для процессов теплопроводности
- •Граничные условия:
- •Теплопроводность в стационарном режиме
- •Многослойная стенка
- •Переменный коэффициент теплопроводности
- •Линейная плотность теплового потока:
- •Критический диаметр цилиндрической стенки
- •Передача теплоты через шаровую стенку
- •Плотность теплового потока
- •Теплопроводность при наличии внутренних источников теплоты
- •Теплопроводность однородной пластины
- •Теплопроводность однородного цилиндрического стержня
- •Нестационарные процессы теплопроводности
- •Аналитическое описание процесса
- •Анализ полученного решения
- •Охлаждение длинного прямоугольного стержня
- •Охлаждение цилиндра конечной длины
- •Приближенные методы решения задач теплопроводности
- •Метод конечных разностей
- •Численные методы решения задач теплопроводности при нестационарном режиме
- •Метод конечных элементов
- •Исследование процессов теплопроводности методом аналогий
- •Электротепловая аналогия
- •Конвективный теплообмен Основные понятия и определения
- •Дифференциальные уравнения конвективного теплообмена
- •Уравнение энергии
- •Уравнения движения
- •Уравнение сплошности
- •Гидродинамический и тепловой пограничные слои
- •Уравнение теплоотдачи.
- •Тепловой пограничный слой
- •Подобие и моделирование процессов конвективного теплообмена
- •Приведение математической формулировки краевой задачи к записи в безразмерной форме
- •Аналогично преобразуем и уравнение движения
- •Граничные условия
- •Условия подобия физических процессов
- •Следствия из условий подобия
- •Обработка результатов опыта
- •Теплоотдача при вынужденном продольном омывании плоской поверхности
- •Интегральные уравнения пограничного слоя
- •Теплоотдача при ламинарном пограничном слое
- •Нагрев диэлектриков в электромагнитном поле. Электромагнитное поле в диэлектрике.
- •Плоская электромагнитная волна в диэлектрике
- •Используя соотношения
- •Решение уравнения (1) имеет вид
- •Виды поляризации
- •Зависимости ε и tg δ от частоты и температуры
- •Температурные характеристики релаксационной поляризации
- •Структура электродугового разряда
- •Характеристика приэлектродных областей и протекающих в них процессов
- •Выделяющаяся на аноде мощность
- •Мощность, выделяющаяся на катоде
- •Электроды дуговых установок
- •Термохимический катод
- •Основные закономерности электродугового столба
- •Особенности дуги переменного тока
- •Устойчивость и регулирование параметров электрической дуги
- •Способы зажигания дуги
- •Процессы переноса в дуговых и плазменных электротехнологических установках
- •Физико-технические основы электронно-лучевого нагрева
- •Основы расчета устройств формирования электронных пучков элу
- •Потери энергии электронного пучка и энергетический баланс элу
- •Основы лазерного нагрева. Основные принципы работы лазеров
- •Типы оптических квантовых генераторов
- •Лазеры твердотельные с оптической накачкой
- •Основы технологии светолучевой обработки
- •Литература
Нестационарные процессы теплопроводности
Общие положения
Это такие процессы теплопроводности, когда поле температуры в теле изменяется не только в пространстве, но и во времени (рис. 11). Они имеют место при нагревании или охлаждении различных энергетических устройств.
Среди практических задач нестационарной теплопроводности большое значение имеют 2 группы процессов: 1) тело стремиться к тепловому равновесию; 2) температура тела претерпевает периодические изменения.
Рис. 11 Характер изменения температуры тела во времени
В первом случае по мере нагрева температура в каждой точке асимптотически приближается к установившейся температуре.
Нестационарные тепловые процессы связаны с изменением внутренней энергии вещества. К примеру, в условиях пере дачи теплоты через стенку при внезапном изменении температуры одного из теплоносителей не вся теплота будет передаваться через стенку: часть уйдет на изменение внутренней энергии самой стенки (её температуры). Только при наступлении стационарного процесса вся теплота будет передаваться через стенку от одной жидкости к другой.
В настоящем разделе рассмотрим наиболее важные задачи, относящиеся к процессам, в которых тело стремиться к тепловому равновесию.
Аналитическое описание процесса
Рассмотрим задачи без внутренних источников тепла. Аналитическое описание процесса теплопроводности включает в себя дифференциальное уравнение и условия однозначности.
Дифференциальное уравнение теплопроводности при отсутствии внутренних источников теплоты имеет вид:
, (1)
где - коэффициент температуропроводности, м2/с (характеризует скорость изменения температуры).
Условия однозначности задаются в виде:
- физических параметров l, с, r;
- формы и геометрических размеров объекта l0, l1,…, ln; (2)
- температуры тела в начальный момент времени: при t = 0 Т0 = f (x, y, z).
Граничные условия могут быть заданы в виде граничных условий 3-го рода:
.
Дифференциальное уравнение совместно с условиями однозначности дает законченную математическую формулировку рассматриваемой задачи. Решение её заключается в отыскании функции
Т = f (x, y, z, t, а, Т0, Тж, l0, l1,… ln),
которая удовлетворяла бы уравнению (1) и условиям (2).
Рассмотрим подробно решение задачи охлаждения плоской однородной стенки. Изучив метод решения задачи, можно понять принцип решения задач и для тел другой геометрической конфигурации.
Охлаждение (нагревание) неограниченной пластины
Дана пластина толщиной 2d. Если толщина пластины мала по сравнению с длиной и шириной, то такую пластину считают неограниченной. Коэффициент теплоотдачи a одинаков для всех точек пластины (рис. 12). Изменение температуры происходит в направлении х. В пространстве задача является одномерной.
Рис. 12 К охлаждению плоской неограниченной пластины. При τ = 0 задано Т0=const и θ0=const.
Начальное распределение температуры задано некоторой функцией
.
Охлаждение происходит в среде с постоянной температурой Тж = const. Отсчет температуры пластины для любого момента времени будем вести от температуры окружающей среды Тж, то есть
.
Дифференциальное уравнение теплопроводности имеет вид:
. (1)
Начальные условия: при t = 0 θ = θ 0 = f (x) – Тж = F(х).
При заданных условиях охлаждения задача становиться симметричной и начало координат удобно поместить на оси пластины.
При этом граничные условия:
а) на оси пластины при х = 0 ;
б) на поверхности пластины при х = d .
Для решения дифференциального уравнения применяем метод разделения переменных. При этом решение дифференциальное уравнения ищем в виде произведения двух функций, из которых одна является функцией только времени t, а другая – только х:
.
После подстановки этого выражения в уравнение (1), получим:
или
.
В этом уравнение легко разделяются переменные
. (2)
Левая часть – функция только t, правая только х. Если зафиксировать аргумент х и менять только t, то при любом его значении левая часть уравнения (2) равна постоянной величине, стоящей в правой части, то есть . Аналогично при фиксации t и изменении х правая часть уравнения для любого значения х должна равняться постоянной левой части, которая зависит только от t, то есть .
Так как равенство (2) должно иметь место при любых значениях х и t, то обе его части должны быть равны одной и той же постоянной величине.
.
Нетривиальное решение для функции y(х) только при e < 0. Положим, что e = - k2 :
.
Система обыкновенных дифференциальных уравнений:
; (3)
. (4)
Постоянная k определяется из граничных условий, а знак «минус» выбирается из физических соображений. Для тепловых процессов, стремящихся к тепловому равновесию, знак может быть только минус.
Уравнению (3) удовлетворяет функция .
Уравнение (4) – функция .
В результате получили общее решение
. (5)
Для того чтобы уравнение (5) было решением поставленной задачи, его нужно подчинить начальным и граничным условиям.
При х = 0 , находим
Þ Þ с2 = 0
Следовательно, частное решение – должно быть отброшено как не удовлетворяющее граничным условиям.
Если учесть, что с2 = 0 и обозначить с1с3 = А, то уравнение (5) примет вид:
.
При х = d Þ ½ умножив и разделив на d, получим:
,
- число Био – безразмерный показатель (характеризует соотношения внутреннего и внешнего тепловых сопротивлений).
Если обозначить kd = m, то:
(1)
Из анализа этого уравнения следует, что при каждом значении Вi существует бесконечное множество решений. Наиболее просто это равнение решается графическим способом.
Рис. 13 К решению уравнения (1).
Обозначим через y1 = сtg m, .
Пересечение котангенсоиды у1 с прямой у2 дает бесконечное множество корней характеристического уравнения m1 < m2 < m3 <…<mn (рис. 13).
Каждому значению Bi отвечает своя совокупность корней уравнения (1). Первые четыре корня такого уравнения приводятся в таблице.
При Вi ® ¥ (внутреннее сопротивление велико по сравнению с внешним) у2 = 0 – совпадает с осью х и корни будут равны:
m1 = p/2; m2 = 3p/2; mn = (2n-1) p/2
При Вi ® 0 (внутреннее сопротивление мало по сравнению с внешним) прямая совпадает с осью ординат и тангенс угла наклона стремиться к ¥ ® корни равны:
m1 = 0; m2 = p,…, mn = (n-1) p.
Следовательно, каждому найденному значению корня m будет соответствовать свое частное распределение температур:
, здесь мы учли, что k =μ/δ.
Путем наложения бесконечного числа таких распределений температур можно получить истинное распределение:
. (а)
Постоянная Аn находится из начальных условий:
.
Это есть разложение четной функции в ряд Фурье. Есть специальные формулы для определения коэффициентов Аn..
Если в начальный момент времени t = 0 температура в любой точке пластины распределена равномерно (Т0 – Тж = θ0 = const), то:
.
Подставляя Аn в выражение (а), получим
. (б)
Уравнение температурного поля (б) целесообразно представить в безразмерной форме. Для этого разделим правую и левую части уравнения на θ0 (начальная разность температур). При этом обозначим: Dn = Аn / θ0. Получим:
, (в)
где Q = θ/θ0 – безразмерная температура; Х = х/d - безразмерная координата; Fo =aτ/δ2 – число Фурье, представляющее собой безразмерное время; Dn = Аn / θ0 – безразмерный коэффициент.
Получим, что температура каждой точки во времени изменяется по экспоненциальному закону. Распределение температуры по координате х (по толщине) – имеет вид косинусоиды с максимумом в центре пластины.