Дифференциальное уравнение теплопроводности.




Изучение любого физического явления сводится к установлению зависимости между величинами, характеризующими это явление. Для сложных физических процессов, в которых определяющие величины могут существенно изменяться в пространстве и времени, установить зависимость между этими величинами достаточно сложно. В таких случаях используют методы математической физики, которые заключаются в том, что ограничивается промежуток времени и из всего пространства рассматривается некоторый элементарный объем. Это позволяет в пределах выбранного объема и данного промежутка времени пренебречь изменениями величин, характеризующих процесс, и существенно упростить зависимость.

Выбранные таким образом элементарный объем dV и элементарный промежуток времени , в пределах которых рассматривается процесс, с математической точки зрения являются величинами бесконечно малыми, а с физической точки зрения – величинами еще достаточно большими, чтобы в их пределах можно было считать среду как сплошную, пренебрегая ее дискретным строением. Полученная таким образом зависимость является общим дифференциальным уравнением процесса. Интегрируя дифференциальные уравнения, можно получить аналитическую зависимость между величинами для всей области интегрирования и всего рассматриваемого промежутка времени.

Для решения задач, связанных с нахождением температурного поля, необходимо иметь дифференциальное уравнение теплопроводности.

Примем следующие допущения:

· тело однородно и изотропно;

· физические параметры постоянны;

· деформация рассматриваемого объема, связанная с изменением температуры, очень мала по сравнению с самим объемом;

· внутренние источники теплоты в теле, распределены равномерно.

В основу вывода дифференциального уравнения теплопроводности положим закон сохранения энергии, который сформулируем так:

Количество теплоты dQ, введенное в элементарный объем dV извне за время dτ вследствие теплопроводности, а также от внутренних источников, равно изменению внутренней энергии или энтальпии вещества, содержащегося в элементарном объеме.

где dQ1 – количество теплоты, введенное в элементарный объем dV путем теплопроводности за время ;

dQ2 – количество теплоты, которое за время выделилось в элементарном объеме dV за счет внутренних источников;

dQ – изменение внутренней энергии (изохорный процесс) или энтальпии вещества (изобарный процесс), содержащегося в элементарном объеме dV за время .

Для получения уравнения рассмотрим элементарный объем в виде кубика со сторонами dx, dy, dz (см. рис.1.2.). Кубик расположен так, чтобы его грани были параллельны соответствующим координатным плоскостям. Количество теплоты, которое подводится к граням элементарного объема за время в направлении осей x, y, z обозначим соответственно dQx, dQy, dQz.

Количество теплоты, которое будет отводиться через противоположные грани в тех же направлениях, обозначим соответственно dQx+dx, dQy+dy, dQz+dz.

Количество теплоты, подведенное к грани dxdy в направлении оси x за время , составляет:

где qx – проекция плотности теплового потока на направление нормали к указанной грани. Соответственно количество теплоты, отведенное через противоположную грань будет:

Разница между количеством теплоты, подведенном к элементарному объему, и количеством теплоты, отведенного от него, представляет собой теплоту:

Функция q является непрерывной в рассматриваемом интервале dx и может быть разложена в ряд Тейлора:

Если ограничиться двумя первыми слагаемыми ряда, то уравнение запишется в виде:

Аналогичным образом можно найти количество теплоты, подводимое к объему в направлении двух других координатных осей y и z.

Количество теплоты dQ, подведенное в результате теплопроводности к рассматриваемому объему, будет равно:

Второе слагаемое определим, обозначив количество теплоты, выделяемое внутренними источниками в единице объема среды в единицу времени qv и назовем его мощностью внутренних источников теплоты [Вт/м3], тогда:

Третья составляющая в нашем уравнении найдется в зависимости от характера ТД процесса изменения системы.

При рассмотрении изохорного процесса вся теплота, подведенная к элементарному объему, уйдет на изменение внутренней энергии вещества, заключенного в этом объеме, т.е. dQ=dU.

Если рассматривать внутреннюю энергию единицы объема u=f(t,v), то можно записать:

, Дж/м3

, Дж/кг

где cv изохорная теплоемкость или единицы объема или единицы массы, [Дж/м3];

ρ – плотность, [кг/м3].

Соберем полученные выражения:

Полученное выражение является дифференциальным уравнением энергии для изохорного процесса переноса теплоты.

Аналогично выводится уравнение для изобарного процесса. Вся теплота, подведенная к объему уйдет на изменение энтальпии вещества, заключенного в объеме.

Полученное соотношение является дифференциальным уравнением энергии для изобарного процесса.

В твердых телах перенос теплоты осуществляется по закону Фурье , значение теплоемкости можно принять . Напомним, что проекция вектора плотности теплового потока на координатные оси определяются выражениями:

Последнее выражение называют дифференциальным уравнением теплопроводности. Оно устанавливает связь между временным и пространственным изменениями температуры в любой точке тела, в котором происходит процесс теплопроводности.

Наиболее общее дифференциальное уравнение теплопроводности в частных производных имеет такую же форму, но в нем величины ρ, l, с являются функциями времени и пространства. Это уравнение описывает большое количество задач теплопроводности, представляющих практический интерес. Если принять теплофизические параметры постоянными, то уравнение будет проще:

Обозначим , тогда:

Коэффициент пропорциональности а2/с] называется коэффициентом температуропроводности и является физическим параметром вещества. Он существенен для нестационарных тепловых процессов характеризует скорость изменения температуры. Если коэффициент теплопроводности характеризует способность тел проводить теплоту, то коэффициент температуропроводности является мерой теплоинерционных свойств тела. Например, жидкости и газы обладают большей тепловой инерционностью и, следовательно, малым коэффициентом температуропроводности, а металлы наоборот имеют малую тепловую инерционность.

Далее, если система тел не содержит внутренних источников теплоты q=0, тогда дифференциальное уравнение принимает форму уравнения Фурье:

Если имеются внутренние источники теплоты, а температурное поле является стационарным, то мы получаем уравнение Пуассона:

Наконец, при стационарной теплопроводности и отсутствии внутренних источников теплоты мы получаем уравнение Лапласа:

Условия однозначности для теплопроводности.

Так как дифференциальное уравнение теплопроводности выведено из общих законов физики, то оно описывает целый класс явлений. Для его решения необходимо задать граничные условия или условия однозначности.

Условия однозначности включают:

1. геометрические условия – характеризуют форму и размеры тела;

2. физические условия – характеризуют физические свойства среды и тела;

3. начальные (временные) условия – характеризуют распределение температур в теле в начальный момент времени, задаются при исследовании нестационарных процессов;

4. граничные условия – характеризуют взаимодействие рассматриваемого тела с окружающей средой.

Граничные условия могут быть заданы несколькими способами.

Граничные условия первого рода. Задается распределение температуры на поверхности тела для каждого момента времени:

tc = f(x,y,z,τ)

где tc – температура на поверхности тела;

x,y,z – координаты поверхности тела.

В частном случае, когда температура на поверхности является постоянной на протяжении всего времени протекания процессов теплообмена, уравнение упрощается:

tc = const

Граничные условия второго рода. Задаются значения теплового потока для каждой точки поверхности тела и любого момента времени. Аналитически выглядит так:

qc = f(x,y,z,τ)

В простейшем случае плотность теплового потока по поверхности тела остается постоянной. Такой случай имеет место при нагревании металлических изделий в высокотемпературных печах.

Граничные условия третьего рода. При этом задаются температура окружающей среды tср и закон теплообмена между поверхностью тела и средой. Для описания процесса теплообмена используется закон Ньютона-Рихмана. Согласно этому закону количество теплоты, отдаваемое или принимаемое единицей поверхности тела в единицу времени, пропорционально разности температур поверхности тела и среды:

где α коэффициент пропорциональности, называется коэффициентом теплоотдачи [Вт/(м2·К)], характеризует интенсивность теплообмена. Численно он равен количеству теплоты, отдаваемому единицей поверхности тела в единицу времени при разности температур равной одному градусу. Согласно закону сохранения энергии количество теплоты, которое отводится окружающей среде, должно равняться теплу, подводимому вследствие теплопроводности из внутренних частей тела, то есть:

Последнее уравнение является граничным условием третьего рода.

Встречаются более сложные технические задачи, когда ни одно из перечисленных условий задать невозможно, и тогда приходится решать задачу методом сопряжения. При решении такой задачи должны выполняться условия равенства температур и тепловых потоков по обе стороны от границы раздела. В общем случае условия сопряженности можно записать:

Решение сопряженной задачи связано с нахождением температурных полей по обе стороны границы раздела.



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2018-01-31 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: