Аналитическое описание процесса

Рассмотрим задачи без внутренних источников тепла. Аналитическое описание процесса теплопроводности включает в себя дифференциальное уравнение и условия однозначности.

Дифференциальное уравнение теплопроводности при отсутствии внутренних источников теплоты имеет вид:

, (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 – безразмерный коэффициент.

Получим, что температура каждой точки во времени изменяется по экспоненциальному закону. Распределение температуры по координате х (по толщине) – имеет вид косинусоиды с максимумом в центре пластины.



php"; ?>