|
x x |
W |
|
||
1 |
|
t |
|
|
t 0
|
|
|
|
x |
|
|
|
|
|
|
x |
|
|
0 |
|
|
0 |
0 |
|
. |
|
|
|
|
|
|||||||
|
|
|
|
l |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Следовательно, первая краевая задача для неоднородного уравнения в общей постановке сведена к изученной ранее задаче вида 83 – 85 о вынужденных колебаниях ограниченной струны с зафиксированными концами.
3.4. Одномерное уравнение теплопроводности
Явление переноса теплоты от более нагретых мест тела к менее нагретым обусловлен изменением температуры u в разных местах тела. Таким образом, рассмотрение данного явления в макроскопической теории, вообще говоря, трансформируется к нахождению нестационарного температурного поля в теле.
Изучим одномерный процесс переноса теплоты теплопроводностью в плоском слое изотропного материала (рис. 8), имея ввиду, что температура u u x, t есть функция
только одного пространственного аргумента x.
Плотность материала, его удельную теплоѐмкость c и коэффициент теплопроводности k, вообще говоря, неоднородной среды примем тоже зависящими лишь от одного пространственного аргумента x.
Описывая математически данное явление, считаем что среда неподвижна, а вариация объѐма материала вследствие вариации температуры является бесконечно малой величиной. Тогда возможно сделать предположение о независимости явления теплопроводности от выполнения механической работы.
В расматриваемом слое материала в качестве некоторой термодинамической системы выделим объѐм V в виде ци-
66
линдра с площадью основания ординатной оси Ox (см. рис. 8).
S
и осью, параллельной ко-
|
|
Рис. 8. Явление теплопереноса |
||||
|
Согласно первомй закону термодинамики, адптирован- |
|||||
ного к рассматриваемому объѐму |
V , |
имеет место соотноше- |
||||
ние |
|
|
|
|
|
|
|
|
dU |
Q |
Q |
, |
99 |
|
|
|
||||
|
|
dt |
1 |
2 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
здесь |
U |
– внутренняя энергия системы, еѐ представляется |
||||
возможным определить, интегрируя объѐмную плотности внутренней энергии x, t по объѐму цилиндра:
|
x2 |
U dV S dx. |
|
V |
x1 |
Тогда вариация внутренней энергии системы за единицу времени
67
dU |
x |
|
|
2 |
|
||
dt |
S |
t |
dx. |
x |
|
||
|
|
|
|
|
1 |
|
|
100
Тепловой поток Q1 сквозь всю замкнутую поверхность
рассматриваемого цилиндра (количество теплоты, выделяемое сквозь данную поверхность за единицу времени), представляется возможным определить интегрированием по поверхности нормальной составляющей плотности теплового
потока |
q. |
нормаль к
Тогда
.
1 |
|
Т |
|
Q |
|
q ndS, |
|
|
|
|
|
здесь
n
– единичная внешняя
Исходя из физического закона Фурье, при теплопереносе теплопроводностью q kgradu. Поскольку в изучаемом
варианте вектор плотности теплового потока q имеет только
одну составляющую
q |
|
k |
u |
, |
|
x |
x |
||||
|
|
|
|||
|
|
|
|
то тепловой поток от рас-
сматриваемого объѐма распространяется только сквозь основания цилиндра, и
|
|
|
u |
|
|
|
|
u |
|
|
x |
|
|
u |
|
|
|
|
|
|
|
|
2 |
|
|
||||||
Q S |
|
k |
|
1 |
S |
|
k |
|
2 |
S |
|
|
k |
|
dx. |
1 |
|
x x |
|
|
x x |
|
|
|
|
||||||
|
|
|
x |
|
|
|
|
x |
|
|
x |
x |
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
101
Внутри рассматриваемого объѐма из-за прохождения эндоили экзотермических реакций, протекания электрического тока, выделения влаги в пористом материале и иных воздействий предсталяется возможным выделение (поглощение) теплота. При условии, что F x, t интерпретируется как объ-
ѐмную плотность или, другими словами, удельная мощность, тепловых источников, за единицу времени в изучаемом объѐ-
ме выделится F 0 или поглотится F 0 количество теплоты
|
|
|
|
x |
|
2 |
|
|
|
2 |
F x, t dx. |
|
FdV S |
|
|||
Q |
|
|
|||
|
|
V |
|
x |
|
|
|
|
|
1 |
|
68
102
После подстановки соотношений
100 |
– |
102 |
в уравнение
99 ,
будем иметь
x |
|
|
|
u |
|
|
2 |
|
|||||
S |
|
t |
|
k |
|
F x, t dx 0. |
x |
|
|
x |
x |
|
|
1 |
|
|
|
|
|
|
103
Поскольку координаты |
x1 |
и |
рутся произвольно, соотношение
x2 оснований цилиндра бе-103 имеет место только в
случае, когда подынтегральная функция нулевая. Следовательно, в рассматриваемом явлении теплопере-
носа локально (в любой точке пространства) обязательно выполнение такого дифференциального выражения:
|
|
|
u |
F x, t . |
t |
k |
|
||
|
x |
x |
|
104
Объѐмная плотность
сжимаемой среды |
|
внутренней энергии
u |
варьируется |
исследуемой в зависимости
неот
|
|
|
|
|
|
d |
~ |
|
|
|
температуры, а |
производная |
du |
c |
интерпретируется как |
||||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
~ |
|
|
|
|
объѐмная |
теплоѐмкость |
c c |
вещества. |
Поэтому |
||||||
|
|
u |
c |
u |
. Исходя из соотношения 104 |
имеет ме- |
||||
t |
u t |
t |
||||||||
|
|
|
|
|
|
|
||||
сто дифференциальное уравнение
c |
u |
|
|
u |
F x, t . |
|
k |
|
|||
|
t |
|
x |
x |
|
105
В случае однородного вещества со свободными от
пературы теплофизическими характеристиками , |
c |
тем- и k
уравнение
105
равносильно
u a2 2u f x, t ,
t x2
106
69
здесь |
a |
2 |
|
k |
– константа, именуемая коэффициентом тем- |
|
|||||
|
c |
||||
|
|
|
|
|
пературопроводности вешества,
Уравнения
f105
x и
, t106
1F x, t .
c
|
есть дифференциальные урав- |
нения в частных производных параболического класса. Они используются в математических интерпретациях явления теплопереноса в неоднородных и однородных телах с одномерным температурным полем. Данные уравнения именуют-
ся уравнениями теплопроводности.
Замечание. Уравнения
105
и
106
моделируют вари-
ацию температурного поля в стержне постоянного поперечного сечения, сделанном из неоднородного или однородного вещества при условии, что боковая поверхность стержня теплоизолирована, площадь поперечного сечения пренебрежимо мала и представляется возможным не принимать во внимание распределение температуры по сечению, предполагая еѐ зависимость лишь от осевой переменной.
Постановка краевых задач.
С целбю использования уравнения теплопроводности для построения математической модели вариации температурного поля в теле, нужно иметь распределение температуры в начальный момент времени, что означает постановку начального условия. Для исследуемого одномерного явления начальное условие
u |
t 0 |
|
|
x |
|
, 0 x l, |
107 |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
определяется исходной функцией x .
Более, необходимо знать тепловой режим на поверхности тела S, что означает постановку граничных условий в каждой точке поверхности тела в произвольный момент вре-
70