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

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

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

В самом общем случае сила сопротивления, проявляющаяся в движущемся потоке, выражается уравнением (4-15), которое мо­ жет быть представлено в следующем виде:

где f — площадь проекции частицы на плоскость, перпендикуляр­

ную направлению

его

движения; р — плотность среды;

w0— ско­

рость

среды вдали

от частицы; С — коэффициент пропорциональ­

ности,

равный общему

коэффициенту сопротивления.

Уравнение

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

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

Р к = /?л. с “Ь ^ т р

(4-35)

то и общий коэффициент сопротивления £ может быть выражен уравнением

£ — £л. с + £тр

(4-36)

где £л. с — коэффициент лобового сопротивления, часто

называе­

мый коэффициентом формы; £тр— коэффициент трения.

потоком

При ламинарном течении частица плавно обтекается

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

Взаимодействие сил на поверхности шарообразной частицы, вы­

ражающееся зависимостью (4-35),

можно представить также урав­

нением

 

 

FK = F, +

F2- F 3

 

 

(4-37)

где

 

 

 

 

 

Л

 

 

 

 

 

 

 

 

 

 

F) =

J

J (— р |г_ £ cos Ѳ) R2 sin Ѳ dQ dq>

 

(4-38)

2Я Я

 

о о

 

 

 

 

 

 

 

dwr

 

 

 

F , - i ;

 

 

 

sin Ѳ

R2 sin Ѳ d0

dq> (4-39)

 

 

 

 

о 0

 

 

 

(ЭѲ

r=R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F3= I

 

JЛ [— (po — Pgz) |r_ ^ cos Ѳ] R2 sin Ѳ de d<p

(4-40)

0

0

 

 

 

 

 

117

В уравнениях (4-38) —(4-40) ро —давление на плоскости z — 0,

проходящей через экватор

шарообразной

частицы,

величина

R2sin ѲіШср — поверхность частицы.

(4-15) и

уравнений

При использовании общей

зависимости

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

р — р0 +

Pgz

 

W

 

Г

Р®о

;

w0

Wr =

 

w0

L =

Г

Re = 2w0Rp

(4-41)

 

R

P

 

Тогда получим:

 

± ._ L

я

 

Г ѳ

dWT-

 

£тр

J L

 

я Re о

оЯ dL

L

дѲ l

 

 

 

 

 

л

 

 

 

 

 

£л.с=—J J

 

(— Р cos Ѳ) sin ѲdQ dtp

(4-43)

 

 

о

о

 

 

 

Таким образом, для того чтобы рассчитать значения £, необхо­ димо знать Р, Wr и Wq как функции обобщенных сферических ко­ ординат L, Ѳ и ф. Для Re <0,1 на стр. ПО показано, как можно применить закон Стокса в условиях ползущих течений.

Для значений Re < 0 ,1 очень трудно установить количествен­ ные зависимости для распределения скоростей и давлений на по­ верхности обтекаемой потоком частицы, но принципиально они мо­ гут быть найдены при решении системы уравнений движения, уравнения неразрывности и соответствующих граничных условий:

1)

при 1 =

1

Wr =WQ = 0

(4-44)

2)

при L =

оо

ЦГ2 =

1

(4-45)

3)

при L =

оо

Р =

0

(4-46)

Отсюда следует, что

Р =

Г (X, Г, Z,

Re)

(4-47)

W =

r ( X , Y , Z ,

Re)

(4-48)

где X, Y, Z — обобщенные координаты.

Таким образом, можно записать обобщенную зависимость:

S = f(Re)

(4-49)

Вид этой функции для обтекания шарообразных частиц диамет­ ром d4 на основе большого числа экспериментальных данных пред­ ставлен на рис. 4-8 *.

Графическую зависимость 4-8 часто называют кривой Редея,

118

По рис. 4-8 видно, как меняется механизм переноса количе­

ства движения.

В области ползущих течений (Re < 0,1) действие силы сопро­ тивления подчиняется закону Стокса и в соответствии с аналити­ ческим решением системы уравнений Навье — Стокса и уравнения неразрывности для шарообразных частиц коэффициент сопротив­ ления рассчитывается по уравнению (4-18) и зависимости 24/Re соответствует прямой участок в логарифмических координатах.

Рис. 4-8. Зависимость коэффициента сопротивления среды от режима обтекания шарообразных частиц.

Следует отметить, что в случае обтекания (или осаждения) шаро­ образных частиц на графике £ = f(Re) переход от ламинарного режима к турбулентному не выражен так отчетливо, как при дви­

жении потоков в трубах.

5

При развитой турбулентности в области

5-10 < К е < . 2 ' Ш

потока сопротивлением трения можно пренебречь, а лобовое со­ противление становится преобладающей силой. При наступлении автомодельности (по числу Рейнольдса) коэффициент сопротивле­ ния можно считать постоянным (£ » 0,44). Этой области функции

£ =

/(Re) соответствует квадратичный

закон Ньютона,

т. е. сила

сопротивления Fк :

: Wo.

в пределах изменения 2 <

<

Промежуточный

режим обтекания

Re < 5 - ІО2 характеризуется меньшей зависимостью

сопротив­

ления от числа Рейнольдса, чем показывает расчет по закону Стокса для Re < 2:

Эта эмпирическая формула часто называется уравнением Аллена.

119

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

Рис. 4-9. Коэффициент сопротивления среды для обтекаемых частиц различной формы:

/ — шар; 2—цилиндр; 3—диск.

Фактор ф для частиц неправильной формы находится как отноношение поверхности шара fm, имеющего такой же объем, как и реальная частица, к поверхности частицы f:

f

v2,t

(4-51)

Ф = -у- =

4,878- у -

Поверхность шара с объемом, равным объему частицы, опре­

деляется по эквивалентному диаметру, равному d3 = V 6 Ѵч/п. От­ куда

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

g = /(Re, ф)

(4-52)

Значения фактора формы для различных твердых тел приве­ дены в табл. 4-1.

Фактор формы для сыпучих тел, используемый обычно при расчете со- противления при течении через насадку или слой сыпучего мелкозернистого материала, равен обратной величине ф:

ф': ± = 0,207 л

y2/s

Для шара ф' = 1, для других форм частиц ф' > 1.

120

 

 

Т а б

л и

ц а

4 - 1

 

З н а ч е н и я

ф а к т о р а ф о р м ы

ф д л я

н е к о т о р ы х

 

 

г е о м е т р и ч е с к и х

т е л

[ 3 ]

 

 

а длина, г —

радиус,

А

высота.

\

Форма частиц

Ф

Форма частиц

 

Ф

Ш а р ...........................................

 

 

1,00

К уб '...............................................

а X

 

0,806

Призма (а X

2 а ) . . . .

0,767

Призма ( а X

2а X

2 а ) . . . .

0,761

Призма ( а Х 2 а Х З а ) . . . .

0,725

Диск г ) ...............................

 

 

0,827

Диск ( h = r ß ) .......................

0,594

Диск =

г / 10).......................

0,323

Диск =

г / 15).......................

0 ,2 2 0

Цилиндр (Л — З л ) .......................

0,860

Цилиндр — Юл) ......................

0,691

Цилиндр = 20л) ...................

0,580

Коэффициент сопротивления для несферических частиц можно определить по зависимости, аналогичной уравнению (4-18) для шаров:

А

( 4 - 5 3 )

Re

где А = /(ф).

Для ламинарного обтекания несферических частиц может быть использовано эмпирическое уравнение

А

_____________2 4

 

 

( 4 - 5 4 )

 

0 , 8 4 3 l g -

Ф

 

0

, 0 6 5

Для шарообразных частиц (ф = 1) уравнение (4-54) дает зна­ чение А = 24, и зависимость (4-53) превращается в уравнение

(4-18).

Для турбулентного обтекания (в соответствии с законом Нью­ тона) коэффициент сопротивления не зависит от числа Re, а только от фактора формы:

£ =

5 , 3 1 — 4 , 8 8 ф

( 4 - 5 5 )

Для шарообразных частиц (ф = 1) получим £ = 0,43.

Общий закон сопротивления среды не зависит от природы сил, вызывающих движение твердых частиц в этой среде.

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

Скорость осаждения

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

121

движения, то не имеет значения, что именно находится в состояний

движения — частица или среда.

Свободное осаждение. Рассмотрим частицу с массой т (рис. 4-10), движущуюся в неподвижной среде под действием внешней силы. Эта внешняя сила может быть силой тяжести или

силой центробежного поля.

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

не уравновесит силу

тяжести.

Затем частица будет продолжать

 

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

 

Эту постоянную скорость и называют скоростью

 

свободного осаждения w0c- Таким образом, при

 

падении частицы имеют место три стадии ее

 

движения: 1) начальный момент падения; 2) дви­

 

жение с увеличением скорости; 3)

равномер­

 

ное движение

(с постоянной скоростью).

Возра­

 

стание СКОРОСТИ

ОТ W = 0

ДО W =

Шкон =

w oc

 

происходит в течение очень короткого промежут­

 

ка времени (например частица пыли диаметром

 

10 мкм с плотностью ртв =

 

2700 кг/м3

дости­

 

гает

постоянной

скорости

 

осаждения

через

 

0,006 с), поэтому для технических расчетов пред­

 

ставляет интерес лишь третья стадия движения

 

тела.

 

 

 

 

 

 

 

Рис. 4-10. Силы,

В условиях динамического

равновесия

прин­

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

цип Д ’Аламбера для движущейся частицы

при­

падающую части-

водит к уравнению:

 

 

 

 

 

 

 

 

FB — Fc — Fn — т

 

 

 

(4-56)

Внешняя сила (сила тяжести) по закону Ньютона может быть

выражена следующим образом:

 

 

 

 

 

 

 

 

 

Рв = та

 

 

 

 

(4-57)

где а — ускорение движения частицы в результате действия

внеш­

ней силы FB.

 

 

 

 

 

 

 

 

 

Сила сопротивления среды по уравнению (4-15):

 

 

 

 

 

Fc

£pw2 f

 

 

 

 

(4-58)

 

 

2

'

 

 

 

 

 

 

 

 

 

 

 

 

где W — скорость частицы относительно среды;

р — плотность

сре­

ды; / — площадь

поперечного сечения частицы;

£ — коэффициент

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

 

 

 

 

 

 

 

 

 

Подъемная, или архимедова, сила пропорциональна массе

среды, вытесненной массой частицы:

 

 

 

 

 

 

 

 

 

п

^

 

 

 

 

(4-59)

Из уравнений

(4-56) — (4-59)

Утв

 

 

 

 

 

 

получил

 

 

 

 

 

 

 

dw

Ра

_ If

 

 

 

 

 

 

 

dx = а

Ртв

5

 

 

 

 

(4-60)

122

Уравнение (4-60) характеризует взаимодействие сил, в поле ко­ торых находится твердая частица или тело. Для его решения не­

обходимо знать

природу

внешней силы

и закон сопротивления.

В случае падения частицы под действием силы тяжести а = g,

где g — ускорение свободного

падения.

Тогда

уравнение (4-60)

примет вид:

 

 

 

 

 

 

 

( , ____ Р \

f

(4-61)

 

d r

® \

pTB) 2m

 

 

Если частица

осаждается

в центробежном

поле, то а = га2,

где

со — угловая скорость, а г — радиус пути частицы. В этом слу­

чае

уравнение (4-60) будет иметь вид:

d w

£рw2

ли2

(4-62)

d r

2 т '

Уравнения (4-61) и (4-62) одинаково важны при решении за­ дач, связанных с разделением неоднородных систем.

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

В соответствии с третьим допущением dw/dx = 0, и тогда при введении массы шарообразной частицы m = (nd2/6)pTB и площади

поперечного сечения частицы f — nd24j 4 в уравнение (4-61) получим:

g

или

g

Откуда

р

\

_

3£рш 2

Ртв

/

 

(4-63)

 

4 d 4p TB

р\ з ;р а 4

(4-64)

р тв /

4с!чр тв

4

(ртв

р) g d ч

(4-65)

з

'

PS

 

Если вынести постоянные величины из-под корня, то получим окончательно следующее выражение для woc (в м/с):

w oc = 0,227 ] / d 4 (pp g ~ Р )"

(4-66)

В приведенных выше уравнениях d4— диаметр осаждающейся частицы, м; ртв и р — плотность частицы и среды соответственно, кг/м3.

Уравнения (4-65) и (4-66) справедливы для ламинарного, пе­ реходного и турбулентного режимов осаждения частиц с различ­ ной сферичностью, причем коэффициент сопротивления £ может быть определен по рис. 4-11.

Для ламинарного режима осаждения (ReOc < 0 , l ) сила сопро­

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

(рис. 4-10), по Стоксу вы­

ражается зависимостью:

 

Рр — ЗяЛуЦф

(4-67)

123

Тогда для равновесия действующих сил и силы инерции т

можно записать:

т dx

m( l

pTB) g

(4-68)

 

Для шарообразных

 

 

nd\

 

 

частиц т ~ g—ртв, следовательно

 

nd\

dw

nd^

 

Р) S 3ndqptfi)

(4-69)

0

Ртв ^

--- 0

(Ртв

или

dw

(ртв р) g

18p®

 

 

(4-70)

 

 

р тв

 

^чРтв

 

 

 

 

Скорость осаждения при dw/dx =

0:

 

 

 

Wnr =

-р ) 8

(4-71)

 

 

 

18р

 

 

 

 

 

 

Зависимость (4-71) носит название формулы Стокса и спра­ ведлива для области ІО'4 < Re0c < 2. В соответствии с верхним и

Р

Рис. 4-11. Зависимость коэффициента сопротивления среды от Reoc и фактора формы (сферичности) ф частиц, осаждающихся под действием силы тяжести

нижним пределами применимости этой зависимости можно рассчи­ тать максимальный диаметр осаждающихся по закону Стокса ча­ стиц, подставив в формулу (4-71) вместо скорости w0с, например, ее выражение из Re0c = wocd4p/[x = 2. При этом в случае осажде­ ния частиц в газовой фазе величиной плотности среды р можно пренебречь.

124

При осаждении в газовой фазе нижний предел применимости формулы Стокса имеет особо важное значение. Когда размер ча­ стицы приближается к длине пути свободного пробега молекулы, скорость осаждения, рассчитанная по уравнению (4-71), будет за­ вышена.

Поправочный коэффициент k при осаждении частиц пыли с

d4 > 10 мкм в воздухе приближенно равен

1. В других случаях k

рассчитывается по эмпирическому уравнению:

k = l + a-£-

(4-72)

где а зависит от свойств среды и меняется от 1,4 до 20 (для воз­ духа а ä: 1,5). Для частиц пыли размером d4 < 0,1 мкм возникает броуновское движение, при котором частицы беспорядочно дви­ жутся во всех направлениях вследствие соударений, не осаждаясь.

Если осаждение проводится в аппарате, диаметр которого Dann сопоставим с размерами осаждающихся частиц, то следует учиты­ вать так называемый пристеночный эффект (стенки аппарата бу­ дут замедлять осаждение). Значения поправочного коэффициента b (на который надо умножить woc, рассчитанную по формуле Сток­ са) в зависимости от отношения <7ч/0 аіт приведены в табл. -4-2.

Т а б л и ц а 4-2

Значения коэффициента 6, учитывающего замедляющее влияние стенки аппарата на осаждение твердых частиц

 

 

 

(Re < 2)

[3]

^ч/^апп

ь

1

йч/-°лпп

ь

0,0

1,000

 

0,5

0,170

0,1

0,792

 

0,6

0,0945

0,2

0,596

 

0,7

0,0468

0,3

0,422

 

0,8

0,0205

0,4

0,279

 

 

 

В области действия закона Ньютона (при турбулентном режиме осаждения)

Для того чтобы определить режим осаждения и, следовательно, выбрать формулу для расчета скорости осаждения, необходимо знать величину критерия Re, в который также входит woc. В связи с этим уравнения (4-65) и (4-71) применимы для расчета woc ме­ тодом последовательных приближений, т. е. на первой ступени расчета приходится задаваться, например, ламинарным режимом осаждения, а затем, определив wос, проверять, лежит ли Reoc в об­ ласти, соответствующей принятому условию. При несовпадении результатов переходят ко второй ступени расчета до получения удовлетворительной сходимости данных.

125

Однако такой трудоемкой процедуры расчета woc можно избе­ жать, преобразовав уравнение (4-65) методом Лященко [4]. Метод

основан на подстановке в уравнение (4-65) выражения woc— Reft-

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

 

а У з

Ртв Р

(4-74)

 

и2

р

 

 

Правая

часть уравнения (4-74)

представляет собой

критерий

Архимеда.

Таким образом

 

 

 

С Re2 = ~

Ar

(4-75)

Зная величину критерия Архимеда (для осаждения частиц задан­ ного размера), по графику Лященко (рис. 4-12) легко рассчитать значение Re, из которого определяется

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

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

В области ламинарного режима осаж­

дения (при

Re <

2,

т. е. в условиях,

характеризующихся

законом

Стокса)

£ = 24/Re

и

уравнения (4-75)

и

(4-71)

примут вид:

 

 

 

 

 

Re =

Ar/18

или

Re = 0,056 Ar

(4-76)

Критическое значение критерия Архиме­ да, соответствующее верхнему пределу числа Re, будет Агкр. і = 36, следователь­

но, существование ламинарного режима осаждения ограничивается

условием Ar ^

36.

 

 

В области действия

закона Ньютона

(в условиях автомодель­

ности критерия

Re) £ =

0,44 и уравнение (4-75) можно предста­

вить так:

 

 

 

 

 

Re =1,74 Ar0,5

(4-77)

В переходной области верхнее предельное значение критерия Архимеда соответствует значению Re = 500 и рассчитывается по уравнению

'

Re = 0,152 Ar0’715

(4-78)

Таким образом, область осаждения в переходном режиме огра­ ничивается изменением критерия Архимеда в пределах 36 < Аг < < 8 ,3 - ІО4,

126

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