книги из ГПНТБ / Романков, П. Г. Гидромеханические процессы химической технологии
.pdfПерепад давления вследствие необратимости процессов, происхо дящих в системе:
Вт_ |
Вж_\ |
(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