![](/user_photo/_userpic.png)
книги из ГПНТБ / Попов В.Д. Основы теории тепло- и массообмена при кристаллизации сахарозы
.pdfкристаллов не смогут обобщить данные по теплоотдаче к утфелю, так как область их действия ограничена рас творами (пунктиры на рис. 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 А |
|
различных |
продуктов, |
про |
|
|||
веденных |
в |
разное |
время |
|
.і |
|
несколькими |
авторами |
и в |
|
|
||
отличающихся |
условиях |
Рис. |
V—31 б. Обобщение дан |
|||
(размеры |
и |
конструкции |
ных |
о теплоотдаче к утфелю |
вакуум-аппаратов, |
режим |
с |
пониженной |
степенью вли |
питания растворами, |
пере |
|
яния числа |
Прандтля. |
|
|
|
||
сыщение и др.). Обобщение |
всех проведенных опытов |
|||
дает сплошную полосу опытных |
точек [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]. Уравнение (V—58), естественно, дает те же чис
ловые |
значения |
аг, что |
и формулы для утфеля, так как |
оно основано на |
одних |
и тех же кривых функции физи |
|
ческих |
параметров |
утфеля. Уточнение уравнения |
(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]. Линии 1—10 соответствуют повы шающейся скорости искусственной циркуляции в трубах от 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