Добавил:
Без скрытых скриптов, криптомайнинга, вирусов и прочего, - чистая литература. 你好,所有那些谁花时间翻译中国 Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Тишин ВБ Новоселов АГ Процессы переноса в технолог аппаратах

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

T T (bt) p ,

(2.58)

н

 

где Tн – начальная температура.

Решим уравнение (2.58) относительно времени t :

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(T T ) p

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

н

 

 

.

 

 

 

 

 

 

(2.59)

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Производная функции (2.58) будет иметь вид

 

 

T

 

 

dT

b p pt p 1

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или с учѐтом уравнения (2.59) она примет вид

 

 

 

dT

 

b p p T

 

 

 

p

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

T

 

p .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

н

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив полученный результат в равенство (2.57), запишем

 

 

 

 

 

 

 

 

 

 

 

 

 

p

1

 

 

 

 

 

c

p

M b p p(T

T )

 

p

 

 

μ

 

 

 

 

 

 

 

 

н

 

 

 

 

 

.

(2.60)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A2n2d 3

Для проверки предлагаемого метода реологических исследований нами были проведены опыты по установлению закономерности перемешивания чистого глицерина в лабораторном смесителе для теста с двумя Z-образными мешалками, вращающимися с различными скоростями. Отношение n1/n2 = 1,5. Эксперименты проводились при следующих числах оборотов мешалок: n2 = 0,4; 0,6; 0,8; 1 об/c; n1 = 1,5 n2; среднее число оборотов nc 1,25 n2 . Диаметр мешалок

одинаков и равен 0,068 м. Начальная температура в опытах колебалась в пределах 19–24 °С. Масса загружаемой в смеситель жидкости

была постоянной, M

V 0,78 кг.

61

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

T

Tн

bt ,

(2.61)

поэтому в равенствах (2.54) – (2.60)

p 1. Поскольку перемешивание

проводилось в ламинарном режиме,

то и показатель степени m 1.

Производная dT / dt b зависит

от

числа оборотов и

физических

свойств жидкости. Например, при Tн = 19,5

n 0,5

и b 2,14 10 4 ;

 

c

 

при Tн = 20,1 nc 0,75 и b 4,1 10 4 . Для неньютоновских сред показатель степени p будет отличаться от единицы.

На основе данных по кинетике разогрева жидкости из уравнений (2.48), (2.49), (2.53) и (2.60) определялись коэффициент A и показатель степени m в уравнении (2.49). В результате обработки опытных данных получено уравнение для расчѐта критерия мощности в исследуемом смесителе:

K N

3450

.

(2.62)

 

 

Re1ц,1

 

В уравнении (2.62) вызывает некоторое сомнение величина показателя степени при критерии Рейнольдса. Известно, что при ламинарном режиме перемешивания должно быть m 1. Отклонение показателя степени от единицы в уравнении (2.62) связано, скорее всего, с ошибками в эксперименте.

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

2.5. Турбулентные течения

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

62

локальной скорости u и локального давления p . Тогда их локальные значения будут равны:

u

 

u';

p p p ,

(2.63)

u

где u' и p' – пульсационные составляющие локальной скорости и ло-

кального давления.

В проекциях на координатные оси (ограничимся осью 0x) получаем ux' , 'x , 'xy , 'xz .

Полагая, что уравнения движения в напряжениях могут быть пригодны для описания турбулентных течений, после введения в уравнения (1.24)–(1.26) осредненных параметров, получим уравнение движения жидкости при турбулентном режиме в проекции на координатную ось 0x :

 

d

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

'

 

 

'

 

 

 

 

 

 

 

 

 

 

 

'

 

 

'

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

'

 

 

F

 

(

 

 

 

 

 

u

u

)

 

 

(

 

 

 

 

u

u

)

 

 

 

 

 

(

 

 

u

u

), (2.64)

 

 

 

 

 

 

x

 

 

 

x

x

 

 

 

 

xy

x

y

 

 

 

 

 

 

xz

x

z

 

dt

x

x

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x'

 

x'

 

 

 

 

x ;

 

 

x'

 

y'

'xy ;

 

 

 

 

 

x'

 

z'

 

 

 

'xz .

 

 

 

 

 

 

(2.65)

 

 

 

 

 

 

 

 

u

u

 

 

 

 

u

u

 

 

 

u

u

 

 

 

 

 

 

 

 

 

Аналогичные уравнения можно записать для осей 0y и 0z. Уравнения (2.64) носят имя О. Рейнольдса.

Наличие в турбулентном потоке пульсаций скорости приводит к сглаживанию профиля скорости по его сечению. Исследования турбулентных течений показали наличие двух зон с различным хаpактеpом изменения осредненной локальной скорости u . У твердой поверхности происходит резкое изменение скорости в пристенном слое толщиной (рис. 2.7, а), значительно меньшем по сравнению с поперечным pазмеpом канала. Считается, что в пределах этого слоя жидкость движется ламинарно.

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

63

а

u

б

 

y

 

 

u( y)

0

x

Рис. 2.7. Распределение скорости по сечению потока при турбулентном течении

Пpенебpегая массовыми силами ( Fx 0 ) и полагая движение равномерным, из уравнения (2.64) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

 

'

 

 

 

 

 

 

 

u

 

 

' '

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

 

 

xy

uxuy )

 

 

 

 

 

 

 

 

uxuy

0.

(2.66)

 

 

 

 

 

 

y

 

y

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В ламинарном пристенном слое толщиной

турбулентные

напряжения t

 

 

x'

 

y'

0 и из уравнения (2.66) следует

 

 

 

 

 

 

u

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После интегpиpования этого выражения имеем

d

 

x / dy

c1.

u

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

y 0

 

x

0,

u

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

c1

0 .

Таким

образом,

касательные

напряжения

в пределах пристенного ламинарного слоя изменяются по закону Ньютона (1.7)

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

ux

0 y

.

(2.67)

 

Величина

64

0

u*

(2.68)

 

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

Из уравнений (2.67) и (2.68) следует

 

 

 

 

 

x u*y

 

 

u

(2.69)

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

u*

ν

 

Обозначим

отношение

 

x / u*

 

– безpазмеpная

скорость,

u

 

а отношение u* y /

– безразмерная координата. В гидромеханике

их называют универсальными координатами. Введя

в уравне-

ние (2.69) и , получим

 

 

 

 

 

 

 

 

 

 

 

 

.

 

(2.70)

Уравнение (2.70) описывает профиль скорости в пристенном ламинарном слое.

Найдем распределение скоростей в турбулентном ядре потока, исходя из предположения, что на границе раздела ламинарного слоя и турбулентного ядра t 0 c1. Так как в турбулентном ядре моле-

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

0 и из уравнения (2.66) следует

 

 

x'

 

y'

0 t .

(2.71)

u

u

Согласно теории турбулентности Л. Прандтля [8, с. 165–185],

 

 

 

' d

 

x

 

 

 

 

 

' d

 

x

 

 

u

'

k l

u

и u

'

k

 

l

u

.

(2.72)

x

 

 

 

 

y

2

 

 

 

 

 

1

 

dy

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

65

Подставляя значения ux' и u'y в уравнение (2.71), получим

 

 

 

 

 

 

 

 

 

 

 

 

2 d

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

'

 

 

u

x

 

 

 

 

 

 

 

 

 

 

 

 

t

k k

2

(l

)

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

k

и k

2

– коэффициенты пpопоpциональности;

l

k k

2

(l' )2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

длина перемешивания.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С учетом значения l

предыдущее равенство примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

2

 

u

x

.

 

 

 

 

(2.73)

 

 

 

 

 

 

 

t

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Значение l 2 (d

 

x / dy)

принято

называть

коэффициентом

 

 

u

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

Подставив в уравнения (2.73) значения ,

получим

 

 

d

 

x

 

 

 

u

 

 

t

 

.

 

(2.74)

 

 

dy

 

 

 

 

 

Согласно теории Прандтля, длина пути перемешивания прямо

пропорциональна расстоянию от стенки,

т. е. l

y. После подста-

новки значения l в уравнение (2.73), с учѐтом равенств (2.68) и (2.74), запишем

 

d

 

x

 

u χy

u

,

 

 

 

*

dy

 

 

 

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

d 1 d .

66

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

 

 

C

1

ln .

 

(2.75)

 

 

Экспериментально установлено, что C

5,5;

0,4. Решая со-

вместно уравнения (2.50) и (2.75), с учетом значений C и , найдем

безpазмеpную толщину пpистенного слоя

11,5 .

 

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

η при η 11,5;

(2.76)

5,5 2,5ln при η 11,5.

Согласно системе уравнений (2.76), безразмерная толщина ламинарного слоя может быть найдена из равенства

u*

11,5.

(2.77)

 

Из уравнений (2.36), (2.44) и (2.68) можно получить

u*

w

 

.

(2.78)

 

 

8

 

 

Из уравнений (2.77) и (2.78) следует, что толщина гидродинамического ламинарного слоя для двухслойной модели

32,5

d

 

.

(2.79)

 

 

 

 

Re

 

В уравнении (2.66) выражение в скобках представляет собой полное касательное напряжение

0 t .

67

Подставляя в это равенство значения 0 и t из уравнений (1.7) и (2.74), запишем

() dduyx . (2.80)

После приведения уравнения (2.80) к безразмерному виду по-

лучим

 

 

 

 

 

1 1

 

 

d

.

(2.81)

 

 

 

 

 

d

 

Из уравнений (2.81) и (2.76) следует:

0 при 11,5,

(2.82)

 

1

при

11,5.

2,5

 

 

 

Таким образом, отметим ещѐ раз, что, согласно двухслойной модели турбулентного потока, турбулентные пульсации в пристенном слое равны нулю. Эта модель пригодна только для решения задач, связанных с переносом количества движения. Решения различных задач, связанных с переносом количества движения, приводятся в курсах «Гидравлика» и «Процессы и аппараты». Попытки применить уравнения (2.82) для решения задач теплообмена и массообмена к успеху не привели и не могли привести. В последующих разделах разберѐмся в причинах этих неудач.

2.6.Гидродинамика газожидкостных потоков

Сдвижением двухфазных газожидкостных потоков приходится иметь дело в различных отраслях пищевой и микробиологической промышленности. Примерами тому могут быть: в пищевой промышленности – газирование диоксидом углерода различного рода безал-

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

68

Во всех приведѐнных примерах имеют место процессы переноса количества движения, теплоты и массы.

2.6.1. Основные понятия и определения гидродинамики многофазных потоков

Приведенные скорости фаз есть отношение их объемных расходов к площади сечения потока. Для двухфазных газожидкостных потоков запишем

w

Qг

;

w

Qж

,

(2.83)

г

Sa

ж

Sa

 

 

 

 

где wг и wж – приведѐнные скорости газа и жидкости;

Qг и Qж

объѐмные расходы газа и жидкости, соответственно; wсм

wг wж

приведенная скорость смеси.

Истинные скорости фаз есть отношение объемного расхода данной фазы к площади сечения, занятой ею, т. е.

w

Qг

;

w

Qж

,

(2.84)

 

 

г.и

Sг

ж.и

Sж

 

 

 

 

где wг.и и wж.и – истинные скорости газа и жидкости;

Sг и Sж

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

Истинные объемные доли фаз. Истинной объемной долей фазы называется отношение площади сечения, занятой данной фазой, к площади сечения аппарата. Доли газа (газосодержание), жидкости и клеток, соответственно, будут равны:

г

Sг

;

ж

Sж

;

к

Sк

,

(2.85)

S

S

S

 

 

 

 

 

 

 

 

 

 

 

где г , ж , к – объѐмные доли газа,

жидкости и клеток в среде;

Sк – площадь поперечного сечения аппарата, занятого клетками.

Так как Sг Sж S , то из уравнения (2.85) для двухфазной

модели

среды получим

г

ж

1, для трѐхфазной модели –

 

 

 

 

 

г

ж

к

1.

 

 

 

 

 

 

 

69

Из равенств (2.84) и (2.85) для двухфазной модели следует:

w

wг

;

w

wж

 

 

wж

.

 

(2.86)

 

 

 

 

 

 

г.и

 

 

ж.и

 

1

 

 

 

 

 

г

 

ж

г

 

 

 

 

 

 

 

 

В дальнейшем при использовании параметра

г

для простоты

 

 

 

 

 

 

 

 

 

 

изложения будем оперировать только понятием «газосодержание». Истинное объемное газосодержание является основным гидро-

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

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

w0

wг

 

 

wж

 

.

(2.87)

 

1

 

 

 

г

 

г

 

 

 

 

 

 

Значения г и w0 зависят от конструктивных особенностей

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

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

тах wж 0 и относительная скорость газа равна истинной.

Величина, называемая объемным расходным газосодержани-

ем, есть отношение

 

 

 

 

 

 

0

 

Qг

 

wг

.

(2.88)

 

 

 

Qг

Qж

 

wсм

 

 

 

 

В случае, если w0 0 , имеет место г

0 . В барботажных

культиваторах, где Qж 0 , 0 1. Расходное газосодержание значительно реже используется в качестве характеристики газожидкостных

70