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

книги из ГПНТБ / Бушмелев, В. А. Процессы и аппараты целлюлозно-бумажного производства учебник

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

ПЕРЕНОС ВЕЩЕСТВА В ТВЕРДЫХ И ПОРИСТЫХ ТЕЛАХ

В твердых телах перенос вещества осуществляется на микроуровне, так как конвективного движения среды в теле нет. Аналогично веще­ ство распространяется и в пористых телах, поскольку жидкости и газы в порах неподвижны или движутся ламинарно. В отличие от молекулярной диффузии перенос вещества в твердых и пористых те­ лах назван м а с с о п р о в о д н о с т ь ю . Уравнение массопровод­ ности аналогично уравнению молекулярной диффузии (12-1). Раз­ ница лишь в том, что вместо коэффициента молекулярной диффузии D в него подставляется коэффициент массопроводности DM, имеющий ту же размерность, что и D. Процесс извлечения вещества из твердой фазы массопроводностыо нестационарен. Сначала удаляется вещество с поверхности твердой частицы, вслед за этим — из более глубоких слоев и только затем — из центра частицы. При этом с течением вре­ мени изменяется не только длина пути потока извлекаемого вещества, но и его концентрация. С поверхности твердой или пористой частицы в окружающую среду вещество переносится вследствие молекуляр­ ной и конвективной диффузии.

РАСЧЕТ КОЭФФИЦИЕНТОВ ДИФФУЗИИ

Наиболее достоверными данными о коэффициентах диффузии яв­ ляются данные эксперимента, которые приводятся в справочниках.

Коэффициент молекулярной диффузии одного газа в другом мо­ жет быть определен по формуле (в м21сек)

 

 

0,43-10

6т 3'2

A l l

(12-3)

 

 

р (»!,3 +

»2,э)2 V

М і М 2 ’

 

 

 

где

р — абсолютное давление, ат\

 

М х и М 2 — молекулярные массы компонентов, кг!кмоль\

 

ѵі

и ѵ 2

— их молярные объемы при температуре кипения, т. е.

 

 

для жидкой фазы (эти величины не следует путать

 

 

с молекулярными объемами газов), см3!моль.

 

Коэффициент молекулярной диффузии газа в многокомпонентной

газовой

смеси

равен

 

 

 

 

 

D f

 

■+

(12-4)

 

 

^1—2 Яг-*

 

 

 

 

 

где

аі> а 2> ап — молярные (объемные) доли компонентов в

 

 

газовой

смеси;

 

 

£>]_2, Пх-з, П і- П— коэффициенты диффузии для бинарных сме­ сей 1—2, 1—3, 1—п, определяемые по фор­ муле (12-3).

Пересчет коэффициента молекулярной диффузии D lt вычисленного

для абсолютной температуры Т г и абсолютного давления р 1г на

усло­

вия с температурой Т 2 и давление р 2,

производится по формуле

 

п

— Г) ( Ъ \3!2.В±

(12-5)

D2

- Di [ t J

р2

 

247

Коэффициент молекулярной диффузии газов и жидкостей в водных растворах равен

 

В = 0,119-10-10- ^ — ■,

(12-6)

где

р, — вязкость воды в сантипуазах

при температуре t\

Т = 273 +

t — абсолютная температура воды;

М— молекулярная масса диффундирующего вещества,

кгікмоль;

V — его молярный объем, см3/моль.

Значения коэффициентов молекулярной диффузии в газах нахо­ дятся в пределах 0,3- ІО-5 -ь 7 - ІО-5 м2/сек\ в жидкостях'— в преде­

лах 0,4 -ІО-9 7 - ІО-9 м21сек. Диффузия в газах проходит примерно в ІО4 раз быстрее, чем в жидкостях.

МАССООТДАЧА

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

М = ßA CF т,

(12-7)

где ДС — движущая сила, равная разности концентраций

вблизи

поверхности раздела фаз и вдали от нее (при потоке веще­ ства от поверхности) или разности концентраций вдали от поверхности и вблизи нее (при потоке к поверхности массо­ отдачи);

ß — коэффициент пропорциональности, называемый

коэффици­

ентом массоотдачи; его размерность [ß] =

[кг/м2 ■сек

(кг/м3) ].

 

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

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

в кг/м3

[р ] =

м/сек. Если движущую

силу выразить в кг/кг, размерность

[ß ] =

кг/сек-м2. Вычисление

коэффициента массоотдачи ß,

как и

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

Основы диффузионного подобия

Первый и второй законы Фика в представленном у нас виде [фор­ мулы (12-1) и (12-2)] для практических расчетов не применяются. Од­ нако они используются для вывода критериев диффузионного подобия.

При совместном решении уравнения массоотдачи (12-3) и первого

248

уравнения

Фика (12-1) получаем равенство —D — = ßAc. Отсюда

- — cdx

dx

_ j ^ сключив знаки д и дифференциалов, отбросив знак

D de J

 

минус и заменив х на стандартное обозначение характерного линейного

размера I, получим критерий Nu' = ^ . Это диффузионный критерий

Нуссельта, характеризующий условия массоотдачи на границе раз­ дела фаз. Критерий Nu' является искомым критерием, так как в него входит определяемая величина коэффициента массоотдачи ß.

Для получения о п р е д е л я ю щ и х к р и т е р и е в исполь­ зуем второе уравнение Фика (12-2). Разделим второй член левой части на правую часть и в полученном безразмерном комплексе отбросим

знаки дифференциалов. После замены

л: на стандартное / получим

— = Ре' — диффузионный критерий

Пекле, который является ме­

рой отношения диффузионных потоков конвекции и молекулярной диф­ фузии и характеризует условия подобия при конвективном массопереносе.

Если в уравнении (12-2) правую часть разделить на первый член

D x

левой части, то после подобных преобразований получим — = Fo' —

диффузионный критерий Фурье, который характеризует условия по­ добия при неустановившемся массообмене. Как и при теплообмене,

критерий Ре' = RePr)., где Р г '= -- -----

диффузионный критерий

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

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

обменного аппарата (например, Гх = у , Г2= у

и т. д.). Крите­

риальное уравнение массоотдачи имеет вид

 

Nu' = А Re“ (Рг')6Г“Г^

(12-8)

где А, а, b, с, d — постоянные, характерные для каждого частного случая массоотдачи.

Равновесие между фазами

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

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

249

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

переход вещества в первую фазу, а через некоторое время снова установится рав­ новесие.

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

 

 

Массопереход

наблюдается лишь

Рис. 12-1. График измерения

в том случае, если система не находится

в равновесии.

Фазовое равновесие

ха­

равновесных

и рабочих кон­

рактеризуется

отсутствием

видимого

центраций:

перехода вещества из одной фазы в дру­

ОА — линия равновесия; ВС — ра­

бочая линия

массообмена

гую. В состоянии

равновесия

при

дан­

 

 

ных температуре

и давлении существует

определенная зависимость между содержаниями вещества в обеих

фа­

зах: любой концентрации вещества х во второй фазе соответствует равновесная концентрация ур в первой фазе. В общем виде эта зави­ симость представлена на рис. 12-1. Кривая ОА — это линия равнове­ сия.

Условия равновесия, определяемые уравнением ур = f (х) при заданных температуре и давлении, дают четкое представление о воз­ можности или невозможности проведения процесса в данном направ­ лении. Если действительные (рабочие) концентрации распределяе­ мого вещества у больше равновесных ур, вещество переходит из пер­ вой фазы во вторую. При у< іур идет обратный процесс. Отсюда можно заключить, что процесс всегда идет в направлении установления равно­ весия.

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

250

I

равновесные концентрации ур и хр находятся в прямой зависимости от соответствующих концентраций вещества в фазах л* и у:

г/р = т|щ X = \ у ,

(12-9)

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

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

Материальный баланс массообмена

Составим уравнение материального баланса массообменного аппа­ рата по распределяемому веществу (рис. 12-2).

Обозначим: G и L — расходы первого и второго массоносителей; у х и х х, у 2 и х г — их начальные и конечные концентрации. Приход

вещества

равен

его

расходу,

т. е.

Gyx +

Lxx =

 

I

= Gyо -|-

Lx,. Отсюда

 

 

 

 

 

Ч*

 

 

 

 

 

х,

 

G(y1— y2) = L(x2— x1) = M,

(12-10)

 

 

 

где М — количество распределяемого вещества.

 

 

В соответствии с этим

уравнением

количество ве­

V

X

щества G (уху 2),

которое

в аппарате передает фаза -

G, равно

количеству

вещества

L (х 2—лу),

которое

 

 

принимает фаза L. По уравнению (12-10) может

 

 

быть определен удельный

расход второго массоноси­

 

 

теля, который равен

 

 

 

 

 

 

 

('S -" )

Уравнение рабочей линии массообмена

В произвольном месте рассечем аппарат горизон­ тальной плоскостью N —N (рис. 12-2) и составим урав­ нение материального баланса по распределяемому веществу для нижней части аппарата: Gyx -+- Lx = X я у — текущие значения концентраций. Отсюда

У = Уі— ^ ( х 2— х).

I Уг

Рис. 12-2. Ма­ териальный ба­ ланс массооб­ менного аппа­ рата

Gy + Lx о, где

( 12- 12)

Это уравнение представляет собой зависимость изменения кон­ центрации у распределяемого вещества в фазе G от его концентрации * в фазе L. Оно называется уравнением рабочей линии процесса. При постоянных у х, х 2, L и G рабочая линия представляет собой прямую

с углом наклона к оси абсцисс, тангенс которого равен / = — .

Рабочая линия процесса массопередачи строится в координатах у—X по двум точкам. Положение этих точек может быть определено

251

по уравнению (12-12) при подстановке в него

^_L___Уі ~ Уг

G

Xо — лу

 

После преобразований получим

 

Уі ~ У =

Уг.

(12-13)

х„ X

Хг, — хх

'

Это отношение, как и зависимость (12-12), представляет уравнение прямой (рис. 12-1), проведенной через точки В и С с координатами со­

ответственно (і/і, х 2)

и (г/2; лу). Таким образом, для построения рабо­

чей линии процесса

массопередачи (линия ВС на рис. 12-1) нужно

знать содержание распределяемого вещества в массоносителях на вхо­ дах в аппарат и выходах из него. Для построения линии равновесия (линия ОА на рис. 12-1) должна быть известна зависимость ур = / (х) в виде расчетных формул, графиков или таблиц. Построение рабочей и равновесной линий является одним из валдаых этапов расчета мас­ сообменного процесса.

Рабочая и равновесная линии могут быть построены и в координа­ тах рX, где р — парциальное давление распределяемого вещества.

Движущая сила массопередачи

Допустим, что рабочая линия процесса ВС расположена выше ли­ нии равновесия ОА (рис. 12-1). Тогда в соответствии с правилом оп­ ределения направления процесса, по которому процесс должен идти в сторону равновесия, вещество будет переходить из первой фазы G во вторую L. Разность рабочей и равновесной концентраций назы­ вается движущей силой массопередачи. Как видно из рис. 12-1, дви­ жущая сила массопередачи по высоте аппарата неодинакова. В верх­ ней части аппарата на выходе фазы G и входе фазы L движущая сила равна АСв — у 2 Урв, гДе урв — равновесная концентрация рас­ пределяемого вещества в верхней части аппарата, определяемая по концентрации х х. В нижней части аппарата на входе фазы G и выходе фазы L величина движущей силы равна АСн = у г — г/ри, где г/рн — равновесная концентрация распределяемого вещества в нижней ча­ сти аппарата, соответствующая концентрации x z. В других сечениях по высоте аппарата движущая сила будет иметь промежуточные зна­ чения. Движущая сила может быть выражена и через концентрации фазы L.

УРАВНЕНИЕ МАССОПЕРЕДАЧИ

Рассмотрим установившийся процесс массообмена между двумя фазами через элементарную поверхность раздела фаз dF (рис. 12-3).

Обозначим: и б2 — толщины граничных пленок первой G и вто­ рой L фаз;

ßx и ß2 — коэффициенты массоотдачи в фазах.

Допускаем, что вдоль всей поверхности массообмена величины ßx и ß2 постоянны.

252

Концентрации распределяемого вещества, усредненные относи­

тельно элементарной поверхности dF, равны: у и х — в ядрах

фаз;

у' и х'

— на границах ядра и пленки; у} и xf — в первой и второй

фа­

зах у поверхности раздела фаз.

В связи с тем что перенос вещества

конвективной диффузией очень

эффективен,

можно принять у'

= у

и х' =

X.

количество

вещества, передаваемое

По

уравнению массоотдачи

в единицу времени из ядра первой фазы к поверхности раздела фаз,

равно М =

р! yf) dF. При

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

процессе

такой же

поток вещества будет двигаться от поверхности

 

 

раздела фаз в ядро второй фазы:

 

 

 

М = ß2 (xf — х) dF.

 

 

 

 

Допустим, что межфазное

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

 

 

сывается

уравнением ур = фд-,

где

ф — коэф­

 

 

фициент пропорциональности. Это выражение

 

 

показывает, что линия равновесия представ­

 

 

ляет собой

прямую со степенью наклона к оси

 

 

X , равной ф. Поскольку фазы на границе их

 

 

раздела находятся в состоянии равновесия,

 

 

величины

yf

и ду можно связать

уравнением

 

 

yf — ф.ху. Заменив в уравнении массоотдачи для

 

 

второй фазы концентрации х и xf соответственно

Рис. 12-3.

К выводу

на X —-и JCt= ^ .получим М — — (уfyD)dF.

общего

уравнения

массопередачи

ф

 

Ф

ф

' ^

 

 

 

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

У— Uf

М

Ff Ур

Л4ф

PidP'

ß2dF ‘

 

 

Сложив левые и правые части и решив полученное выражение от­

носительно М, получим

 

M = K a ( y - y P)dF,

(12-14)

где Ка — общий коэффициент массопередачи, отнесенный к фазе G.

Он равен

 

* • - 7 7 1 •

(12‘І5)

Р .+ Р>

 

Уравнение (12-14) называется общим уравнением массопередачи. Оно выведено для какой-то элементарной поверхности межфазового контакта dF, которой соответствует движущая сила уур. Всей по­ верхности массообменного аппарата F соответствует средняя движущая

сила АСс. Тогда общее уравнение массопередачи примет вид

 

М = K gF&Cq.

(12-16)

Это уравнение по форме аналогично общему уравнению теплопе­ редачи.

253

Размерность [/<ßl = [кг/сек-м2- (кг/м3) ]. Следовательно, общий коэффициент массопередачи показывает, какое количество вещества переходит в единицу времени из первой фазы во вторую через единицу

поверхности при движущей силе,

равной единице. Для [ДС0 ]

=

[кг/м3\

сокращенная

размерность [/CG]

= [м/сек ], а при [ АСа I

=

[кг/кг ]

размерность

t/Co 1 — [кг/сек-м2].

 

 

Для аппаратов со свободным объемом общее уравнение массопере­ дачи целесообразно выразить через рабочий объем V и объемный ко­

эффициент

массопередачи

К:

 

 

 

М = КѴЬСв.

(12-17)

Объемный

коэффициент массопередачн К — fKa, где

/ — удельная

поверхность межфазового

контакта, м2/м3. Размерность

1 зависит

от размерности [/С0].

Если равновесная концентрация выражена через концентрацию фазы L и равна хр, то движущая сила для элементарной поверхности dF межфазного контакта равна х—л-р, а для всей поверхности ACL.

Уравнение массопередачи примет вид

 

 

M = K lFACl .

(12-18)

Коэффициент массопередачи, отнесенный к фазе L, равен

 

 

 

(12-19)

# і

ßa

 

При сравнении уравнений (12-16) и (12-18) получим соотношение

( 12-20)

*l АСа

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

РАСЧЕТ МАССООБМЕННЫХ ПРОЦЕССОВ

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

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

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

254

чительно больше поверхности насадки. К четвертой группе относятся тарельчатые аппараты со ступенчатым взаимодействием фаз. Поверх­ ность межфазного контакта в них создается при взаимодействии фаз путем распределения легкой фазы G в фазе L в виде струй или пузырь­ ков.

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

Расчет поверхностных, насадочных и объемных массообменных аппаратов

При прямых линиях равновесия рабочие поверхности и объемы аппаратов определяются по общим уравнениям массопередачи (12-16) и (12-17), где количество передаваемого вещества М рассчитывается по уравнению материального баланса (12-10), а коэффициент массо­ передачи по формуле (12-15).

Средняя движущая сила процесса определяется по формуле средне­

логарифмической разности

 

АС = АС|' ~ ДС|(,

(12-21)

2 ,3 lg ^ s

 

АСК

 

где ДСИ и ДСК— наибольшая и наименьшая движущие силы на входе и выходе массоносителей в аппарате.

Если соблюдается условие 2 > - ^ ^ ^ 0 ,0 5 , средняя движущая

АСК

сила может быть определена, как среднеарифметическая, т. е. АС =

= 0,5

С„ + АСк). Определение

ДСН и ДСК рассмотрено выше

(стр. 252).

равновесия расчет, выполненный

При

большой кривизне линий

по рассмотренной методике, приводит к неточностям. В таких случаях массообменные аппараты рассматриваемого типа рассчитывают по числу единиц переноса (ЧЕП) и высоте единицы переноса (ВЕП). Для выяснения сущности этих характеристик рассмотрим массопере­

дачу через элементарную

поверхность межфазного

контакта dF =

= SfdH, где S — сечение

аппарата; f — удельная

поверхность ра­

бочего объема; dH — элементарная высота рабочего объема. Допустим, что движущая сила процесса при этом равна у ур, а концентрация распределяемого вещества на пути dH изменяется на величину dy.

При коэффициенте массопередачи К а количество передаваемого

ве­

щества по уравнению массопередачи равно

dM = K GS[ (уур)

dH.

По уравнению материального баланса это

количество равно

dM =

= — Gdy, где G — расход массоносителя. Знак минус здесь

указы­

вает на то, что концентрация вещества в фазе в процессе массообмена уменьшается.

255

Поскольку левые части в этих двух уравнениях равны, то равны и правые. Приравняв правые части и решив уравнение относительно элементарной высоты аппарата, получим

dH =

dy

а

У — Ур

i<Gsf '

 

Проинтегрировав это выражение от 0 до Я и от у 1до у 2 (концентрации на входе и выходе), получим общую высоту рабочего объема аппарата:

Я = - ^ — Г

( 12- 22)

K aSf,) у — ур

л

Рис. 12-4. Определение числа единиц переноса методом графического ин­ тегрирования

Величина Г —

= п называется

J У — Ур

л

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

п =

Оно характеризует из­

менение рабочих концентраций на единицу движущей силы.

Обозначив

G

= h , из уравнения

Н

------

(12-22) получаем h = — .

 

KGSf

 

n

Эта величина, равная отношению рабочей высоты аппарата к числу единиц переноса, называется высотой единицы переноса.

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

Высота рабочего объема аппарата равна произведению числа еди­

ниц переноса на высоту единицы переноса:

 

Н = nh.

(12-23)

Число единиц переноса обычно определяют с помощью графиче­ ского интегрирования (рис. 12-4). Для этого по оси абсцисс отклады­ вают значения рабочих концентраций у в фазе G в пределах от у г до

цл, а по оси ординат — соответствующие значения —-— . Проведя

УУр

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

256

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