книги из ГПНТБ / Зубченко А.В. Новое в кинетике кристаллизации сахара
.pdfРазность химических потенциалов для пересыщенно го раствора определяют соотношением
(j-i — fi2 = k Т In |
С / С |
(25) |
тде С — концентрация пересыщенного сахарного раствора; |
||
С со — концентрация насыщенного сахарного |
раствора. |
|
Обозначим |
|
|
= |
— р . |
|
k Г
тогда
2/3 |
(26) |
|
Зависимость АФ/kT от /г для сахарных растворов с коэффициентом пересыщения а от 1,2 до 1,4 при темпе ратуре 7.0° С показана на рис. 4. По уравнению (26) наш-
|
|
1^, -s—" |
-~ |
1 |
|
|
|
|
3D |
/ |
|
|
Рис. |
4. |
Зависи |
||
|
2 |
|
||||||
|
|
мость |
АФ/kT |
от п |
||||
|
|
|
||||||
.} |
|
|
|
для сахарных |
ра- |
|||
|
|
|
стеоров при |
|
темпе |
|||
3D |
//А |
|
|
|
||||
N |
|
ратуре |
сыщения |
а: |
||||
|
|
|
70°С |
и |
ко |
|||
20 |
|
|
|
эффициентах |
пере |
|||
|
|
|
|
|
|
|
|
|
10 |
|
|
|
/ — 1.20; |
2— |
1,25; |
3 — |
|
|
|
|
|
1,30; |
А — 1.40, |
|||
О г100 200 Ж Ш |
|
|
||||||
500 |
п • |
|
|
|
|
ли значение количества молекул" пк, составляющих кри тический зародыш при данном пересыщении раствора...
Радиусы зародышей критических размеров определили по уравнению
4те |
, |
« к ^ 7 — d . |
(27) |
20
Полученные результаты приведены в табл. 3.
|
|
Т а б л и ц а |
3 |
Значения /г„ и г„ в зависимости от пересыщения |
|||
сахарного |
раствора |
при температуре |
70°С |
К о э ф ф и ц и е нт |
"к |
Гц, мкм |
|
пересыщения |
|
||
1,20 |
675 |
3,86 |
|
1,25 |
372 |
3,16 |
|
1,30 |
236 |
2,72 |
|
1,40 |
ПО |
2,14 |
|
Как видно из табл. 3, с увеличением пересыщения ра створа количество молекул в зародыше критического раз м е р а и его радиус резко уменьшаются.
. Для «пересыщенного раствора функция распределе ния / сначала убывает (рис. 5), проходя через минимум, а затем возрастает. Одновременно с присоединением от-
Рис. 5. Зависи мость функции рас пределения / от /1
при коэффициенте пересыщения а=
дельных молекул к образовавшимся зародышам проис ходит их отрыв до тех пор, пока они не достигнут крити ческого размера и только тогда начинается свободный рост зародышей. Следовательно, скорость образования таких центров пропорциональна минимальному значению функции распределения /к , которое можно вычислить по
уравнению |
|
. |
• • „. |
.Л, = Л в - 1 ' 3 * я Я в 1 |
(2S> |
2L
но так как
Ф = |
н |
е я 2 ' 3 |
= 4г. г- <з, |
|
kT |
|
|
|
|
ТО |
|
|
|
|
|
. |
_ |
° 4 |
|
U |
= Ae |
3 , 1 |
т . |
(29) |
Для нахождения функциональной |
зависимости / от п |
II определения минимального значения fK сделаны расче ты для сахарных растворов при различных коэффициен тах 'Пересыщения и температуре 70° С (табл. 4).
|
|
|
|
|
Т а б л и ц а |
4 |
|
|
|
Зависимость функции распределения зародышей / |
|||||
|
|
от количества |
молекул п для сахарного |
раствора |
|||
|
с коэффициентом |
пересыщения а=1,2 при температуре |
70°С |
||||
|
|
|
2/ |
|
—ДФ/АГ |
|
|
п |
|
|
<? Л " |
|
|
||
|
|
е |
|
|
|||
|
|
|
|
|
|||
•450 |
- 83,0 |
58,6 |
141 |
5S.0 |
9,54-10-м |
95,4 ' |
|
470 |
- 86, 5 |
60,4 |
145 |
58,5 |
8,45-10-2« |
84,5 |
|
520 |
- 95,7 |
64.5 |
155 |
59,3 |
4,77-10-20 |
47,7 |
|
570 |
-105.0 |
68.5 |
165 |
60,0 |
3,02-10-м |
30,2 |
|
€00 |
—110.0 |
71,1 |
171 |
60,8 • |
8,13-10-=' |
8,1 |
|
620 |
— 114,0 |
72,6 |
175 |
01,0 |
7,78-10-2' |
7,7 |
|
650 |
-120,0 |
• 75.0 |
181 |
61,1 |
6,99-10"'' |
6,9 |
|
700 |
— 129,0 |
79,0 |
191 |
62,0 |
3,01 10-2- |
3,0 |
|
720 |
—132,5 |
80,4 |
194 |
61,5 |
6,02-10-27 |
6,0 |
|
800 |
-147,0 |
86,0 |
207 |
60,0 |
3,01-10-2u |
30,1 |
|
850 |
—156,5 |
89,7 |
216 |
59,5 |
4,77-10-2в |
47,7 |
|
•900 |
—165,6 |
93,0 |
224 |
58.5 |
7,78-10-2« |
77,8 |
|
|
|
|
|
|
Т а б л и ц а |
5 |
|
|
Зависимость функции распределения /,; от коэффициента |
||||||
|
пересыщения а растворов сахарозы |
при температуре 70° С |
|||||
а |
С, % масс. |
"к |
о П |
о п 3 |
|
is |
/ к |
3 |
|
|
|||||
|
|
|
|
|
|
|
|
1,20 |
7S.58 |
675 |
62,0 |
1,19-10-2' . |
1,64-Ю-о -5,8125 |
||
1,22 |
79,85 |
522 |
52,0 |
2,61-Ю--"3 |
3,65.10-' • -1,4370 |
||
1,23 |
79,98 |
470 |
48,5 |
7.82-10-32 |
1,095 |
0,0390 |
|
1,24 |
80,10 |
412 |
44,6 |
4,27-10-2» |
5.98-10 |
1.7767 |
|
1,245 |
80,16 |
394 |
43,3 |
1,57- Ю-'» |
2,10-102 |
2,3423. |
|
1,25 |
80,23 |
372 |
41.5 |
9,48-10"19 |
13,30-102 |
3,1239 |
|
1,255 |
S0.30 |
350 |
40,0 |
4,24-10-18 |
5,98-10я |
3,7769 • |
|
1,26 |
80,36 |
339 . |
39,0 |
1,15-10-1' |
1,62-10' |
4,2065 |
22
Зная минимальное значение fK npH различном пересы щении раствора, можно определить границу метастабиль ной зоны сахарных растворов, выше которой при дайной температуре возможно 'Самопроизвольное образование , центров кристаллизации. Данные для определения зна чений / к приведены в табл.' 5, а зависимость lg /„ от ко эффициента пересыщения раствора показана на рис-.-6.
|
|
|
•igh |
|
|
|
k |
|
|
|
3 |
|
|
|
2 |
Рис. 6. |
Зависимость lg/ K |
^ |
|
от коэффициента пересы- |
_ ^ |
||
щения |
при |
температуре: |
|
J - 50°С; |
2 - |
70°С; 3 — 90°С. |
|
|
|
|
-2 |
|
|
|
-3 |
|
|
|
-4 |
|
|
|
-5 |
Как видно из табл. 5 и рис. 6, при a=d,-20 1,23 lg / к отрицателен. Следовательно, в этом интервале пересыще ний пои температуре 70° С вероятность достижения обра зовавшимися зародышами критических размеров нич тожно 'мала. При а>,1,йЗ lg f,(- становится положитель ным. Таким образом, границу метастабильности сахар ных растворов при температуре 70° С, выше которой тео ретически возможно спонтанное образование центров кристаллизации, необходимо .признать при а>1,23.
Изменение температуры раствора перемещает грани цу метастабильной области. Результаты определения функции распределения зародышей в пересыщенных ра створ ак при температуре 60 и 90° С показали, что пони жение температуры (при постоянном пересыщении) вы зывает расширение метастабильной зоны. При темпера туре 50° С ее граница соответствует коэффициенту пересыщения а—11,00-Н1,25 (см. рис. 6). И наоборот, с по вышением температуры образование центров кристалли зации возможно при меньшем коэффициенте пересыще ния. При температуре 90° С предел метастабильной зоны соответствует а=>1,205.
23
ОБРАЗОВАНИЕ ЦЕНТРОВ НОВОЙ ФАЗЫ ПО ТЕОРИИ ФОЛЬМЕРА —ФРЕНКЕЛЯ
Термодинамические закономерности образования нозой фазы Binар вые развиты Гиббсом (98], а затем исполь зованы Фольмером |jl'06] для разработки теории зарож дения (центров новой фазы при стационарном состоянии системы.
Классическая флуктуационная теория фазовых пре вращений разработана и экспериментально проверена для процесса конденсации пара, В пересыщенном паре давлением р статистические флуктуации способствуют об разованию зародышей, содержащих различное количе ство молекул «.
При образовании зародыша новой фазы в исходной фазе возникает поверхность раздела двух фаз. На ее соз дание затрачивается свободная энергия системы. Поэто му переход метаетабильной фазы в стабильную на первой зародышевой стадии (вопреки второму началу термоди намики) сопровождается увеличением свободной энергии системы, которую Гиббс [98] принимал за меру устойчи вости метаетабильной системы.
Химический потенциал одиночной молекулы в паре
.определяется по фор муле
н = |xto + k Tin z\ , |
(30) |
а для каждой молекулы в массе жидкости
(лэ = (ijO + k F i n / „ , |
(31) |
где ('«> — концентрация молекул в насыщенном паре.
Следовательно, i |
молекул в массе -жидкости |
имеют |
свободную энергию i |
\i2. Если небольшую каплю отделить |
|
от большой массы жидкости, то свободная энергия |
моле |
кул этой капли увеличится на величину поверхностной-
энергии капли с F i . Тогда |
свободная |
энергия образова |
|
ния капли жидкости из пересыщенного пара |
равна |
||
' д Ф = / (^ |
- в ) + о |
• |
(32) |
Принимая форму зародыша сферическую (соответст вует минимуму свободной поверхностной энергии), полу чим з н ачение i к апли
24
где Ti — радиус капли;
V2 — объем молекулы в жидком состоянии.
Подставляй значение |
I |
в уравнение (32) |
и заменяя |
Fi—4 лгг2, получим |
|
|
|
4я |
|
(1,) г/> + 4к гС- а. |
(34) |
Д Ф = — (,u3 |
- |
Приравнивая производную от АФ по г нулю и определяя гк , получим
2oV2 |
|
2а V3 |
(35) |
|
Н — F-2 |
kT |
In |
||
|
||||
|
|
|
Вторая производная отрицательная, следовательно, кри вая АФ (г) имеет максимум.
Как видно из рис*. 7, образование зародышей с очень маленькими радиусами сопровождается возрастанием свободной анергии системы, или термодинамического по
тенциала. |
Выше некоторого значения радиуса |
гк, назы |
ваемого |
критическим, термодинамический |
потенциал |
|
4ФА |
|
Рис. 7. Изменение термодинамиче ского потенциала системы при об разовании зародышей новой фазы.
уменьшается, а следовательно, центры новой фазы, до стигшие г к , 'будут свободно расти. [Потенциал при' г к так же называется критическим. Его величину можно найти, подставляя в уравнение (34) значение (д-i—цг) из урав нения: (35),
Д ф к = — |
ъаг1, |
(36) |
Подставляя в уравнение |
(36) значение гк из уравне |
||
ния (35) ..получим |
|
|
|
|
1 6 * а з у а 2 |
(37) |
|
3 £ 2 |
^(lnP/P™)-' |
||
|
Из уравнений (35) и '(37) видно, что с увеличением коэффициента пересыщения гк и ДФК уменьшаются, что облегчает образование зародышей. Чем выше 'коэффици ент пересыщения* тем меньшие то величине зародыши становятся способными к дальнейшему росту.
Рассматривая,.процесс образования щ зародышей в единице объема паровой фазы равновесным, применили закон действия масс, согласно которому
л; = г,ехр — (Л Ф,-/ЛГ): |
(3S) |
где га,-— число зародышей, состоящих из / одиночных молекул.
Скорость образования центров новой фазы равна про
изведению |
концентрации |
зародышей |
критического |
||||
размера пк |
на |
частоту |
gK, |
с которой они |
присоединяют |
||
одну молекулу, |
|
|
|
\ |
|
|
|
|
I |
= |
к «к = |
g* |
ехр - |
(Л Фк7е Г). |
. (39) |
Подставляя |
в уравнение |
(39) |
значение |
ДФК из урав |
|||
нения (37), получим |
|
|
|
|
|||
|
|
|
|
|
16- о3 |
i y |
|
|
|
|
|
З*-' Г 3 (In |
Р,Р"„ |
|
|
|
|
/ = |
ЯГк «1 е. |
|
• |
(40) |
Зависимость / от пересыщенного пара Р/Л» представ лена на рис. 8. Из рис. 8 видно, что скорость образования центров новой фазы незначительно мала, пока пересыще ние достигает критического значения — метастабилкного предела ло Оствальду, а затем резко увеличивается.
Основные положения рассмотренной теории образо вания центров новой фазы применены и для пересыщен ных растворов. Однако при этом необходимо учитывать ряд особенностей конденсированной среды. Пересыщенкый раствор отличается от газообразной фазы присут ствием молекул растворителя, которые затрудняют пере в о д молекул растворенного вещества в кристаллическую
фазу. Присутствие в системе молекул второго коМпонен-
26
та и гидратация молекул растворенного вещества значи тельно, усложняют обмен мономерньтми молекулами на поверхности зародышей, уменьшая их частоту обмена и, следовательно, скорость обра зования центров кристаллиза ции.
По |
Френкелю [89], вероят |
||||
ность |
перехода |
молекул |
в |
||
кристаллический |
зародыш |
из |
|||
пересыщенного |
раствора |
про |
|||
порциональна |
|
выражению |
|||
гхр—(&U/RT), |
в |
котором |
|||
\0—некоторая |
энергия |
акти |
|||
вации, |
являющаяся |
функцией |
|||
вязкости раствора |
по уравне |
||||
нию |
|
|
|
|
|
. i U'RT
т] = А е
где 11 —вязкость раствора; А —постоянная^величина; R— газовая постоянная.
Тогда произведение gKii представить выражением
Р/Р-
(41)
Рис. 8. Зависимость скоро сти образования центров но вой фазы от пересыщения пара.
= В\ IB уравнении (40) можно
В = S i ехр ( - Д UJRT) , |
(42) |
где Si — постоянный коэффициент. |
|
Подставляя |
значение В в уравнение (40) |
и за'.меняя |
||||||
R |
„ |
М |
• |
получим |
|
|
|
|
k = -тг- |
и V? = |
|
|
|
|
|||
Nn |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16» Лг0 At- с3 |
|
|
I |
= В, |
ё-4 |
U'RT |
3R3 |
Т" р= (in С Coo)1 |
(43) |
|
|
|
|
где М — молекулярная масса кристалла;
. р—плотность кристаллического вещества; Л'о — число Авогадро.
Однако дальнейшие исследования показали, что для пересыщенных, растворов уравнение (43) не описывает экспериментальную 'кривую скорости образования кри сталлов, так как оно 'Справедливо только при условии стационарного распределения зародышей.
Процесс возникновения центров кристаллизации про текает во времени. Под действием тепловых флуктуации'
27
в каждый момент времени изменяется распределение за родышей по размерам. Поэтому 'не все зародыши, достиг шие величины пк-и станут .центрами кристаллизации. Часть из них, потеряв .по одной или две молекулы, .перей дет в более низший класс. Поэтому фактическая скорость образования центров [кристаллизации будет меныше, чем вычисленная по уравнению (43). Эти вопросы рассматри ваются ниже. '
НЕСТАЦИОНАРНЫЙ РЕЖИМ ОБРАЗОВАНИЯ ЗАРОДЫШЕЙ
.НОВОЙ ФАЗЫ ПО ТЕОРИИ ФРЕНКЕЛЯ — ЗЕЛЬДОВИЧА
Френкель [89], Зельдович [19], Баккер и Доринг [96], развивая дальше теорию фазовых переходов, исследовгли процессы 'возникновения зародышей новой фазы в усло виях нестационарного режима, при котором скорость об разования центров кристаллизации -зависит от времени.
Схематично образование зародышей новой фазы три стационарном и нестационарном режимах, показано на рис. 9. Пусть в начальный момент времени система со
держит |
определенное |
число |
л,- зародышей |
одинакового |
||||||||
|
2 |
|
|
размера. |
Функция |
рас |
||||||
f |
|
|
пределения / |
равна |
нулю |
|||||||
|
|
|
|
везде, |
кроме |
определен |
||||||
|
|
|
|
ной |
точки |
(прямая |
/ ) . П о |
|||||
|
|
|
|
макроскопической |
теории |
|||||||
|
|
|
|
каждый |
зародыш • изме |
|||||||
|
|
|
|
няется |
во времени |
вполне |
||||||
|
|
|
|
определенным |
образом |
в |
||||||
|
|
|
|
зависимости |
от |
|
своих |
|||||
|
|
|
|
размеров |
и внешних |
ус |
||||||
|
|
г |
п |
ловий, |
следовательно, |
че |
||||||
|
|
рез |
определенный |
|
про |
|||||||
Рис. 9. Изменение функции'"распре |
межуток |
времени |
t |
число |
||||||||
деления зародышей во времени при |
зародышей rii останется |
|||||||||||
стационарном {1, 2) и |
иестацио- |
|||||||||||
|
нарном (3) |
режимах. |
неизменным |
и |
тождест |
|||||||
|
|
|
|
венным |
первоначальному, |
|||||||
т. е. прямая / переместится на другое место (прямая |
2). |
|||||||||||
По теории," учитывающей флуктуации, |
прямая |
/ |
через |
|||||||||
время t |
превратится |
в |
кривую 3. |
Вследствие |
флуктуа |
ции зародыши, которые первоначально не отличались друг от друга, будут различны и распределение размо ется. Наряду с общим увеличением п, за счет присое-
28
динения к зародышам щ-i по одной молекуле образу ются и такие зародыши, которые по своим размерам сдвинулись влево и уменьшили значение щ путем отде ления от них одной молекулы. Эти процессы особенно интенсивно протекают, по выражению Зельдовича, в критической области, для которой значение свободной энергии системы отличается от максимального не более чем на величину kT.
Обозначим через ,|3г среднее число -молекул, которые отрываются в единицу времени с единицы поверхности зародыша яг*2. Тогда произведение Fi^dt пред ставляет собой вероятность отрыва молекулы с поверхно сти зародыша за время Ш. Аналогичным образом обозна чим через gidi вероятность присоединения очередной мо лекулы на единице поверхности зародыша. Если между обеими фазами системы имеется устойчивое равновесие, то отрыв и присоединение молекул протекает с одинако вой скоростью и будет справедливо следующее соотно шение:
= / ? + 1 f c + i . |
(44) |
где ft" — функция равновесного распределения зародышей по числу молекул 1.
Подставим в соотношение (44) значение f° из уравне ния (23), тогда-
|
|
|
Д Ф , - + 1 - Д Ф . |
|
||
|
gi=*h+ie |
к |
Т |
. |
(45) |
|
Если исходное распределение |
зародышей по их раз |
|||||
мерамнеравновесно, то оно с течением времени |
должно |
|||||
меняться. Число |
зародышей |
,+ь переходящих |
за еди |
|||
ницу времени из класса |
I в класс t-Ы, составит |
|
||||
|
h, i+i = Sifi |
— Pi+i fi+i |
(46) |
|||
Перепишем уравнение!(44) в виде |
|
|||||
|
gih |
= gifi°-J^ |
|
н |
|
|
|
|
|
/1 |
|
|
|
о г |
й _ |
,о A ± i _ _ а |
f о Il±l., |
|
||
P,-+i ft+i |
P=t+i |
/i+i^o |
= |
gift0-л |
|
|
|
|
|
/1+1 |
|
'l+l |
|
29