Постановка задач ТМО

Имеем объем, на который воздействуют тепловые нагрузки необходимо определить численное значение q V и распределение ее по объему.

Рис.2-Внешние и внутренние источники трения

1. Определить геометрию исследуемого объема в любой выбранной системе координат.

2. Определить физические характеристики исследуемого объема.

3. Определить условия, инициирующие процесс ТМО.

4. Уточнить законы, определяющие перенос тепла в исследуемом объеме.

5. Определить начальное тепловое состояние в исследуемом объеме.

Задачи, решаемые при анализе ТМО:

1.«Прямые» задачи ТМО

Дано: 1,2,3,4,5

Определить: распределение температур в пространстве и во времени (далее 6).

2.«Обратные» задачи ТМО (инверсные):

а) обратные граничные задачи

Дано: 1,2,4,5,6

Определить: 3;

б) обратные коэффициенты задачи

Дано: 1,3,4,5,6

Определить: 2;

в) обратная ретроспективная задача

Дано: 1,2,3,4,6

Определить: 5.

3.«Индуктивные» задачи ТМО

Дано: 1,2,3,5,6

Определить: 4.

ФОРМЫ ПЕРЕНОСА ТЕПЛА И ТЕПЛОВЫЕ ПРОЦЕССЫ

Различают 3 формы переноса тепла:

1) теплопроводность в твердых телах (определяется микрочастицами, а в металлах свободными электронами);

2) конвекция (определяется макрочастицами подвижной среды);

3) тепловое излучение (определяется электромагнитными волнами).

Теплопроводность твердых тел

Общие понятия

Поле температур – это совокупность значений температуры в исследуемом объеме, взятая в некоторый момент времени.

t(x, y, z, τ) - функция, определяющая поле температур.

Различают стационарное и нестационарное поле температур:

стационарное - t(x,y,z);

нестационарное - t(x, y, z, τ) .

Условием стационарности является:

Возьмем некое тело и соединим точки с равными температурами

Рис.3-Градиент температур и тепловой поток

grad t - градиент температуры;

с другой стороны: .

Закон Фурье ‑ тепловой поток в твердых телах пропорционален градиенту температуры, поверхности, через которую он проходит и рассматриваемому интервалу времени.

Коэффициент пропорциональности называется коэффициентом теплопроводности λ , Вт/м·К.

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

;

Для бесконечно малой поверхности и промежутка времени:

Уравнение теплопроводности (уравнение Фурье)

Рассмотрим бесконечно малый объем: dv =dx ·dy ·dz

Рис.4-Тепловое состояние бесконечно малого объёма

Имеем ряд Тейлора:

Аналогично:

; ; .

В общем случае имеем в кубике q V . В основе вывода лежит обобщенный закон сохранения энергии:

.

В соответствии с законом Фурье:

; ; .

После преобразований имеем:

.

Для стационарного процесса:

Пространственная мерность задач определяется количеством направлений, в которых происходит перенос тепла.

Одномерная задача: ;

для стационарного процесса: ;

для :

для : ;

a - коэффициент температуропроводности, .декартова система;

k = 1 , ξ = x - цилиндрическая система;

k = 2 , ξ = x - сферическая система.

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

Условие однозначности – это условия, позволяющее выделить из множества допустимых решений одно-единственное, соответствующее поставленной задаче.

Страница 4

. (2.24)

Уравнение (2.24) называется дифференциальным уравнением теплопроводности (или дифференциальным уравнением Фурье) для трехмерного нестационарного температурного поля при отсутствии внутренних источников теплоты. Оно является основным при изучении вопросов нагревания и охлаждения тел в процессе передачи теплоты теплопроводностью и устанавливает связь между временным и пространственным изменениям температуры в любой точке поля. Отоларингология лазер применение лазеров.

Температуропроводность является физическим параметром вещества и имеет единицу м2/c. В нестационарных тепловых процессах a характеризует скорость изменения температуры.

Из уравнения (2.24) следует, что изменение температуры во времени для любой точки тела пропорционально величине a. Поэтому при одинаковых условиях быстрее увеличивается температура у того тела, которое имеет большую температуропроводность.

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

, (2.25)

где qV - удельная мощность источника, то есть количество выделяемой теплоты в единице объёма вещества в единицу времени.

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

, (2.26)

где r - радиус-вектор в цилиндрической системе координат;

Полярный угол.

2.5 Краевые условия

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

· геометрическая форма и размеры тела,

· физические параметры среды и тела,

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

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

Обычно начальные условия распределения температуры задаются для момента времени t = 0.

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

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

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

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

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

2.6 Теплопроводность через шаровую стенку

С учётом описанной в разделах 2.1 - 2.5 терминологии задачу данной курсовой работы можно сформулировать так. Постоянный тепловой поток направлен через шаровую стенку, причем источником теплоты является внутренняя сфера радиусом R1. Мощность источника P постоянна. Среда между граничными сферами изотропна, поэтому её теплопроводность c является функцией одной переменной - расстояния от центра сфер (радиуса) r. По условию задачи . Вследствие этого температура среды тоже является в данном случае функцией одной переменной - радиуса r: T = T(r), а изотермические поверхности это концентрические сферы. Таким образом искомое температурное поле - стационарное и одномерное, а граничные условия являются условиями первого рода: T(R1) = T1, T(R2) = T2.

Из одномерности температурного поля следует, что плотность теплового потока j так же, как теплопроводность и температура, являются в данном случае функциями одной переменной - радиуса r. Неизвестные функции j(r) и T(r) можно определить одним из двух способов: или решать дифференциальное уравнение Фурье (2.25), или использовать закон Фурье (2.11). В данной работе избран второй способ. Закон Фурье для исследуемого одномерного сферически симметричного температурного поля имеет вид:1 4

z
x
ЛЕКЦИЯ 4
Задачи теплопроводности в различных системах координат.
Декартова система координат
T
T
T
q
i
j
k
T T x, y , z , t
y
x
x
y
T
T T T
c
qV
t x x y y z z
c
T T
qV
t x x
(1)
(2)
(3)
В практике часто встречаются такие условия, которые приводят к необходимости записи уравнения
теплопроводности в иной форме, более удобной для представления решения и его физической
трактовки.
Зависимость вида уравнения
от используемой системы
координат можно исключить,
используя операторную запись
1 T
q
T V
a t
2
x
2
2
y
2
2
z 2
a c
T
c
div gradT qV
t
или
c
T
T qV
t
(4)
Слагаемые, выражающие тепловыделение и аккумулирование энергии, инвариантны относительно
системы координат (т.е., неизменны); но слагаемые, выражающие результирующий кондуктивный
тепловой поток, зависят от геометрии и, следовательно, от системы координат.

Цилиндрическая система координат
z
c
dr
r
dz
r , z
z
x
T
div q q
t
q T
x r cos
y
r , z
(5)
y r sin
(6)
1 1 2
2
r 2 2 2
r r r r
z
d
y
dr
d
dy
dx
z
qr
(7)
1 T 1 T 1 2T 2T qV
r 2 2 2
a t r r r r
z
x
1 T 1 T
r
qV
a t r r r
T
1 T
T
; q
; qz
r
r
z
a
(9)
T Ts
c
(8)

r ,
Сферическая система координат
z
dr
r ,
r
d
x
1 T
div q q
a t
q T
y
1 2
1
1
2
2 r
2
sin
2
r sin 2
r r r r sin
T
1 T
1 T
; q
; q
r
r
r sin
(10)
1 T 1 2 T
1
T
1
2T qV
2 r
2
sin 2
2
a t r r r r sin
r sin
(11)
d
qr
1 T 1 2 T qV
2 r
a t r r r
x r sin cos
y r sin sin
z
(12)
z r cos
y
x

Уравнения теплопроводности для тел канонической формы
Запись уравнений в различных системах координат особенно удобна,
когда нужно найти распределение температуры в телах канонической
формы – в цилиндре или шаре. В этих случаях уравнения существенно
упрощаются при задании особых условий, когда поле температуры
зависит только от одной координаты.
параллелепипед
пластина
цилиндр
сфера
c
T T T T
qV
t x x y y z z
1 T 2T qV
2
a t x
qe
1 T 1 T qV
r
a t r r r
1 T 1 2 T qV
r
2
a t r r r
T Ts
z
y
x

1 T 1 n T qV
r
n
a t r r r
Три последних
уравнения вместе:
n 0
n 2
n 1 цилиндр
плоскость
T T0
T* T0
t
t*
(13)
сфера
r
r*
1 1 n
qV
n
Fo
На доске
Число Фурье
at*
Fo 2
r*
qV 1:
at*
at
1: 2
2
r*
r*
(14)
qV r*2
qV
T* T0
q
T* T0 V r*2
1 n
1
n
Fo

Стационарные задачи теплопроводности в различных системах координат
Цилиндрическая стенка: стационарный процесс теплопроводности в
цилиндрической стенке (трубе) с внутренним радиусом r1;
d 1 2r1
r1
1 T 1 T 1 2T 2T qV
r
a t r r r r 2 2 z 2
r2
Te1
2
1
T1
d1
T2
Te 2
dT
u
dr
du 1
u 0
dr r
T C1 ln r C2
q
d2
(17)
dT
C
1 (18)
dr
r
d 2T
1 dT
0
2 r dr
dr
(15)
ln u ln r ln C1
(16)
Удельный тепловой поток не
постоянен по толщине и убывает по
направлению к внешней поверхности
В стационарных условиях постоянным должен быть полный тепловой поток проходящий через
участок цилиндрической трубы длиной l и равный
Q q F q 2 rl
Удельный тепловой поток
убывает с радиусом
!!!
(19)
Площадь поверхности
увеличивается с радиусом
Температура по толщине трубы изменяется нелинейно даже при постоянном
коэффициенте теплопроводности
Постоянные интегрирования могут быть найдены из граничных условий.

r r1: T T1; r r2: T T2
T1 C1 ln r1 C2 ,
Система линейных
уравнений
T2 C1 ln r2 C2 ,
T ln r2 r T2 ln r r1
T 1
;
ln r2 r1
q
Q
Погонный тепловой поток
qп
(20)
dT
C
1
dr
r
dT
T
l 2 r
2 l ,
dr
ln r2 r1
Вт
Q
2
T , T T1 T2
l ln r2 r1
(21)
(22)

(температуры стенок – неизвестны)
T C1 ln r C2
Можем поступить аналогично:
r r1:
Поступим иначе:
(23)
T
T
1e T Te1 ; r r2:
2e Te2 T
r
r
Конвективный тепловой поток на единицу длины
трубы должен равняться погонному тепловому потоку
вследствие теплопроводности:
qп 1e Te1 T1 2 r1
2
T1 T2
qп
ln r2 r1
qп Kc Te1 Te2
1
Kc
, Вт/(М·К)
1
1 r
1
ln 2
2 1e r1 2 r1 2 2e r2
qп 2e T2 Te2 2 r2
Коэффициент теплопередачи для
цилиндрической стенки
Rc
1
1
1 r
1
ln 2
Kc 2 1er1 2 r1 2 2er2
плоская стенка
R
1 L 1
1 2
1 L 1
K
1
2
1
Вт/(М2·К)
Из системы уравнений (23) мы можем найти
и температуры стенок и подставить в (20)
Полное термическое
сопротивление трубы
(24)
(25)
(26)
Размерность
отличается от
размерности К для
плоской стенки!
T ln r2 r T2 ln r r1
T 1
;
ln r2 r1
Можно
на доске

В безразмерных переменных
r1
d 2
d
r2
2
1 d
0
d
(27)
d
Bi
d
(28)
r1 r2:
Te1
2
1
d1
d2
Задание
на дом:
1:
T Te 2
r
; r* r2
Te1 Te 2
r2
d
Bi 1
d
(29)
2er2 1e
Bi
2e
C1 ln C2
Te 2
C1
Bi C1 ln C2
C1 Bi C2 1
(30)
А) Перейти аккуратно к безразмерным переменным
Б) Найти постоянные интегрирования из системы (30)
В) Построить для разных значений параметров

10.

Принципы
последовательного
и
параллельного
соединений термических сопротивлений в цепь,
справедливые для плоской стенки в прямоугольной
системе координат, можно применить и для задачи о
теплопроводности в полом цилиндре.
Электрическая аналогия
2
Q
1
Q
T0
r3
r2
r1
T1
T2
Ts
RT
ln r2 r1
2 l
Жидкость течет в трубе, R 1 1
0
F 2 r1l
покрытой изоляционным
материалом
dT
T
l 2 r
2 l,
dr
ln r2 r1
T
Q
,
ln r2 r1 2 l
В форме
закона Ома
Термическое сопротивление
полого цилиндра
Конвективное термическое
сопротивление жидкости
Имеем последовательное соединение конвективного сопротивления жидкости с двумя
кондуктивными термическими сопротивлениями. Если задана температура жидкости и температура
внешней поверхности:
T0 Ts
T
Q
А)
R
full
r
r
1
1
1
ln 2
ln 3
2 1r1l 2 l 1 r1 2 l 2 r2
(31)
Сопротивление
изоляции
Если заданы температуры внутренней и внешней поверхностей
Б)
T
Q
R full
T1 Ts
r
r
1
1
ln 2
ln 3
2 l 1 r1 2 l 2 r2
(32)

11.

Пример
1 185
В алюминиевой трубе, имеющей теплопроводность
Вт/(м К), течет водяной пар
◦
при температуре 110 С. Внутренний диаметр трубы – 10 см, наружный диаметр – 12
Te
см. Труба расположена в помещении с температурой
30◦С; коэффициент
e
конвективной теплоотдачи от трубы
к воздуху
равен 15 Вт/(м2К). 1) Требуется
найти тепловой поток на единицу длины трубы, если труба не теплоизолирована.
2)Чтобы снизить тепловые потери от трубы, она была покрыта слоем теплоизоляции
(2 0 ,2 Вт/(м К)) толщиной 5 см. Найти тепловой поток на единицу длины от
тепоизолированной трубы. Предположить, что конвективное термическое
сопротивление пара пренебрежимо мало.
Решение. Для трубы без теплоизоляции наиболее существенными являются
кондуктивное термическое сопротивление самой трубы и конвективное термическое
сопротивление комнатного воздуха. Поскольку конвективным термическим
сопротивлением пара можно пренебречь, температура внутренней поверхности
трубы равна температуре пара. Тепловой поток на единицу длины трубы следует из
соотношения T T
110 30
80
q
0
e
ln r2 r1
1
2 1
2 r2 e
ln 6 5
1
2 185 2 0 ,06 15
1,57 10
4
0 ,177
452 Вт/м.
Для трубы с теплоизоляцией нужно добавить термическое сопротивление
теплоизоляции, и соотношение для теплового потока примет вид
q
T0 Te
80
138
ln r3 r2 1,57 10 4 0 ,096 0 ,482
ln r2 r1
1
2 1
2 r3 e
2 2
Вт/м.

12.

Многослойная цилиндрическая стенка
qc
Tn T1 1
n
d
1
ln i 1
2 i
di
, d i 2r1
qc
i 1
Остается справедливым понятие
эквивалентного коэффициента
теплопроводности
экв
ln d n 1 d1
n
i 1
T1
T2
1
(33)
T3
2
(34)
1 d i 1
ln
i di
r1 d1 2
... ...
Tn 1
n 1
Tn
n
Tn 1
r2 d2 2
Температура Ti 1
Ti 1 Ti
2 экв T1 Tn 1
ln d n 1 d1
на границе между i-м и i+1-слоями
qc 1 d 2 1 d3
1 d
ln ln ... ln i 1
2 1 d1 2 d 2
i
di
(35)
Коэффициент теплопередачи:
Kc
1
1
1d1
n
i 1
1 di 1
1
ln
2 i di 2 d 2
(36)

13.

r1
Радиальный тепловой поток в трубе обратно пропорционален логарифму
наружного радиуса (возрастает сопротивление радиальной проводимости);
r2
Рассеяние тепла от наружной поверхности прямо пропорционально этому
радиусу (увеличивается площадь охлаждающей поверхности)
qc K c Te1 Te 2
Kc
1
,
1
1 r2
1
ln
2 1r1 2 r1 2 2 r2
Следовательно, существует определенный радиус, при
котором потери тепла максимальны!
Если при фиксированном (малом) внутреннем радиусе увеличивать
толщину стенки трубы (т.е., увеличивать внешний радиус r2), то действие
логарифма в формуле для термического сопротивления окажется более
сильным, чем при большем внутреннем радиусе

14.

Критический диаметр теплоизоляции
qc Kc Te1 Te2
Kc
1
,
1
1 r2
1
ln
2 1r1 2 r1 2 2 r2
dqc
0
dr2
Условие экстремума:
дает
r2 * 1
2
Критический радиус
Частный случай нулевого внутреннего сопротивления, 1 1 0
y
q
2 Te1 Te 2
1
r
, x 2,
ln x x
r1
2r1
(38)
0 Внешнее сопротивление также равно нулю
r1 r2
Толщина стенки равна 0
1: x 2r2
Для заданного внутреннего радиуса величина критического
внешнего радиуса увеличивается, если увеличивается
теплопроводность трубы или если понижается коэффициент
теплоотдачи на внешней поверхности
(37)
Bi 1

15.

изоляция
Существование критического наружного радиуса приводит к тому, что при
некоторых реальных условиях, вопреки привычным представлениям,
потеря тепла изолированной трубой фактически может быть снижена
путем уменьшения толщины изоляции
d1
d2
Полное термическое сопротивление для двухслойной трубы, сечение которой
изображено на рисунке, определяется по формуле
d3
Rc
1 2
труба
Условие
экстремума:
d2 d3 *
d3 d 2
(39)
- толщина изоляции
Термическое сопротивление теплопроводности изоляции (I) растет с увеличением
толщины изоляционного покрытия; термическое сопротивление теплоотдачи изоляции
(II) – падает (так как увеличивается поверхность теплоотдачи)
dRc
1
1
0
dd3 2 2 d3 2 d 32
Rc
d2 d3 *
1
1
1 d2
1 d3
1
ln
ln
K c 1d1 2 1 d1 2 2 d 2 2 d3
II
(I)
d 3 *
22
8 32
0
d3 * 2 2
2
не зависит от
d2
(40)
(т.е., не зависит от диаметра самого трубопровода)
В критической точке полное термическое
сопротивление – минимально!
увеличение толщины изоляции уменьшает теплоотдачу
нанесение выбранного покрытия первоначально приведет к возрастанию
теплоотдачи, и лишь при достижении критического диаметра тепловой поток будет
убывать; затем достигнет той величины, которая была без изоляции и лишь потом
приведет к желаемому эффекту

16.

Задача для полого шара
(шаровая стенка)
d 2T
dr
2
2 dT
0
r dr
(41)
Рассматриваем пространственно одномерную стационарную
задачу теплопроводности в шаровой стенке с заданными
радиусами внутренней и внешней поверхностей. Одномерность
задачи означает, что распределение температуры в стенке
зависит только от радиуса
С помощью замены
переменных
r1
dT
u
dr
du
2u
Общее решение
dr
r
C
C
dT C1
ln u 2 ln r ln C1; u 21 ; T r 1 C2 ;
2
r
dr r
r
r2
Граничные условия первого рода
r r1: T T1
C1
C2
r1
T 1 r 1 r2 T2 1 r1 1 r
T r 1
1 r1 1 r2
r r2: T T2
(42)
Плотность потока тепла
Полный тепловой поток
Q
T1
T2
C1
C2
r2
(43)
(44)
dT
r2
T1 T2
q
2 C1
dr
1 r1 1 r2
r
(45)
dT
4
T1 T2
4 r 2 4 C1
dr
1 r1 1 r2
(46)

17.

Граничные условия третьего рода
T r
Общее решение
не изменяется
C1
C2
r
T
r r1: -
1 T Te1
r
T
r r2: -
2 Te2 T
r
(47)
2r2 C1 2r22C2 2r22Te2
C1
1r1
1r12
2 r22
2 r2
r1
r2
1r1 C1 1r12C2 1r12Te1
1r12 Te 2 Te1
dT C1
2
dr r
C2
(48)
Полный тепловой поток Q не
зависит от текущего радиуса
1r1 T 1r12 T
2 r2 e 2 2 r22 e1
1r1 1r12
2 r2 2 r22
(49)
В пределе при идеальном теплообмене сред с заданными температурами и
шаровой стенки (т.е., при бесконечных коэффициентах теплоотдачи) решение задачи с
граничными условиями третьего рода переходит в решение задачи с граничными
условиями первого рода.
4
Q
T T
1 1 1 2
r1 r 2
=
тепловому потоку,
4 r1 2 1 Te1 T
приходящему к
внутренней стенке
=
тепловому потоку,
4 r 2 2 2 T Te 2
покидающему
внешнюю стенку

18.

Распределение температуры в шаровой стенке
для граничных условий третьего рода
Дома:
воспроизвести все
решение
1 1
1 1
T1 T 2
r r
r1 r
2
T r
1 1
r1 r 2
Температуры стенок:
T1
r12 1Te1 s Te 2
2 Te1
r2 2
r12 1
s 1 2 r12 1
r
2 2
r12 1
r12 1
s Te 2 2 Te1
r2 2
2
r1 1 2
s 1 2 r1 1
r
2 2
r12 1Te 2
T2
Проводимость шаровой стенки:
s
1 1
r1 r 2
r1r 2
r 2 r1

19.

Решения простейших задач в безразмерной форме
Соберем решения стационарных задач для тел канонической формы с
граничными условиями первого рода вместе
T p T1 T1 T 2
r
r2
Дома: воспроизвести!
Tc
1 1
1 1
T1 T 2
r r
r1 r
2
Ts
1 1
r1 r 2
T1 ln r 2 r T 2 ln r r1
l n r 2 r1
T T2
T1 T 2
r
r2
0,8
p 1
ln
ln
1 1
1
1
1 1
c
p
0 1
0,6
r2
1
r1
2
0,2
0,0
0,0
В плоской стенке качественное распределение
температуры (линейное) не зависит от ее
толщины. А вот в цилиндрической и шаровой –
нелинейно меняется с радиусом;
характер
распределения (кривизна кривой) зависит от
соотношения внешнего и внутреннего радиусов.
1
3
0,4
0,2
0,4
0,6
0,8
Распределение температуры в плоской
(1), цилиндрической (2) и шаровой (3)
стенке. Сплошные линии
;
10
пунктирные линии - . 5

20.

В случае граничных условий третьего рода решения простейших задач
зависят от параметров, характеризующих теплообмен.
Для одинаковых коэффициентов теплоотдачи.
T Te 2
Te1 Te 2
r
r2
1 2
0,8
для пластины
1
p 1 1 2
1 1
2 Bi
2
1
2 Bi
для цилиндра:
0,6
3
0,4
3
1
2
0,2
1 2 ln 2 ln
ln
1 1
2
1 Bi ln
1 Bi ln
c
для сферы:
s
1
1 1 1 2
1
1 Bi 1
1 1 Bi
2
Bi
r1
1
1 1 Bi
0,0
0,2
0,4
0,6
1
0,8
2
Распределение температуры
вдоль координаты в плоской (1),
цилиндрической (2) и сферической
(3) стенках в условиях
конвективного теплообмена.
Сплошные линии - Bi 2 ;
пунктирные - Bi 1 0

21.

Примеры: сосуд Дьюара (Dewar bottle)
Частица металла, покрытая пленкой окисла
Задание на дом:
1.Сформулировать задачу о распределении температуры в двухслойной
шаровой оболочке при ее конвективном охлаждении, пользуясь материалом
лекции. Тепловой контакт между слоями считать идеальным. Привести
задачу к безразмерной форме. Построить точное аналитическое решение
этой задачи.
2.*Рассчитать температуры внутренней и внешней поверхностей шаровой
оболочки в задаче 1, а также температуру на контакте; определить полный
тепловой поток, уходящий с поверхности шара, принимая, что температуры
среды внутри оболочки – 175 С, температура окружающей среды – 25 С;
коэффициенты теплоотдачи одинаковы и равны – 28,8 ккал/(м2·час·град);
внутренний, и внешний радиусы оболочки – 3 см и 5 см, толщина
внутренней оболочки – 25 мм. Внутренняя оболочка изготовлена из
материала с теплопроводностью 1,45 ккал/(м час град); внешняя из
материала с коэффициентом теплопроводности 0,137 ккал/(м·час·град). Как
будет изменяться тепловой поток при изменении толщины внешней
оболочки в пределах от 25 мм до 300 мм?

22.

d 2T
Te 2
2
T1
Te1
T2
1
xmax
qV
0;
2
dx
Г.у. первого рода: r r1:
qV const
T T1;
(1)
r r2:
T T2 (2)
Г.у. третьего рода:
r r1:
-
T
1 T Te1 ;
r
r r2:
-
T
2 Te2 T
r
Первый «способ» решения:
Решается задача элементарным интегрированием:
qV x 2
T x
C1x C2
2
dT
q
V x C1;
dx
(4)
Подставляя общее решение в г.у., найдем постоянные интегрирования.
Максимум находится на некотором расстоянии от поверхностей.
Положение максимума можно найти из условия (условие экстремума)
dT
q x
V C1 0
dx
(5)
dT
0
dx
(3)

23.

Задачи с внутренними источниками тепла
ТЕПЛОПРОВОДЯЩАЯ ПЛОСКАЯ СТЕНКА С ОБЪЕМНЫМ ТЕПЛОВЫДЕЛЕНИЕМ
Te 2
2
T1
Te1
1
2
1
Поступим немного иначе. (Второй «способ»
решения)
qV x 2
T x
C1x C2
общее
решение
2
(4)
Поместим начало координат в точку, где
температура максимальна
T2
1; 2
- расстояния от максимума до краев пластины
0
C1 0
Граничное условие справа перепишем следующим образом:
x 2:
dT
dx
2
2 T Te 2
2
2
q
V
2
2 C2
Te 2 qV 2
2
(6)
Так как плоскость x=0 можно считать теплоизолированной, все тепло, выделившееся в
пластине справа в единицу времени, должно быть отведено в окружающую среду
посредством теплоотдачи с правой стенки. В противном случае будет нарушено условие
стационарности
qV 2 - количество тепла, выделяющееся в объеме пластины толщиной =1 в единицу времени
Слева – выражение для потока теплоотдачи с единицы площади поверхности пластины

24.

Аналогичные рассуждения для левого слоя пластины толщиной
1 2
приводят к выражению
2
q
V
2
1 C2
Te1 qV 2
2
(7)
С помощью равенств (6), (7) находим положение
максимума
2
2 1 2 Te1 Te 2 qV 2 1 2
2qV 1 2 1 2
(8)
Определяя постоянную С2, (подходит любое из равенств) , находим общее решение.
Наиболее простой вид оно принимает, если
1 2 ;Te1 Te2 Te
1 2 2
тогда
qV qV 2
C2
Te
2
8
и
2
q
qV
2
T x
x V Te
2 2
2
(9)
(10)
qV 2 qV
Тем ниже, чем выше теплопроводность пластины
Tmax T x 0
Te
8
2
q
Температура стенки Ts T1 T2 V Te растет с ухудшением теплоотдачи
2

25.

Граничные условия первого рода
T1
2
1
T2
0
qV 22
C2 T2
2
dT
dx
2 T1 T2
2 1
2
qV 2
(11)
qV 2 2
C2 T1
2
2
qV 2 T1 T2
2
T x T2
x
1
2
2 2
qV
При очень больших значениях
x 2:
qV x 2
T x
C1x C , C1 0 (4)
2
2
Граничные условия третьего рода переходят в граничные
условия первого рода. Поэтому такое же решение мы
получим, используя предыдущее решение
2 T Te 2
2
(12)
T x T2 T2e
2
(13)
Следовательно, из симметричной задачи с граничными условиями третьего рода (10) находим
2
qV
2
T x
x Ts
2 2
Tmax T x 0
q
V Ts
8
2
Температура
стенок
(14)
Это же равенство следует из предыдущего решения при условии равенства температур стенок

26.

Рассмотрим бесконечный сплошной цилиндр, равномерно нагреваемый (или
охлаждаемый) с боковой поверхности. В объеме цилиндра находится источник тепла
постоянной интенсивности. Требуется найти распределение температуры для
установившегося режима.
d 2T 1 dT q
dr
u dT dr
2
r dr
q r
du
r
u V 0
dr
V
или
0
(1)
d ru qV r
0
dr
qV r 2
ru
C1
2
q r C
dT
V 1
dr
2
r
Общее решение
Первый
интеграл
(3)
qV r 2
T
C1 ln r C2
4
Условие в центре для
сплошного цилиндра
dT dr 0; r 0
(2)
(4)
C1 0

27.

Цилиндр с объемным тепловыделением
dT
T Te
r R
dr
qV 2
qV R
2
qV R qV R 2
T
R
r
Te
C2
Te
4
2
2
4
q
q R
q R
Tmax V R 2 V Te
Ts V Te
4
2
2
Внешнее условие:
плотность теплового потока на поверхности цилиндра:
полный тепловой поток с поверхности цилиндра:
q Ts Te
Q qF
(5)
(6)
(7)
qV R
2
qV R
2 Rl qV R 2l
2
Задача об охлаждении цилиндра с объемным тепловыделением представляет, в
частности, интерес для нахождения распределения температуры в катодах,
используемых в плазмотронах для генерации потоков ионов. В практическом
приложении эта задача может быть переформулирована так: найти мощность
источника, достаточную для распыления катода при условии, что для этого нужно
достичь температуру плавления материала катода
Используя общее решение (4), можно найти распределение температуры по толщине
стенки полого цилиндра или по толщине цилиндра, покрытого защитным слоем
(рассмотрим далее). В первом случае нужно задать условия на внутренней поверхности
цилиндра. Во втором случае потребуется дополнительное условие на границе раздела
двух материалов с разными свойствами, т.е. граничное условие четвертого рода.

28.

Шар с объемным тепловыделением
qV r 2 C1
Дома: покажите,
T
C2 (2)
(1)
что общее решение
6
r1
dr 2
(1) имеет вид (2)
dT
Условия:
dT dr 0; r 0 и dr T Te ; r R
q
q
дают С1 0 и
C2 Te V R V R 2
3
6
2
qV
qV 2 r (3)
T Te
R
R 1
3
6
R
q
q
Tmax Te V R V R 2 (4)
максимальная температура
3
6
q
q
Температура поверхности
Ts Te V R V R 2 (5)
3
6
R 2 dT
1
Полный поток тепла через поверхность
Q
R 3qV
4 dr r R 3
шара
qV R
qV 2 qV R
T
Te
Tmax
R
Te
цилиндр
s
2
4
2
Сравните
d 2T
2 dT qV
0
r dr
Плоский слой Tmax
qV qV 2
Te
2
8
q
T s V Te
2
с (4), (5)

29.

Пример 1. Найти максимальную силу тока, который можно пропускать по
алюминиевой проволоке (λ=204 Вт/(м·К)) диаметром 1 мм, чтобы ее
температура не превышала 200 С. Проволока подвешена в воздухе с
температурой 25 С. Коэффициент конвективной теплоотдачи от проволоки к
воздуху равен 10 Вт/(м2·К). Электрическое сопротивление Re/l на единицу
длины проволоки есть 0,037 Ом/м.
Решение. Воспользуемся формулой (66), из которой следует
qV
Re I 2
R 2l
Tm ax
qV R R
I 2 Re
Te
1
Te
2
2
2 R l
R
1 2
Подставляем заданные значения физических величин:
200 25
I
2
2 1 0 3
Отсюда находим силу тока:
1 0 3 2 1 0
0,0 3 7 1
2 204
2 10
I 12,2 A

30.

Провод с изоляцией
Строгая математическая постановка задачи:
d 2T1
dr
2
d 2T2
Первое условие есть условие симметрии;
второе говорит о том, что тепловой
контакт между проводом и изоляцией –
идеальный, а третье соответствует
конвективному теплообмену провода с
изоляцией с окружающей средой.
dr
2
1 dT2
0
r dr
r 0: dT dr 0
r R: 1
r R
(1)
R r R
(2)
(3)
dT1
dT
2 2 ; T1 T 2
dr
dr
r R: 2
Общее решение задачи:
1 dT1 qV
0
r dr
1
dT2
T2 Te
dr
qV r 2
T1
C1 l n r C 2
4 1
T2 C3 l n r C 4
(4)
(5)
Дома: покажите
справедливость

31.

Провод с изоляцией
qV r 2
T1
C1 l n r C 2
4 1
Общее решение задачи:
T2 C3 l n r C 4
Из условия (3) имеем:
C1 0
q R
C
1 V 2 3
R
2 1
Условия (4) дают:
qV R 2
C3
2 2
qV R 2
qV R 2
C2
l nR C 4
4 1
2 2
Из условия (5) следует:
qV R 2
C3 2 qV R 2
2
ln R C 4 Te
R
R 2 2
2 2
Находим:
qV R 2
q R
C 4 Te
l n R V
2 2
2
qV R 2 2 1 qV R 2 R
C 2 Te
ln
1
4 1 R 2 2
R

32.

Следовательно, распределение температуры в проводе с изоляцией
описывается формулами
qV R 2 2 1 qV R 2 R qV r 2
T1 Te
ln
1
4 1 R 2 2
R 4 1
и
qV R 2 2 qV R 2 R
T 2 Te
ln
2 2 R
2 2
r
Окончательное решение представим в виде:
T Te
i i
T Te
qV R 2
T Te
1
r
R
1
Bi K
2
1 1 2
ln 1
4
K 2
4
2
K K 1
ln
2Bi
2
Определим поток тепла с поверхности
проводника
q T2 R Te
Q R2l T2 R2 Te
K Bi 1
K Bi 1
Перейдите дома к
безразмерным переменным
0 1
Bi
1 1
K
Q
R2 2 l T* Te
1
2
R
2
K
Bi
- изоляция не отводит тепло от проводника с током
- возможно остывание проводника за счет потерь тепла в
окружающую среду
R

33.

Пример 2. Пусть по длинной алюминиевой проволоке диаметром 1 см
течет электрический ток силой тока 1000 А. Проволока покрыта слоем
резиновой изоляции толщиной 3 мм (λ2=0,15 Вт/(м·К)). Температура
наружной поверхности изоляции 30 С. Найти температуру внутренней
поверхности изоляции. Омическое сопротивление проволоки на единицу
длины 3,7·10-4 Ом/м.
Решение. Для решения этой задачи воспользуемся второй формулой для Т2
рассмотренной сопряженной задачи. С учетом того, что задана температура
2
внешней поверхности изоляции, т.е.
Re I 2
Re I 2
R
T2 r R Te
ln
qV
2
l
2
R
R l
2
2
1000
0 . 005 0 . 003
273 30 3 . 7 10 4
ln
477 . 6
2 3 . 14 0 . 15
0 . 005
Используя значение коэффициента теплопроводности алюминиевой проволоки
1 232 Вт/(м·К) и формулу для Т, можем рассчитать температуру в центре
1
провода. В рассматриваемых условиях имеем
2
Re I 2
Re I 2
R Re I
T1 r R Te
ln
T2 r R
l 2 2 R l 4 1
l 4 1
3 . 7 10 4 1000
477 . 6
477 . 7
4 3 . 14 232
2

34.

Задание на дом.
1.Ток силой I=200А пропускается через проволоку из нержавеющей стали
диаметром 2 мм и длиной 1 м. Электрическое сопротивление проволоки –
0.125 Ом, коэффициент теплопроводности 17Вт/(м·К). Температура
поверхности проволоки 150 С. Требуется рассчитать температуру на оси
проволоки.
2.Предположить в этой же задаче, что проволока покрыта слоем изоляции
(коэффициент теплопроводности изоляции 0,15 Вт/(м·К)), а коэффициент
теплоотдачи на поверхности изоляции равен 60 Вт/(м2К). Как нужно
изменить силу тока (увеличить или уменьшить), чтобы температура
поверхности проволоки осталась равной 150 С.

35.

Эффективные (эквивалентные) теплофизические свойства
Реально используемые в машиностроении и окружающие нас материалы
являются многокомпонентными и многофазными. Это относится к сталям,
сплавам, интерметаллидным композитам, спеченным материалам,
волокнистым композитам, композитам на полимерной основе, смесям,
растворам и т.д.
Если для исходных компонентов (из которых композиты синтезируют в
разных технологиях) или дано используемых материалов со свойствам все
более или менее ясно, то для вновь разрабатываемых материалов
определение свойств представляет собой серьезную проблему.
Стандартные экспериментальные методы могут не работать или становятся
дорогими или трудоемкими
Для расчета необходимо знать свойства составляющих, структуру и взаимное
влияние физических явлений друг на друга.
Без данных о физических свойствах невозможен ни один научный
или инженерный расчет
Дульнев Г.Н., Заринчак Ю.П. Теплопроводность смесей и композиционных
материалов

36.

Модели для расчета свойств:
корпускулярные (молекулярные), континуальные и комбинированные
В корпускулярных моделях изучают свойства на основе знаний о природе,
строении и характере взаимодействия частиц. Расчет физических свойств в
этом случае возможен лишь с использованием данных о других свойствах.
Классификация гетерогенных структур:
Дульнев, стр.10-52 (открыть)
Композиты: стр.106-130

37.

Известны многочисленные способы расчета эффективных коэффициентов
теплопроводности гетерогенных и пористых материалов
В простейшем приближении для процесса теплопроводности в отдельной
микрообласти (которую рассматривают как представительный объем)
справедливы физические уравнения
JT ,k k grad Tk , div JT ,k 0
Граничные условия на поверхностях раздела областей с идеальным
тепловым контактом имеют вид:
T
T
k k k 1 k 1 ; Tk Tk 1
n
n
Для определения эффективной теплопроводности материала (состоящего из
различных фаз) необходимо определить распределения физических полей во
всех микрообластях, а потом уже перейти к квазигомогенной среде, для
которой справедливы соотношения
JT * T
1
J k dV ;
V
1
Tk d
T
V
V
Установление вида этой
Эффективный коэффициент: f k , k ;
зависимости и является
основной задачей
- доли фаз
различных теории.
JT
T

38.

Двухфазная система
1
J
J1dV1 J 2dV2 1 1 T1 2 2 T2
V
V2
V1
1 V1 V , 2 V2 V
(1)
1 1 1 2 2 2 ;
k
T 1 T1
2 T2
Tk T
T
2
1 1 2 2 1
Следует из
предыдущего
, k 1,2
- средний по объему градиент
Система двух уравнений (1) содержит три неизвестных. Для е замыкания
требуется дополнительная информация, например, сведения о структуре
гетерогенной системы, данные специально поставленного эксперимента.
Решение проблемы замыкания таких систем и привело к появлению всего
разнообразия методов определения коэффициентов переноса (не только
коэффициента теплопроводности), которое известно в литературе

39.

1. В случае простейшей структуры, представляющей собой систему
неограниченных пластин, параллельных потоку J
1 2 1
и
1 1 2 2
2. Если слои - перпендикулярны потоку
1 T1 2 T2 ;
1 2 2 1
1 2
1 2
1
Типы структур неоднородных сред весьма разнообразны. Так, в случае
двухфазных сред, к которых фазы (микрообласти, содержащие разные фазы)
могут быть распределены в пространстве как хаотически, так и упорядоченно,
можно выделить структуры, содержащие одну из фаз в виде изолированных
изомерных (1) или анизотропно ориентированных (2) включений в
непрерывной другой фазе, зернистые системы с непрерывным каркасом (3) и
порами (4), волокнистые системы из волокон (5) и пор (6), статистически
неоднородные (микронеоднородные) системы из близких по размерам
компонентов (7), слоистые системы из параллельных (8) и перпендикулярных
(9) потоку слоев. Можно представит себе системы, состоящие из отдельных
подсистем с различными структурами описанного типа. Дополнительно
каждая из фаз, входящих в структуры может быть как многокомпонентной, так
и однокомпонентной. В любом случае требуется расчет свойств каждой из фаз
или их экспериментальное определение.

40.

Уравнение Кондорского
3 1 1 3 2 1 2
3 1 1 3 2 1
Оделевского (метод
1
эффективной среды)
4
16
2
2 1
1 V1 V , 2 V2 V
13
2 1
1 2
Интегральный метод
Двусторонние оценки (оценки
Хашина-Штрихмана)
Шермергор:
1 2
1
2
1
1
2 1
1
1
1 3
1 3
1 2 1 2
1 2 1 2
1 2
1 2
1 2 1
1 1 2
Индекс 1 относится к матрице, а «2» - к включениям
Несмотря на упрощенные модели сред, некоторые из известных формул
позволяют проводить вполне достоверные оценки, хотя число формул для
различных частных случаев сред быстро возрастает с увеличением числа фаз.

41.

Дома:
Имеется композит. Матрица - сплав на основе фольфрама (считаем его
коэффициент теплопроводности равным теплопроводности вольфрама).
Частицы (включения) карбид титана.
Используя выписанные выше формулы рассчитать зависимости
эффективных коэффициентов теплопроводности композита от доли
включений (ξ= от 0 до 0,75). Построить на одном графике.
Какой вывод можно сделать?

42.

Свойства зернистых и пористых материалов
На эффективную теплопроводность пористых материалов при прочих равных
условиях оказывает влияние теплопроводность твердой фазы. При этом для
одних пористых материалов (на основе А12О3, BeO, MgO и др.) коэффициент
теплопроводности с ростом температуры уменьшается, в то время, как для
других, изготовленных на основе SiO2, ZrO2, - увеличивается. Решающее
влияние на эффективную теплопроводность оказывает пористость, поскольку
сами поры вследствие низкой проводимости газа являются эффективным
барьером на пути распространения тепла. Однако здесь имеются иные
механизмы теплопереноса (конвекция, излучение).
Самые простые модели основаны на представлении пористого или
дисперсного материала в виде плоских чередующихся слоев, составленных и
твердого каркаса (остова) и воздуха.
1
1
2
2
1
1 1 2
- доля пор; пористость
- теплопроводность воздуха или другого вещества, заполняющего
пористое пространство

43.

Модели, представленные на рис в центре связывают с именами
Maxwell–Eucken (Максвелла-Эйкена). Результат имеет вид
1
2
2 1 2 2 1 2
2 1 2 2 1 2
2 2 1 2 2 1 1
2 2 1 2 2 1 1
1 1
2
0
1 2
2 2
непрерывным является твердый каркас
непрерывным является пористое
пространство
модель теории эффективной среды

Распространение теплоты теплопроводностью в плоской и цилиндрической стенках при стационарном режиме (граничные условия первого рода)

Однородная однослойная плоская стенка. Рассмотрим распространение теплоты теплопроводностью в однородной однослойной плоской стенке толщиной 8 при ее неограниченной ширине и длине.

Ось х направим перпендикулярно стенке (рис. 7.4). По обеим поверхностям стенки как в направлении оси у, так и в направлении оси г благодаря равномерному подводу и отводу теплоты температуры распределены равномерно.

Поскольку стенка в направлении этих осей имеет бесконечно большие размеры, то соответствующие температурные градиенты Ж/йу = (к/(к = = 0, и, таким образом, влияние на процесс теплопроводности торцевых поверхностей стенки отсутствует. При этих упрощающих задачу условиях стационарное температурное поле является функцией только координаты х, т.е. рассматривается одномерная задача. Применительно к данному случаю дифференциальное уравнение теплопроводности примет вид (при д^дх = 0)

Даны граничные условия первого рода:

Рис. 7.4.

Найдем уравнение температурного ноля и определим тепловой поток Ф, проходящий через участок стенки площадью А (на рис. 1Л стенка не обозначена, поскольку располагается в плоскости, перпендикулярной плоскости рисунка). Первое интегрирование дает

т.е. температурный градиент является величиной постоянной по всей толщине стенки.

После второго интегрирования получим искомое уравнение температурного поля

где а и Ь - постоянные интегрирования.

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

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

Так как? > ? СТ2 , то проекция градиента на ось х отрицательна, как

это и следовало ожидать при выбранном направлении оси, совпадающем с направлением вектора поверхностной плотности теплового потока.

Подставляя значение постоянных в (7.24), получим окончательное выражение для температурного ноля

Линия а-Ь на рис. 7.4, так называемая температурная кривая , показывает изменение температуры но толщине стенки.

Зная температурный градиент, можно, пользуясь уравнением Фурье (7.10), найти количество теплоты 8{), проходящей за время т через элемент площади поверхности??4, перпендикулярной оси т.

и для участка поверхности площадью А

Формула (7.28) для теплового потока и поверхностной плотности теплового потока примет вид

Рассмотрим распространение теплоты теплопроводностью в многослойной плоской стенке, состоящей из нескольких (например, трех) плотно прилегающих друг к другу слоев (см. рис. 7.5).

Рис. 7.5.

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

Для первого слоя

для второго и третьего слоев

где Х 2 , А 3 - теплопроводности слоев; 8 1? 8 2 , 8 3 - толщина слоев.

На наружных границах трехслойной стенки считаются известными температуры? Ст1 и? СТ4 . По плоскостям раздела слоев устанавливаются температуры? СТ2 и? СТз, которые рассматриваются как неизвестные. Уравнения (7.31)-(7.33) решим относительно разностей температур:

а затем почленно сложим и тем самым исключим неизвестные промежуточные температуры:

Обобщая (7.36) для гг-слойной стенки, получим

Для определения промежуточных температур? СТ2 , ? СТз по плоскостям разделов слоев используем формулы (7.34):

Наконец, обобщая вывод на и-слойную стенку, получим формулу для температуры на границе г-го и (г + 1)-го слоя:

Иногда пользуются понятием эквивалентной теплопроводности Я экв. Для поверхностной плотности теплового потока, проходящего сквозь плоскую многослойную стенку,

где - суммарная толщина всех слоев многослойной стенки. Сравнивая выражения (7.37) и (7.40), заключаем, что

На рис. 7.5 в виде ломаной линии изображен график изменения температуры по толщине многослойной стенки. В пределах слоя, как было доказано выше, изменение температуры следует линейному закону. Тангенс угла наклона ср, температурной прямой к горизонтали

т.е. равен абсолютному значению температурного градиента ^1"ас1 Таким образом, по наклону прямых аЬ, Ьс и с

Следовательно,

т.е. температурные градиенты для отдельных слоев многослойной плоской стенки обратно пропорциональны теплопроводностям этих слоев.

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

Однородная однослойная цилиндрическая стенка. Найдем для стационарного режима теплопроводности температурное поле и поверхностную плотность теплового потока для однородной однослойной цилиндрической стенки (рис. 7.6). Для решения поставленной задачи используем дифференциальное уравнение теплопроводности в цилиндрических координатах.

Ось 2 направим по оси трубы. Примем, что длина трубы по сравнению с диаметром бесконечно велика. В этом случае можно пренебречь влиянием торцов трубы на распределение температур вдоль оси 2. Будем считать, что в связи с равномерным подводом и отводом теплоты температура на внутренней поверхности повсеместно равна? СТ1 , а на наружной поверхности - ? СТ2 (граничные условия первого рода). При этих упрощениях (к/ = 0, а ввиду симметрии температурного поля относительно любого диаметра?/?/?Лр = 0. Изотермическими поверхностями в этом случае будут поверхности цилиндров, соосные с осью трубы. Таким образом, задача сводится к определению одномерного поля температур? = / (г), где г - текущий радиус цилиндрической стенки.

Рис. 7.6.

Дифференциальное уравнение теплопроводности (7.19) при условии dt/d т = 0 примет вид

Введем новую переменную

которая является градиентом температур (grad ?).

Подставляя переменную и в (7.43), получим дифференциальное уравнение первого порядка с разделяющимися переменными

или

Интегрируя, получаем

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

Подставляя значение и из (7.44) в (7.45), получаем и

где ап Ь - постоянные интегрирования.

Следовательно, кривая распределения температур по толщине стенки является логарифмической кривой (кривая а-Ь на рис. 7.6).

Определим постоянные а и Ь, входящие в уравнение температурного поля, исходя из граничных условий первого рода. Внутренний радиус поверхности обозначим г х, наружный - г 2 . Соответствующие диаметры обозначим (1 Л и (1 2 . Тогда имеем систему уравнений

Решая данную систему уравнений, получаем

Уравнение температурного ноля примет вид Температурный градиент определяем но формуле (7.45):

Так как? СТ1 > ? СТ2 , а г, г 2 , то проекция grad? на радиус-вектор имеет отрицательное значение.

Последнее показывает, что для данного случая тепловой поток направлен от центра к периферии.

Для определения теплового потока, проходящего через участок цилиндрической поверхности длиной Ь, воспользуемся уравнением

Из (7.46) следует, что тепловой поток, проходящий сквозь цилиндрическую поверхность, зависит от соотношения наружного и внутреннего радиусов г 2 / г х (или диаметров с1 2 / (1 {), а не от толщины стенки.

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

Из (7.47)-(7.49) следует

Многослойная цилиндрическая стенка. Рассмотрим распространение теплоты теплопроводностью в трехслойной цилиндрической стенке (трубе) длиной А (рис. 7.7) с внутренним диаметром с1 х и наружным диаметром (1 Л. Промежуточные диаметры отдельных слоев - с1 2 и Х 2 , Х 3 .

Рис. 7.7.

Известными считаются температура? СТ) внутренней и температура? СТ4 наружной поверхности. Подлежит определению тепловой поток Ф и температуры? СТ2 и? СТз на границах слоев. Составим для каждого слоя уравнение вида (7.46):

Решая (7.51)-(7.53) относительно разностей температур, а затем почленно складывая, получим

Из (7.54) имеем расчетное выражение для определения теплового потока для трехслойной стенки:

Обобщим формулу (7.55) на и-слойную стенку трубы:
где i - порядковый номер слоя.

Из (7.51)-(7.53) находим выражение для определения температуры на границах промежуточных слоев:

Температуру? ст. +) на границе?-го и (г + 1)-го слоя можно определить по аналогичной формуле

В литературе приведены решения дифференциального уравнения теплопроводности для полого шара при граничных условиях первого рода, а также решения для всех рассмотренных тел при граничных условиях третьего рода. Мы эти проблемы не рассматриваем. За рамками нашего курса остались также вопросы стационарной теплопроводности в стержнях (ребрах) постоянного и переменного поперечных сечений, а также вопросы нестационарной теплопроводности.

где с р , Дж/(кг×К) – изобарная теплоемкость; r , кг/м 3 – плотность; l , Вт/(м×К) – коэффициент теплопроводности; w х, w y , w z – проекции вектора скорости движения жидкости; q v , Вт/м 3 – объемная плотность внутреннего тепловыделения жидкости.

Уравнение (1.12) записано для случая l=const .

Дифференциальное для твердых тел называется дифференциальным уравнением теплопроводности и может быть получено из (1.12) при условии w х = w y = w z = 0, с р = с v =с:

,

где - коэффициент температуропроводности, характеризует скорость изменения температуры в теле. Значения а = f (t) для различных тел приводятся в справочниках.

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

(1.13)

описывает нестационарное температурное поле твердых тел с внутренним тепловыделением (с внутренними источниками тепла). Такими источниками тепла могут быть: джоулева теплота, выделяемая при прохождении электрического тока по проводникам; теплота, выделяемая ТВЭЛами ядерных реакторов и т.д.

Дифференциальное уравнение теплопроводности (1.13), записанное в декартовых координатах, можно представить в цилиндрических (r , z , φ) и сферических (r , φ , ψ).

В частности, в цилиндрических координатах (r – радиус; φ – полярный угол; z - аппликата) дифференциальное уравнение теплопроводности имеет вид

(1.14)

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

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

· геометрические условия , характеризующие форму и размеры тела;

· физические условия , характеризующие свойства участвующих в теплообмене тел;

· граничные условия , характеризующие условия протекания процесса на границе тела;

· начальные условия , характеризующие начальное состояние системы при нестационарных процессах .

При решении задач теплопроводности различают:

· граничные условия первого рода , когда задается распределение температуры на поверхности тела:

t c = f (x, y, z, τ) или t c =const ;

· граничные условия второго рода , когда задается плотность теплового потока на поверхности тела:

q c = f (x, y, z, τ) или q c =const ;

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

В соответствии с законом Ньютона-Рихмана тепловой поток, передаваемый с 1м 2 поверхности в среду с температурой t ж ,

В то же время этот тепловой поток подводится к 1м 2 поверхности из глубинных слоев тела теплопроводностью

Тогда уравнение теплового баланса для поверхности тела запишется в виде

(1.15)

Уравнение (1.15) является математической формулировкой граничных условий третьего рода.

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

Контрольные вопросы и задания

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

2. Почему в правой части уравнения (1.3) стоит минус?

3. Проанализируйте с помощью справочной литературы зависимость λ(t) для металлов, сплавов, теплоизоляционных материалов, газов, жидкостей и ответьте на вопрос: как изменяется коэффициент теплопроводности с изменением температуры для этих материалов?

4. Как определяется тепловой поток (Q , Вт) при конвективной теплоотдаче, теплопроводности, тепловом излучении?

5. Запишите дифференциальное уравнение теплопроводности в декартовых координатах, описывающее трехмерное стационарное температурное поле без внутренних источников теплоты.

6. Запишите дифференциальное уравнение температурного поля проволоки, которая длительное время находится под напряжением при постоянной электрической нагрузке.

2. ТЕПЛОПРОВОДНОСТЬ И ТЕПЛОПЕРЕДАЧА
ПРИ СТАЦИОНАРНОМ РЕЖИМЕ

2.1. Теплопроводность плоской стенки

Дано: плоская однородная стенка толщиной δ (рис. 2.1) с постоянным коэффициентом теплопроводности λ и постоянными температурами t 1 и t 2 на поверхностях.

Определить: уравнение температурного поля t=f (x) и плотность теплового потока q , Вт/м 2 .

Температурное поле стенки описывается дифференциальным уравнением теплопроводности (1.3) при следующих условиях:

· т. к. режим стационарный;

· т.к. отсутствуют внутренние источники теплоты;

· т.к. температуры t 1 и t 2 на поверхностях стенки постоянны.

Температура стенки является функцией только одной координаты х и уравнение (1.13) принимает вид

Выражения (2.1), (2.2), (2.3) являются математической постановкой задачи, решение которой позволит получить искомое уравнение температурного поля t= f (x) .

Интегрирование уравнения (2.1) дает

При повторном интегрировании получим решение дифференциального уравнения в виде

Зависимость t= f (x) , согласно (2.5) – прямая линия (рис. 2.1), что справедливо при λ=const .

Для определения плотности теплового потока, проходящего через стенку, воспользуемся законом Фурье

С учетом получим расчетную формулу для плотности теплового потока, передаваемого через плоскую стенку,

Формулу (2.6) можно записать в виде

где

Величина называется термическим сопротивлением теплопроводности плоской стенки.

На основании уравнения

q R=t 1 – t 2

можно сделать вывод о том, что термическое сопротивление стенки прямо пропорционально перепаду температур по толщине стенки.

Учесть зависимость коэффициента теплопроводности от температуры, λ(t) , можно, если в уравнения (2.6) и (2.7) подставить значения λ ср для интервала температур t 1 –t 2 .

Рассмотрим теплопроводность многослойной плоской стенки , состоящей, например, из трех слоев
(рис. 2.2).

Дано: δ 1 , δ 2 , δ 3 , λ 1 , λ 2 , λ 3 , t 1 =const , t 4 =const .

Определить: q , Вт/м 2 ; t 2 , t 3 .

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

Температуры на границах слоев t 2 и t 3 можно рассчитать по уравнениям (2.8) – (2.10) после того, как найдена плотность теплового потока (q ) по (2.12).

Общий вид уравнения (2.12) для многослойной плоской стенки, состоящей из п однородных слоев с постоянными температурами на наружных поверхностях и , имеет вид

2.2. Теплопроводность цилиндрической стенки
при граничных условиях первого рода

Дано: Однородная цилиндрическая стенка (стенка трубы) с внутренним радиусом r 1 , наружным – r 2 , длиной , с постоянным коэффициентом теплопроводности λ , с постоянными температурами на поверхностях t 1 и t 2 .
(рис. 2.3).

Определить: уравнение температурного поля
t = f (r) , тепловой поток, передаваемый через стенку
Q , Вт.

Дифференциальное уравнение теплопроводности в цилиндрических координатах (1.14) для условий данной задачи:

принимает вид

Порядок решения системы уравнений (2.15) – (2.17) тот же, что и в случае плоской стенки: находится общий интеграл дифференциального уравнения второго порядка (2.15), который содержит две константы интегрирования
с 1 и с 2 . Последние определяются с помощью граничных условий (2.16) и (2.17) и после подстановки их значений в решение дифференциального уравнения (общий интеграл) получаем уравнение температурного поля цилиндрической стенки t = f (r) в виде

Если взять производную от правой части уравнения (2.18) и подставить в (2.19), получим расчетную формулу для теплового потока цилиндрической стенки

(2.20)

В технических расчетах часто тепловой поток вычисляется для 1 м длины трубы:

и называется линейной плотностью теплового потока .

Запишем уравнение (2.20) в виде

где – термическое сопротивление теплопроводности цилиндрической стенки .

Для трехслойной цилиндрической стенки (трубы, покрытой двумя слоями тепловой изоляции) с известными постоянными температурами поверхностей (t 1 и t 4 ), с известными геометрическими размерами (r 1 , r 2 , r 3 , r 4 , ) и коэффициентами теплопроводности слоев (λ 1 , λ 2 , λ 3 ) (рис. 2.4) можно записать следующие уравнения для теплового потока Q :

Температуры на границах слоев (t 2 , t 3) можно рассчитать по уравнениям (2.21).

Для многослойной цилиндрической стенки , состоящей из п слоев, формулу (2.22) можно записать в общем виде

(2.23)

Эффективный коэффициент теплопроводности для многослойной цилиндрической стенки, как и для многослойной плоской стенки, определяется из равенства суммы термических сопротивлений многослойной стенки термическому сопротивлению однородной стенки той же толщины, что и многослойная. Так, для двухслойной тепловой изоляции трубы
(рис. 2.4) эффективный коэффициент теплопроводности (λ эф) определ ится из равенства

2.3. Теплопроводность плоской и цилиндрической стенок
при граничных условиях третьего рода (теплопередача)

Граничные условия третьего рода состоят в задании температуры жидкости (t ж) и коэффициента теплоотдачи () между поверхностью стенки и жидкостью.

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

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

Теплообмен между поверхностью и средой (теплоносителем) может быть конвективным , если теплоноситель – жидкость (вода, нефть и т.д.) или радиационно-конвективным , когда теплота передается путем конвективного теплообмена и излучением, если теплоноситель – газ (дымовые газы, воздух и т.д.).

Рассмотрим теплопередачу через плоскую и цилиндрическую стенки при условии только конвективного теплообмена на поверхностях. Теплопередача с радиационно-конвективным теплообменом (сложным теплообменом) на поверхностях будет рассмотрена позже.Вт/м 2 теплопередачи (Q

если a 1 и a 2 соизмеримы.

Теплопередача через многослойную цилиндрическую стенку рассчитывается по формуле

(2.35)

где F 1 и F 2 – площади внутренней и наружной поверхностей многослойной цилиндрической стенки.