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

книги из ГПНТБ / Попов В.Д. Основы теории тепло- и массообмена при кристаллизации сахарозы

.pdf
Скачиваний:
9
Добавлен:
24.10.2023
Размер:
11.64 Mб
Скачать

Рис. VI17. Коэффициент диффузии сахарозы в насыщенных растворах (П. В. Головин и И. А. Бе_локонь).

12J T7S

100

и

10

&

0,1

t

0,01 4

<?0Ш " 0,001 0,01

аг5 Of Ц751 т и т

г щ

~

l,/in

Рис. VI—18. Расчет сред­

него

размера кристаллов;

/ — II

продукт: 2—1 продукт;

 

3 — рафинад.

Рис. VI—19. Разность тем­ ператур в пограничном слое кристаллов:

=84%;

2 —СВ,

3

Сй =86%,

г'

~~я

а

'У* у ,«ЇІ4І

V

1

0,1

100

1000

Рис. VI—20. Критериальное обобщение опытных данных по кри­ сталлизации при охлаждении (И, А. Белоконь).

сит от теплоты кристаллизации

г к р , массы

выпавшей

твердой фазы AG, коэффициента теплопроводности ра­

створа Я, его толщины б и других

факторов. Величину

At мы находили по рис. VI—19

в

зависимости от кон­

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

косвенно

учитывающей

часть из этих факторов.

 

 

 

 

Учесть все многообразие факторов

одним

приемом

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

Обобщающая

прямая (рис. VI20)

соответствует

уравнению:

 

 

 

 

 

 

р о ' = 5 8 2 - 1 0 2 р 7 ^ Ь - 4 -

 

 

( V I - 4 5 )

Область применения формулы (VI—45):

 

Fo' = 0,2 - н 600;

П а = 1 5 + 1 0 0 ;

Рг' =

(0,5 -н 100) 10»;

Gr = 9- 1 0 - ю -г-0,1;

Дбм

= 6 7 , 4

-,100.

 

После подстановки физических величин в критерии

подобия уравнение

(VI—45)

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

в размер­

ную расчетную формулу для определения

времени т ро­

ста единичного кристалла с фиксированным

начальным

размером /„:

 

 

 

 

 

 

 

ДО V 0 , 3

 

AGv0,3

 

 

* = 3 ' 3

4 ,3 ,,0,6 4

= С 1 73 — 0764

-

( V I - 4 6 )

 

'кр а і

 

'кр гкр

 

 

где / К р = / н + Д/ — конечный размер кристалла, м; Д/ — приращение его линейного размера за время т.

К фиксированным

условиям

относятся

постоянный

начальный линейный

размер кристалла

/ н = 1 мм (так

как в критерии подобия F0' должен быть

масштабный

раз,мер, не изменяющийся в начальных

условиях) и по­

стоянная начальная

расчетная

толщина

пограничного

слоя 6 = 0,3 мм.

 

 

 

 

Как показывают опыты и расчеты, изменение этих величин отражается лишь на постоянном множителе уравнения (VI—45), не влияя на показатели степени при критериях. Так, изменение начального размера кри-

14* 211

сталлов в большую сторону сдвигает обобщающую пря­ мую вправо при сохранении ее наклона (рис. VI—21).

Начальный размер кристаллов и Д б у изменяют кон­ станту расчетного уравнения

(VI—46), что можно учесть с по- е>^\ мощью графика на рис. VI—22.

А»1

wo ев

гв

 

 

 

 

 

 

ю

 

 

 

 

 

 

в

 

 

 

 

 

 

 

¥

 

 

 

 

 

 

 

г

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

am am ojагat

г г

ч

т го

to швzooевв

70 7J /В

ff ЗО SfJOff

 

 

 

 

 

 

 

ч щ

Рис. VI—21. Рост

 

кристаллов

при

Рис. VI—22.

Константы

различных

начальных

размерах:

расчетного

уравнения

/ — 0,4 мм;

2 — 1,12

мм;

3 — 3,1

мм

для времени

кристалли­

(И.

А.

Белоконь).

 

зации.

Таким образом, проникновение в область микрокине­ тики кристаллизации также привело к обнадеживающим результатам, которые со временем дадут «сопряжение» уравнения (VI45) с кинетической моделью (VI15).

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

КИНЕТИКА ТЕПЛООБМЕНА ПРИ КРИСТАЛЛИЗАЦИИ

Кинетика теплообмена в вакуум-аппаратах

Основное критериальное уравнение теплоотдачи к ут­ фелю (см. главу V)

Nu = CKnf P r - m exp (— я к , ) (VI—47)

не содержит в явном виде времени процесса т.

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

Получив обобщенные кинетические уравнения кри­ сталлизации типа (VI15), пригодные в качестве ма-

тематической модели для вакуум-аппаратов и кристалли­ заторов в области массообмена, естественно искать та­ кие же уравнения для теплообмена. Процессы эти вза­ имно обусловлены и связаны общими закономерностями. Разница между нимд в том, что при кристаллизации ос­ новная выходная величина — Кр растет по экспоненци­ альному закону, асимптотически приближаясь к мас­ штабному пределу Крт, а теплоотдача к утфелю падает соответственно росту кристаллов, начиная с некоторой масштабной величины Nun. Это падение отражено мно­ жителем ехр(—яК р), который обращается в единицу при Кр = 0, превращая уравнение (VI—47) в уравнение теп­ лоотдачи к пересыщенному раствору без твердой фазы. Учитывая экспоненциальную зависимость содержания кристаллов от времени процесса, естественно предполо­ жить экспоненциальный закон падения теплоотдачи во времени, по крайней мере в течение основного периода уваривания — от заводки кристаллов до их предельного содержания /Ср = 50-=-55%.

В качестве начальной масштабной величины следует принять

Nu„ = CKnh Pv^m ,

( V I - 4 8 )

где индексом «О» обозначены величины, относящиеся к моменту времени, когда сформировалось состояние меж­ кристального раствора. Это состояние соответствует при­ мерно моменту времени то.5 от начала видимого роста кристаллов (см. рис. VI—3) и во всяком случае нахо­ дится в интервале Ат. Вести отсчет от максимального пересыщения в момент заводки кристаллов нет смысла, так как в этот момент теплоотдача не достигает макси­ мального масштабного значения. Некоторое незначи­ тельное содержание кристаллов в интервале Ат (по­ рядка 5—7%) не может заметно повлиять на числовое значение Nu0 в силу характера зависимости пь-р от /Ср.

В качестве

конечной предельной

величины

следует

принять такое

наименьшее значение

критерия Нуссельта

Nu M H H , которое

соответствует переломной

точке

кривой

N U = / ( T ) , возникающей в определенных

условиях при

соприкосновении кристаллов. Теоретическое критическое объемное содержание сферических тел одинакового раз­ мера, вступающих в свободное соприкосновение, состав­ ляет около 53%, чему соответствует в утфеле I продук-

та общая концентрация сухих веществ около 90—91%. В полидисперсной системе предельное содержание кри­ сталлов колеблется в соответствии с фракционным со­ ставом.

Рис. VI23. Схема измене­ ния интенсивности теплоот­ дачи во времени уваривания.

Рис. VI24. Изменение относи­ тельного теплового потока во времени цикла варки I продук­ та (Н. А. Шевандин).

yswoo

— "А

 

Рис. VI25. Изменение

ко­

 

 

 

Г*

 

эффициента

теплопередачи

1 500

1

 

в

процессе нормальной

вар­

Г

п

 

ки

утфеля

I продукта

гоо

 

 

 

(Н. А. Шевандин).

 

50

/50

гоо тмим

 

 

 

 

 

 

 

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

Nu =

NuM „„ + (Nu„ — Nu„„H ) exp

т

(VI—49)

Т

 

 

 

где т — время

от начала подкачки или, приближенно, время от за­

водки

кристаллов;

 

 

G — постоянная времени затухающего процесса.

 

Формфактор и ^ 1 , так как на кривых в пределах основного периода варки не ожидается точек перегиба (рис. VI—23). В первом приближении для I продукта примем п—\. Из-за сложной взаимосвязи тепло- и массообмена значения постоянной времени 0 в уравнениях (VI—49) и (VI—15) могут не совпадать.

Поиск в указанном направлении начат Н. А. Шевандиным [268]. Обработка его данных, опытов М. Л. Вайс­ мана и А. П. Соколова для I продукта подтвердила эк­ споненциальный характер падения тепловых потоков

(рис. VI24) и коэффициента теплопередачи (рис. V I — 25) в течение второго периода варки при непрерывных подкачках и нормативных условиях варки. Спрямление усредняющих линий в полулогарифмических координа­ тах позволяет найти приближенные формулы для фикси­ рованных условий. Так, при Дб„ = 91-^-92, В = const и At — const коэффициент теплопередачи для утфеля I про­ дукта

 

 

k = k 0 e x p ( - ~ y

 

( V I - 5 0 )

где 0 = 230 мин;

 

 

 

 

 

ko — коэффициент теплопередачи,

определяемый по формулам и

графикам

главы V, для начала второго периода

варки.

Поскольку

в

вакуум-аппаратах

ai~k,

уравнение

(VI—50) приближается к обобщению вида

(VI49);

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

значение 8 экстраполирует зависимость

к нулю вместо &МИН (соответственно NUMHH),

Ч Т О , впро­

чем, не отражается на расчете из-за

небольшого накло­

на обобщающей

прямой на рис. VI25.

 

Экспоненциальное падение

коэффициента

теплопере­

дачи обусловлено в основном

двумя

факторами: ростом

кристаллов и накоплением несахаров в межкристальном растворе.

Кинетика теплообмена в кристаллизаторах с охлаждением

Основные данные по кинетике теплообмена

в мешал­

ках-кристаллизаторах

получены и обобщены

И. С. Гу-

лым

[47], частично совместно с другими

авторами [52,

53,

55, 56, 57, 58, 59 и др . ] . В принципе

исходные рас­

суждения

те же, что и для вакуум-аппаратов. Основное

уравнение

экспоненты

 

 

 

 

 

* =

* ш ш +

( V — £ м и н ) е х р ^ — — ^

 

(VI—51)

из-за недостатка

данных заменяется

приближенным

уравнением

 

 

 

 

 

 

 

* = * 0 е х р ( - - g - ) -

 

( V I - 5 2 )

Постоянную времени 0 находят при обобщении опыт­ ных данных в полулогарифмических координатах; вели­ чину k0 рассчитывают по полученным ранее уравнениям

(см. главу V ) . Для расчетов по формуле (VI—52) реко­ мендуется среднее значение 8Ч ас находить по эмпириче­ ской формуле:

6 = 7,6 . 10 -3 ехр ( 0 , 1 С В М 0 ) .

( V I - 5 3 )

Область применения формул (VI—52) и (VI—53): начальная концентрация маточного раствора СВм0= = 80-^-90%, окружная скорость вращения мешалки в кристаллизаторах ш = 100ч-110 м/ч.

 

Рис. VI—26. Обобще­

-^x-Qsa-L

ние

опытных

данных

по кинетике

теплооб­

 

 

мена

в мешалках с ох­

 

лаждением (И. С. Гу-

 

лый,

Н . А . Шевандин)

о.)

иг

O,J

до

asi

 

Зависимость

0 = f(CSM O ),

разумеется,

упрощена; в

действительности

должно

быть 0 = / ( С 5 м О ,

Дбм0, Кр, /у,

П. . .), что предстоит выяснить в будущем.

 

Характер обобщения (VI—52) иллюстрируется рис. VI—26.

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

КИНЕТИКА НАКОПЛЕНИЯ НЕСАХАРОВ

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

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

процесса. Кинетический подход обнаруживается в рабо­ тах последних лет (5, 13, 43, 44, 59, 89, 98, 144, 192 и др.]. .

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

детально

исследована

в

основополагающих

работах

проф. С.

3. Иванова

[89].

В КТИППе также

были по­

ставлены работы по изучению разложения сахарозы при кристаллизации, в частности работа Л. И. Требина [236].

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

створа

в данный

момент времени — его

концентрации,

пересыщения и доброкачественности.

 

 

 

Первым приближением является установление рас­

четных

значений

Дб

маточного раствора

в процессе

его

насыщения несахарами. Обобщив

имеющиеся данные,

Н. А. Шевандин

[267]

нашел, что

в

вакуум-аппаратах

 

 

Дбм

=

Дб0 ехр

 

,

( V I - 5 4 )

где Дбм

— текущее

значение доброкачественности

межкристального

Дб0

раствора в момент времени т;

 

 

 

 

— начальная доброкачественность

межкристального

ра­

 

створа;

Дбо^Дбу.

 

 

 

 

Значения постоянной времени 9і в мин: для рафинад­ ного утфеля — 5900; I продукта—1500; I I продукта 1300.

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

Сухие вещества утфеля складываются из сахара и не­ сахаров: СВу = Сху + Нсу.

В процессе уваривания сухие вещества накопляются в объеме вакуум-аппарата периодического действия за счет добавления сахара Сх2 второго потока (подкачек)

к сахару Сх\ и несахару Нс{ первого потока (начального набора). В течение периода сгущения ті происходят тер­ мическое разложение и химические превращения части сахара первого потока; во времени кристаллизации тК р этот процесс продолжается, и в конце уваривания при­ ращение несахаров составляет АНс к введенным с ис­ ходным сиропом и подкачками.

Рис. VI27. Схема кинетики накопления несахаров при уваривании утфеля.

Поэтому

Нсу = Нсг + Неї + АНс.

За счет накопления несахаров количество сахара в утфелях убывает:

Сху = Cxi + Схъ — АНс;

при этом накопившиеся несахара удерживают в раство­

ре дополнительное количество сахара АСх,

за счет чего

снижается выход сахара в виде кристаллов

С*к р .

Оба процесса концентрации

несахаров

(вследствие

подпитки и кристаллизации и

вследствие

разложения

сахара) определяют закономерность текущего нараста­

ния

несахаров, начиная от

масштабного значения Нсй

для

начала кристаллизации;

эта закономерность харак­

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

нает падать доброкачественность межкристального ра­ створа по закону (VI—54).

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

о

і

г

з

*

 

 

 

 

г

Рис VI—28. Обобщение данных по

накоплению

несахаров при уваривании I

продукта

(В. Д. По­

пов и Е. А. Недужко).

 

 

В связи с тем что в

продуктах с

малым содержа­

нием несахаров постоянные 6 велики и определяются неточно, нами с Е. А. Недужко были предприняты по­ пытки перехода к другим обобщенным координатам, по

аналогии

с кинетическим

уравнением кристаллизации

(VI—27).

Для обобщения

была принята кинетическая

модель с формфактором я = 1 :

Ї - ( Т ) . - [ ( Ї ) . - ( Ї ) . 1 Ч Ч ) '

— соответственно отношения несахара к

воде маточного раствора текущее, мак­ симальное (соответствующее Крт), и на­ чальное (в момент заводки кристаллов).

Результат обобщения показан на рис. VI—28 со сред­ ней постоянной времени 02 = 650 мин.

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