Материал: Уравнения математической физики. Пашуева И.М., Ускова Н.Б

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

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

Источник: https://studfile.net/preview/16563675/