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

8) ТЕПЛОВЫЕ ПРОЦЕССЫ И АППАРАТЫ

.pdf
Скачиваний:
13
Добавлен:
10.10.2020
Размер:
1.26 Mб
Скачать

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

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

1 ТЕПЛООБМЕН

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

1.1 Кондуктивный теплообмен в плоской стенке

 

 

 

Рассмотрим

теплообмен

в

 

 

 

неподвижной

плоской

стенке

из

Т1

 

 

однородного

 

 

материала,

 

 

теплофизические

свойства которого

 

 

 

 

 

 

постоянны (ср,

λ, ρ = const)

 

 

 

 

Т2

( рис. 1.1).

Общее

уравнение

 

 

нестационарной

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

 

 

 

Фурье имеет вид:

 

 

 

 

λ

 

T a(

2T2

2T2 2T2 )

 

 

 

 

(1.1)

 

δ

 

 

 

 

t

x

y

z

 

 

 

 

 

 

 

 

 

у

х

Рис.1.1. Распределение температуры в плоской стенке

1

Процесс стационарный, тогда

T

0 . Считаем, что высота и длина гораздо

 

t

 

больше толщины стенки δ, следовательно, отсутствует, тогда температура изменяется тогда:

T Ty z

теплообмен по этим направлениям лишь вдоль одной координаты х,

0 .

Поскольку а 0 , имеем:

d 2T 0 dx2

Очевидным решением этого уравнения является:

 

 

 

 

 

 

 

dT

c

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T c1x c2

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при x 0

 

 

T T1 ,

 

 

 

 

 

 

 

 

 

 

 

при x

 

 

T T2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Находим с

 

T

и T

c T

,

c

 

 

 

(T2

T1 )

.

2

 

 

 

 

 

 

1

2

1

 

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dT

 

 

 

T2

 

T1

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

Распределение T по толщине :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

x

(T

 

 

T ) T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.2)

(1.3)

(1.4)

(1.5)

Из полученного уравнения (1.5) видно, что в плоской стенке распределение Т является прямолинейным.

Поток тепла за счет теплопроводности определяется по закону Фурье:

2

 

 

qм

dT

(

T2 T1

)

 

(1.6)

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

qм (

T1 T2

)

T

 

 

 

Т

 

 

(1.7)

 

 

 

 

/

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь

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

- термическое

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

 

 

 

 

 

 

 

 

 

 

 

 

Для многослойной стенки термическое сопротивление отдельных стенок

необходимо суммировать:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q м

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.8)

 

 

 

 

n

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

i

 

 

Определим количество теплоты, передаваемое за время t через площадь F:

 

 

 

Qм qм F t

 

(1.9)

Тогда расход тепла определяется как:

 

qм F

 

Qм

(1.10)

Здесь F – поверхность пластины, t – время.

Однако, приведенные расчетные формулы не всегда достаточны для практического использования. Как, например, учесть термическое сопротивление стенки при теплопередачи. Большей частью бывает, что температуры поверхностей Т1 и Т2 заранее неизвестны, но зато определены температуры Тср1 и Тср2 обеих сред, омывающих стенку, и, кроме того, соответствующие коэффициенты теплоотдачи α1 и α2. Для случая теплопередачи расход тепла запишется:

 

F

(Òñð1 Òñð2 )

Qì Ê F T

 

 

 

 

 

 

 

(

1

 

 

 

1

)

 

 

 

 

 

 

 

1

 

 

2

 

1.2. Кондуктивный теплообмен в цилиндрической стенке

Исходное уравнение в цилиндрической системе координат r, , z , имеет вид:

3

T

a(

1

 

T

)

1 2T

 

2T

 

 

 

 

 

(r

 

 

 

 

 

)

 

t

r

r

r

r 2

2

z2

(1.11)

T1

r φ

T2

R2

R1

Рис.1.2 Распределение температуры в цилиндрической стенке.

Считаем, что процесс теплообмена стационарный и длина цилиндра

достаточно велика для того, чтобы пренебречь потоком тепла к его торцам вдоль

оси

z , процесс осесимметричный. При этих

условиях температура является

функцией только одной координаты – радиуса r

(рис. 4.2):

 

 

 

1 d

 

(r

dT

) 0

или

 

 

 

 

 

 

 

 

 

 

 

 

 

r dr

 

 

 

dr

 

 

 

d 2T

 

 

1

 

dT

0

(1.12)

 

 

 

dr2

 

 

 

 

 

 

 

r dr

 

Написав уравнение (1.12) в виде:

drd ( dTdr ) 1r ( dTdr ) 0

и разделив переменные, получим:

4

 

d (

dT

)

 

 

 

 

 

 

 

dr

 

 

 

dr

 

 

 

 

 

 

 

dT

r

 

 

 

dr

 

 

 

 

Выполняя интегрирование, находим:

ln( dTdr ) ln r C

Положив, что С=lnC1 ,где С1 некоторая новая постоянная, получим:

dT C1 dr r

Вторичное интегрирование дает:

dT C1 drr

T C1 ln r C2

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

при

r R1

T T1

при

r R2

T T2

T1 C1 ln R1 C2

T2 C1 ln R2 C2

Отсюда

C1

 

T2 T1

 

ln

R2

 

 

 

 

 

 

 

R1

Окончательно:

 

 

T T1

C2 T1 C1 ln R1

 

T2

T1

ln

r

 

 

 

 

ln

 

R2

R

 

 

R1

1

 

 

 

 

 

 

(1.13)

(1.14)

Как видно из уравнения (1.14) имеет место логарифмический закон распределения температуры по радиусу цилиндра.

Градиент температуры на внутренней поверхности цилиндра равен:

5

 

dT

 

T2 T1

(

 

)r R1

 

 

 

dr

(ln

R2

)R

 

 

 

 

 

 

 

 

 

 

R1

1

 

 

 

 

 

 

В правой части уравнения для любого r в знаменателе вместо R1 брать r.

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

qM

 

dT

 

T1 T2

dr

(ln

R2

) r

 

 

 

 

 

 

 

 

 

 

 

 

 

R1

необходимо

(1.15)

Как видно из уравнения (1.15) тепловой поток зависит от координаты r, уменьшаясь с возрастанием r.

Количество теплоты находим как:

QM qM F(r) t

(1.16)

Здесь F=2π rL – внутренняя поверхность цилиндра, t – время, L – высота цилиндра.

Расход тепла определяется как:

Если труба

получим:

(1.18)

Здесь T T1

.

 

 

 

 

 

Q qМ

F (r)

(1.17)

многослойная и состоит из n слоев,тогда для потока тепла

qM

 

 

 

 

T

 

 

 

 

 

 

 

 

 

n

 

1

 

 

Ri 1

 

 

 

 

 

 

 

 

 

 

ln

r

 

 

 

 

 

 

 

 

 

 

 

Ri

 

1

 

i

 

Tn - общая разница температуры.

Зависимость qм и F от радиуса r не позволяет использовать традиционную форму уравнения теплопередачи для цилиндрической стенки. В этом случае

используется коэффициент теплопередачи KTL отнесенный к единице длины:

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Я T Я )

 

Q K

TL

L (T

,

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

KTL

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

ln

R2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

1

 

 

 

R

 

R

2

 

 

 

 

 

1

 

 

 

 

 

1

 

2

 

 

 

6

Я

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

поверхность.

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

1.3 Конвективный теплообмен

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

Конвективный перенос теплоты описывается уравнением Фурье-Кирхгофа. Закономерности течения среды описываются уравнениями Навье-Стокса (ламинарный режим) и Рейнольдса (турбулентный режим), а также уравнением неразрывности.

Исследование закономерностей конвективного теплообмена можно провести

визотермической и неизотермической постановке.

Визотермической постановке сначала решаются уравнения Навье-Стокса и неразрывности, затем полученные значения скоростей используются для решения уравнения Фурье-Кирхгофа. Полученные таким способом значения коэффициентов теплоотдачи впоследствии уточняются, корректируются.

Внеизотермической постановке уравнения Навье-Стокса, неразрывности и Фурье-Кирхгофа решаются совместно, с учетом зависимости теплофизических свойств среды от температуры. Как показывают экспериментальные данные,

зависимости ср(Т), λ(Т) и ρ(Т) слабые, а µ(Т) – очень сильная. Поэтому обычно учитывается только зависимость µ(Т). Она, эта зависимость, может быть представлена в виде зависимости Аррениуса или, проще, в виде алгебраического уравнения. Таким образом, возникают так называемые сопряженные задачи.

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

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

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

TC f (x, y, z, t)

(1.19)

7

простейший случай, когда ТС=const;

- граничные условия второго рода – задается распределение теплового потока на стенке:

qC

 

dT

 

 

 

f (x, y, z, t)

(1.20)

dn

n 0

 

 

 

 

 

 

 

 

 

 

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

 

dT

 

 

 

(T Я Т Г )

 

dn

(1.21)

 

 

 

 

 

 

n 0

 

 

 

 

 

 

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

1.3.1 Гидродродинамический и тепловой пограничные слои на плоской пластине

Рассмотрим поток, обладающий неизменными теплофизическими характеристиками (ρ, μ, λ, cp = const), совершающий вынужденное движение вдоль плоской полубесконечной тонкой пластины и обменивающейся с ней теплом. Предположим, что неограниченный поток со скоростью wx0 и температурой Т0 набегает на полубесконечную пластину, совпадающую с плоскостью х – z и имеющую температуру Тст = const.

Выделим гидродинамический и тепловой пограничные слои с толщиной δг и δт соответственно (область 99% изменение скорости wx и температуры T). В ядре потока wx0 и Т0 постоянны.

Ламинарные пограничные слои (рис. 1.3)

Проанализируем уравнения неразрывности и Навье-Стокса. Задача

двумерная, поскольку

wz,

 

=0. По экспериментальным данным известно, что в

 

 

 

z

 

 

 

 

 

 

 

 

 

гидродинамическом пограничном слое

p

0 . В ядре потока wx0=const, поэтому,

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

 

p

согласно уравнению Бернулли x

0 , в погран. слое то-же самое:

x 0

 

 

 

 

2 w

X

 

2 w

X

 

 

Как известно х»δг

поэтому

 

 

 

x2 .

 

y 2

 

 

 

Следовательно, имеем:

8

 

 

wX

wY

0

,

 

 

(1.22)

 

 

 

x

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

wX

w

X

wy

w

x

 

2 w

X

(1.23)

 

 

 

 

 

x

y

y

2

 

 

 

 

 

 

 

 

уУ

w

0

 

U

0

 

 

 

 

 

 

 

 

 

 

x

, ,T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Wx

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

w Wx

x

 

 

 

 

(T-TСТ)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ Г (х)

δ Т (х)

ТСТ

0 0

х X

Рис. 1.3. Гидродинамический и тепловой ламинарные пограничные слои на плоской пластине.

Записывать аналогичные уравнения для оси у не имеет смысла, так как wy может быть найдена из уравнения неразрывности (1.22).

Используя аналогичные процедуры можно упростить и уравнение ФурьеКирхгофа:

wX

T

wy

T

a

2T

(1.24).

x

y

y

2

 

 

 

 

 

Система дифференциальных уравнений (1.22)-(1.24) составляет изотермическую математическую модель плоского стационарного теплового ламинарного пограничного слоя. Сформулируем граничные условия на границе с

9

пластиной, т.е при у=0: при любом х скорость wX=0 (условие прилипания). На границе и вне гидродинамического погранслоя, т.е. при у≥δг (х), а также при х=0 для любого у: wX=wX0. Для поля температуры аналогичные рассуждения.

Итак, граничные условия:

wx (x, 0) = 0, x > 0;

wx (x, ∞) = wx0,

wx (0, y) = wxo ;

(1.25)

T (x, 0) = TCT, x > 0;

T (x, ∞) = T0,

T (0, y) = T0

(1.26)

Точное решение задачи в виде бесконечных рядов было получено Блазиусом. Имеются более простые приблежённые решения: методом интегральных соотношений (Юдаев) и по теореме импульсов (Шлихтинг). А. И. Разиновым задача была решена методом сопряжённого физического и математического моделирования.

Были получены профили скоростей wx (x, y), wy (x, y) и температур Т, а также толщины пограничных слоев δг (x) и δт (х):

Г (x) А

 

x

 

(1.27)

 

 

 

 

w

0

 

 

 

 

X

 

 

 

 

 

T (x) Г (x) Pr 1/ 3 , Pr ≥ 1

(1.28)

Pr = ν/a

У Разинова А=5.83; Юдаева А=4.64; Блаузиуса =4; Шлихтинга А=5.0. Примерный вид найденных зависимостей приведён на рисунке 4.3.

Как известно, для газов Pr ≈ 1, капельных жидкостей Pr > 1.

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

теплоотдачи. Локальные значения (x)

и NuГ,x:

 

 

 

 

 

 

 

 

 

 

 

 

(x)

 

2

 

0,343

w0

 

 

 

 

 

 

x

,

 

 

(x)

x

 

 

 

 

 

 

 

 

 

Nuг,x =

x

0.343 Re0.5X ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Rex

w0

x

 

 

 

 

x

 

 

 

(1.29)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( ) и Nu Г ,

по участку длиной :

Усреднённые значения

 

 

1

(x)dx

 

4

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

o

 

 

( )

 

10