
Белоглазов Фильтрование технологических пульп
.pdfобусловленный формой |
mrm<rm.rQ,„ |
поверхности частиц и пороз- |
||||
ностью слоя. |
' |
остоянием |
||||
Установлена следующая зависимость: |
|
|||||
|
|
Д, |
= Y R p -/~ 2 |
(1.89) |
||
|
|
|
|
|
' |
|
потока. J |
к° эффициенты, |
обусловленные |
режимом движения |
|||
АО 2 (ри?Ф1.7)И0Что п о 7 тв% ™ ется“ e |
S |
значения от 1 |
||||
44], в которых предложено пепепя л |
! ^ |
Р ^ о м |
исследований [5, 6, |
|||
виде суммы двух слагаемых |
в о ш |
о ^ ™ В СЛ° е пРеАставлять в |
||||
потока в первой степени, а во второе - |
с к о р о Т т Г п о т ^ |
В своих работах [5, 6] ЭрГУН,„ Допуская, что зернистый слои представляет собой со вокупность извилистых ш е роховатых трубок (причем все зернистые слои имеют при мерно одинаковые характе ристики шероховатости), для монодисперсного слоя сфери
ческих частиц предлагает со отношение
|
|
|
31 сл Ржи |
|
1 - |
е |
(1.90) |
||
105 |
Re |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
м н й 17 \ 3аВИСИМОСТЬ показателя сте- |
Экспериментальные |
дан |
|||||||
данным ЛеваРИТеРИЯ РеЙН0ЛЬ^ |
по |
ные |
позволяют |
|
определить |
||||
|
|
значение |
коэффициента |
со |
|||||
|
|
противления: |
31 |
сл |
/4 |
= |
1 7S |
||
|
|
|-р |
|
|
|
|
|
it ГО. |
|
|
|
Тогда |
|
|
|
|
|
|
|
ОС |
|
е3 |
' |
|
|
|
|
(1-91) |
|
м ера% Г ^ им Т и Г нГ у“ |
|
|
*' |
м |
10 • Оказыва- |
||||
ется, что пои п п п и «п . 1 т „ |
|
ж |
Аавлений в зернистом слое м ож еГ бы ^Г н а^ен и ^ о о т Х е н и я ^
Ар = |
1,75рж[/2 1 - в |
ОС |
(1.92) |
40 |
|
которое называется уравнением Эргуна. При малых значениях параметра ржШ /[цж(1 — в)] оно переходит в уравнение Козени—
Кармана, а при больших — в уравнение Бурке—Пламмера.
В ряде работ [46, 80] при расчете сопротивления слоя в условиях турбулентного режима предлагается учитывать проницаемость слоя
и шероховатость частиц: |
|
|
|
|
|
|
— |
= ^ |
и |
+ 2рж— U2, |
(1.93) |
|
1*ОС |
к п р |
d 3 |
|
|
где кпр — |
приведенный коэффициент проницаемости слоя, м2; Д' — |
||||
относительная шероховатость твердых частиц. |
|
||||
Этот |
подход наиболее |
обстоятелен, посколькувходящий в |
|||
уравнение параметр кпр = |
рж(/Лос/Др позволяет учитывать реальные |
свойства пористых сред комплексно, в то время как сама пористая среда считается псевдогомогенной системой.
Для засыпок, образованных монодисперсными частицами не правильной формы, коэффициент проницаемости может быть опре делен по уравнению [51]
к |
= - |
— |
d 2. |
(1.94) |
"Р |
9 |
к„ |
3 |
|
|
|
'■к |
|
|
Для определения коэффициента проницаемости засыпок, обра зованных полидисперсными частицами с известной функцией рас пределения пор по эффективным диаметрам, используется выра жение
g |
^max |
(1.95) |
кпр = — |
j d 3f(d3)d(d3)2 . |
min
В критериальном виде для расчета гидравлического сопротив ления слоя зернистых материалов при турбулентном режиме при менимы зависимости: в переходной области (60 < Re < 7000) —
Eu = 7,6Re-°'2r, |
(1.96) |
в области развитого турбулентного течения (Re > 7000) —
Eu = 1,ЗГ. |
(1.97) |
При расчете гидравлического сопротивления среды важно знать границу перехода ламинарного режима движения жидкости в тур булентный, характеризуемую величиной критического числа Рей нольдса ReKp. Сложность заключается в том, что для различных
пористых сред и засыпок из частиц неправильной формы величина ReKp не является постоянной и зависит от порозности, проницае
41
мости и относительной шероховатости поровых каналов. Практи чески неизменными значения Rекр остаются только для засыпок с
схожими структурами порового пространства и одинаковой шерохо ватостью частиц.
Средние значения критических чисел Рейнольдса для засыпок из различных материалов по данным [72] приведены в табл. 1.12.
Таблица 1.12
Средние значения Reip для засыпок из различных материалов
Виды засыпок |
Размер частиц, мм |
Порозность |
Reкр |
Активированный уголь, глинозем, |
1,50-10,0 |
0,48-0,52 |
0,29 |
силикагель |
|
|
|
Антрацит, керамзит |
• 0,25-10,5 |
0,38-0,48 |
0,16 |
Гравий |
2,50-16,0 |
0,30-0,43 |
0,55 |
Керамическая крошка |
0,20-7,5 |
— |
0,11 |
Металлургический шлак |
0,20-1,5 |
— |
0,12 |
Однородные сферические частицы |
0,25-30,0 |
0,35-0,52 |
0,59 |
Однородный песок |
0,25-5,5 |
0,35-0,44 |
0,27 |
Щебень |
3,00-11,0 |
0,40-0,45 |
0,39 |
Значения Rexp можно найти и через приведенный коэффициент проницаемости слоя £пр, используя зависимость [115]
^■ехр = P*U ijknp Дгж. |
(1.98) |
Результаты графической обработки опытных данных [5] в виде зависимости ReKp = f(knp/e2/3) представлены на рис. 1.8.
Рассмотренные выше уравнения Козени—Кармана, Бурке— Пламмера и Эргуна относят к группе так называемых капиллярных моделей, сопротивление слоя в которых зависит от извилистости каналов lnop/h oc, пористости б и удельной поверхности пор / .
Однако кроме них в литературе предложено множество иных иде ализированных моделей фильтрации [6, 41, 44, 46, 80, 90], обзор которых дан в работе [48].
Согласно струйным моделям фильтрации, основным фактором, определяющим характер движения потока в пористой среде, явля ется эффект сжатия и расширения струй, возникающий под вли янием конфигурации отдельных ячеек слоя [6, 44, 90]. Для системы со сложной структурой порового пространства уравнения струйной модели малопригодны, поскольку требуют громоздких расчетов эф фекта сжатия и расширения струй.
В моделях внешнего обтекания частиц дисперсной фазы жид костью, фильтруемой через слой, течение рассматривается как со вокупность последовательных обтеканий отдельных частиц либо в
Рис. 1.8. Аппроксимации экспериментальных зависимостей ReKp = f(knp/e2'3)
для различных материалов:
1 — засыпки из полированных сферических частиц; 2 — пористые материа лы из округлых гладких частиц; 3 — крупнозернистые засыпки из шерохо ватых частиц неправильной формы; 4 — мелкозернистые засыпки
предположении квазистационарности поля скоростей вокруг от дельной частицы, либо с учетом влияния пульсаций скорости в условиях стесненного обтекания частиц в абсолютно однородных слоях. Эти модели также предполагают определение формы и коэф фициентов гидравлического сопротивления одиночных частиц, что является непростой задачей для частиц неправильной формы.
1.5. Тепломассообменные процессы в зернистом слое
Теоретического описания закономерностей тепло- и массообмена в пористых средах на данный момент пока не предложено, что можно объяснить сложностью математических расчетов и большим числом влияющих факторов, учесть которые не представ ляется возможным даже в идеализированных случаях. Поэтому в основу расчета интенсивности тепломассообмена в ходе фильтра ции пульпы на фильтрах положены идеализированные мате матические описания, базирующиеся на критериальных регресси онных уравнениях и позволяющие проводить технические расчеты
42 |
43 |
в строго ограниченном диапазоне значений критериев Re, Рг, Ре и других. Применяемые уравнения носят эмпирический характер, что потребовало большого объема экспериментальных исследо ваний.
Отсутствие общепризнанных и нормативно закрепленных методик анализа тепло- и массообменных процессов в пористых средах затрудняет исследования в этой области и препятст вует созданию комплексной теории тепломассообмена в пористых слоях.
Обзор литературы показывает, что при обработке опытных данных по тепломассопереносу при фильтрации различные ис следователи в качестве определяющих параметров процесса используют разные величины. Приводимые в литературе регрес сионные уравнения отражают ход процесса в различных диапа зонах определяющих критериев, поскольку каждый из исследо вателей разбивал области Re, Рг, Ре, Gr, Sh и других критериев посвоему. Вдобавок, удовлетворительная воспроизводимость резуль татов экспериментов зачастую достигается лишь в той пористой
среде, для |
которой |
получено данное критериальное |
уравнение |
|
[37, |
38]. |
|
|
|
|
Большой |
вклад в |
обобщение экспериментальных |
данных по |
тепломассообмену в пористых средах внесли Г. А. Витков, Л. П. Хлопанов и С. Н. Шерстнев [45 —48], которые графическим путем сопо ставили приводимые в литературе критериальные уравнения и по лучили усредненные регрессионные зависимости.
Поиск аналогичных обобщений в литературе оказался безус пешным, поэтому сопоставить результаты работ [45-48] оказалось не с чем. В связи с этим ниже приводится краткое изложение результатов исследований Г. А. Виткова, Л. П. Хлопанова и
С.Н. Шерстнева с некоторыми дополнениями и замечаниями.
Вцелях достижения максимальной воспроизводимости ре зультатов эксперимента и упрощения вычислений на практике при анализе процессов тепломассообмена в зернистом слое исполь зуются идеализированные модели. Основные упрощения касаются формы и размеров частиц дисперсного слоя (засыпки). Обычно применяют засыпки, состоящие из частиц правильной геометри ческой формы (шар, цилиндр и т. д.), имеющих одинаковые размеры или определенный диапазон размеров.
Наибольшее число критериальных уравнений тепломассопереноса получено для засыпок, состоящих из сферических частиц, более всего приближенных по форме к дисперсным частицам реаль ных технологических пульп. Зависимости для конвективного тепло- и массопереноса в таких слоях носят сходный характер.
Как известно, тепломассоперенос в дисперсных слоях напрямую связан с гидродинамическими особенностями обтекания элементов структуры слоя. Модель стесненного обтекания сферических час тиц, изложенная в работах Д. Хаппеля [118], Р. Пфеффера [17, 18], К. Касика [13], Д. Карберри [3] и других авторов, легла в основу
теоретического решения задачи тепломассообмена в зернистом слое, образованном из шаров. Однако полученные ими описания весьма приближенно характеризуют процесс. Точного решения этой задачи до настоящего времени не найдено.
Приближенные аналитические решения, предложенные авто рами названных выше работ для описания тепломассообмена в условиях стесненного обтекания частиц, в критериальной форме излагаются ниже. Каждой из приводимых зависимостей соответ ствует свой диапазон изменения эквивалентного критерия Рей нольдса. Включенный в состав регрессионных зависимостей кри терий Прандтля, наряду с гидродинамикой обтекания, позволяет учесть еще и температурный фактор. Ограниченность примени мости приводимых уравнений узкими диапазонами чисел Re, Рг и Ре является следствием сложного влияния большого числа факторов, определяющих закономерности гидродинамики и тепломассопереноса в пористой среде и не учитывающихся полностью и в явном виде в критериальных уравнениях.
При значениях 102 < Res < 103, когда течение жидкости по
каналам соответствует развитому турбулентном режиму, величину критерия Nu, характеризующего интенсивность конвективного
теплообмена, можно определить из уравнений: |
|
|||||
Nu |
= 0,93[е - |
0,75(1 |
- е)(е - |
0,2)]~1/2Pr1/3Re1/2; |
(1.99) |
|
|
|
Nu « Pr1/4Re1/2; |
|
(1.100) |
||
Nu |
= 1.26------ |
[1 - |
(1 - |
е ) 1/3 ] Рг1/3 Re1/3 |
(U 0 1 ) |
|
1 |
e) |
-------- Г75------------- 2 ' |
||||
|
2 - 3(1 - |
+ 3(1 - e) |
- 2(1 - s)2 |
|
Можно заметить, что в области развитого турбулентного тече ния большое влияние на интенсивность теплообмена оказывает показатель порозности слоя, соответствующий доле объема осадка, не занятого твердой фазой. Объясняется это тем, что с умень шением эквивалентного диаметра каналов существенно возрастают скорость и турбулентность потока жидкости.
При меньших скоростях движения жидкости в уравнения необходимо вводить поправку на испарение жидкости с поверх ности непористых частиц, растворение элементов зернистого слоя в потоке различных жидкостей, сублимацию зерен в газовом потоке, нагревание и охлаждение потока или слоя, в то время как пока затель порозности теряет свое первостепенное значение [32]. Для таких условий можно использовать обобщенные выражения, основанные на обширном экспериментальном материале.
В приводимых уравнениях в качестве характерного размера использован эквивалентный диаметр пор йэ. За скорость газа U
принято ее усредненное значение по сечению слоя. При 10“2 < Re < 30
45
44
|
Nu3 = l,2P r1/3 Re3/3; |
( 1. 102) |
|
при 0,6 |
< Рг < 10 и 1 < Re3 < 2 |
|
|
|
Nu3 = 0,515 Рг1/3 Re°'85; |
(1.103) |
|
при 2 < Re3 < 30 |
|
|
|
|
Nu3 = 0,725 P r1/3 Re°'47; |
(1.104) |
|
при 0,6 |
< Рг < 6404 и 0 < Re3 < 5105 |
|
|
|
Nu |
= 0,395Pr1/3Re:j; |
(1.105) |
при 10-2 < Re3 < 2 и 10-2 < Рг-Агэ < 5407 |
|
||
|
Nu3 |
= 0,32(РгАгэ)1/4, |
(1.106) |
где Агэ — эквивалентное значение критерия Аррениуса.
Теоретические решения, полученные Д. Хаппелем, Р. Пфеффером, К. Касиком и Д. Карберри для частиц сферической формы, дают сопоставимые результаты, удовлетворительно согласующиеся сэкспериментом в области переходных чисел Рейнольдса, и подтверждают существование аналогии между процессами переноса
тепла или вещества в слое при Re -» 0, Рг » 1 и Рг » |
1. |
Узость диапазонов Re, Рг и Ре в уравнениях (1.99) |
—(1.106) и при |
менимость только для сферических частиц, слагающих слой, затруд няют использование этих формул при проведении численных расче тов для реальных технологических пульп [87]. В связи с этим очевидна необходимость поиска модели, которая устранила бы указанные недостатки. Попытка создания таковой предпринята Г. А. Витковым, Л. П. Хлопановым и С. Н. Шерстневым [45 —48]. Разработанная ими модель тепломассопереноса в пористых средах построена на безразмерных комплексных величинах, учитывающих относительную значимость известных критериев Архимеда, Грасгофа, Прандтля, Рейнольдса, Шмидта и Шервуда.
Предположим, что среда является псевдогомогенной, в качестве определяющих параметров примем ^&пр (кпр — константа проница
емости, м2), среднюю скорость фильтрации U и особенности внут ренней структуры пористой среды, учитывающиеся коэффици ентом дополнительных сопротивлений турбулентной фильтрации п.
Критериальное уравнение межфазного переноса может быть записано в виде
46
/ |
|
|
|
|
|
I 3/2„ |
п |
t V n p |
t F,k пр |
рак.пр |
(1.107) |
J - * £ - Re 1 + |
Re+ |
i=i______ , |
i=i |
р и |
|
ц |
л[к |
р[/ |
р[/ |
|
где J — отнесенная к единице объема сила, которая для случая переноса импульса равна Др /1 и вызывает перенос тепла или массы в пористой среде; р — вязкость жидкой фазы, кг/м 3; р — дина мический коэффициент вязкости жидкости, Паю; а — удельная поверхность слоя, м2/м 3.
Обозначим
|
I |
пр Z FiKp |
ракпр |
(1.108) |
А = Re 1 + —— Re+ i=i |
1 = 1 __________ |
р U |
||
yfk |
р[/ |
р[/ |
|
|
|
, 3 / 2 |
|
|
|
|
S = J |
^ . |
|
(1-109) |
|
|
Ц |
|
|
Безразмерный комплекс (1.108) имеет смысл обобщенного критерия гидродинамического подобия для процессов конвектив ного переноса в пористых средах, инвариантного по отношению к виду и особенностям структуры пористой среды. В частных случаях (при п = 0) он переходит в критерии подобия Рейнольдса, Ар химеда, Грасгофа и др.
С учетом принятых обозначений на основании (1.107) справед ливо равенство
S = А.
Сила, вызвавшая перенос массы, будет равна
з 3 tcliD |
(1.110) |
J = ( p7/ D) 3 — |
V
где р — усредненный по объему слоя коэффициент массообмена, м/с; 1 — толщина слоя, м; D — коэффициент молекулярной диф фузии, м2/с; V — объем жидкости, находящейся в порах осадка, м3.
Выделим в объеме пористого слоя некоторый единичный объем в виде куба со стороной, равной толщине слоя 1 = ^ к пр , внутри которого осуществляется элементарный перенос массы. Объем жид-
47
кой фазы, находящейся внутри этого элементарного объема, с уче
том коэффициента порозности е будет равен V = в* 372 . Тогда, под ставив в выражение (1.109) значение J из уравнения (1.110), получим
|
Зк |
Dp |
|
|
S = |
I |
( 1. 111) |
|
D |
||
|
И |
|
|
где и/Ар = рг |
диффузионный |
критерий Прандтля |
(Шмидта); |
|
/ D - Nu? — диффузионное число Нуссельта (Шервуда).
В кРитеРиальном виде формула (1,111) может быть записана
следующим образом:
|
|
|
|
|
S = 371 |
Рг |
|
|
( 1. 112) |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
т е ч е ™ |
™ ™ ^ 1^ |
- 5^ ' 107’ И УЧ" ТЫВ“ |
' " го |
|
’Трбуденшог» |
||||
Г.ГЧОЛО |
величина |
к |
- |
32, для процесса массопереноса в пористой |
|||||
среде окончательно получим |
|
|
|
|
|||||
|
1/3 |
|
|
|
п |
пр |
|
пр |
ракпр |
Nu = |
|
Рг |
|
Re |
i=i |
||||
|
|
(=1 |
|||||||
<3пJ |
|
1 + - — Re+ |
|
||||||
9 |
|
|
|
л[к |
HU |
н и |
|
н и |
|
|
|
|
|
|
|
|
|
|
(1.113) |
Ввиду подобия процессов тепло- и массопереноса, крите риальное уравнение теплопереноса может быть записано анало гичным образом, с заменой диффузионных критериев Nu и Рг со
ответствующими им критериями подобия для процессов тепло обменаСогласно этому обобщенное критериальное уравнение меж-
фазного тепломассопереноса в пористых средах можно записать в ЬИДс
|
- ) |
1/3 |
|
NUTW^ Ргга(?) |
.1 /3 |
(1.114) |
|
|
|||
|
.3п ) |
|
|
где Л функция пористости зернистого слоя.
п послеАних слагаемых уравнения (1.113) отвечают наличию внутренних течении в порах осадка, возникающих вследствие
48
каких-то внешних воздействий. Если пренебречь воздействием внешних факторов, то эти слагаемые можно опустить. С учетом этого допущения уравнение (1.113) примет вид
|
|
/ |
Е Л"3 |
|
|
л 1/3 |
Nu |
|
= Рг.1/3 |
Зп |
Re+ |
-Re |
(1.115) |
|
т(д) |
щ(ег)' |
|
5,65 |
|
Выражение (1.115) по форме отвечает тепломассопереносу в изолированной системе в условиях турбулентного течения жидко сти или газа по каналам осадка. Для ламинарного режима (Re < Rexp)
в уравнении (1.115) исчезает поправка на турбулентность потока и уравнение упрощается:
ту |
1/3 |
' £ V/3 •Re 1/3 |
(1.116) |
ш{дг) ~ |
m(g)' |
37Г |
|
Для удобства сопоставления уравнений (1.103) и (1.116) пере пишем уравнение (1.103) в виде
1 |
,1 /3 |
|
•Re173. |
(1.117) |
|
Nura(ff) = U P r |
|
кК
\к /
Нетрудно заметить, что уравнения (1.116) и (1.117) внешне сход ны. Отличие состоит в наличии коэффициента пропорциональности и характерных параметрах в критериях Нуссельта и Рейнольдса. Поскольку величина константы Козени—Кармана кк зависит от
формы, сечения и извилистости пор, а для уравнения (1.103) в работе [32] не приведены диапазон значений пористости и виды засыпок, для сравнения коэффициентов в выражениях (1.116) и (1.117) можно принять средние значения кк = 4,65 и е = 0,39, что
справедливо для стационарных зернистых слоев с неупорядоченной укладкой частиц. В этом случае значения коэффициентов в урав нениях будут равны 0,34 и 0,29 соответственно. Заметим, что раз ница между ними составляет 15 %.
Анализ конвективного теплообмена в слое частиц включает два аспекта: определение коэффициентов теплообмена между теплоно сителями и установление закономерностей распределения темпе ратур в объеме слоя. Здесь основное внимание уделяется оценке интенсивности теплопереноса; вопросы расчета теплоносителей подробно рассматриваются ниже.
Характерной особенностью конвективного теплообмена в слое частиц является высокая интенсивность переноса теплоты между теплоносителями. Она объясняется ранней и довольно сильной
49
турбулизацией газовой фазы из-за сложной формы каналов для прохода газов со случайным чередованием сужений и расширений, поворотов и т. д. В целом, теплообмен в слое — сложное явление, определяющееся многими факторами [84], главными из которых являются следующие:
1) температура на поверхности и в глубине дисперсного слоя определяется не только передачей теплоты от газа к частицам материала (внешний теплообмен), но и отводом теплоты внутрь его (внутренний теплообмен);
2)внешний теплообмен в плотном слое осадка включает пере дачу теплоты от частицы к частице за счет теплопроводности и конвекции;
3)внутренний теплообмен определяется размером и формой частиц, их теплопроводностью и интенсивностью развития тепло вых явлений на их поверхности.
При расчетах теплообмена в неподвижном слое рекомендуется использовать формулы В. Н. Тимофеева, которые позволяют вычис лить средний для всего слоя коэффициент теплоотдачи:
NuM = 0,106Reci при 20 < ReM< 200; |
(1.118) |
NuM= 0,61 R e ^ 7 при ReM > 200. |
(1.119) |
Формулы В. Н. Тимофеева дают хорошие результаты для тер мически тонких частиц правильной сферической или близкой к ней формы [43]. В них отсутствует величина порозности слоя е. Влияние £ на интенсивность теплопереноса от газа к поверхности частиц слоя можно выявить из формулы Р. Ешара:
Nu3 = 2 б / |
(1 - |
s) + Re°'5+ 0,005Re3, |
(1.120) |
где Nu3 = (<хДАг)б/(1 |
- |
в); Re3 = ( U d /v ) /( 1 - |
е); aF — коэф |
фициент конвекции близ поверхности частиц, Вт/(м2 К); Хг— коэф фициент теплопроводности газа (воздуха), Вт/(М'К); vr — кинема
тический коэффициент вязкости газа, м/с. Формула (1.120) спра ведлива для однородного по размерам частиц слоя при Рг = 0,7 и 100 < ReM < 4000. Число Рейнольдса рассчитывается по скорости
фильтрации, теплофизические свойства соответствуют средней тем пературе системы.
Для полидисперсного неподвижного слоя, включающего час тицы фракций d v d2, ..., dn, формула Ешара принимает вид
Nu3 = [2 S cm/ ( 1 - e J ]/(У/У, d/d,.) + Re05 + 0,005Re3. (1.121)
Здесь S CM — средняя порозность смеси частиц (слоя); У. и У — объем соответственно фракции и слоя в целом:
50
НУ, / |
V, |
d1 / d i) = |
|
|
|
|
У, / |
V + (У2 / |
y )(d t / d2)2+...+(Vn / |
V)(d, / |
d j 2 |
(U 22) |
|
Vi / |
У + (V2 |
/ y)dj / d2+.. .+(У„ / |
У )dt / |
dn |
|
п эквивалентный диаметр частиц слоя, входящий в числа Ren и Nun,
d, = [У/lV d,) + V2/(Vd2) + ...]->.
Строгое решение задачи теплообмена в слое при неравно мерном газораспределении может быть обеспечено лишь на основе достоверной информации о поле скоростей газа в агрегате. При этом для вычисления локальных коэффициентов теплоотдачи аРл
(значений aF в данной точке) в условиях неподвижного слоя может быть рекомендована формула В. Б. Щербатского
Ыид = [0,31 - 0,3ехр(-0,0076 Re^)] - Re^+0’66ехр(-° 0175Re-)!, (1.123)
и которой в качестве определяющих параметров при расчете ло кальных чисел Nu3 и Rea приняты эквивалентный диаметр частиц и
локальная скорость газа в слое. Для слоя шероховатых частиц неправильной формы коэффициент п = 0,67, а для слоя шаро образных частиц п = 0,625.
Все указанные формулы характеризуют конвективный тепло обмен в слое термически тонких частиц. В тех случаях, когда на результирующую интенсивность теплопереноса заметно влияет вну тренний теплообмен, следует использовать суммарный коэффи
циент теплопередачи k v определяемый по формуле [103] |
|
kv = [ а '1 + d 2 / [А(1 - s)XH]]~\ |
(1.124) |
и которой коэффициент А = 60 для частиц сферической формы и 75 — для частиц неправильной формы, а км — теплопроводность
вещества частиц.
1.6. Общие закономерности процессов фильтрования
Ранее уже отмечалось, что движение сплошной среды сквозь слой дисперсного материала — достаточно распространенное явле ние в промышленных технологиях, и в первую очередь это отно сится к процессу фильтрования.
Под фильтрованием понимают процесс разделения неодно родных систем (суспензий и аэрозолей) при помощи пористых перегородок, пропускающих сплошную (дисперсионную) среду и задерживающих дисперсную твердую фазу.
51
Аппарат для проведения этого процесса называется фильтром, осевший на перегородке слой твердых частиц с некоторым содер жанием жидкости между ними — осадком, а прошедшая через него жидкость — фильтратом. На практике могут представлять ценность либо оба продукта фильтрования (осадок и фильтрат), либо один из них.
Как указывалось, движущей силой процесса является разность давлений в разделяемой суспензии и за фильтрующей перего родкой. Она расходуется на преодоление гидравлических сопротив
лений слоя осадка и перегородки. |
процесса |
фильтрования, |
где |
||||
На рис. 1.9 представлена |
схема |
||||||
р, > р2 и движущая сила Др = |
р, — р2. |
|
|
|
|
|
|
|
|
|
Фильтрование и филь |
||||
|
|
тры можно |
классифици |
||||
|
|
ровать |
по |
нескольким |
|||
Суспензия |
признакам [106]: |
|
|||||
|
|
|
1) |
по |
действующей |
||
|
|
силе: а) гидростатический |
|||||
Осадок |
напор |
(дренирование), |
|||||
Фильтрующая |
б) |
повышенное |
давление |
||||
перегородка |
перед |
перегородкой |
или |
||||
|
|
вакуум позади нее, в) цен |
|||||
Фильтрат |
тробежная сила; |
|
|
||||
|
|
|
2) |
по целенаправлен |
|||
|
|
ности: получение а) сухо |
|||||
Рис. 1.9. Схема процесса фильтрования |
го |
осадка, |
б) |
чистого |
|||
|
|
фильтрата, в) сухого осад |
|||||
|
|
ка |
и |
чистого |
фильтрата |
||
|
|
одновременно; |
|
|
3)по режиму фильтрования: а) непрерывное, б) периодическое (последнее, в свою очередь, может протекать с постоянной ско ростью, при постоянном давлении или при изменении обеих величин);
4)по механизму процесса фильтрования: а) с образованием осадка на поверхности фильтрующей перегородки, б) с закупоркой пор фильтрующей перегородки;
5)по характеристикам осадка: а) сжимаемый и б) несжи маемый.
Фильтрование используется для разделения суспензий, содер жащ их твердые частицы, крупность которых не превышает размеров пор фильтрующей перегородки. Современное фильтра ционное оборудование способно отделять частицы крупностью от единиц микрометров до нескольких миллиметров.
Фильтрование сопровождается химическим и физико-хими ческим взаимодействием между жидкой средой, дисперсной фазой, перегородкой и корпусом фильтра. В подавляющем большинстве случаев это взаимодействие приводит к нежелательным последст-
52
пням, например быстрому росту гидравлического сопротивления, разрушению перегородки, нарушению целостности фильтра и пр. Поэтому структура осадка и перегородки, а также свойства жид кости играют важную роль при фильтровании. От этих факторов зависит в итоге и режим процесса [4, 58].
Общая теория фильтрования описывает лишь те закономер ности, которые свойственны всем перечисленным в классификации типам и режимам протекания процесса. Ее основой является эм пирический закон Дарси, открытый на основе изучения фильтрации жидкости сквозь естественные грунты еще в 1852—1855 гг. Со гласно этому закону объем фильтрата, проходящего через единицу поверхности фильтра за единицу времени, прямо пропорционален разности давлений и обратно пропорционален общему сопротив лению осадка и фильтрующей перегородки. Дарси установил, что при постоянном сопротивлении грунта расход жидкости, прохо дящей через поперечное сечение слоя, прямо пропорционален гидравлическому уклону, равному потере напора по высоте слоя, и обратно пропорционален сопротивлению грунта:
Q = |
FAH |
(1.125) |
---------= Fi / RCJl, |
Л А
где О — объемный расход фильтрата, м3/с; АН — потери напора в слое, м; Лсл — высота слоя, м; Ra — сопротивление слоя грунта,
м^е/м3; F — площадь фильтрации, м2; i — гидравлический уклон. При рассмотрении закономерностей общей теории фильтро
вания вводится допущение, согласно которому поры осадка имеют цилиндрическую форму, причем диаметры всех цилиндров одина ковы, а их длина равна толщине фильтрующей перегородки (так называемый идеальный фильтр). При переходе к реальным случаям вводят поправки (например, при помощи эмпирических коэффи циентов фильтрования).
На основании закона Дарси и в предположении, что значение морозности фильтрующей перегородки постоянно по сечению, уравнение кинетики фильтрования, согласно (1.83), можно пред ставить в виде [57]
d V |
е3 |
Др |
Fdx |
k K{1 - е ) 2 |
(1.126) |
цжЛ0 |
где V — объем фильтрата, м3; т — продолжительность фильт рования, с.
В действительности уравнение кинетики фильтрования имеет б< >лее сложную форму, поскольку величина г может изменяться как но сечению слоя, так и во времени, что обусловлено непостоянством скорости движения фильтрата и зависимостью порозности от
53
состояния поверхности, на которой происходит образование осадка. Влияние скорости фильтрования особенно значимо на начальном этапе формирования осадка, поскольку в этот период мелкие час тицы суспензии проникают в поры перегородки, вызывая ее за купоривание. Забивка пор фильтрующей перегородки может вы зывать резкий скачок ее сопротивления, что негативно отражается на всех последующих стадиях разделения и снижает регенера тивную способность фильтроткани.
Характеристики осадка сказываются на эффективности филь трования в неменьшей степени. Как уже отмечалось, основным влияющим фактором является показатель сжимаемости осадка. Со противление несжимаемых осадков неизменно во времени вследст вие постоянства порозности и незначительно зависит от перепада давлений и скорости осаждения. У сжимаемых осадков по мере роста давления фильтрования происходит уменьшение порозности и, как следствие, наблюдается рост сопротивления движению фильтрата. Отметим также, что у сжимаемых осадков, в отличие от несжимаемых, зависимость между порозностью и давлением филь трования носит нелинейный характер, а у сильно сжимающихся осадков (гидроксидов железа и меди, волокнистых, высокопористых материалов и др.) увеличение Ар сверх некоторого критического значения приводит к уменьшению скорости фильтрования (рис. 1.10). Рисунок демонстрирует наличие некой оптимальной величины движущей силы, при которой производительность фильтра О = AV/Дт в случае сжимаемого осадка максимальна.
При фильтровании с закупориванием пор образования осадка не происходит, а дисперсные частицы закупоривают поры фильтру ющей перегородки. Для реализации этого режима необходимо, чтобы
крупность |
твердых частиц |
не превышала |
диаметра |
межзерновых |
|||||||
|
|
|
каналов. |
|
Следствием |
этого |
|||||
|
|
|
будет |
проскок некоторой |
|||||||
|
|
|
части |
дисперсной |
фазы |
в |
|||||
|
|
|
фильтрат, |
что |
приведет |
к |
|||||
|
|
|
повышению его мутности. |
|
|||||||
|
|
|
|
Предельную |
задерж и |
||||||
|
|
|
вающую |
способность пере |
|||||||
|
|
|
городка |
приобретает |
после |
||||||
|
|
|
непродолжительной |
работы |
|||||||
|
|
|
в |
результате |
уменьшения |
||||||
|
|
|
эффективного |
размера пор |
|||||||
|
|
|
осевшими |
в |
них |
мелкими |
|||||
|
|
|
частицами или из-за обра |
||||||||
|
|
|
зования |
сводиков |
над вхо |
||||||
Рис. 1.10. Влияние давления фильтро- |
дом в поры. Поэтому мут- |
||||||||||
ность |
фильтрата |
в |
начале |
||||||||
вания на |
производительность |
фильтра |
процесса еще не свидетель- |
||||||||
для несжимаемого 1 и сжимаемого 2 |
ствуют о непригодности ис- |
||||||||||
осадков |
|
|
пытываемой перегородки. |
|
Если подобная картина наблюдается при фильтровании с образованием осадка, то для того, чтобы исключить появление мутного фильтрата в начале процесса, при снятии осадка оставляют его тонкий слой на перегородке или начинают фильтрование при низком перепаде давлений.
Изменение давления при фильтровании не всегда целесооб разно, так как это требует использования дополнительных регу ляторов, в широких пределах варьирующих мощность насосных установок. В связи с этим наиболее часто используемым режимом протекания процесса является фильтрование с образованием осадка, проводимое при постоянном перепаде давлений.
Для определения типа процесса фильтрования при постоянном давлении устанавливают линейные зависимости между различными переменными, которые могут быть измерены в реальном процессе. Ниже приведены линейные зависимости для различных типов филь трования [81]:
с полным закупориванием пор — q—и;
с постепенным закупориванием пор — х—x/q; промежуточный — т— 1/и;
с образованием осадка — V—т/V,
здесь q — удельная производительность фильтра, т. е. объем филь трата, получаемый с единицы поверхности фильтрующей перего родки, м3/м 2; и — объем осадка, отложившегося на фильтрующей пе регородке при прохождении единицы объема фильтрата, м3/м 3, т — продолжительность фильтрования, с; V — объем фильтрата, м3.
Во всех случаях процесс фильтрования стремятся организовать гак, чтобы он протекал с образованием осадка и по возможности исключалась закупорка пор фильтрующей перегородки. Для этого в последнее время широко используются специальные вспомога тельные вещества, образующие на перегородке осадок, улавлива ющий мелкие частицы, что особенно эффективно при фильтрова нии разбавленных суспензий.
Сопротивление фильтрования R в общем случае можно пред ставить в виде суммы сопротивлений осадка Roc и фильтрующей
перегородки Rn : R = Roc + Rn. |
|
Сопротивление осадка пропорционально его толщине h |
: |
Кс = Ч с |
° С (!-127) |
где Roc — сопротивление осадка при фильтровании под действием
разности давлений, Н'С/м3; г — удельное сопротивление осадка, 11-с/м4, т. в. сопротивление осадка высотой 1 м, отложенного на площади 1 м2.
Объем осадка, отложившегося на фильтре, можно выразить как произведение площади F фильтра на толщину hoc осадка или через
объем V прошедшего фильтрата. Тогда объем осадка будет равен uV. ( Следовательно,
54 |
55 |
Fh^ = uV, |
(1.128) |
откуда толщина осадочного слоя составит |
|
ho c = u — = uq. |
(1.129) |
F |
|
Подставив значение Лос в уравнение (1.127), найдем сопро |
|
тивление осадка |
|
Я с = ruq |
(1.130) |
и получим следующее выражение длясопротивления фильт рованию:
R = |
Яос+ |
Rn = ruq+ Rn. |
(1.131) |
|
Отсюда скорость движения жидкости в расчете на свободное |
||||
пространство аппарата составит |
|
|
|
|
гг |
v |
Я |
Ар |
(1.132) |
U= — = -2- |
= - ^ . |
|||
|
Fx |
х |
R |
|
Как видно из уравнений (1.131) и (1.132), сопротивление R по мере образования осадка и увеличения его толщины возрастает, а скорость фильтрования уменьшается. Перепишем уравнение (1.132)
в дифференциальной форме с учетом (1.131). Тогда |
|
|
dq |
Ар |
(1.133) |
U = |
= ---£---- , |
|
dx |
ruq + Rn |
|
или |
|
|
ruq + R„ |
|
|
dx |
= ------- d q . |
|
Ap
Интегрируя это выражение в пределах от 0 до т и от 0 до q, определим продолжительность фильтрования:
Ч
x = J ШЯ , К ' dq.
о Ч Ар Ар
После интегрирования получим
1 = п |
+ M |
(1,134) |
2Ар Ар
56
Решив уравнение относительно q, определим удельную произ водительность фильтра:
н |
2Ар т |
ч |
(1.135) |
Vги ) |
ги |
уГи |
|
Из уравнений (1.134) и (1.135) следует, что сопротивление, встречаемое потоком фильтрата, по мере накопления осадка воз растает, поэтому постоянство этого потока во времени (следова тельно, и максимальная производительность фильтра) может поддерживаться лишь при непрерывном увеличении разности давлений.
Такой рабочий режим обеспечивается путем нагнетания сус пензии поршневым насосом. При использовании сжатого газа и вакуумирования Ар = const и с ростом высоты слоя скорость дви жения фильтрата снижается, а производительность фильтра падает. Если суспензия подается центробежным насосом, то в пределах его рабочей характеристики по мере нарастания слоя осадка проис ходит увеличение Ар, которое сопровождается уменьшением потока фильтрата.
1.7. Фильтрование с образованием несжимаемого осадка
Как уже отмечалось, основным свойством несжимаемых осадков является постоянство их пористости и сопротивления потоку движущейся по каналам жидкости в течение всего процесса фильтрования. Данное утверждение справедливо и для несжима емых фильтрующих перегородок (например металлических, порис тых керамических или стеклянных), но численное значение их удельного сопротивления может существенно отличаться от тако вого для осадка. Раздельный же учет Rx и Rn сильно усложняет
модель. Для упрощения расчетов допускают, что перегородка имеет небольшую толщину и ее сопротивление потоку фильтрата пре небрежимо мало [94].
В таком случае в уравнении (1.126) можно объединить посто янные величины, характеризующие свойства данного осадка, и с
учетом R * Roc получить |
|
|
R = М 1з~2£)2 • |
I1' 136» |
|
s d |
Ч |
|
После соответствующего преобразования уравнение (1.126) при мет вид
57
d V |
А р |
(1.137) |
Fdx |
= -----*— . |
|
R px hoc |
|
Выразим толщину отложившегося осадка hoc через объем про
шедшего фильтрата, учитывая пропорциональность объемов осадка и фильтрата.
Связь между массовой долей дисперсной фазы в исходной сус пензии ам и толщиной осадка описывается уравнением ам = hocFpm-
•(1 — s). Массовое отношение Т : Ж можно выразить как
ам _ KcFpmi1-*
1 ~ ам РЖ(У + *РКс)
где Урж — масса фильтрата; hocFzрж — масса жидкости, задержав
шейся в порах осадка. Используя последнее уравнение, определим высоту осадка:
h = . °мУРж
F[ ( l - E ) ( l - a J p m - а мерж]
иего объем после прохождения единицы объема фильтрата:
и= --------------Р£«---------------.
(1 - е )(1 - х )р ш - а мерж
Окончательно связь между толщиной осадочного слоя и объ емом фильтрата выражается зависимостью
Л= Vu/F.
ОС
На основании полученного соотношения уравнение кинетики фильтрования приведем к виду
^ = |
(1.138) |
dx |
Rpx uV |
Дифференциальное уравнение (1.138) можно проинтегрировать для двух режимов фильтрования (при постоянном давлении и при постоянной скорости).
В результате интегрирования (1.138) при Ар = const получим
т |
_ |
КРЖЦ у |
(1.139) |
V |
|
2АpF 2 |
|
Таким образом, фильтрованию при постоянном перепаде дав лений отвечает линейная зависимость между отношением х /7 и V.
Итогом интегрирования выражения (1.137) при постоянной скорости будет
тДржп
— = |
, V. |
(1.140) |
V |
ApF |
|
Последнее уравнение свидетельствует о наличии прямо про порциональной зависимости между Ар и V.
Если фильтрующая перегородка не отвечает условию, приня тому в начале раздела, и ее сопротивлением потоку фильтрата пренебречь нельзя, выражение (1.126) необходимо переписать в общем виде с раздельным учетом сопротивлений осадка и пере городки:
= ____ *£_____ . |
[1.141) |
№+ « ,)
Следует учитывать, что под сопротивлением перегородки по нимаются не только помехи, создаваемые ее материалом, но и препятствие со стороны тонкого слоя осадка, оставшегося на ней после предыдущего цикла. При фильтровании суспензии с неиз менными во времени свойствами на конкретной фильтрующей перегородке в условиях устоявшегося процесса сопротивление Rn
можно считать постоянным. Для того чтобы проинтегрировать уравнение (1.137) и выявить закономерности кинетики процесса с учетом Rn, необходимо видоизменить его, учтя толщину фильтру
ющей перегородки Лп и допустив приблизительное равенство удельных сопротивлений осадка и перегородки:
d V |
Ар |
(1.142)
Fdx гцж(Л0С+Л п)
В соответствии с уравнением (1.138) запишем
d V |
ApF |
ApF2 |
|
Vu и |
(1.143) |
|
um |
|
|
— + h |
|
|
v F |
/ |
Полученное выражение называется основным уравнением процесса фильтрования. Путем его интегрирования определяется связь между объемом получающегося фильтрата и продолжи тельностью процесса, учитывающая характеристики суспензии (рж),
58 |
59 |