
книги из ГПНТБ / Попов В.Д. Основы теории тепло- и массообмена при кристаллизации сахарозы
.pdfиТиванйём. Математическая модель такбго реактора с гомогенным типом процесса записывается в виде:
|
^ , 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] для на чального и конечного потоков тепла.
По аналогии с уравнением (VII—12) в искомый кри терий надо включить коэффициент теплопередачи либо тепловой поток. Если исходить из аналогии с критерием Нуссельта, то искомым будет критерий Кирпичева:
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-го рода. |
|
||
Необходимой предпосылкой |
обобщения |
в |
форме |
(VII—16) является рассмотрение определенного сочета ния процессов теплоотдачи со стороны пара и со стороны утфеля, которое однозначно при фиксированной (подоб ной) структуре потока утфеля, подобном распределении
тепловой нагрузки по высоте трубы и оптимальных |
запа |
||||||
сах водосодержания [233]. При этих условиях |
предпола |
||||||
гаемая критериальная |
система, |
может выглядеть |
так: |
||||
|
Ki = / ( K / , |
Pa, Arn, |
N , |
. . . ), |
(VII—18) |
||
где дополнительно Am=E/RT— |
|
критерий |
Аррениуса; |
|
|
||
Л' — мольная |
доля сахарозы |
в растворе; |
|
|
|||
Е —энергия |
активации; |
|
|
|
|
|
|
R — универсальная газовая |
постоянная. |
|
|
В этом уравнении критерий Ра заменяет критерий твердой фазы яК р, а критерии Агп'и N появляются вме сто критерия Прандтля, что обосновано способом обоб щения в мольной системе координат, описанным в гла ве I I .
Уравнение (VII—18) отражает стационарные условия теплообмена при квантовании процесса по времени.
Критерии гидродинамического подобия
При обработке опытных данных применимы извест ные критерии гидродинамического подобия Рейнольдса
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. Существенным