Дифференциальное уравнение теплопроводности
При решении задач, связанных с нахождением температурного поля, необходимо иметь дифференциальное уравнение теплопроводности.
Температурное поле – совокупность значений температур во всех точках рассматриваемого пространства для каждого момента времени .
Для упрощения вывода этого дифференциального уравнения сделаны следующие допущения:
– физические параметры постоянны;
– деформация рассматриваемого объема, связанная с изменением температуры, является очень малой величиной по сравнению с самим объемом;
– внутренние источники теплоты в теле распределены равномерно.
В основу вывода дифференциального уравнения теплопроводности положен закон сохранения энергии в формулировке:
количество теплоты dQ, введенное в элементарный объем извне за время dτ теплопроводностью, а также от внутренних источников, равно изменению внутренней энергии или энтальпии вещества (в зависимости от рассмотрения изохорного или изобарного процесса), содержащегося в элементарном объеме.
(*)
где dQ1 – количество теплоты, Дж, введенное в элементарный объем теплопроводностью за время dτ;
dQ2 – количество теплоты, Дж, которое за время dτ выделилось в элементарном объем за счет внутренних источников;
dQ – изменение внутренней энергии или энтальпии вещества, содержащегося в элементарном объеме , за время dτ.
Для нахождения составляющих выделим в теле элементарный параллелепипед со сторонами dx, dy, dz. Параллелепипед расположен так, чтобы его грани были параллельны соответствующим координатным плоскостям.
Количество теплоты, которое подводится к граням элементарного объема за время dτ в направлении осей Оx, Оy, Оz обозначим соответственно dQx, dQy, dQz.
Количество теплоты, которое будет отводиться через противоположные грани в тех же направлениях, обозначим соответственно dQx+dx, dQy+dy, dQz+dz.
Количество теплоты, подведенное к грани dydz=dF в направлении оси Ох за время dτ, составляет ,
где qx – проекция плотности теплового потока на направление нормали к указанной грани.
Количество теплоты, отведенное через противоположную грань элементарного параллелепипеда в направлении оси Ох
.
Разница количеств теплоты, подведенного к элементарному параллелепипеду и отведенного от него за время dτ в направлении оси Ох
Функция является непрерывной в рассматриваемом интервале dx и может быть разложена в ряд Тейлора
Если ограничиться двумя первыми членами ряда:
Аналогично можно найти количество теплоты, подводимое к элементарному объему в направлениях двух других координатных осей Oy и Oz.
Количество теплоты dQ, подводимое теплопроводностью к рассматриваемому объему, будет равно
Обозначим через , Вт/м 3 , количество теплоты, выделяемое внутренними источниками в единице объема в единицу времени.
Тогда
Третья составляющая уравнения (*) найдется в зависимости от характера термодинамического процесса изменения системы.
В случае рассмотрения изохорного процесса вся теплота, подведенная к элементарному объему, уйдет на изменения внутренней энергии вещества, заключенного в этом объеме, т.е.
где – изохорная теплоемкость единицы массы, Дж/(кг·К);
ρ – плотность вещества, кг/м 3 .
Подставляя полученные выражения в уравнение (*), получим
,
Проекции вектора плотности теплового потока на координатные оси Ох, Оу, Оz определяются законом Фурье:
; ; .
где λ – коэффициент теплопроводности (физический параметр вещества, характеризующий способность проводить теплоту), Вт/(м∙°С).
Подставляя полученные выражения проекций вектора плотности теплового потока в уравнение (*), опуская индекс при с, ипринимая теплофизические характеристики постоянными, получим
(***)
Выражение (***) называется дифференциальным уравнением теплопроводности. Оно устанавливает связь между временнЫм и пространственным изменением температуры в любой точке тела.
и
Тогда выражение (***) имеет вид:
Выражение (***) в цилиндрической системе координат:
где r – радиус-вектор;
φ – полярный угол;
Коэффициент пропорциональности а, м 2 /с, называется коэффициентом температуропроводности и является физическим параметром вещества.
Он характеризует скорость изменения температуры, т.е. является мерой теплоинерционных свойств тела. Поэтому при прочих равных условиях выравнивание температур во всех точках пространства будет происходить быстрее в том теле, которое обладает бόльшим коэффициентом температуропроводности.
Коэффициент температуропроводности зависит от природы вещества.
Например, жидкости и газы обладают большой тепловой инерционностью и, следовательно, малым коэффициентом температуропроводности.
Металлы обладают малой тепловой инерционностью, т.к. они имеют большой коэффициент температуропроводности.
Если система тел не содержит внутренних источников теплоты (qυ=0), то
Если имеются внутренние источники теплоты, но температурное поле соответствует стационарному состоянию, т.е. , то
При рассмотрении изобарного процесса вся теплота, подведенная к объему, уйдет на изменение энтальпии вещества, заключенного в этом объеме:
(**)
Если рассматривать энтальпию единицы объема как , то
где сp – изобарная теплоемкость единицы массы, Дж/(кг·К).
В итоге (**) имеет вид:
| | следующая лекция ==> | |
СПОСОБЫ ПЕРЕНОСА ТЕПЛОТЫ | | | Условия однозначности для процессов теплопроводности. Дифференциальное уравнение теплопроводности описывает явление теплопроводности в самом общем виде |
Дата добавления: 2016-02-09 ; просмотров: 4555 ; ЗАКАЗАТЬ НАПИСАНИЕ РАБОТЫ
Лекция 4. Вывод уравнения теплопроводности
При построении математической модели распространения тепла в стержне сделаем следующие предположения:
1) стержень сделан из однородного проводящего материала с плотностью ρ;
2) боковая поверхность стержня теплоизолирована, то есть тепло может распространяться только вдоль оси ОХ;
3) стержень тонкий — это значит, что температура во всех точках любого поперечного сечения стержня одна и та же.
Рассмотрим часть стержня на отрезке [х, х + ∆х] (см. рис. 6) и воспользуемся законом сохранения количества тепла:
Общее количество тепла на отрезке [х, х + ∆х] = полному количеству тепла, прошедшему через границы + полное количество тепла, образованного внутренними источниками.
Общее количество тепла, которое необходимо сообщить участку стержня, чтобы повысить его температуру на ∆U, вычисляется по формуле: ∆Q= CρS∆x∆U, где С — удельная теплоемкость материала ( = количеству тепла, которое нужно сообщить 1 кг вещества, чтобы поднять его температуру на 1°), S — площадь поперечного сечения.
Количество тепла, прошедшее через левый конец участка стержня за время ∆t (тепловой поток) вычисляется по формуле: Q1 = -kSUx(x, t)∆t, где k — коэффициент теплопроводности материала ( = количеству тепла, протекающего в секунду через стержень единичной длины и единичной площади поперечного сечения при разности температур на противоположных концах, равной 1°). В этой формуле особого пояснения требует знак минус. Дело в том, что поток считается положительным, если он направлен в сторону увеличения х, а это, в свою очередь, означает, что слева от точки х температура больше, чем справа, то есть Ux CpS∆x∆U = kSUx(x + ∆х, t) ∆t — kSUx(x, t)∆t.
Если это равенство поделить на S∆x∆t и устремить ∆х и ∆t к нулю, то будем иметь:
Отсюда уравнение теплопроводности имеет вид
Ut = a 2 Uxx,
где — коэффициент температуропроводности.
В случае, когда внутри стержня имеются источники тепла, непрерывно распределенные с плотностью q(x,t), получится неоднородное уравнение теплопроводности
Начальные условия и граничные условия.
Для уравнения теплопроводности задается только одно начальное условие U|t=0 = φ(х) (или в другой записи U(x,0) = φ(х)) и физически оно означает, что начальное распределение температуры стержня имеет вид φ(х). Для уравнений теплопроводности на плоскости или в пространстве начальное условие имеет такой же вид, только функция φ будет зависеть, соответственно, от двух или трех переменных.
Граничные условия в случае уравнения теплопроводности имеют такой же вид, как и для волнового уравнения, но физический смысл их уже иной. Условия первого рода (5) означают, что на концах стержня задана температура. Если она не изменяется со временем, то g1(t) ≡ Т1 и g2(t) ≡ Т2, где Т1 и Т2 — постоянные. Если концы поддерживаются все время при нулевой температуре, то Т1= Т2 = 0 и условия будут однородными. Граничные условия второго рода (6) определяют тепловой поток на концах стержня. В частности, если g1(t) = g2(t) = 0, то условия становятся однородными. Физически они означают, что через концы не происходит теплообмен с внешней средой (эти условия еще называют условиями теплоизоляции концов). Наконец, граничные условия третьего рода (7) соответствуют случаю, когда через концы стержня происходит теплообмен с окружающей средой по закону Ньютона (напомним, что при выводе уравнения теплопроводности мы считали боковую поверхность теплоизолированной). Правда, в случае уравнения теплопроводности условия (7) записываются немного по-другому:
Физический закон теплообмена со средой (закон Ньютона) состоит в том, что поток тепла через единицу поверхности в единицу времени пропорционален разности температур тела и окружающей среды. Таким образом, для левого конца стержня он равен Здесь h1 > 0 — коэффициент теплообмена с окружающей средой, g1(t) — температура окружающей среды на левом конце. Знак минус поставлен в формуле по той же причине, что и при выводе уравнения теплопроводности. С другой стороны, в силу теплопроводности материала поток тепла через этот же конец равен Применив закон сохранения количества тепла, получим:
Аналогично получается условие (14) на правом конце стержня, только постоянная λ2 может быть другой, так как, вообще говоря, среды, окружающие левый и правый конец, бывают разные.
Граничные условия (14) являются более общими по сравнению с условиями первого и второго рода. Если предположить, что через какой-либо конец не происходит теплообмена со средой (то есть коэффициент теплообмена равен нулю), то получится условие второго рода. В другом случае предположим, что коэффициент теплообмена, например h1, очень большой.
Перепишем условие (14) при х = 0 в виде и устремим . В результате будем иметь условие первого рода:
Аналогично формулируются граничные условия и для большего числа переменных. Для задачи о распространении тепла в плоской пластине условие означает, что температура на ее краях поддерживается нулевой. Точно так же, условия и внешне очень похожи, но в первом случае оно означает, что рассматривается плоская пластина и края ее теплоизолированы, а во втором случае оно означает, что рассматривается задача о распространении тепла в теле и поверхность его теплоизолирована.
Решение первой начально-краевой задачи для уравнения теплопроводности.
Рассмотрим однородную первую начально-краевую задачу для уравнения теплопроводности:
Найти решение уравнения
удолетворяющее граничным условиям
и начальному условию
Решим эту задачу методом Фурье.
Шаг 1. Будем искать решения уравнения (15) в виде U(x,t) = X(x)T(t).
Найдем частные производные:
Подставим эти производные в уравнение и разделим переменные:
По основной лемме получим
Теперь можно решить каждое из этих обыкновенных дифференциальных уравнений. Обратим внимание на то, что используя граничные условия (16), можно искать не общее решение уравнения б), а частные решения, удолетворяющие соответствующим граничным условиям:
Шаг 2. Решим задачу Штурма-Лиувилля
Эта задача совпадает с задачей Штурма-Лиувилля, рассмотренной в лекции 3. Напомним, что собственные значения и собственные функции этой задачи существуют только при λ>0.
Собственные значения равны
Собственные функции равны (См. решение задачи)
Шаг 3. Подставим собственные значения в уравнение а) и решим его:
Шаг 4. Выпишем частные решения уравнения (15):
В силу линейности и однородности уравнения (15) их линейная комбинация
Шаг 5. Определим коэффициенты An в (19), используя начальное условие (17):
Приходим к тому, что начальная функция φ(x) разлагается в ряд Фурье по собственным функциям задачи Штурма-Лиувилля. По теореме Стеклова такое разложение возможно для функций, удовлетворяющих граничным условиям и имеющих непрерывные производные второго порядка. Коэффициенты Фурье находятся по формулам
Вычислив эти коэффициенты для конкретной начальной функции φ(x) и подставив их значения в формулу (19), мы тем самым получим решение задачи (15), (16), (17).
Замечание. Используя формулу (19), можно также, как в лекции 3, получить решение первой начально-краевой задачи для уравнения Ut = a 2 Uxx. Оно будет иметь вид
где
Учебное пособие: Дифференциальное уравнение теплопроводимости
Название: Дифференциальное уравнение теплопроводимости Раздел: Промышленность, производство Тип: учебное пособие Добавлен 00:27:14 02 февраля 2009 Похожие работы Просмотров: 3441 Комментариев: 22 Оценило: 3 человек Средний балл: 5 Оценка: неизвестно Скачать |