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

Курс лекций по ТМО_ Абузова Ф[1].Ф

..pdf
Скачиваний:
48
Добавлен:
17.03.2015
Размер:
3.1 Mб
Скачать

Линейная плотность теплового потока для диаметра d (d1 ≤d ≤d2 ):

ql

= π k d (tж1 tж2 );

(3.38)

q

 

=

Q

=

π (tж1 tж2 )

.

(3.39)

l

l

Rц

 

 

 

 

 

 

 

 

 

 

kΣ

 

 

Критический диаметр цилиндрической стенки

Термическое сопротивление теплопередаче однослойной цилиндри-

ческой стенки

Rцk

=

1

 

+

 

1

ln

d2

+

1

.

(3.40)

α1 d1

 

 

d1

 

 

 

2

λ

 

α2 d2

 

 

 

 

 

 

 

 

 

 

Rα

 

 

 

R λ

 

 

Rα

2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

Если d2 возрастает, то Rα2 уменьшается, а Rλ – возрастает.

Примем, что α1 , d1, α2 и λ – постоянные величины, и найдем экстре-

мум:

d

(Rцk )=

1

 

1

1

=0 .

(3.41)

 

 

 

2

dd2

2 λ d2

α2 d2

 

 

Решая (3.41), получаем критический диаметр цилиндрической стенки:

d

 

=d

 

=

2 λ

.

(3.42)

2

кр

 

 

 

 

α2

 

 

 

 

 

 

 

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

мини-

мальной величине Rцk .

Для цилиндрической стенки существует минимальное значение Rцk и

максимальное значение ql .

33

Rцk

Rц

Rц

α2

 

λ

Rцα1

Рис. 3.3. Зависимость термических сопротивлений от d2

С учётом слоя изоляции имеем

RцkΣ

=

1

+

1

ln

d2

+

1

ln

dиз

+

1

.

(3.43)

 

2 λст

d1

2 λиз

 

 

 

 

α1 d1

 

 

 

d2

α2 dиз

 

Примем, что d1, d2 , λст , α1 , α2 – постоянные величины, а диаметр изоляции dиз изменяется.

Рис. 3.4. Зависимость линейной плотности теплового потока ql от dиз

Из рис. 3.4 видно, что при возрастании dиз от d2 до dкр наложение изоляции увеличивает потери. Линейная плотность теплового потока ql у

голой трубы ( dиз =d2 ) и трубы с dиз = dиз.эф одинакова. При возрастании

34

диаметра изоляции от dкр и дальше значение ql уменьшается. Для рацио-

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

ловия

dкр <d2 .

(3.44)

Таким образом, соответствующее значение d2 расположено на рис. 3.4 правее dкр . В случае несоблюдения условия (3.44) выбирают материал

изоляции так, чтобы λ =

α2 d2

.

 

 

 

 

(3.45)

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В промышленности λ

 

 

≤0,2

 

Вт

 

. Если принять,α

 

= 4

Вт

то

 

 

 

 

 

 

 

м2 град

 

 

 

 

из

 

 

 

м град

2

 

 

 

 

 

 

 

 

 

 

2 0,2

 

Вт

 

 

 

 

 

критический диаметр d

 

=

2

λ

=

м град

=0,1 м.

 

 

 

 

 

 

 

 

 

 

 

 

кр

 

 

 

 

 

 

 

 

 

 

 

 

 

α2

 

 

 

 

 

Вт

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

м2 град

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как правило, для большинства диаметров в промышленности усло-

вие (3.44) выполняется и оптимальная изоляция выбирается исходя из минимума приведённых годовых расходов:

Рпр =Э+E К,

(3.46)

где Е – нормативный коэффициент сравнительной экономической эф-

фективности;

ΔК – капитальные вложения, [ К]= руб.;

Э – эксплуатационные затраты, [Э]= руб. .

 

 

год

 

Коэффициент Е находится по формуле

 

E =

1

,

(3.47)

 

τсл

где τсл – срок службы, [τсл ]= лет.

Заметим, что для плоской стенки оптимальной изоляции нет.

35

Тепловая изоляция

Вкачестве тепловой изоляции используют волокнистые, порошковые

ипористые материалы, заполненные воздухом. Основа препятствует кон-

векции, воздух создаёт термическое сопротивление.

Коэффициент теплопроводности в значительной степени зависит от

увлажнения, т.к. λ20°

=0,0296

Вт

, λ20°

=0,599

Вт

, поры заполня-

м град

м град

возд

 

вода

 

 

 

 

 

 

 

ются водой, есть капиллярный эффект.

Тепловую изоляцию рассчитывают по СНиП 2.04.14–88 « Тепловая

изоляция оборудования и трубопроводов», Госстрой, 1989 г.

Особенно низкий коэффициент теплопроводности имеет вакуумно-

−4

Вт

 

пористая изоляция при высоких температурах: λ =10

 

.

При этом

м град

следует учитывать предельные температуры, при которых изоляция раз-

рушается (по СНиП).

1.3.3. Нестационарная теплопроводность

Физические представления о процессах нагревания и охлаждения тел

Рис. 3.5. Нагрев тела в печи

При нагревании и охлаждении тел температурное поле зависит ещё и от времени – это нестационарные процессы теплопроводности.

Различают 2 группы процессов:

36

тело стремится к тепловому равновесию со средой, и на бесконечности эти температуры выравниваются, но на начальных стадиях температу-

ра оси отстает от температуры поверхности;

температура периодически изменяется – это процессы с периодиче-

ским нагреванием и охлаждением (регенеративные теплообменные ап-

параты).

Лекция 4

Охлаждение (нагревание) неограниченной пластины

Рис. 4.1. Охлаждение неограниченной пластины

Задача симметрична относительно оси, qv =0 . По (1.39)

∂t = a 2t . ∂τ ∂x2

Введём избыточную температуру

ϑ = t tж .

Из (4.1) и (4.2) можно получить

ϑ =a 2ϑ . ∂τ ∂x2

(4.1)

(4.2)

(4.3)

37

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

1.Геометрические: ширина и длина пластины намного больше толщины.

2.Физические: λ, с, ρ , a – const.

3.

Временные:

τ =0 ;

t(x,0) = f (x) ; J0 = f (x) −tж .

4.

Граничные:

пусть тело остывает.

 

 

x =0 ;

ϑ =0 ;

 

 

 

∂x

 

 

 

 

x = δ ;

ϑ =−

α

ϑ;

 

 

 

 

 

 

∂x

λ

−λ

∂ϑ (tc t

ж ).

 

x

 

Цель – получить зависимость вида: ϑ = f (x,τ) .

Введём безразмерные координаты и числа подобия:

относительное превышение температуры

Θ =

t tж

=

ϑ

.

(4.4)

to tж

 

 

 

ϑo

 

относительные координаты

Χ =

x

;

(4.5)

 

 

δ

 

число Био – число краевого подобия, характеризует связь между по-

лем температур твёрдого тела и условиями теплоотдачи на его поверх-

ности:

Bi =

α δ

;

(4.6)

 

 

λ

 

число Фурье – безразмерное время, характеризует связь между скоро-

стью изменения температурного поля и физическими характеристиками и размерами тела:

Fo =

α τ

.

(4.7)

 

 

δ2

 

Цель: получить безразмерное превышение температуры над темпе-

ратурой окружающей среды, т.е. перейти от ϑ = f (x,τ) к θ = f (Χ, Bi, Fo) .

38

Необходимо решить уравнение (4.3). Решение ищем в виде

 

 

ϑ =ϕ(τ) ψ(x) .

 

(4.8)

Продифференцируем (4.8) и подставим в (4.3):

 

 

 

 

ϕ(τ) ψ(x) =a 2ψ(x) ϕ(τ) .

 

(4.9)

 

 

∂τ

 

 

∂x2

 

 

 

 

ϕ (τ)

=a

ψ′′(x)

 

 

 

 

ϕ(τ)

 

.

 

(4.10)

 

 

ψ(x)

 

Уравнение (4.10) интегрируют с учётом условий однозначности, по-

сле преобразований, перейдя к безразмерным переменным, имеем.

 

2 sin µn

 

 

 

 

 

 

Θ =

 

2

 

 

 

 

 

cos(µn Χ) exp(−µn

Fo) ,

(4.11)

 

 

 

n=1 µn

+sin µn cos µn

 

 

где µn – корень характеристического уравнения:

 

ctg µn =

µn

.

(4.12)

 

 

 

 

Bi

 

 

При Fo ≥0,3 берут 1-ый член ряда (µn 1 ), получают

 

на оси:

 

 

ΘΧ=0 =N(Bi) exp(−µ12 Fo);

(4.13)

на поверхности:

 

 

ΘΧ=1 =P(Bi) exp(−µ12 Fo).

(4.14)

Функции N(Bi) и P(Bi) затабулированы и для зависимостей (4.13) и (4.14) имеются графики.

Количество теплоты, отданное пластиной в процессе охлаждения

Количество внутренней энергии тела относительно среды (в началь-

ный момент времени):

Qп = 2 δ f ρ c (tо tж ).

(4.15)

Рассмотрим процесс от τ =0 до τ = τ1 .

К моменту времени τ = τ1 внутренней энергии осталось:

39

Q = 2 δ f ρ c (tо tж ) (1−M(Bi) exp(−µ12 Fo)).

 

 

 

 

 

Q1 =2 δ f ρ c (

 

 

 

 

 

 

 

 

 

(4.16)

 

 

 

 

 

t1 tж ),

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1 – средняя температура тела.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Изменение внутренней энергии тела:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

 

ж )

1

ж

п (

1)

 

 

 

Q = Q

п

1

t

о

t

 

 

 

 

 

, (4.17)

 

 

 

 

 

 

 

 

 

Q

=2 δ f ρ c

 

 

1−

 

 

 

 

= Q 1−Θ

 

 

 

 

 

 

 

 

 

 

 

tо t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ж

 

 

 

 

где Θ1 – средняя безразмерная избыточная температура к моменту вре-

мени τ = τ1 .

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t1 tж

 

;

 

 

 

(4.18)

 

 

 

 

 

 

 

 

Θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tо tж

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

Χ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

Θ dΧ .

 

 

(4.19)

 

 

 

 

 

 

 

 

Θ

 

 

 

 

 

 

 

 

 

 

 

 

Χ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим (4.11) в (4.19) и проинтегрируем:

 

 

 

 

 

 

 

 

2 sin2 µn

 

 

2

 

 

 

 

 

 

Θ =

 

 

 

 

 

 

 

 

 

2

n sin µn cos µn

 

exp(−µn

Fo) .

 

(4.20)

 

 

 

 

 

n=1 µn

 

 

 

 

 

 

При Fo ≥0,3 берут 1-й член ряда (µn 1 ), получают

 

 

 

 

 

 

2 sin2 µ1

 

 

 

 

 

 

 

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Θ =

 

exp(−µ1 Fo)

=M(Bi) exp(−µ1

Fo) .

(4.21)

µ12 1 sin µ1 cos µ1

Значения функции М(Bi) затабулированы.

(4.22)

Аналогично решены задачи для бесконечного цилиндра и шара [1].

Влияние числа Вi на процессы нестационарной теплопроводности

1) Bi →∞ (Bi >100 ).

40

Рис. 4.2. Влияние числа Bi на процессы нестационарной теплопроводности Bi →∞

Здесь Fo0 <Fo1 <Fo2 <Fo3 .

Bi = α δ = 1 δ = Rλ . λ 1 α λ Rα

(4.23)

Если Rλ Rα Bi →∞.

Значение коэффициента теплоотдачи большое, охлаждение (нагре-

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

2) Bi →0 (Bi <0,1).

Рис. 4.3. Влияние числа Bi на процессы нестационарной теплопроводности Bi →0

41

Здесь Fo0 <Fo1 <Fo2 <Fo3 .

Если Rα Rλ Bi →0 .

Процесс охлаждения (нагрева) определяется интенсивностью тепло-

отдачи на поверхности пластины.

3) 0,1<Bi <100 .

Рис. 4.4. Влияние числа Bi на процессы нестационарной теплопроводности 0,1<Bi <100

Здесь Fo0 <Fo1 <Fo2 <Fo3 .

Если Rα ≈Rλ 0,1<Bi <100 .

Охлаждение (нагрев) определяется как внутренним, так и внешним термическим сопротивлением.

Охлаждение тел конечных размеров

Эти расчёты основываются на теореме перемножения решений.

Решением задачи представляется произведение безразмерных тем-

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

торых образовалось рассматриваемое тело.

Задачи нестационарной теплопроводности сведены в табл. 1.

42