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

книги из ГПНТБ / Пирумов, А. И. Обеспыливание воздуха

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

Определению суммарного коэффициента захвата пылевых частиц единицей длины волокна в фильтре посвящено много ■серьезных исследований1. В последнее время все большее ■распространение получает метод расчета улавливания очень мелкодисперсных частиц (^ < 5 мкм), разработанный Н. А. Фуксом и его сотрудниками [47]. для высокоэффективных фильтров при скоростях фильтрации менее 1 м/с. Метод осно­ вывается на определении коэффициента захвата Ks «веерной» модели фильтра, представляющего собой систему решеток из параллельных волокон, повернутых на произвольные углы отно­

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

В элементарном слое фильтра толщиной dh, образованного волокнами диаметром D, уложенными с одинаковой плотностью

,

4 a dh „

упаковки та, общая длина волокон

равна'— —-.П ри протека-

 

я и 1

нии со скоростью w через этой слой аэрозоля, концентрация частиц в котором равна Nо, каждой единицей длины волокна улавливается KsDN0w/t(\—та) частиц данного размера.

Убыль частиц в потоке с переменной концентрацией N со­ ставит — wdN, или

d N

 

4 a

(1.12)

N

=

я й ( 1 - а )

 

Учитывая, что величина

а

очень мала, после интегрирова­

ния получим:

1,72 a hKs

lg ( 1 - У , ) = 2 -

я D

Практический интерес представляет закономерность распре­ деления осаждающейся пыли по толщине фильтра. Для выяс­ нения этой закономерности запишем выражение (1.12) в виде

d N

_ dh

N

Я ’

где Я — постоянная фильтрации. После интегрирования получим:

(i.i3)

м

 

откуда следует, что .при Л=Я

коэффициент проскока равен

£-1 = 0,368.

В соответствии с этим постоянную фильтрации можно опре­

делить как толщину слоя фильтра, в котором

задерживается

63,2% всех частиц:

 

Я= -----Ц г - .

(114)

2-31§-д^

 

1 Полный обзор развития и современного состояния теории фильтрации приведен в работе \[107].

21

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

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

тивности можно получить с помощью

приводимых

теоретичес­

ких зависимостей.

действия электрических

Отделение пылевых частиц в поле

сил обычно протекает в две стадии:

сначала они

заряжаются

в поле коронного разряда, а затем осаждаются под действием электрических сил1.

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

ные и отрицательные ионы.

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

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

1 Принцип электрической очистки газов от пыли детально рассмотрен в работе [106].

22

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

Закономерности инерционной сепарации пыли в криволинейных потоках

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

Движение частицы в криволинейном потоке при отсутствии внешних воздействий описывается уравнением

dV

1

(1.15)

dt

 

 

В работе [75] давалось

приближенное решение

уравнения

.(1.15), при этом движение частицы рассматривалось в подвиж­ ной системе координат. Скорость частицы принималась равной v = w + v c, а ее ускорение определялось согласно теореме Кориолиса.

Уравнение (1.15) при этом может быть представлено в виде

dvc + т o)X(toxR) + tn

dm

d \iB vc. (1.16)

X R + 2 m ( » X v c) = —3 я

dt

dt

 

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

представ­

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

и изменить ее скорость.

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

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

23

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

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

скорость пылевых час­

тиц, равная до этого скорости потока

(v0= w 0), изменяется

сравнительно медленнее. В силу этого в начальный момент вре­ мени t —О можно принять vc= w 0w(R0).

В соответствии с указанным, в правой части потока, пред­ ставленного на рис. 1.2, скорость vc направлена к его пери­ ферии. Легко убедиться, что вектор mdvc/dt направлен в ту же сторону, что и вектор центробежной силы, хотя не всегда пол­ ностью совпадает с ним. Слева от точки М эта сила обращена к центру потока.

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

Четвертый член уравнения представляет собой силу Корио­ лиса. В рассматриваемом случае она направлена также к периферии потока.

Сепарация частиц на начальном участке движения. Расчет сепарационного движения частиц в криволинейном канале мо­ жет быть произведен с помощью ЭВЦМ1.

Схема потока представлена на рис. 1.2. В неподвижной сис­

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

w —k/R

и R = y r

х2-\-у2. Тогда уравнение

(1.15) в проекциях на оси координат имеет вид

 

 

 

 

rf2 X

 

1

 

dx

+

k

.

y

 

0-

 

 

 

 

dt2

 

Г

 

dt

 

 

 

 

 

 

 

 

 

T

x2 + y2

 

 

 

(1.17)

 

 

 

 

1

 

dy

 

k

 

 

 

 

 

 

 

d 2 u

 

 

 

X

— 0.

 

 

 

 

dt2

+

т ' dt

 

T

X 2 + у2

 

 

 

 

 

 

 

 

 

 

Для возможности

 

расчета

на ЭВЦМ

по методу Рунге —

Кутта уравнения

(1.17)

преобразуются

в

систему

уравнений

первого порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

dy

-

 

 

 

 

du

 

_1_

_k_

 

 

У .

 

dq

1

, _fe_

x

 

dt

~~

x

x

 

x2 + y2

 

dt

=

x

q ~r x ‘

x2 + y2

'

Начальные условия имеют вид

 

 

 

 

 

 

 

t =

0; x0 = R0;

(/о =

0;

w0 = 0;

д0 = ш0.

 

 

При

расчете принимали

 

т=3,06-10~3с

(d=20

мкм),

х =

= 27,5-10-3

с { d = 60

мкм)

и

т =

172,1 -10—3 с (с?=150 мкм).

1 В работе [75] уравнение (1.15) решалось приближенно с помощью разложения функции в степенной ряд и введения некоторых упрощений.

24

Расчет проводили для значений R q, равных 0,1; 0,15;

0,2; 0,3

и

0,4 м. Средняя скорость входа воздуха w0= 15 м/с;

=0,005

м;

•/?2=0,5 м; й— 2,93 м2/с. Для каждого

интервала

 

времени

At определялись координаты частицы x(t),

y(t), R — У

х2-{-у2,

скорости dxfdt и dy/dt и угол <р*.

 

 

траек­

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

тории частиц на неограниченном расстоянии.

 

 

 

Рис. 1.4. Траектории в криволинейном потоке

/ — частиц крупностью 20 мкм; 2 — то же, 60 мкм; 3 — то же, 150 мкм; 4 — то же, 60 мкм по упрощенному р асч ет

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

окружности очень большого радиуса. В области

траек-

* Расчет проведен на ЭВЦМ «Мир-2» вычислительного отдела

ЦНИИ-

Промзда.ний.

 

25

тории частиц искривлены меньше, чем в области ау>ш0. Каи и следовало ожидать, чем крупнее частица, тем меньше искрив­ лена ее траектория. График позволяет оценить возможности циклонной сепарации. Частицы крупностью 20 мкм и меньше,, входящие в циклон на расстояние более Vs радиуса от его на­ ружной стенки, должны совершить более двух оборотов в сво­ ем спиральном нисходящем движении, чтобы достигнуть этой стенки. Между тем, продолжительность пребывания частицы в циклоне может не удовлетворять этому условию.

Анализ уравнений (1.17) приводит к ряду интересных выво­ дов. Прежде всего перейдем от декартовых координат к поляр­ ным, введя в уравнения значения первых и вторых производ­ ных от координат x=Rcosq> и r/=Psincp. После преобразований получим:

 

d2 R

l

d ф \ 2

1 d R

 

 

(1.18)

 

dt2 - * ( Л + , ■

- o;

 

 

 

 

d R

d ф

 

d2 ф

1

d ф

k

(1.19)

2 ——

dt

R + R2 d ! - + r Ri

dt

x - ° '

dt

 

откуда

 

 

 

 

 

 

 

 

 

R2

dt ~ k + c*e~ t/X■

 

(1.20)

 

 

 

 

Применительно к воздушному

потоку Rdcpldt = w. В дан-

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

 

k +

c9e ~ tlx

 

v<t>~

R

Постоянная

интегрирования определяется из условия R = Rq

п и9 = v0 = w0

при t=0. Тогда c0= w 0Ro — k и

Из выражения (1.21) видно,

что с течением времени танген­

циальная составляющая скорости движения частицы стремится к тангенциальной составляющей скорости воздушного потока, что в пределе приводит к квазистацйонарному состоянию дви­ жения. Оно достигается тем скорее, чем меньше время релаксаот нуля, является участком нестационарного движения частицы, в пределах которой значение е ~*1х еще существенно отличается от нуля, является участком нестационарного движения частицы, т. е. начальным участком. Длина этого участка различна для

частиц разного размера.

Для практических

расчетов можно

принять, что Уф = ш при

7т, т. е. при e~tlx

^10~3.

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

2S

Отметим, что из выражений (1.18) и (1.20) имеем:

d2R

,

1

 

dR

[k + (w0R0 — k ) e - t/xf

==0-

( 1. 22)

dt2

^

х

'

dt

R3

 

Сепарация частиц за пределами начального участка движе­ ния. В этом случае время t достаточно велико и уравнение

(1.22) принимает вид [96, 119]

d2R

J _

dR

k2

dt2 +

x '

dt

(1.23)

~ R3

Полученное уравнение

типа Эмдена — Фаулера не может

быть решено в элементарных функциях. Принимая обозначение dR/dt— vR, приведем уравнение (1.22) к следующему виду:

dvR

[ k + (w 0R 0 - k ) e - V * ] 2

1

 

т г

- ---------- v ? ------------- - ■

<r24

В фазовой плоскости (R, vR) уравнение (1.24) представляет собой изменение углов наклона касательных к кривой vR в за­ висимости от R (рис. 1.5). Очевидно, что значение dvR/dR = 0 соответствует максимуму скорости vR. Из уравнения (1.24) в этом случае имеем:

 

 

 

VR макс

+ с0 е т Т

(1.25)

 

 

 

R*3

 

 

 

 

где

R* — значение R,

соответст­

 

вующее VR Макс-

 

 

 

 

Для частиц, входящих на по­

 

ворот в точке М ( с о =

0 ) ,

а так­

 

же

для

очень

мелких

частиц

 

(е~4/т—>0) независимо ют их на­

 

чального положения, т.

е. для

 

всех частиц, у траекторий

кото­

 

рых

нет

начального

участка,

 

уравнение

(1.25)

имеет вид

 

 

 

 

k2

 

(1.26) Рис. 1.5. Изменение окозости

 

 

 

R*3

 

 

 

 

 

стицы в фазовой плоскости

 

Уравнение (1.26) представляет собой кубическую гиперболу, на которой расположена точка, соответствующая максимально­ му значению vR, т. е. точка максимума функции vR= f(R).

Отбрасывая в порядке первого приближения первый член в

уравнении (1.23), получим равенство

Ь2

» , р = т —

,

(1-27)

совпадающее с уравнением (Ы 0), поскольку vcvR. Очевидно, что на рис. 1.5 уравнение (1.27) представляет собой кривую 1, а

27

ее участок в области между значениями Ям и R2 дает величины скорости vr, принимаемые при пользовании формулами типа (1.10) для расчетов сепарации. Действительная величина ско­ рости определяется кривой 2, а погрешность расчета оценивает­ ся площадью, заключенной между кривыми 1 и 2.

Интегрируя уравнение (1.27),

получим при t= 0, Я=Яо и

Cl= #0

 

 

Д4 = Я« + 4т&.

(1.28>

Пользуясь этим уравнением,

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

траектории

частиц для рассмотренного выше течения. Положение частицы

на траектории определяется

величинами Я и <р по формулам

(1.20) и (1.28), а именно: RdR=xkd(f, и после

интегрирования

2xk(f=R2jrc2. При ср=0 и Я = Я 0 постоянная

интегрирования

с2= — Я о и

 

 

Ф =

Я1 -

(1.29>

2 х k

 

Вычисленная для частицы размером 60 мкм при Яо— 0,2 м траектория нанесена на рис. 1.4 (кривая 4) для возможности сопоставления ее с приведенными на том же графике более точными траекториями частиц. В то время как истинные тра­ ектории на всем протяжении отклоняются от прямой первона­ чального движения в направлении течения, т. е. частицы как бы «сдуваются» потоком, траектории, построенные по упрощенной формуле, отклоняются в противоположную сторону. Такой характер движения не может быть объяснен действием какихлибо физических факторов. Он, несомненно, является следстви­ ем того, что величина членов дифференциального уравнения, отброшенных при составлении упрощенной формулы, соизмери­ ма с величиной оставленных членов, т. е. настолько велика, что их игнорирование приводит не только к количественным, «о и к качественным погрешностям.

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

х № 3 t 3 fe4

vi r

Аналогично можно получить последующие приближения. Расчеты показали, что для рассмотренного примера вблизи значения Я2 второе и третье приближения достаточно точно оп­ ределяют величину vR, однако в области между значениями Ям и Я* истинные значения скорости таким путем получить нельзя.

28

На рис. 1.5 показано изменение скорости движения части­ цы, вошедшей в поток в равновесной точке М. Для других час­ тиц фазовая картина изменения vR отличается расположением гипербол 1. Из рис. 1.5 следует, что выражениями типа (1.10) можно пользоваться для оценки величины радиальной состав­ ляющей скорости движения частиц за пределами начального участка, а для мелких частиц — вблизи наружных стенок. Рас­ считывать с их помощью траекторию частиц и соответственно «минимальные диаметры» сепарирующихся частиц нельзя, так как при этом вместо истинной скорости vR, изменяющейся по кривой АВ, принимаются значения скорости по кривой А'В'У намного большие на участке Rm R*.

Сепарация пыли в пограничных областях турбулентных потоков

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

Рассмотрим движение частицы, находящейся в потоке, вблизи выделенной площадки. Очевидно, что такая частица бу­ дет вовлечена в пульсирующее поперечное движение, при этом в тем большей мере, чем меньше ее размеры и масса. Степень увлечения частицы размером 1 мм (при плотности, равной 1) составляет всего 2%, но частица размером 100 мкм увлека­ ется на 70%, а частица размером менее 30 мкм — практически полностью [ПО].

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

сил турбулентное перемешивание не может нарушить этого распределения. При неравномерном начальном распределе­ нии частиц перемешивание приводит к их распространению по всему сечению потока.

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

дне.

Экспериментальные исследования распределения пульсаций скоростей позволили подойти к количественной оценке явления взвешивания частиц потоками [16]. В частности, исследования Е. М. Минского [61] показали, что среднеквадратичные пуль­

сации скоростей w 'г изменяются по следующей закономер­ ности:

V w 'z = 1,05 вуср ] / о , 1 2 5 - ^ - 1 + 1,72

1,81

(1.30)

где wср — средняя скорость течения в канале; X — коэффициент трения; ki — коэффициент корреляции; I — расстояние от оси до рассматриваемой точки; — гидравлический радиус канала.

Как видно из формулы (1.30), наименьшее значение пульса­ ций скоростей соответствует оси потока, когда 1=0, а наиболь­ шее их значение — его периферической части, например вблизи дна канала (для пограничного слоя потока формула неприме­ нима). Подобная закономерность наблюдается в распределении пульсаций скоростей и в циклонах [120].

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

оседает до прежнего уровня.

 

 

__

Пользуясь выражением

(1.30), из условия vs(d)=-f/rw'*полу­

чим для граничного размера частиц:

 

 

 

Цв

-I Г

х Г

( i у -8'

-----

У 0,125

1 + 1,72 —

.(1.31)

p g

r

fei L

V RT /

J

30

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