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

книги / Физико-химические основы технологии дисперсных систем и материалов

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
9.33 Mб
Скачать

ше, представлением о «дырках» в укладке частиц. Движение такой «дырки» происходит со скоростью Vh, которую можно оп­ ределить из уравнения:

Зя

гк г\а!гD20

F,

( 1. 112)

8 '

2h

В левой части этого уравнения записано выражение для силы

тяжести, действующей на частицу,

а в правой — сумма сил

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

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

со стороны ее соседей 2К/2. Здесь использована формула [20], характеризующая силу сопротивления относительному сближе­ нию сфер со скоростью vh в предположении, что при переходе частицы в нижний вакантный узел такая сила действует на нее со стороны половины ее соседей (гк/2).

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

фмакс ф

3/t

(P — Po) gha

фмакс

vh ' Dn■Vh' 3zK

П

 

 

( 1 .1 1 3 )

где h — среднее расстояние между поверхностями.

При записи соотношения (1.113) была использована формула (1.112), которая позволяет оценить скорость движения «дырки»:

(Р— Po) g D \h

Vh'

Згк

Т|

( 1 . П 4 )

 

 

При этом было учтено, что среднее расстояние между поверх­ ностями частиц составляет:

Ц) [ 1— (ф/фмакс)1/3].

(1.115)

Формулы (1.114) и (1.115) описывают седиментацию суспензии при Fg^> UJhK.

Далее рассмотрим процесс седиментации суспензии при условии Fg^tiUJhK. В этом случае необходимо учитывать воз­ можность объединения частиц в конечные агрегаты диаметром

£> = (Ш0

(где

— параметр

структуры,

который характеризует

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

частиц в суспензии).

Для описания

седиментации

суспензии

при

указанном

условии

можно

вос­

пользоваться

формулами

(1.111) — (1.115), полагая

в них za, фа,

фа,макс»

ра»

ра,0» гц» D и

d вместо гк, <р, фмакс,

р, po, t],

и А

(где ра

и

p3i о— плотности

агрегатов

и

суспензии

вокруг

них,

Фа, Фмакс и 2а — плотность

упаковки

агрегатов,

ее

предельное

значение и координационное

число).

Величина

ра

и ра>0 выра­

жаются через объемное содержание частиц твердой фазы в аг­ регате фа и в окружающей суспензии ф:

Ра = РФа + Po (1 фа) J Ра,о = РФ+ Po (1 — ф) •

41

Для качественных оценок можно принять

Р — Ро Ра — Ра,0« Р ~ Ро ~ Ра ~ Ра,о « 1 Кг/м3,

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

Найдем среднее расстояние между поверхностями агрегатов d , полагая, что d - ^ D . При движении агрегатов относительно друг друга со скоростью, равной скорости движения «дырки» Vh, наибольшее давление среды на их поверхности можно оце­ нить по формуле [20]:

 

р ^3/8^vhD/d\

 

(1.116)

Величину d можно определить

из

сравнения р с прочно­

стью агрегата

 

гк

 

 

 

Ркр

4Е/а

 

А*

(1.117)

я/IKD2о

2

6JT.D0/IV

 

 

Подстановка полученного выражения для d в (1.113) дает оцен­ ку скорости оседания агрегатов:

8

(Ра-ра,о)3£3Д4

(1.118)

27z3a

Pi Лер

 

Для концентрированной

суспензии характеристикой седи­

ментации является скорость роста объемного содержания твер­ дой фазы ф в нижней части сосуда. Зная эту величину, можно определить скорость роста высоты осадка osi. В концентри­ рованной дезагрегированной суспензии ysi~Os, а в сильно аг­

регированной nsi~H sa. Сравнивая (1.113)

и (1.118),

видим, что

в дезагрегированной суспензии osi пропорциональна

а в силь­

но агрегированной usi

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

g 3. Рост

отношения

vslfg с увеличением g

в последнем случае отражает факт уско­

рения седиментации при разрушении агрегатов. При dKD зависимость меняется — дальнейшее уменьшение агрегатов уже замедляет седиментацию. Это можно учесть, если рассмотреть седиментацию отдельного агрегата. Максимальное сдвиговое напряжение на его поверхности при движении агрегата со ско­ ростью

Г8а

1

 

(Ра — Ра,о) gPa

18

'

ТЦ

 

оценивается по формуле [21]:

3

p ^ — ^vJD .

(1.119)

(1.120)

Подставляя (1.119) в (1.120) и сравнивая Р с РкР, можно оценить размер агрегата D. Используя (1.119), получим окон-

42

Рис. 1.4. Зависимость безразмер­ ного параметра, характеризующе­ го скорость седиментации в цент­ рифуге суспензии полистирольного латекса, от относительного ус­ корения.

Точки — экспериментальная

зависи­

мость;

прямая — расчет

по

формулам

(1.111),

(1.118) и (1.121)

при

3«100

нательное выражение для оценки скорости оседания агрегатов:

vsa------- Г 8^

- -

(М21)

 

'll (Ра Ра,о) ё

 

 

В [41] экспериментально

исследована

зависимость

скорости

оседания агрегированной

суспензии

полистирольного

латекса

(D0« l мкм, ф«22% ) от ускорения g,

обусловленного центро­

бежной силой инерции центрифуги. Величину g изменяли в пре­ делах от 1 до 20go (где g0~9,8 м/с2 — ускорение свободного падения в поле тяжести Земли). По данным реологических из­ мерений [41], для этой системы гцлФ.ОЗ Па-с, Ркр?»0,5 Н/м2.

На рис. 1.4 представлена зависимость безразмерного пара­ метра gvs(go)lg0vs(g), характеризующего скорость седимента­ ции данной модельной системы, от относительного ее ускорения g/go в центрифуге. При полном разрушении агрегатов (при больших значениях g) этот безразмерный параметр сохраняет постоянное значение, что соответствует (1.111). При меньших значениях g зависимость проходит через минимум, чему соот­ ветствует наибольшая скорость седиментации. В этом случае агрегаты еще достаточно велики и способствуют ускорению се­ диментации, но уже малы и не мешают взаимному движению. Наконец, при дальнейшем снижении g седиментация замедля­ ется, так как уменьшается d, т. е. структура «сдерживает» се­ диментацию.

Как отмечалось в [41], существует тесная связь седиментационных и реологических характеристик агрегированных си* стем. Подробнее эта связь будет проанализирована ниже при сопоставлении данных, приведенных на рис. 1.4, с реологиче­ ской кривой для исследуемой системы. Здесь же дается каче­ ственный анализ влияния сдвигового напряжения на седимен­ тацию агрегированных систем.

В сдвиговом потоке суспензии размер агрегатов определяет­

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

43

седиментации с е. Для этого следует воспользоваться результа­ том, приведенном в разд. II.1:

D / d t t ] / r гкР Кр/(гце) .

(1.122)

Рассуждая далее так же, как при выводе формулы (1.118) и ис­ пользуя (1.122), получим:

(Ра

Ра,о) £ 3 Eg

(1.123)

 

 

кр

Аналогично можно показать, что

vsa ~ D3 (Ра — Ра,о) g P 2

(1.124)

%р2кр

 

где Р — напряжение сдвига.

Соотношения (1.123) и (1.124) показывают, как с ростом градиента скорости сдвига и напряжения сдвига по мере углуб­ ления разрушения структуры агрегированной системы ускоря­ ется седиментация. Это происходит, однако, пока d<^D, т. е.

при Б<^БКр1= 2кРкр/т11. При дальнейшем увеличении градиен­

та б разрушение агрегатов приводит уже к замедлению седи­

ментации— зависимость vsa от б

имеет максимум, аналогично

рассмотренной выше зависимости

osa от g. Поэтому наиболее

интенсивная седиментация происходит при значениях градиента,

близких к критическому б Кр,ь при котором происходит раз­ рушение агрегатов. Полученные выше выводы можно исполь­ зовать также для качественного анализа седиментации в полидисперсных системах. Пусть система содержит частицы разме­ ром D0 в пределах Di<t.D0<^D2 . При течении системы со ско­

ростью сдвига г наиболее ускорится процесс оседания

частиц

размером DKР. Величина DKP определяется из уравнения:

 

екр,1 (Д ф )

~ гкр кр (£*Кр)/f\i.

(1.125)

Такие частицы объединены

в агрегаты размером D ^ $ D o ^D o -

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

данном градиенте б не смогут объединяться в агрегаты, и по­ этому скорость их оседания будет меньше.

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

зовавшемся через 1 ч

после начала течения с градиентом б =

= 10 с-1. Эксперимент

проводили в ротационном вискозиметре

44

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

Скорость сдвига 8—10 с~!, продолжительность те­ чения 1 ч; с — массовое содержание фракции

40 —

30 -

 

•ч«5>

 

 

 

 

СЗ'"-

 

 

 

с коаксиальными цилиндрами,

tS

20

 

 

че­

 

 

 

рез который пропускали исследуе­

 

 

 

мую суспензию с массовым содер^

10

 

 

жанием твердой фазы «67% . Из

 

I

_____I______ L

графика видно, что по сравнению с

 

исходной суспензией в осадке по­

0

во

то ZW

вышено содержание фракций

ча­

 

 

Ла,мим

стиц размером 25-^90 мкм. Именно частицы этих промежуточных фракций особенно ускоренно се-

диментируют в сдвиговом потоке. В то же время содержание частиц фракций 90-У250 мкм в осадке понижено.

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

работку таких систем. Это обстоятельство

следует

учитывать

при проведении различных химико-технологических

процессов

в дисперсных системах.

процессов, наоборот, целесообразно

Вместе с тем в ряде

ускорить седиментацию

(«высаживание»)

частиц из полидис-

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

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

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

модели безразлично, колеблется ли весь

сосуд или

какая-то

его стенка, можно полагать, что вибрация

вызывает

хаотиче­

ское движение частиц. В результате устанавливается

динами­

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

45

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

Предположим, что функция распределения частиц по скоро­ стям имеет вид:

f (vp) ~ г’2р/”2р ■ехр ( — у%/~2р ),

(1.126)

где vp и vp — скорость движения частицы и ее среднее значение*.

Это предположение дает основание использовать величину Q=Mvp2/2 в качестве эффективной температуры, так как она играет ту же роль, что величина k^T в кинетической теории га­ зов. Условно величину 0 можно назвать вибротемпературой.

Определим число частиц в единице объема n(z) на

высоте

z, отсчитываемой от дна сосуда, в котором находится система.

К данной системе, рассматриваемой как некоторый газ,

приме­

нимо уравнение состояния

Ван-дер-Ваальса

[43, 44]:

Р +

.

«0

(1.127)

Ап-> = - -----— .

^

 

1 — Вп

 

Параметры А и В для данной системы частиц определяются по формулам:

оо

А = л j

U (х) х2 dx ~ Ua D2о

В — я7)3о/(6срмаКС) ,

(1.128)

где U(х) — энергия

взаимодействия частиц как функция расстояния х между

их центрами.

 

 

 

Формула для определения параметра А применима при 03> ^>Ua. Как известно, в газе, находящемся в гравитационном по­ ле, устанавливается распределение частиц по высоте n(z). Для идеального газа (Л—Л), В—»-0) оно имеет вид [44]:

п (г) = пх ехр ( — Mgz/Q),

(1.129)

где tti — нормировочная постоянная.

Чтобы найти n(z) для реального газа, запишем химический потенциал газа:

Р (г) = Ро +

0 In (р + Ап2) + - — —— ^----- ------

— Ап + Mgz,

(1.130)

 

(TinD3о/6фМакс) — 1

 

 

где р0 — постоянная.

 

 

 

* Усреднение при вычислении

vp можно проводить как по времени, так

и по фазовому

пространству (т.

е. по скоростям

и координатам

частиц).

Для системы твердых сфер, которая здесь рассматривается, оба способа ус­ реднения эквивалентны в соответствии с эргодической теоремой [42].

46

Из условия постоянства по высоте z получим:

п

1 — я п х D 30/6cp

 

1

1п

1 — я л D 30/6cp _

1 — JtttD30/6cpMaKc

пх

 

1

Mgz

п пг

1

Я7Ц 7)3/6фмакс

0

( 1 . 1 3 1 )

0

Учитывая, что 0 М /а и принимая

во внимание (1.128), послед­

ним слагаемым в (1.131) можно пренебречь. Тогда приближен­

ное решение (1.131)

будет иметь вид:

 

 

 

 

 

 

 

 

 

 

 

п = пг при

z<^Ls

 

 

 

(1.132)

 

 

 

п ях nj exp [ — Mg (z Ls) /0]

при z > L s

(1.133)

Величину Ls можно

найти

в результате

сравнения

второго и

третьего членов (1.131):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

[ 1

 

7)3о/бфмакс]

1 >

 

 

 

 

 

£ * ~ [ 0 / № ) ] П - я п ^ У б ф м а к с ]- !.

 

(1.134)

Величина П\ определяется из условия нормировки:

 

 

 

 

j

п (г) dz ^

пг Ls +

 

0nj/(Mg)

= N,

 

( 1 . 1 3 5 )

где

N — число

частиц,

приходящихся

на единицу

площади горизонтальной

плоскости.

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим через Lsi= n/6D03N высоту

слоя

частиц «газа»

при

максимально

плотной

упаковке,

т. е. при

«конденсации».

Из

(1.134)

и

(1.135)

получим:

 

 

 

 

 

 

 

 

/

я ^ Д 3,,

у

_

/

MgLsl \

т г Д30

 

MgLsl =

Q

 

\

бф м акс

/

 

\

9

/

бф м акс

 

0

 

( 1 . 1 3 6 )

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

яях Р 3р

_

 

Ф

_j

44gLsl

__

[

j

, / 7VfpLSi

( 1 . 1 3 7 )

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

бф м акс

 

фмакс

2 0

 

 

 

\

2 0

/

Подставляя в (1.127) вместо Вп величину nrtiDo3/6(pMaKc, полу­ чим:

( 1 . 1 3 8 )

Формула (1.138) справедлива как при MgLsi^>0, так и при

MgLsi <&в.

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

47

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

цессы в системах, содержащих

твердую

фазу

как

в жидкой,

так и газообразной среде [15]. В последнем случае

(порошки)

применение вибрации наиболее эффективно,

под

действием

вибрации порошки переходят

в состояние,

называемое вибро­

ожижением. В этом состоянии порошок ведет себя подобно жид­ кости или плотному газу. При более интенсивных воздействиях порошок переходит в состояние виброкипения [45, 46].

В рамках изложенной выше модели виброожиженный и виб­

рокипящий слой порошка можно рассматривать

как плотный

газ,

находящийся

под действием

гравитационного поля,

кото­

рый

описывается

уравнением газового

состояния

(1.138),

при­

чем объемное содержание ф твердой

фазы в нижней плотной

части слоя отвечает соотношению (1.137). Резкое

снижение ф,

т. е. переход в состояние виброкипения, происходит при

 

 

 

йр » ] / " M^ S1

= VgLTi.

(I • 139)

Такое условие предлагалось использовать как критерий пере­ хода системы к виброкипению в [47].

Рассмотренные выше методы могут быть использованы и для описания динамического состояния при'любом другом способе передачи импульса частицам. Все отличие сводится лишь к способу вычисления vp. Например, для слоя порошка, псевдо­ ожиженного потоком газа [48]:

Vp ~

kf (U/Uf — l)!/2,

(1.140)

где и — скорость передачи газа;

щ — значение и, соответствующее

переходу

к псевдоожижению; k; — эмпирическая постоянная.

 

Оценку vp при вибрационном воздействии можно провести, рассматривая движение отдельной выделенной частицы между двумя соседними, которые совершают гармонические колебания с круговой частотой со и амплитудой а. При этом имеют места преимущественно «лобовые» столкновения частиц и движение частиц прямолинейно. При очередном столкновении с соседней частицей изменение кинетической энергии составляет:

АЕк = М [ (vp (1 — е) -(- асо cos соt)2— у2р ]

(I . 141)

где е — коэффициент неупругости удара.

Усредняя (1.141) по времени и приравнивая полученную ве­ личину энергии Д-Едис, диссипируемой при столкновении двух

48

частиц, определим vP. Величину Л£ДИс можно оценить, исполь­ зуя (1.11):

 

h

 

 

 

 

Л£дис~2 | F0(A) dJi «

3/4itr|D»0Tp In (А/Аннн).

(1.142)

Для

^МИИ

 

 

 

получим:

 

 

 

 

9т) "In (Й/Амии) +

 

9т) ■In (А/Аииа).

(1.143)

 

v \ 2ар£>0

2a

2apD0

 

где a = 1—(1—e)2

Хотя ф и, следовательно, h зависят от скорости частиц, од­ нако эта зависимость в (1.143) логарифмически слабая, и мож­ но принять h t t D 0. Тогда представляется возможным выделить два способа определения vp в зависимости от асо:

- ^

о р Д ,

(а м )2

пр и

у - м ^

ti In

( А / А ш « )

>> a a ffl

( , _ 1 44)

 

Т1

 

 

 

 

 

 

 

 

 

Vv ~

47=

ПРИ

Vaaco>—~

In (Л/Акр)

 

(1.145)

 

 

У

а

 

ры0

 

 

 

В (1.144)

учтено,

что вычисление

по

(1.142)

и все

последу­

ющие расчеты выполнимы при

ц/(рА>) 'п (7Г/Лмин) С vp

Ясно, что (1.144) имеет смысл лишь при а<10~2. Для порошка, содержащего частицы размером Д )~ 100 мкм в воздушной сре­ де ( п » Ы 0 -5 Па-с) скорость частиц определяется по формуле (1.145). На рис. 1.6 показан результат эксперимента, проведен­

ного

с порошком электрокорунда:

П0^120 мкм, высота слоя

4 см.

Ширина сосуда, в котором

проводились испытания, со­

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

Рис. 1.6. Зависимость скорости вибрации асо, соответствующей переходу порошка электрокорунда (£>о=120 мкм) в состояние вибро­ кипения, от частоты вибрации.

Точки — экспериментальная

зависи­

мость;

прямая — расчет по

формуле

(1.145)

при а~0,25

 

С О .Г и ,

4— 1493

49

 

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

действие на непрерывную

структуру,

как

будет показано в

разд. II. 1, может привести к возникновению

в ней

разрывов

сплошности. Интенсивная

вибрация,

под действием

которой

порошок переходит в состояние виброкипения, если ускорение асо2 больше g, сосуд достаточно широк и его дно перфорирова­ но, приведет к «подбрасыванию» слоя, что нарушит его одно­ родность. Аналогично при псевдоожижении газом слоя порошка в нем возникают пузыри, если скорость газа достаточно велика (в 2—3 раза больше минимальной скорости псевдоожижения). Следует отметить, что при анализе не учитывалась возможность проявления разнообразных волновых процессов под воздейст­ вием вибрации, которые возникают, если длина звуковой волны с!со — скорость звука) сравнима с размером емкости, в кото­ рой находится порошок.

ГЛАВА II

ТЕОРИЯ ТЕЧЕНИЯ СТРУКТУРИРОВАННЫХ ДИСПЕРСНЫХ СИСТЕМ

11.1. АНАЛИЗ ПОЛНОЙ РЕОЛОГИЧЕСКОЙ КРИВОЙ КОНЦЕНТРИРОВАННЫХ ДИСПЕРСНЫХ СИСТЕМ

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

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

50

Соседние файлы в папке книги