Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теплопередача учебное пособие

..pdf
Скачиваний:
38
Добавлен:
15.11.2022
Размер:
4.77 Mб
Скачать

Требуется найти стационарное распределение температуры в массиве при заданных температурах на поверхностях труб tтр и на плоской поверхности массива tF .

Для решения задачи применим рассмотренные выше метод источников и принцип наложения.

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

В соответствии с выражением (1.77) для произвольной точки P на поверхности трубы (см. рис. 1.13) можно написать уравнения для перепадов температур относительно всех независимо действующих пар источников и стоков:

 

 

 

q

 

r′′

 

Θ

1=

 

l

 

ln

1

;

 

2πλ

 

 

 

 

 

 

r1

 

 

 

 

 

ql

 

 

r2′′

 

Θ

2=

 

 

 

ln

 

 

 

;

(1.96)

 

2πλ

 

r2

 

 

 

 

 

 

 

 

 

 

q

 

 

r′′

 

Θ

3=

 

l

ln

3

,

 

 

2πλ

 

r3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где r1, r1′′ , r2, r2′′ , r3, r3′′ – расстояния от рассматриваемой точки

трубы P до источников и стоков.

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

 

 

 

q

 

r′′

 

r′′

 

r′′

 

Θ 1+ Θ

+2 Θ = Θ3

(=x, y)

l

+ln

1

 

+ln

2

 

ln

3

.

(1.97)

 

r1

r2

 

 

 

 

2πλ

 

 

r3

 

Откуда

находим

термическое

 

 

сопротивление

труб

в массиве:

51

 

1

r′′

r′′

r′′

 

R =

 

ln

1

 

 

2

 

3

.

(1.98)

 

 

r2

 

 

2πλ

r1

 

 

r3

 

При достаточно больших отношениях h и s прибли- d d

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

 

1

 

 

h

 

 

 

 

h

 

2

 

 

 

 

 

 

 

 

 

 

R =

ln

4

1

+

 

2

 

 

.

(1.99)

 

d

 

 

 

2πλ

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Первый множитель произведения под знаком логарифма

4 h соответствует термическому сопротивлению массива при d

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

Термическое сопротивление массива увеличивается с глубиной заложения труб hd и уменьшается с увеличением расстояния между трубами sd .

Удельная мощность тепловыделения определится по фор-

муле

ql =

 

 

tтр tF

 

 

 

 

 

 

 

.

(1.100)

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

h

 

 

 

h

 

2

 

 

ln

4

1

+

2

 

 

 

 

 

 

 

 

 

 

 

 

2πλ

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Температурное поле для трех труб в массиве определяется по формуле

52

 

 

 

 

 

 

t(x, y) tF

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tтр tF

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

+ ( y0 + y)2

 

( x s)2 +

( y0 + y)2

 

( x + s)2

+ ( y0

+ y)2

 

 

ln

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.101)

 

+ ( y0 y)2

( x s)2 +

( y0 y)2

( x + s)2

+ ( y0

y)2

 

x2

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

d

 

 

 

h

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ln

4

1+ 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теплопередачу трех труб в массиве можно приближенно определить и в том случае, если задается температура окружающей среды t0 и известен коэффициент теплоотдачи на поверхно-

сти αF . Для этого используется тот же подход, что был предложен в предыдущем подразделе.

1.4.11. Теплопроводность тел с внутренним источником тепла

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

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

точника теплоты qV , Вт/м3.

53

Рис. 1.14. К расчету теплопроводности однородного цилиндрического стержня

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

2t+

qV=

0 .

(1.102)

 

λ

 

 

Теплопроводность однородного цилиндрического стерж-

ня. Рассмотрим задачу стационарной теплопроводности однородного цилиндрического стержня бесконечной длины с внутренним источником теплоты(рис. 1.14) [1, 2].

На поверхности тела задано граничное условие третьего рода: температура окружающей среды tж ; коэффициент теплоотдачи α. Задана мощность внутреннего источника теплоты qV . При

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

 

 

 

d 2t

+

1 dt

+

q

= 0

 

 

 

 

 

 

 

 

 

 

 

 

V

.

(1.103)

 

 

 

dr 2

r dr

λ

 

 

 

 

 

 

 

 

 

Граничные условия:

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

= 0

dt

 

= −

α

(tс tж ) .

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

dr r =0

 

dr r =r0

 

 

 

Необходимо найти распределение температуры по радиусу и тепловой поток на поверхности цилиндра, а также значение температур на оси t0 и на поверхности tс .

54

В уравнении (1.103)

 

произведем замену

переменной

dt dr = u , тогда уравнение примет вид

 

 

du

+

u

+

qV

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

r

λ

 

или

 

 

 

 

 

 

 

 

 

 

 

 

d (ur )

+

qV

 

rdr = 0 .

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В результате интегрирования получим

 

 

ur +

 

q r 2

 

= C1 ,

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

q r

2

 

 

 

 

 

 

r +

V

 

 

= C1 .

(1.104)

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

После преобразования и повторного интегрирования выражения (1.104) запишем

t = −

q r

2

+ C1 ln r + C2 ,

 

V

 

(1.105)

 

 

 

 

 

где C1 и C2 – постоянные интегрирования.

При r = 0 получаем, что C1 = 0 . При r = r0 , с учетом того, что C1 = 0 , получим

dt

= −

qV r0

 

 

 

 

 

 

.

(1.106)

 

 

dr r =r0

 

 

 

Подставив последнее выражение в граничное условие на поверхности цилиндра, получим

qV r0

=

α

(tс tж ) ,

(1.107)

 

 

λ

 

 

 

 

 

55

откуда

t

с

=

qV r0

+ t

ж

.

(1.108)

 

 

 

 

 

С другой стороны, температуру tс определяем по форму-

ле (1.105):

 

 

 

= −

q r 2

+ C2 .

 

 

 

 

 

 

 

 

tс

 

 

V

0

 

 

 

 

 

 

 

(1.109)

 

 

 

4λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда C2 определяем по формуле

 

 

С2

= tж +

q r

+

 

q r

2

 

 

 

V 0

 

V

0

 

 

 

.

 

(1.110)

2α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4λ

 

 

 

 

 

 

 

Подставив C1 и C2 в выражение (1.105), получим

 

 

 

q r

 

 

 

q r 2

 

 

 

 

 

r2

 

 

t = tж

+

 

V 0

 

 

+

v

0

1

 

 

 

 

.

(1.111)

 

2α

 

 

4λ

 

2

 

 

 

 

 

 

 

 

 

 

 

 

r0

 

 

При r = 0 найдем температуру на оси цилиндра

 

 

= tж +

 

q r

+

q r2

 

 

t0

 

 

V

0

 

 

V

0

.

 

(1.112)

 

 

2α

 

4λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Плотность теплового потока на поверхности цилиндра

q = α

t

t

=

 

qV r0

.

 

(1.113)

 

 

 

 

 

 

( с

 

 

 

ж )

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Количество тепловой энергии, которую отдает поверхность цилиндра длиной l в окружающее пространство, определяем по формуле

Q = qF =

qV r0

2π r0l= qVπ r02l .

(1.114)

 

2

 

 

56

Рассмотрим граничные условия первого рода, т.е. на поверхности цилиндра задана температура tс. Эти условия соответствуют частному случаю граничного условия третьего рода,

если полагать, что коэффициент теплоотдачи α → ∞

, а tс = tж .

В этом случае уравнение (1.111) примет вид

 

 

q r2

 

 

 

 

r2

 

 

t = tс +

V 0

 

1

 

.

(1.115)

4λ

 

2

 

 

 

 

 

r0

 

 

Температура на оси цилиндра при r = 0

 

t0 = tс +

q r2

 

 

 

 

V

0

.

 

(1.116)

 

4λ

 

 

 

 

 

 

 

 

 

 

Теплопроводность цилиндрической стенки. Рассмотрим задачу стационарной теплопроводности однородной цилиндрической стенки (трубки) бесконечной длины с равномерно распреде-

ленным внутренним источником теплоты qV

(рис. 1.15) [1, 2].

Внутренний радиус трубки равен r1 , внешний –

r2 . Коэффициент

теплопроводности λ – постоянный.

 

 

Примем, что градиенты тем-

 

 

ператур по продольной и окруж-

 

 

ной координатам

равны

нулю,

 

 

а изменение температуры

будет

 

 

происходить только по радиусу.

 

 

Уравнение

теплопроводности за-

 

 

пишем в виде

 

 

 

 

 

 

d 2t

 

1 dt

q

 

 

 

 

 

 

+

 

 

 

+

V

= 0 .

 

Рис. 1.15. Отвод теплоты

 

dr 2

r

 

dr

λ

 

Решением этого

дифферен-

через наружнуюповерхность

цилиндрической стенки при

циального уравнения является вы-

наличии внутренних источ-

ражение (1.105):

 

 

 

ников теплоты

 

 

 

 

 

 

 

 

 

 

 

57

t = −

q r

2

+ C1 ln r + C2 .

 

V

 

(1.117)

 

 

 

 

 

Постоянные интегрирования C1 и C2 в уравнении (1.117)

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

 

q = 0 , или

dt

 

 

 

 

 

= 0 ;

 

 

 

 

при r = r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr r

=

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

при r = r

dt

 

 

 

 

= −

α

(t

 

 

t

 

 

 

) .

 

 

 

 

 

 

 

 

 

 

2

 

ж

2

 

 

 

2

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

r =r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнения (1.104) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

= −

qV r

+

C1

.

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

dt

 

 

 

 

= −

qV r1

+

 

 

C1

= 0 , откуда C1

 

qV r12

 

При r = r1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

.

 

 

 

 

 

 

 

r1

 

 

dr r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При r = r2

с учетом найденного выражения для C1 получим

 

 

 

 

 

 

 

 

q r2

 

 

 

q r2

 

 

 

 

 

 

 

 

 

 

t2

= −

 

 

V

2

 

 

 

+

 

 

 

V

1

 

ln r2 + C2 .

 

 

(а)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

= −

 

α

 

 

(t2 tж2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

dr r

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t2 = tж2 +

 

q r

 

 

 

q r 2

 

 

(б)

 

 

 

 

 

 

V

2

 

 

 

V 1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r2

 

 

 

58

При равенстве левых частей в уравнениях (а) и (б), при-

равнивая правые части, находим

 

 

 

 

 

 

 

 

 

 

 

 

 

q r

 

q r 2

 

q

 

r 2

q r

2

 

C

 

= t

ж2

+

V 2

+

V 2

V

 

1

V 1

ln r .

 

 

 

 

 

2

 

 

 

 

r2

 

2

После подстановки постоянных интегрирования C1 и C2

в уравнение (1.117) запишем уравнение для распределения температуры:

 

q r

 

r

2

 

q r2

 

r

2

r

 

r

 

2

 

t = tж2 +

V 2

1

1

 

 

+

V 2

1

+

1

 

2 ln

 

 

 

 

. (1.118)

r2

r2

r2

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Температуру на внешней поверхности определяем из выражения

 

q r

 

r

2

 

t2 = tж2 +

V 2

1

1

 

.

(1.119)

r2

 

 

 

 

 

 

 

 

 

 

 

 

Плотность теплового потока на внешней поверхности цилиндрической стенки

q = α(t2 tж2 ) =

q r

 

r

2

 

V

2

1

1

 

.

(1.120)

2

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Температура на внутренней поверхности трубы может быть вычислена с помощью выражения (1.118) при условии, что r = r1 :

 

q r

 

r

2

 

q r 2

 

r

2

 

r

r

2

t1 = tж2 +

V 2

1

1

 

 

+

V 2

1

+

1

 

2 ln

1

1

 

.

r2

r2

r2

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим частный случай, когда на внешней поверхности трубы задано граничное условие первого рода, т.е. задана температура t2 :

59

 

q r

2

 

r

2

r

 

r

2

 

t = t2 +

V 2

1

+

1

 

ln

 

 

 

.

(1.121)

 

 

 

 

 

 

 

 

 

r2

 

 

r2

 

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

Перепад температур на внутренней и внешней поверхностях трубы определяем по формуле

 

q r

2

r

2

 

r

 

t1 t2 =

V 1

 

2

 

2 ln

2

1 .

 

r1

r2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1.5. НЕСТАЦИОНАРНЫЕ ПРОЦЕССЫ ТЕПЛОПРОВОДНОСТИ

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

Скорость теплового процесса при нестационарном режиме определяется коэффициентом температуропроводности а = λсρ ,

который здесь имеет такое же важное значение, как и коэффициент теплопроводности при стационарном режиме.

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

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

60