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

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

.pdf
Скачиваний:
53
Добавлен:
17.03.2015
Размер:
3 Мб
Скачать
☆

Таким образом, получаем: Qmz =1,14 (при прочих равных условиях).

Qmп

Лекция 13

2.3. Дополнение к теплообмену при фазовых превращениях. Кипение

Кипение - процесс интенсивного парообразования во всем объеме жидкости, перегретой относительно температуры насыщения, с образова-

нием паровых пузырей.

На p–t диаграмме кипение – во всей области от тройной точки F до критической точки K.

Рис. 13.1. Диаграмма фазового перехода вблизи точки F Различают зарождение паровых пузырей на твёрдой поверхности, к

которой подводится тепло, и в объёме.

На твёрдой поверхности кипение происходит в центрах парообразо-

вания; объёмное кипение наблюдается при значительном перегреве t (относительно ts) и при быстром сбросе давления. В современной энерге-

тике парообразование чаще встречается на твёрдой поверхности.

Механизм теплообмена при кипении отличается от механизма теп-

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

дачи (через Nu – в табл. 2) при кипении не могут быть использованы, т.к.

нарушается условие сплошности (неразрывности) потока (5.16).

103

∂ρ

+

∂(ρ wx )

+

∂(ρ wy )

+

∂(ρ wz )

=0 .

(13.1)

∂τ

∂x

∂y

∂z

 

 

 

 

 

Различают 2 основных режима кипения: пузырьковый и плёночный.

Рис. 13.2. Пузырьковый (А) и плёночный (Б) режим кипения Теплота, подводимая в единицу времени для поддержания процесса

кипения:

Q =r M′′,

(13.2)

где r – удельная теплота парообразования, [r]= кДж ;

кг

M′′– расход пара, M′′ = кг .

с

По мере увеличения Q на поверхности образуется слой пара, усло-

вия отвода теплоты ухудшаются λ

 

λ

 

: 0,017

Вт

0,68

Вт

.

пара

воды

 

 

 

 

 

м град

 

м град

 

 

 

 

 

 

Коэффициент теплоотдачи α изолирующего слоя мал, происходит пережог поверхности.

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

ной жидкости возникновение пузырьков возможно при перегреве.

Рис. 13.3. Давление на пузырёк

104

Давление на пузырёк

p

 

= p

 

+

2 σ

,

(13.3)

п

ж

R

 

 

 

 

 

 

 

 

 

 

 

 

где σ – сила поверхностного натяжения, [σ]= Н ;

м

R – радиус пузырька.

Давление жидкости

pж = ρ g h +po .

(13.4)

Рис. 13.4. Перегрев t

Перегрев t = tп −ts зависит от физических свойств и чистоты жид-

кости – чем чище жидкость, тем больше перегрев. Сильно перегретая жид-

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

ческих примесей придает кипению спокойный характер, так как все приме-

си – центры парообразования.

Если процесс парообразования идёт на твёрдой поверхности, она может быть шероховатой и на ней могут быть абсорбированы частички га-

за – всё это уменьшает молекулярное сцепление жидкости с твёрдой по-

верхностью и уменьшает t . При прочих равных условиях на полирован-

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

Отдельно рассмотрим зависимость q и t : t = f (q) – кривая кипе-

ния.

105

Рис. 13.5. Кипение воды в большом объёме при p ≈ pатм:

1 – свободная конвекция; 2 – неустойчивая область закипания; 3 – пузырьковый режим; 4 – переходный режим; 5 – устойчивое плёночное кипение; 6 – лучистый перенос тепла

При кипении происходит турбулизация пограничного слоя пузырька-

ми пара от поверхности и теплообмен сначала увеличивается с возраста-

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

танным образованием паровых зародышей в объёме жидкости. На участке от qкр1 до qкр2 t возрастает, α снижается. На несмачиваемой поверхности

плёночный режим может наступить при малых перегревах (DА) – t

уменьшается, а значение коэффициента теплоотдачи α увеличивается,

q

 

= q

 

=1,2 106

Вт

,

q

 

=q

 

= 2,5 104

Вт

. Необходимо уменьшать q

кр1

max

 

кр2

min

 

 

 

 

м

2

 

 

 

 

м

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

до значений, меньших, чем qкр2 . Момент перехода пузырькового режима в пленочный сопровождается разрушением поверхности теплообмена. Это –

106

кризисы кипения с коренным изменением механизма и интенсивности

теплообмена α

 

=

qкр1

= 4 104

Вт

.

 

 

 

 

кр1

 

t

м2 град

В области 5 лучистый теплообмен незначителен, а в области 6 – зна-

чителен за счёт теплового излучения.

qкр1 и qкр2 – смена режимов теплообмена и его интенсивности – это

кризисы I рода, а кризисы II рода – ухудшение теплоотдачи в момент высыхания кольцевой плёнки жидкости на стенке трубы при переходе стержневого режима кипения в дисперсный.

Теплоотдача при пузырьковом режиме в большом объёме

Существует много эмпирических формул для расчета α.

 

 

α =

 

 

q

 

,

(13.5)

t

 

−t

 

 

 

 

c

s

 

 

 

 

 

 

 

 

 

=B(Δt)k pl ,

(13.6)

α

где p – абсолютное гидростатическое давление (с учётом po); t, p, q – относятся к режимным параметрам;

B, k, l – зависят от формы, режима течения, природы жидкости, вре-

мени, степени обработки поверхности.

Более подробно вопрос изложен в следующем источнике:

Кутателадзе С.С., Боришанский В.М. Справочник по теплопередаче.

М-Л.: Госэнергоиздат, 1959. – 414 с.

С увеличением времени α уменьшается, так как стабилизируется процесс кипения.

Нефть и нефтепродукты при p =const кипят в интервале температур начала и конца кипения от tнк до tкк . Интервал гостируется условиями сме-

сеобразования в двигателях.

107

Структура двухфазного потока при течении внутри труб

кипящей жидкости

∙ Вертикальная труба

Поверхностью нагрева является стенка трубы. В трубе при движении

жидкости постоянно увеличивается паровая и уменьшается жидкая фаза.

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

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

Рис. 13.6. Кипение воды в вертикальной трубе:

1 – однофазная жидкость; 2 – пристенное кипение; 3 – пузырьковый (эмульсионный) режим; 4 – пробковый режим; 5 – стержневой режим;

6 – влажный пар

В вертикальной трубе различают 3 основные структуры потока:

oI – область подогрева (экономайзерный режим) – до сечения с tc = ts ;

108

oII – область кипения (испарительный участок) – от tc = ts до tc > ts ,

при этом энтальпия смеси стремится к энтальпии пара hсм → hпара ;

oIII – область подсыхания влажного пара – наблюдается в основ-

ном в длинных трубах.

Если скорость w увеличивается при заданных q, tвх и lтр, то длина I

области возрастает, а III – уменьшается.

Если q увеличивается при заданных w, tвх и lтр, то длина I области

уменьшается, а III – увеличивается.

∙ Горизонтальная (с небольшим наклоном) труба

Рис. 13.7. Кипение воды в горизонтальной трубе:

А – нормальное течение; Б – течение с увеличенной скоростью

Рис. 13.8. Пробковый (А) и кольцевой (Б) режим кипения

При кольцевом режиме кипения жидкость течёт по кольцу, внутри ко-

торого находится пар.

Объёмное и массовое паросодержание определяется

β =

Vп

,

(13.7)

 

 

Vсм

 

109

x =

Mп

,

(13.8)

 

Mсм

где Vп, Mп, Vсм, Мсм – массовый и объёмный расходы пара и смеси.

Теплоотдача при вынужденном движении кипящей жидкости в трубах

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

когда жидкость нагрета до температуры насыщения, коэффициент тепло-

отдачи может быть подсчитан по следующим формулам из [2]:

∙

при

αк

<0,5 ;

α =αw ;

αw

 

 

 

 

 

 

 

 

 

 

∙

при

αк

 

≥2 ;

α = α

;

 

 

 

 

 

αw

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

∙

при

0,5 ≤

αк

<2 ;

α

=

4 αw +αк

,

 

 

 

 

 

 

 

 

αw

αw 5 αw −αк

где αw – коэффициент теплоотдачи при вынужденном движении кипящей жидкости, определяемый по формулам однофазной жидкости;

αк – коэффициент теплоотдачи при развитом (турбулентном) пузырь-

ковом кипении в большом объёме, определяемый по формулам:

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

Nu =c Rem Pr

3

,

(13.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∙

при

Re ≥10−2 ;

 

 

 

 

 

c =0,125 ;

m =0,65 ;

∙

при

Re <10−2 ;

 

 

 

 

 

c =0,0625 ;

m =0,5 .

 

 

Re

 

=

 

 

g l

;

 

(13.10)

 

 

r

ρ′′ ν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Nu =

αк l

 

;

 

(13.11)

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l =

c

p

ρ′

 

 

,

(13.12)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

)

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

r

ρ′′

 

 

 

 

 

 

где r – теплота парообразования;

σ – поверхностное натяжение;

110

ρ′, ρ′′ – плотность жидкости и пара;

Ts – абсолютная температура насыщения.

Все теплофизические параметры взяты при температуре насыщения.

Величина l измеряется в метрах.

Формулы справедливы при:

0,86 ≤Pr ≤7,6 ;

10−5 ≤Re ≤104 ;

p = 4,5 103 ÷1,75 107 Па.

Для воды значения l и

l

в зависимости от температуры приве-

r ρ′′ ν

 

 

 

 

дены в [ 2, табл. 9-1]. Если перейти от q к

t (q = α Δt ), то приведённые

формулы для Nu примут следующий вид:

 

∙

при

λ Δt

Pr

r ρ′′ ν

 

 

 

∙

при

λ Δt

Pr

r ρ′′ ν

 

 

 

1

3

1

3

 

 

−3

 

 

λ Δt

 

1,86

 

0,952

 

 

 

 

 

 

 

 

 

≥1,6

Nu

=2,63 10

 

 

 

 

 

 

 

Pr

 

;

 

′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r ρ

 

ν

 

 

 

 

 

 

 

λ Δt

 

 

2

 

 

 

<1,6

 

=3,91 10−3

 

 

 

Pr

 

;

 

 

Nu

 

 

3

 

 

 

 

 

 

 

 

 

 

 

r ρ′′ ν

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Приведённые выше две расчётные формулы для Nu применимы

 

 

λ Δt

1

 

при условии [2]:

0,5 ≤

Pr

3

 

≤200 .

r ρ′′ ν

 

 

 

 

 

 

Лекция 14*

2.4. Приближение пограничного слоя

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

Рассмотрим продольное стационарное обтекание плоской поверхно-

сти тела безграничным потоком жидкости при ламинарном режиме. Ско-

рость и температура потока постоянны и равны соответственно wж и tж.

* Как правило, материал лекций 14 и 15 разбирается на семинаре.

111

Рис. 14.1. Изменение скорости в гидродинамическом пограничном слое

Как было показано ранее, при омывании тела поток жидкости как бы разделяется на две части – пограничный слой и внешний поток.

Теория гидродинамического пограничного слоя впервые дана Л. Прандтлем (1904 г.).

Для течения жидкости внутри пограничного слоя справедливо усло-

вие

∂wx ≠0 , т.е. скорость меняется поперёк потока. Вне пограничного

 

∂y

 

 

слоя и на его внешней границе

∂wx =0 и wx = w

æ . Во внешнем потоке

 

 

∂y

 

преобладают силы инерции, вязкостные силы здесь не проявляются. На-

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

римы. Из уравнения (5.21) для стационарного поля скоростей при омыва-

нии плоской пластины, бесконечной в направлении оси OZ, получим сле-

дующие уравнения движения (gx =0 ):

 

 

 

∂wx

 

 

 

∂wx

 

2

wx

 

2

 

 

 

 

1

 

∂p

 

 

 

 

 

 

 

 

= ν

∂

 

∂

wx

 

 

 

 

wx

 

+wy

 

 

 

2

+

 

 

2

 

−

 

 

 

;

(14.1)

 

 

 

 

 

 

 

 

 

 

 

∂x

 

 

 

∂y

 

∂x

 

 

∂y

 

 

 

ρ ∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂wy

 

 

 

∂wy

 

2

wy

 

2

 

 

 

 

1

 

∂p

 

 

 

 

 

 

 

 

 

∂

 

∂

wy

 

 

 

 

w

x

 

 

+w

y

 

 

= ν

 

 

 

+

 

 

 

 

−

 

 

 

.

(14.2)

 

 

 

 

2

 

 

2

 

 

 

 

∂x

 

 

∂y

 

∂x

 

∂y

 

 

ρ ∂y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из выражения (5.16) имеем уравнение сплошности для плоской за-

дачи:

112