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

книги из ГПНТБ / Романков, П. Г. Гидромеханические процессы химической технологии

.pdf
Скачиваний:
29
Добавлен:
22.10.2023
Размер:
10.97 Mб
Скачать

Перепад давления вследствие необратимости процессов, происхо­ дящих в системе:

Вт_

Вж_\

(5-177)

wr

wж J

 

Перепад давлений как мера работы на поверхности раздела, вы­ полненной вследствие разных скоростей газовой и жидкой фаз:

тж—ГПж-

dps = —

- J U

(5-178)

wx I

 

Возможны различные интерпретации приведенных выше соотно­ шений [63] с оценкой составляющих общего перепада давления в системе в соответствии с режимом движения потока. Подробный расчет Дрсм в зависимости от гидродинамических параметров си­ стемы Ф и /, являющихся функциями от L/G, ргІРж, Рж/Рг и а, при­ веден в монографиях Сийрде, Кафарова [60] и других.

Газосодержание

Во многих случаях для теоретического рассмотрения гидроди­ намических характеристик двухфазных потоков большое значение имеют локальные значения объемного газосодержания системы газ — жидкость [60, 61, 63]: ф = тг/т, где тг — время, в течение ко­ торого окружение данной точки заполняется газом, х — интервал времени, достаточно длинный для того, чтобы исключить влияние случайной флуктуации.

Среднее газосодержание в данном объеме V:

Ф =

-

^ | ф dV

(5-179)

 

 

V

 

Профиль скорости газа можно выразить зависимостью

 

 

 

I

 

- w-

-== (R r ) m

(5-180)

^макс

 

 

а объемное газосодержание:

 

і_

 

 

 

 

_ ^ _

=

( £ _ r ) n

(5-181)

Фмакс

 

 

 

где R — радиус трубы; г — текущий радиус; m и п — показатели степени.

Отношение средних скоростей газа и жидкости можно выразить зависимостью:

Wr _ 1—ф

(5-182)

® я _ К — <Р

 

где К ~~-параметр потока, изменяющийся в пределах 0,6 < К < 1,0 и зависящий от профилей скорости и газосодержания. Для случая, когда справедливы уравнения (5-180) и (5-181):

2 (от + n +

mit) (m +

n + 2mn)

(5-183)

{ n + 1)(2 n +

1 ) { m +

1) { 2 m + 1)

 

237

Если учесть скорость скольжения (относительную) потока как раз­ ность между скоростями газа и жидкости

®отн = WT — шж

(5-184)

го для параболических профилей скорости жидкости и газосодержания получим:

Wr __ ^QTH I

' — Ф

(5-185)

wm ~ П)ж

/С —Ф

 

Таким образом, выражая средние скорости в уравнении (5-185) через объемные расходы фаз G и L, получим уравнение для сред­ него объемного газосодержания:

Ф= - - - -р— - - - - - - -

G + L

“ +

и)отн (1 — ф) F

 

 

 

или паросодержание) потока <р как

Среднее газосодержание л(— ф

 

 

долю сечения, занятую

паром или газом, можно рассчитать, поль­

зуясь моделью равных

скоростных

напоров [64]

PCMw2cu в

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

Ф = Гі +

Рж

kn +

(1 — k) п

+

kn V* (1 + kn)'12

(5-187)

L

 

Рж

\ Рж

I

 

где k — отношение массы воды, содержащейся в гомогенной смеси

(в «стержне») к общей массе воды в системе; п = ------1; фт — массовое паросодержание, которое можно рассчитать из уравнения

Уравнение (5-187) хорошо описывает опытные данные ф при k — 0,4 независимо от режима двухфазного прямоточного потока, давления, скорости и изменения энтальпии.

Все приведенные выше зависимости относятся к двухфазным системам газ — жидкость и частично к системам жидкость — жид­ кость.

Взвешенный слой как двухфазная система

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

238

щенной оценки. Значительный интерес в этой области представляют работы Буевича и других исследователей [47, 48]. При исследова­ нии гидродинамики химических реакторов, работающих со взве­ шенным слоем мелкодисперсного катализатора, используется боль­ шое число физических моделей. В одной из них Черепановым и Гупало [48] взвешенный слой рассматривается как двухфазный идеальный жидко-упругий пластичный континуум. Если обозна­ чить жидкую или газовую фазу G, а твердую S, то уравнения дви­ жения газовой (жидкой) фазы и твердых частиц получим, исходя из баланса действующих сил:

для идеального жидкого потока

p W j W { j + P t = рF{ (г.) - X t

(5-188)

для твердых частиц

РтвО — е) U j U i j + <Tf, / , / = ( ! — е ) ( р тв7'г (тв.) - рТНг-)) +

(5-189)

где р и р — плотность и давление газа соответственно; оц и щ — ло­ кальные значения скоростей газа и твердых частиц соответственно;

Fi — компоненты

внешних массовых сил, действующих

на

поток;

Х{ — компоненты

силы взаимодействия

между потоком

и

части­

цами твердой фазы; оі — напряжения,

действующие на

частицу;

е —локальное значение порозности.

 

 

 

Уравнение (5-188) и (5-189) должны решаться совместно с урав­ нением неразрывности потока газа

(ги>і)і = 0

(5-190а)

и твердой фазы

 

[(1 — е) ui)i = 0

(5-1906)

Сила взаимодействия Хі равна нулю, если причина, ее вызы­ вающая (например, взаимное перемещение фаз), исчезает. Вели­ чина Хі может быть выражена зависимостью:

X t = у ( о ц

щ) (1 + %I w — а I )

(5-191)

где р — вязкость потока; k

и К— некоторые коэффициенты, зави­

сящие от порозности е.

 

 

Значения коэффициентов k и к в зависимости от порозности молено определить с учетом свойств потока, размеров твердых частиц и других параметров.

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

 

111 ~ iilnln- 2pnul+n

(5-192)

где

I — определяющий линейный размер; ие — \w — и \ — локаль­

ная

скорость в перемещающихся вместе с частицами координатах.

239

По Козени (см. стр. 174), I определяется как гидравлический диаметр межзернового канала: /==Ѵб^чб(1— е)-1, где d4 — экви­ валентный диаметр частицы. Тогда

1 X 1 » y } - nd4n- h n~ 2(1 - г)2~ п ßnul+n

(5-193)

Если (ічривц-1<С 1, то [X] не зависит от р и, следовательно,

п — 0.

С другой стороны, если сЦшеЦ-1 ^

1,

то

|Х|

не

зависит

от ц и

в этом случае п — 1.

Сравнивая с зависимостью

(5-191),

можно

сделать вывод, что

 

 

 

 

 

 

 

 

д

 

 

..

r

Р (1 -

е )

 

(5-194)

 

 

k

 

йче

 

k

а У

 

 

 

 

Откуда

,9

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d~4e

 

5

P rf4e

 

 

(5-195)

R

а (1 - е)2

 

1 - е

 

 

■ “

 

 

 

где k', k", а' и а" — константы.

Например, Эргун нашел экспериментально, исследуя движение потока в зернистом слое (см. стр. 176), что k' = 150; k'jk" — 87,5; a' = 6,67ІО“2; а" = 1,14■ 10~2.

Фаза S, состоящая из твердых частиц, характеризуется как упругая пластичная среда, в которой действуют напряжения, удов­

летворяющие следующим условиям:

 

 

 

а2<стк;

<гз<ак

1

(5-196)

1т„ | < /С + 1

I tg Ѳ

I

 

где an и тп — нормальное и касательное напряжения на плоскости

с нормальным вектором п

(для их определения

можно применить

круги Мора); К -— коэффициент

адгезии; ак — константа, опре­

деляемая силами когезии

между

частицами

при напряжении;

0 — угол внутреннего трения.

Величины К, Ѳ и ок зависят от формы и размеров частиц, порозности слоя е и от предыстории образования зернистой твердой фазы. Можно предположить, что стк = 0, а К и 0 не зависят от предыстории образования слоя в случае, например, если слой раз­ вивается при условии однородности перехода во взвешенное со­ стояние (гистерезисная петля BCD на рис. 5-31 отсутствует).

Связь между напряжениями оі, о2, Оз и т. д. в пространстве

слоя может быть найдена на основе общей теории упругости —

см. специальную литературу, например работы Седова [65].

Для

практически важного случая, когда внешняя сила — сила

тяжести,

т. е. FTB — Fr = —gj, где j — единичный вектор вдоль

оси у, для начала образования взвешенного слоя можно записать:

w = w vj; eWy = V ceK

(5-197)

где Vсек — секундный расход потока газа или жидкости.

240

Гидравлическое

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

слоя

выразится зависимостью

dp

 

 

 

pV^p., de

(5-198)

- ^ = - P T B f f ( l - e ) - e g p + - i r — ^

 

a связь порозности с расходом потока уравнением:

 

^

§d4e ( р тв

р) __^

а

р d4Vсек

(5-199)

 

Р Ѵ сек (1

в)

Р ( 1 — в)

 

 

Для однородного взвешенного слоя, если d4, ртв, р и р не зави­ сят от высоты слоя у, изменение по высоте dejdy равно нулю:

Пневмотранспорт

Смесь мелкодисперсного (порошкообразного материала) с воз­ духом (или другим газом) также представляет собой однородную двухфазную систему.

Экспериментально установлено, что скорость движения жидко­ сти в цилиндрической трубе (или скорость истечения жидкости из сосуда под давлением) не может превысить скорость звука С в данной жидкости. Это ограничение является решающим, если рас­ четные скорости превышают скорость звука [66—68].

Скорость звука Ссм в двухфазном потоке при пневмотранспор­ тировании можно определить по приведенному модулю упругости системы Есм:

Здесь рем — плотность двухфазного потока:

Рем = ерв + (1 — е) рхв

(5-201)

где рв и ртв — плотность воздуха и твердого зернистого материала соответственно.

Для изотермического потока

(5-202)

где р — давление; R — газовая постоянная; Т — абсолютная тем­ пература.

Таким образом, плотность смеси воздуха и транспортируемого порошкообразного материала можно выразить через скорость зву­ ка в воздухе:

Рем = 8 ~~~2 Г (1

е ) Ртв

(5-203)

Cr

 

 

Соотношение между изменением первоначального объема двух­ фазной системы и изменением давления воздуха в смеси dp можно записать, при условии, что скорость фильтрации воздуха прене­ брежимо мала (по сравнению со скоростью звука):

dV = dVa + dVrB = ZiLrfp +

I lL d p =

y ( - ^ - +

± j Z ± \ dp (5-204)

Ьв

hjB

\

-^тв J

241

где Ев = рв- ^ ---- модуль

объемной

деформации воздуха;

UOq

 

 

Е

Егв — 5Ті— 5 ~~\ — модуль объемной деформации материала твердых

частиц — модуль Юнга, — коэффициент Пуассона). Допуская, что контакт воздуха с частицами материала обеспе­

чивает изотермичность, и используя уравнение состояния р/рв = = const, получим Ев — р.

Тогда модуль упругости двухфазного потока:

 

Е(

 

 

 

Р

(5-205)

 

 

 

,

1‘■ --

 

 

 

8

8

 

 

 

 

8 +

—ß--- Р

 

 

 

 

 

ТВ

Подставив

уравнения (5-203) и (5-205) в уравнение (5-200), по­

лучим:

 

_______________________________

 

 

 

 

(1 —в)р

 

(5-206)

 

 

е2

 

 

Ртв

 

 

 

8£'тв

 

Р

 

 

 

 

 

Поскольку практически порозность двухфазного потока при пневмотранспорте в > 0,5, а давление не превышает 10 ат (т. е. р <С 1 - ІО5 кгс/м2), при модуле Еув « 0,5-10‘° кгс/м2 без заметной

£7

7.0

 

 

 

 

 

■X

0,9

 

 

 

 

 

 

0,8 ДЛ

 

 

 

 

 

0.7 h

i t 2

1

_!______ I______ I-----------1—

0,6 I

I

 

50

1G0

150

200

250

300

 

 

 

Сем, м/с

 

 

Рис. 5-37. Зависимость скорости распространения

звука

С в двухфазном

потоке

при

пневмотранспорте

апати­

 

 

тового

концентрата:

 

 

/ —при 1 ат; 2 при 5 ат.

погрешности можно принять член в квадратных скобках в уравне­ нии (5-206) равным единице. Окончательно скорость распростра­ нения звука в двухфазном потоке можно рассчитать по уравнению:

где Св — 290 м/с при Т = 293 К.

На рис. 5-37 показана зависимость скорости распространения звука в смеси при пневмотранспорте апатитового концентрата.

Можно предположить [66], что звуковая скорость является верх­ ним пределом скорости пузырей воздуха, проходящих через взве­

242

шенный слой порошкообразного материала, причем

эта

ско­

рость может быть относительно малой. Так, например,

Грек

[68],

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

Сем = 4 м/с.

Исследования гидродинамики гидро- и пневмотранспорта в плотном слое, проведенные Гаспаряном с сотрудниками [69] в лабо­ ратории двухфазных потоков АН Арм. ССР, внесли значительный вклад в развитие теории и практики взвешенного слоя.

Для двухфазных систем среди других параметров большое зна­ чение имеет величина уноса жидкости газовым (паровым) потоком или твердых мелкодисперсных частиц потоком газа или жидкости. Работы Кутепова с сотрудниками [70] по исследованию струйной сепарации капельного уноса позволили выявить область опти­ мальных нагрузок по жидкости (минимальный унос).

Обобщенное уравнение для расчета уноса и (в кг жидкости/кг пара) в случае сепарации двухфазных потоков газ — жидкость пленкой жидкости имеет вид:

 

и = 3,93. IO-7We0,88Re4'26/iCp 4

(5-208)

где

We — критерий Вебера; Re — критерий Рейнольдса;

КР =

=

' L- р- —- — критерий давления.

 

 

}Мѵж-ѵп)

 

 

Трехфазные системы

 

Одно из основных направлений организации однородного взве­ шенного слоя успешно развивающихся в последние годы, связано с созданием взвешенного слоя твердых частиц при одновременном протекании газа и жидкости [71]. Трехфазные системы широко при­ меняются в химической промышленности при экстрагировании из твердых материалов, сушке и обезвоживании, растворении и в дру­ гих процессах.

На кафедре процессов и аппаратов ЛТИ им. Ленсовета разра­ ботана конструкция аппарата (см. рис. 6-22) со щелевым подво­ дом ожижающего агента. Такой аппарат может быть использован и при одновременной обработке твердого зернистого материала га­ зом и жидкостью [72].

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

Представляют

интерес

данные, полученные в диапазоне Re =

= ®ос^чРж/рж >

500 (где

и>ос — скорость осаждения частиц, рав­

ная скорости витания; рж— плотность

жидкости; рж— вязкость

жидкости; d4— диаметр твердых частиц)

по определению порозно-

сти трехфазной системы [73]. Порозность слоя е, рассчитанная как

243

доля объема слоя, занятая газом и жидкостью, может быть опре­ делена по зависимости:

(5-209)

где С — 2,65 для движения трехфазной системы в трубе круглого сечения и С = 2,85 для кольцевого сечения; о — поверхностное на­ тяжение.

ПЕРЕМЕШИВАНИЕ ЖИДКИХ С Р Е Д *

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

Перемешивание представляет собой процесс многократного пе­ ремещения частиц текучей среды друг относительно друга во всем объеме аппарата, происходящий под действием импульса, переда­ ваемого среде механической мешалкой, струей жидкости или газа. Таким образом, перемешивание — это типичный пример смешан­ ной задачи гидродинамики.

Перемешивание осуществляют в целях: 1) обеспечения равно­ мерного распределения твердых частиц в объеме жидкости, а также равномерного распределения и дробления до заданной дисперсности жидкости в жидкости или газа в жидкости; 2) интен­ сификации нагревания или охлаждения обрабатываемых масс; 3) интенсификации массообмена в перемешиваемой системе. Иначе говоря, перемешивание — процесс, при котором градиенты темпера­ тур и концентраций в среде, заполняющей аппарат, стремятся к ми­ нимальному значению.

В связи с тем,

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

в объеме аппарата

с мешалкой, а также распределение времени

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

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

По производственному признаку возможно разделять все типы мешалок на две группы: 1) быстроходные (пропеллерные, турбин­

ные

и другие, у которых обычно окружная скорость концов ло­

*

Раздел написан при участии И. С. Павлушенко,

244

пастей порядка 10 м/с, а отношение D/dM> 3); 2) тихоходные (якорные, рамные и иные, у которых, как правило, окружная ско­ рость порядка 1 м/с, а D/dM> 2).

Физическая модель движения жидкости в аппарате с вращающейся мешалкой

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

Для вращательного движения жидкости систему уравнений Навье — Стокса можно записать в виде:

др

 

2

 

 

wT

 

дг =

Р

г

 

{ 1 . д (штг) \

(5-210)

І г ’

дг

)

 

др

 

 

 

dz ■= — 98

 

В случае плоского вращательного движения вокруг

оси z

(wp — 0 и wa = 0) система (5-210)

имеет общее решение

 

ют ~ С , г + - ^

(5-211)

На оси вращения, т. е. при г = 0, будет шт = 0, следовательно,

С2 = 0.

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

wT = сог

(5-212)

Иначе говоря, вдоль оси вращения безграничного объема жидкости в области 0 < г < гв существует бесконечный цилиндрический вихрь.

В области вне цилиндрического вихря, как следует из уравне­ ния (5-211)

wT — C2/r

(5-213)

а так как на границе обеих областей должны быть справедливы

оба уравнения (5-212) и (5-213),

то С2 =

и для периферийной

области гв < г < оо

 

 

и)т =

юл^г

(5-214)

Такое плоское течение безграничной массы жидкости назы­ вается комбинированным вихрем Ренкина [75]. Преграды, имею­ щиеся в аппарате (стенки и дно сосуда, собственно мешалка, свободная поверхность жидкости), будут оказывать значительное

245

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

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

статических

вихря,

один

из которых расположен над

зоной, оме­

ѵѵт

 

 

 

 

 

 

таемой мешалкой, а другой под

 

 

 

 

 

 

нею [76].

Диаметры

верхнего

и

 

 

 

 

 

 

 

нижнего

 

вихрей

практически

не

 

 

 

 

 

 

 

различаются и с достаточной точ­

 

 

 

 

 

 

 

ностью

определяются

эмпириче­

 

 

 

 

 

 

 

ской зависимостью [77]:

 

 

 

 

 

 

 

 

1

 

 

da

 

Reu

 

 

(5-215)

 

 

н Л

 

 

du

103 -t-l,6Reu

 

 

 

 

 

 

 

 

 

Далее, сопоставление теоре­

/

і

!

а

! ш

 

 

тической

 

и

экспериментальных

/

 

1

 

!

 

 

эпюр тангенциальных

скоростей

/

 

!

І

!

 

 

жидкости

в

аппарате

с

 

вращаю­

 

 

г ь

Гм Г„

/

г *

щейся мешалкой показывает, что

 

 

 

 

 

 

 

Рис. 5-38. Сопоставление теорети-

между

областью

центрального

ческой (1) и экспериментальной (2)

цилиндрического вихря и перифе­

эпюр

тангенциальных

скоростей

рийной областью существует пе­

 

 

жидкости:

 

вихрь;

реходная

область,

в

 

которой

/ —центральный

цилиндрический

окружная

скорость течения уже

// —переходная

область; /// —периферий­

 

 

ная область.

 

 

не пропорциональна

расстоянию

не может

быть описана

законом

от оси вращения

(5-212), но еще

площадей

(5-214). Ширина этой

переходной области при турбулентном режиме перемешивания сравнима с размерами центральной и периферийной областей (рис. 5-38), а следовательно, пренебрегать ее существованием при анализе реального движения жидкости в аппарате с мешалкой нельзя [78].

Вращательное движение жидкости, как показывают экспери­ ментальные данные [78], вызывается тем, что лопасти мешалки оказывают давление на жидкость преимущественно той частью по­ верхности, которая находится в переходной области.

Под действием центробежной силы, возникающей при вращении любого типа мешалки с достаточно большой частотой (Иец > ІО2), жидкость будет стекать с лопастей мешалки в радиальном направ­ лении. Этот поток, двигаясь в плоскости вращения мешалки, дойдет до стенки сосуда, где он разделится на две части: одну, текущую вдоль стенки аппарата вниз ко дну сосуда, и другую, текущую вверх к свободной поверхности жидкости.

Возникновение радиального течения приводит к тому, что в по­ лосе переходной области, которая ометается мешалкой, создается зона пониженного давления, куда и устремляется жидкость, теку­

246

Соседние файлы в папке книги из ГПНТБ