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

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

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

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

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

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

пределения Гаусса. Поэтому, имея в виду

связь СВУ и

Кр,

содержание кристаллов следует вводить в формулы

в виде фактора твердой фазы:

 

 

 

 

 

 

х = я 0 я 2 р ,

 

 

( V - 5 2 )

где

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

константа,

усредняющая

влияние фрак­

 

ционного состава

кристаллов

и

другие факторы взаимосвя­

 

зи

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

с суммарным

эффектом

действия твер­

 

дой

фазы.

 

 

 

 

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

а =

0 ехр (—х), .

(V—53)

где а и а 0 — коэффициенты

теплоотдачи к сиропу

при содержании

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

пересыщении.

 

 

Такой же должна быть связь и между

безразмерны­

ми коэффициентами теплоотдачи:

 

Nu =

Nu0 ехр (— х).

(V—54)

Для статических условий в факторе х будут фигури­ ровать значения Кр для данного момента времени; вре­

мя

же

т

в

явном виде не проявляется из-за усреднения

Кр

на

отрезке конечной

разности

Ат. При

желании

можно

 

выразить Kp = f(x)

и превратить

уравнение

(V—54)

в

кинетическое

уравнение

теплообмена такого

вида:

 

 

 

 

 

 

 

 

 

 

 

Nu = Nu0

ехр | -

(y)*]

.

( V - 5 5 )

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

п — формфактор кинетической кривой теплоотдачи.

гг
' и * з р я

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

Рис. V—31 а. Обобщение результатов опытов с применением исходной формулы В. И.Толубинского.

водства оказались недоста-

го

обобщения

в форме

(V—54),

и обобщение про­

 

изведено

иным

способом

 

(см.

ниже).

 

Ю

у '

^ —

г

На

рис.

V—31а

по­

8

VU

казан

результат обобщения

 

в системе критериев (V—48)

 

У У J V

лишь

для

нескольких

варок

 

Аи 5 А

различных

продуктов,

про­

 

веденных

в

разное

время

 

несколькими

авторами

и в

 

 

отличающихся

условиях

Рис.

V31 б. Обобщение дан­

(размеры

и

конструкции

ных

о теплоотдаче к утфелю

вакуум-аппаратов,

режим

с

пониженной

степенью вли­

питания растворами,

пере­

 

яния числа

Прандтля.

 

 

 

сыщение и др.). Обобщение

всех проведенных опытов

дает сплошную полосу опытных

точек [144].

На рис. V—316 обобщены: опыты И. А. Кухаренко по варке рафинада на Дерюгинском сахарном заводе в

змеевиковом аппарате при высоких тепловых потоках (1929 г.); опыты автора на конической поверхности на­ грева ( I I рафинад, Одесса, 1950 г.) и в промышленном аппарате системы Маликова (1951); промышленные опыты ЦИНСа в аппаратах системы «Проток» и с труб­ ной подвесной камерой (М. Л. Вайсмана, 1952) и др.

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

Nu Рг° - 3 ехр ( 0 , 2 5 л 2 р ) = 54К^'6

(V—56)

(при d0 /=280 м/ч для воды при атмосферном давлении);

N u P r ° ' 2 e x p ( 0 , 2 5 n 2 p ) = 3 5 , 7 К ° > 6

( V - 5 7 )

(при d0f = 560 м/ч при тех же условиях).

Новые данные о степени влияния числа Прандтля и

частоты отрыва паровых пузырей [105] вызвали

попыт­

ку дальнейшего изменения

констант

уравнения

(V—57)

с целью уточнения расчета

аг к кипящему утфелю [133,

134].

 

 

 

Обработка имеющихся

и новых

опытных

данных

привела нас к следующей формуле:

 

 

N u P r 0 , I e x p ( 0 , 2 5 я 2 р ) = 2 1 К ^ 6 .

(V—58)

Формула получена в основном по данным промыш­ ленных испытаний вакуум-аппаратов I продукта, охва­ тывает диапазон чисел Прандтля маточного раствора от 250 до 2100, содержание кристаллов до 54%, учитывает исправленные значения теплопроводности раствора [121, 224, 264] и уменьшенную степень влияния вязкости [105]. Уравнение (V58), естественно, дает те же чис­

ловые

значения

аг, что

и формулы для утфеля, так как

оно основано на

одних

и тех же кривых функции физи­

ческих

параметров

утфеля. Уточнение уравнения

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

Аналогичны результаты обобщения тех же данных

в других

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

(V—31) и (V—38) с до­

бавлением яК р. Как

видно

из формулы

(V—31),

кри­

терий

Рг

выделен

из исходного уравнения (V—28) с

целью

удобства расчета и

сопоставления

показателей

степени. Обобщение

в общем

удовлетворительное

[1-56,

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

С таким же результатом получено обобщение и в си­ стеме (V—38) [163], что приводит к мысли о возмож­ ности обобщения данных по теплоотдаче к утфелю в простейшей системе безразмерных координат вследствие относительного постоянства пересыщения при варке дан­ ного продукта и определенной закономерности накопле­ ния несахаров в маточном растворе. Системе (V—38) соответствует уравнение

N u P r ° ' 3 5 e x p ( 0 , 1 5 я 2 р ) = 8 , 2 - 1 0 - з К ° - 6 .

(V—59)

Одновременно с обобщением данных прямых опытов проверяли выведенные ранее эмпирические формулы Д Л Я общего коэффициента теплопередачи путем выделения значений аг для типовых условий варки и нанесения рас­ четных точек на совместный график с опытными точ­ ками. Расчетные точки оказались в непосредственной близости к усредняющим прямым на всех графиках (например, на рис. V—31). Этого следовало ожидать, так как эмпирические формулы получены из тех же опытных данных, что и критериальные формулы, путем их усреднения и обобщения другим способом. Можно утверждать, что эти формулы в своих областях не менее надежны, чем критериальные формулы, хотя и не обла­ дают широтой диапазона применения.

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

Первая попытка, основанная на предпосылках А. И. Рычкова и Б. Г. Попова, доведена В. Н. Горохом

(ВНИИСП)

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

пригодного

для утфеля I рафинада при температурах

насыщения

60—85° С:

 

 

 

Е, = 0 , 2 5 9 К 7 - 0 , 3 8 К ° ' 6 е х р ( - 0 , 2 я 2 р ) ,

( V - 6 0 )

где Es =

—эбулиоскопический критерий;

 

Лф

 

Кд = ~zr— температурный фактор (критерий).

 

Уравнение (V—60) проверено на стендовом вакуумаппарате непрерывного действия.

Вторая попытка обобщения произведена Е. А. Недужко [133] с использованием данных Цой у Сека по критическим тепловым нагрузкам qKV для чистых сахар-

ко

 

 

 

 

 

 

 

 

 

 

 

 

Ко,0 а

 

D

г,

 

 

-

 

 

 

 

 

OS

 

0

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

о

 

 

d

о

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

мо" о

ла•п

f o

%

 

 

 

 

0.3

 

 

 

 

 

?оPry

 

 

 

 

 

 

 

 

 

 

 

о

 

 

^ о о

 

 

Э

 

 

 

 

 

 

 

 

 

**

 

 

 

 

 

 

с

 

 

 

У

 

 

 

 

 

 

 

 

 

>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

>?Ч

 

 

 

 

 

 

 

 

 

 

 

 

£<

 

 

к ,

 

 

 

 

 

 

 

 

 

" Л о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О

 

A3

0,5

«7

 

/,0

 

 

20

3,0 • -

5,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

• л

 

Рис. V—32. Обобщение данных о теплоотдаче при естественном кипе­ нии утфелей рафинадных, I и I I продуктов и при перемешивании (Е. А. Недужко).

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

 

Е.ехр [ 9 , 6 (

""де*6"

) 4 ] - 0 , 0 6 3 К ° - 6

^ Q 7 0 ' 4 3

К р ' 3 6

X

 

 

Х е х р ( - 0 , 4 5 я 2 р ) ,

 

 

 

( V - 6 1 )

где

Qs =q/qKr> — безразмерная

удельная

тепловая

нагрузка.

 

 

Обобщение проверено обработкой данных на ЭЦВМ

«Урал-4». Результаты представлены

на

рис. V—32. По

вертикальной

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

комп­

лекса Ко:

 

 

г

/ 100 Д б м \

*1

 

 

 

 

 

 

 

Es exp ( 0 , 4 5 я к р > х р | 9 , б (

 

д б м

)

\

( V - 6 2 )

 

К о = -

 

 

К 0 . 6 4 К 0 , 3 6

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая большую разницу в физических парамет­

рах утфелей

различных

 

продуктов,

обобщение в полосе

с

отклонением ± 2 0 %

 

следует

признать,

успешным,

а

цель — достигнутой.

 

Основное

значение

формулы

(V—61) в единстве программы расчета на ЭЦВМ для решения последующих задач оптимизации процесса не­ зависимо от вида продукта. Основное преимущество формул типа (V—60) и (V—61) в отсутствии вязкости, как менее надежного расчетного параметра. Недостатки формулы — непригодность для описания микрокинетики и наличие эмпирического множителя, учитывающего ка­ чество продукта (Дб).

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

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

утфельных

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

(см.

главу

I V ) . Однако мы не отрицаем другой

формы

обоб­

щения,

например,

с помощью

эффективных

параметров

утфеля, который рассматриваем как условную сплошную среду. Полагаем, что этот вопрос будет рассмотрен после детального изучения реологических свойств утфелей. Кроме того, возможно использование критериев подобия иного типа, рассмотренных ниже, в главе V I I , в работах' Р. Я. Ладиева [119] и др.

Влияние перемешивания утфеля на теплоотдачу при его кипении

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

Рассмотрим результаты лабораторных и полупро­ мышленных опытов, проведенных в КТИППе по иссле­ дованию влияния перемешивания на теплоотдачу при уваривании утфелей [133, 136].

Как известно из первых опытов Л. С. Стермана и Н. Г. Стюшина, проведенных с водой, и последующих исследований ряда авторов с растворами, общая картина процесса кипения при включении побудителя циркуля­

ции

(насоса или

мешалки) выглядит так, как показано

на

рис. V—33.

Левый график заимствован из работы

Р. Я. Ладиева [119]. Линии 110 соответствуют повы­ шающейся скорости искусственной циркуляции в трубах от 0,6 до 4 м/с, линия 11 — естественной циркуляции. Правый график, построенный по нашим данным [144], соответствует кипению сиропа на горизонтальной плите без перемешивания (наклонная линия) и при включе-

 

 

 

4

5

6

7 8

9K)3q

Рис. V—33. Влияние эффекта

перемешивания

на

теплоотдачу

к

воде

(Р. Я. Ладиев)

и сиропу (В.

Д .

Попов).

 

 

 

 

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

Опыты Е. А. Недужко были поставлены при переме­ шивании утфеля, кипящего на горизонтальной плите, и при установке пропеллерной мешалки в циркуляцион­ ной трубе вакуум-аппарата. Результаты обработаны в критериальной форме и нанесены черными точками на общий график (см. рис. V—32); по вертикальной оси отложены значения комплекса

Е , е х р ( 0 , 4 5 я * р ) ехр

9 , 6 ^

100 -Дбм

у

 

Кох

 

 

Дбм

 

(V—63)

 

 

53

 

W ° , 6 4 K - 0 , 3 6

 

 

 

 

 

J

 

 

 

V

Re„

 

 

 

где Re3KB = Ren + 0,025 Re4 — эквивалентный

критерий

 

Рейнольдса,

выражающий

суммарный

гидродина-

мический эффект от естественной и ис­ кусственной циркуляции;

 

Р е и

 

 

 

 

 

 

 

 

ReH = - j ^ -

—критерий

Рейнольдса

для

процесса

 

 

кипения,

впервые введенный М. А. Ки-

 

nd2

чигиным;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R e 4 = —

— центробежный

критерий

Рейнольдса

 

 

для мешалки

с числом

оборотов п и

 

 

диаметром

лопастей

d;

эмпирический

 

 

множитель

0,025 найден

при

обработ­

 

 

ке опытов.

 

 

 

 

 

Комплекс

Koi отличается от комплекса

Ко наличием

отношения

/ ? І 2 к в \

отражающего

степень

влияния

 

V ReH /

 

 

 

 

 

 

 

перемешивания на теплоотдачу.

 

 

 

 

 

 

То же отношение,

введенное

 

в

уравнение

(V—58),

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

дующих уравнений:

 

 

 

 

для всех

продуктов

 

 

 

 

Е 5 е х р ( 0 , 4 5 л к

2 р ) е х р [ 9 , б ( Ш

~ ^ Ы

) ' ] = 0 , 0 6 3 К ° - 6 4

К » ' 3 6 X

 

 

Re4

 

 

 

 

для утфеля

І и I I рафинада

и I

продукта

 

N u P r 0 ' 1

ехр ( 0 , 2 5 я 2

р ) = г і К » ' 6

^ - ^ ) 0 , 5 3

. ( V - 6 5 )

Обе формулы требуют проверки на крупных вакуумаппаратах в промышленных условиях.

Критические тепловые потоки <7кр Для сахарных рас­ творов приведены ниже.

СВ, %

1 0 6 ? кр-

СВ, %

1 0 6 "кр>

 

Вт/м2

 

Вт/м2

60

0,99

80

0,69

65

0,91

85

0,60

70

0,84

90

0,53

750,76

Вработе Л . Г. Белостоцкого (ВНИИСП) [18] отме­ чено увеличение интенсивности теплопередачи в лабо-

12 В. Д. Попов

177

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

КО Э Ф Ф И Ц И Е Н ТЫ Т Е П Л О П Е Р Е Д А Ч И

ВМЕШАЛКАХ-КРИСТАЛЛИЗАТОРАХ

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

В связи с медленным изменением теплофизических параметров утфеля в процессе его охлаждения в после­ довательно соединенных мешалках-кристаллизаторах их можно рассматривать как аппараты непрерывного дей­ ствия с установившимся процессом. Вместе с тем за время охлаждения параметры процесса существенно из­ меняются; так, при охлаждении межкристального рас­ твора от 80 до 40° С коэффициент кинематической вяз­ кости увеличивается на 1800%- Поэтому так же, как и для вакуум-аппаратов, при расчете мешалок необ­ ходимо применять метод конечных разностей и «кванто­ вать» процесс во времени для интервальных расчетов.

Первое систематическое исследование мешалок-кри­ сталлизаторов с охлаждением как тепло- и массообменных аппаратов выполнил И. С. Гулый [47]. Детально результаты этой работы изложены в ряде публикаций [48, 52, 53, 55, 56, 57, 59 и др . ] . Итогом ее является ком­ плексный расчет кристаллизатора, связывающий охлаж­ дение утфеля с фазовым превращением сахарозы во времени процесса кристаллизации. Для расчета разра­ ботан ряд формул и рекомендаций. Выше были рассмот­ рены температурные поля в кристаллизаторе (глава I V ) . Здесь приведем основную расчетную формулу для ко­ эффициента теплоотдачи.

При чистой поверхности нагрева для кристаллизато­ ров всех типов I , I I и I I I продуктов коэффициент тепло­ отдачи от утфеля к холодной стенке рассчитывается по уравнению

( V - 6 6 )

Физические параметры отнесены к маточному рас­ твору; обобщение произведено в пределах К е ц = 0,24-0,4;

Рг== (8.5-154) 103 (рис. V—34). Поправка на влияние температурного поля в объеме утфеля принята по дан­ ным М. А. Михеева [129]. Изменение Nu во времени обусловлено закономерностью истощения патоки и на­ растания несахаров; благодаря этому в первом прибли-

Рис. V—34. Теплоотдача в мешалках-кристаллиза­ торах с охлаждением (И. С. Гулый).

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

Уравнение (V—66) приводится к размерной формуле вида

 

 

 

a 1 = / l 1 ! B 0 ' , S M l S ,

 

 

(V—67)

где

Ai — функция

физических

параметров

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

 

створа; Ai—f

( С В М . ty,

Дбм);

 

 

 

 

w — окружная скорость мешалкц, м/ч;

 

 

 

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

 

 

 

Формула для Ai

имеет вид:

 

 

 

 

 

Л . = 0,89 •

v ° ' 1 8 J i ° T ' 2 5

с 0 ' 6 3 р ° . 6 3 Я 0 - 3 7

( V - 6 8 )

 

„ 0 , 25 л 0 , 2 5 л , 0 , 2 5

j0,55

 

 

 

 

 

 

 

Р—'

V„„'

а

 

 

 

 

 

"СТ гст

 

 

 

 

где d — диаметр лопасти мешалки.

 

 

 

 

 

Все величины

и

ai = k — в

технической системе еди­

ниц,

сохраненной

в

данном

случае

в

связи с размер-

12* 179

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