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

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

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

иТиванйём. Математическая модель такбго реактора с гомогенным типом процесса записывается в виде:

 

^ , i ^ . ( C „ -

c „ , ;

- f _ . ,

< V H - „

где

Усек секундная

подача

(расход)

потока, м3 /с;

 

V — объем аппарата, м3 ;

 

 

Свх,

Свах и С — концентрации вещества на

входе,

выходе и в

любой точке аппарата, кг/м3 ;

I- —длина реакционного пространства, м.

Модель идеального вытеснения соответствует порш­ невому течению потока без продольного перемешивания при одинаковом времени пребывания всех частиц в ап­ парате x=V/VCeK- При таком течении для элемента по­ тока и для всего потока справедливо уравнение

дх дх oL

где w — линейная скорость потока в осевом направлении проточной части аппарата х, м/с.

Этой модели соответствуют реальные трубчатые аппа­ раты проточного типа с большой длиной труб при тур­ булентном течении потока.

Степень приближения системы к предельным моде­ лям характеризуется критерием Боденштейна:

 

 

wL

 

іде

w — средняя линейная скорость потока, м/с;

 

DL

— коэффициент продольного перемешивания, м2 /с.

Для систем полного вытеснения ,DL-^0

и Во - ^оо; си­

стемы

полного смешения характеризуются значениями

0L-^oo

и

Во->-0. Для промежуточных

типов систем

0<DL<<x>

и о о > В о > 0 .

 

Диффузионная модель может быть однопараметрической (с продольным перемешиванием) и двухпараметрической [с продольным и поперечным (радиальным) пе­ ремешиванием]. Она ближе к реальному процессу в ап­ паратах типа, описанного в работе [269]. Основой для однопараметрической модели служит модель вытеснения, осложненная обратным перемешиванием; характерным параметром является коэффициент продольного переме­ шивания DL м2/с. Уравнение модели:

дС

 

д*С

дС

—— = D,

w .

(VII—3)

дх

1

дх*

дх

'

При отсутствии продольного перемешивания DL 0 и уравнение (VII—3) превращается в уравнение (VII—2) для модели идеального вытеснения. При опи­

сании

двухпараметрической

модели вводят

дополни­

тельно

 

коэффициент

радиального

перемешивания D R .

Для

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

типа с циркуляцией в_

каждой

секции, описанных

в

работах [25,

26, 27, 84

и др.],

применима

я ч е е ч н а я

м о д е л ь ,

состоящая

из последовательно соединенных ячеек с идеальным пе­

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

между

ячейками. Характерным параметром

модели

является

число ячеек т; если

т—\, ячеечная

модель переходит

в модель идеального

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

при т-^оо ячееч­

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

 

 

1

dC

VCPK

 

 

 

 

• 7 7

= ~1Г

<с'-1

- с<>•

<vlI-4>

 

 

т

ах

V

 

 

 

где

г ' = 1 , 2, 3,

т.

 

 

 

 

 

 

Параметры ячеечной и диффузионной модели связа­

ны

отношением,

достаточно точным при т > 1 0 :

 

 

 

 

Lw

Во

 

 

 

Кроме упомянутых,

возможны

различные

комбиниро­

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

потоков — основного

и

питающего, что отразится на ки­

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

уравнения

превращения

типа

(VI—15). Это было показано в главе V I при переходе

от

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

 

— к координате [l/Le)

для

не-

 

 

 

6

 

 

 

I

 

 

 

 

 

 

 

прерывно действующих

аппаратов,

в которых т =

— ;

 

 

 

 

 

 

wk

Lo—wkQ

— линейная

постоянная;

w — скорость цирку­

ляции;

k — кратность

обмена раствора между

секциями

аппарата.

 

 

 

 

 

 

16 В. Д. Попов

241

Для одной секции приближенное значение /~г|)0 , где гро= отношение объема секции к ее поверхности

нагрева.

Таким образом, в числе критериев, характеризующих подобие процессов превращения при кристаллизации,

будут критерии типа Во, В о д = тг"» т/8, I/LQ.

Другими кинетическими критериями могут быть [69]: критерий Дамкелера

 

 

 

их

 

 

 

 

 

Da =

,

(VII—6)

где

и — объемная скорость процесса

превращения,

к г / ( м 3 - с ) ;

 

 

критерий

Маргулиса

 

 

 

 

 

 

Ma' ~

,

(VII—7)

 

 

 

w.

 

 

 

где

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

м/с;

 

 

 

критерий равновесности

Дьяконова

 

 

 

 

Ра = J —

,

(VII—8)

 

 

1 X

 

 

где

х— степень фазового перехода в долях единицы; в нашем

слу­

 

чае

х=Кр/Крт.

 

 

 

 

 

Для описания процессов

на

гранях кристаллов

[13,

69, 98, 102, 111, 144 и др.] используются критерии диф­

фузионного

подобия:

 

 

6/

wl

v

Dx

N u , = J ^ .

Ре' = — ;

P r ' = — ; Fo' =

— и др . ,

где D — коэффициент диффузии;

 

/ — определяющий размер.

 

В отдельных случаях

может оказаться

удобным кри­

терий, включающий общую скорость процесса и, проте­ кающего в диффузионной области [132]:

 

и/2

 

 

K D = ~ -

( V . I - 9 )

Кроме упомянутых, применимы и другие критерии

подобия

[69, 257, 272].

 

Среди различных вариантов задач по расчету про­

цессов

превращения основной задачей,

по-видимому,

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

Ко = / ( К і , К„ Кз, . . . )•

(VII—10)

В нашем случае таким критерием явился относитель­ ный выход кристаллов Кр/Кр?, выраженный в функции

безразмерного времени процесса — .

8

Критерии подобия процессов теплообмена

Эти критерии общеизвестны (они перечислены в гла­ ве V ) . Это различные критерии теплового напряжения, содержащие тепловой поток:

v

У

п

<1<1р

п

ЧСІи

К / =

„ , , ;

Р е и =

— —

;

 

Re„ - — — и др.

'

rp"d0f

 

гр"а

 

 

Ф v

Физические параметры среды отражены рядом крите­ риев: Рг, Кр и т. п.

Для установившегося процесса теплоотдачи искомым является критерий Nu, выражаемый в функции от про­

чих критериев

подобия. Специфичным

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

к утфелю является

новый

критерий

твердой

фазы

Кр/(100Кр),

связывающий

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

с про­

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

Кр=!(т).

Одной из форм связи может быть уравнение кинети­

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

вида:

 

 

 

 

 

Nu =

N u 0 e x p ^ — - ^ - j .

 

(VII—11)

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

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

16* 243

Предложение выдвинул Г. М. Щеголев [274] для описа­ ния теплопередачи в испарителях с паровым обогревом. Учитывая влияние на тепловой поток факторов, дейст­ вующих по обе стороны стенки, используя уравнение Коши и анализ размерностей, он предложил критери­ альную зависимость в следующем виде:

<7макс<* , „ Я Ж < Ш pd2 I At Vn \

где

^макс максимальное значение теплового потока;

 

d и I — диаметр

и длина кипятильной

трубы;

 

"Уж

И Я Ж вязкость

и

теплопроводность

кипящей

жидкости;

Vn

и Уж — количество

пара на выходе

и

жидкости

на входе

 

в трубу;

 

 

 

 

 

k>KdAt

соответствующий

критерию Кутателадзе.

 

критерии,

rp"Vm

Прочие обозначения общеизвестны (см. главу V ) . Аналогичные выражения даны в работе [274] для на­ чального и конечного потоков тепла.

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

kl

(VII—13)

Ki = — ,

л э к в

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

Хэкв эквивалентный коэффициент теплопроводности;

/— определяющий размер системы.

Вотличие от критерия Нуссельта для случая тепло­ отдачи, в котором определяющим размером служил от­ рывной диаметр или критический радиус [91] парового пузыря, в критерии Ki должен быть масштабный размер,

характерный

для реакционного

объема,

например, / =

= — м3 2

— удельный объем

утфеля

на 1 м2 поверх-

F

пости нагрева — величина, обратная удельной поверхно-

F

сти нагрева ф 0 = — .

Для определенности следует принять неизменный объем начального сиропа до подпитки V0, т. е. l=Vo/F.

Под величиной Яэкв понимают эквивалентный коэф­ фициент теплопроводности условной трехслойной стенки толщиной

б ^ б і + бст + ба,

(VII—14)

где б і = Х с т / а і — р а с с т о я н и е

от стенки

до направляющей

точки со

стороны пара (приграничных

условиях 3-го рода);

бг = Яст / а2 — расстояние

от стенки

до направляющей,

точки со

стороны утфеля;

 

 

 

Хот и* бет теплопроводность и толщина

материала

поверх­

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

 

 

 

При указанных условиях

бі + бет + 62

ЭК В '

Si

SCT

б 2

 

^

Л к

Л с т

Лу

1

1

 

 

После подстановки и преобразований уравнение для эквивалентного коэффициента теплопроводности прини­ мает вид:

 

Лет , s

^ст|

 

1

, J ,

1

 

 

я

« і

 

«а

_

В и

 

B i 3

(VII—15)

^•ст .

Ост

, \ т

 

1

. 1

,

1

 

ОІІАК

Я с т

с%2^у

 

А к Вії

\'СТ

 

 

 

 

аі&с т

 

 

 

 

 

 

 

 

где Віі =

критерий

Био для поверхности

стенки

со стороны

 

 

пара;

 

 

 

 

 

 

В12

= ~т

— то же, со стороны утфеля;

 

 

 

Х,(

и Я У коэффициенты

теплопроводности

конденсата и ут­

 

 

феля.

 

 

 

 

 

 

Д ля конкретных условий варки уравнение (VII—15) упрощается. Следовательно, обработку опытов можно предсказать в форме

 

kl

ql

 

я » , . . . ) •

(VII—16)

К і

=

= / ( Я і .

 

^экв

^-экв

 

 

 

Откуда искомый тепловой поток

 

 

 

< 7 > - ^ у 2 - А / / ( я 1 ,

я , , . .

. ) ,

(VII—17)

где rti, Я2 — определяющие

критерии

процесса.

 

Целесообразность представления результатов опытов в такой форме подтверждается построением расчетных диаграмм теплообмена даже в том случае, когда Ігфаг [159], Преимущество формул типа (VII—16) по сравне-

нию с эмпирическими формулами

(см. главу

V)

состоит

в более полном и явном учете действующих

факторов,

отраженных граничными условиями 3-го рода.

 

Необходимой предпосылкой

обобщения

в

форме

(VII16) является рассмотрение определенного сочета­ ния процессов теплоотдачи со стороны пара и со стороны утфеля, которое однозначно при фиксированной (подоб­ ной) структуре потока утфеля, подобном распределении

тепловой нагрузки по высоте трубы и оптимальных

запа­

сах водосодержания [233]. При этих условиях

предпола­

гаемая критериальная

система,

может выглядеть

так:

 

Ki = / ( K / ,

Pa, Arn,

N ,

. . . ),

(VII—18)

где дополнительно Am=E/RT

 

критерий

Аррениуса;

 

 

Л' мольная

доля сахарозы

в растворе;

 

 

Е —энергия

активации;

 

 

 

 

 

 

R — универсальная газовая

постоянная.

 

 

В этом уравнении критерий Ра заменяет критерий твердой фазы яК р, а критерии Агп'и N появляются вме­ сто критерия Прандтля, что обосновано способом обоб­ щения в мольной системе координат, описанным в гла­ ве I I .

Уравнение (VII18) отражает стационарные условия теплообмена при квантовании процесса по времени.

Критерии гидродинамического подобия

При обработке опытных данных применимы извест­ ные критерии гидродинамического подобия Рейнольдса

Re= — , Фруда F r = —

, Галилея Ga= — , Архимеда

V

gl

V2

Ar = Ga — , Эйлера

Eu=

и др. Специфичными для

Р

 

рвд2

движения утфеля будут обобщенный критерий Рейнольд­ са, критерий Бинґама, критерий твердой фазы (см. гла­ ву IV) и некоторые другие.

Постановка задачи [162, 164, 165, 166]

заключается

в совместном

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

следующих

уравнений.

а) Уравнения

теплообмена

для статики

и кинетики

* процесса:

Nu=7_(Ki, Кз. . . .): N u = N u » е х Р ( — у - ) - (VII—19)

б) Аналогичные уравнения массообмена (скорости кристаллизации) имеют вид:

 

 

 

Г т

1

N u '

= / e ( K | , К'2

, • • . ) ; АГр =

/ С р г е х р

j"MF (т)

dx

 

 

 

L

о

(VII—20)

 

 

 

 

 

в)

Уравнения

гидравлических

сопротивлений

конту­

ра, при которых обеспечивается оптимальный пьезомет­

рический

уровень:

 

Ей = /з ( К " , К 2 '

, .

. . ) ; Eu = Eu0cp ( K J ' , К'2', . . . , т) ,

 

 

 

(VII—21)

где K i , Кг,

К р К

2 ,

... — определяющие критерии различных про­

 

 

 

цессов.

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

Укажем на неиспользованные возможности аналити­

ческого расчета

оптимального пьезометрического

уровня

и опрокинутой

циркуляции, заложенные в трудах

акаде­

мика В. И. Толубинского [233].

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

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

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

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

расчеты на каком-то этапе неизбежно осложняются не­

обходимостью

согласования

кинетических

коэффициен­

тов со скоростями превращения.

 

 

При дальнейшем

развитии достигнутых

результатов

в стройную теорию рабочих процессов

кристаллизации

необходимо учесть следующие

соображения.

 

А. В. Лыков

[122,

132 и др.] давно

доказал, что про­

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

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

h = L 1 1 X 1 + L 1 2 X 2 + Ll3X3

+ . . . ,

( V I I - 2 2 )

где X,, Х2, ... •

термодинамические силы;

 

коэф­

Ьи-

диагональный

(с одинаковыми индексами)

 

фициент, отражающий линейную зависимость дан­

 

ного потока от основной термодинамической

силы

 

(например, от

градиента

температур

для

тепло­

 

вого потока);

 

 

 

 

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

Для другого вида энергии или вещества

/2 =

LuXi -\- L 2 2 X 2 -f- L23X3 т •

• •

(VII—23)

При наличии

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

сил

составляется

і уравнений Онзагера с і членами в правой части; это будут термодинамические уравнения движения.

Согласно второму принципу взаимности

перекрестные

коэффициенты равны: L 1 2 = L 2 i; L 1 3 — L 3 1 и

т. д. Это оз­

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

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

Кроме того, большие перспективы совместного рас­ смотрения процессов переноса тепла, вещества и коли­ чества движения при фазовых переходах содержатся в обобщенной теории переноса, созданной Е. В. Толубинским [234].

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

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

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

ренной

воды, "выход утфеля

и т. п. Тогда все

расходные

величины будут измерены одним масштабом

(например,

G

G

G

G

 

 

 

 

— и т. п.).

Определенные затруднения

G H

ОН ОН ОН

 

 

вызовет выбор масштаба для оси времени.

 

 

Относительное время т/тц применимо только для ва­

рок

данного

продукта, мало отличающихся

активным

временем цикла тц . В других случаях следует перейти к безразмерному времени превращения т/9. Существенным

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