
книги из ГПНТБ / Вигдорович, В. Н. Совершенствование зонной перекристаллизации
.pdfпечивает беспрепятственного протекания процесса по схеме (I), становится энергетически возможным, а кинетически выгодным так называемый процесс бездиффузионной кристаллизации, при котором образующаяся твердая фаза по химическому составу соответствует жидкой фазе на всех этапах кристаллизации. Очевидно, что ни раз деления, ни очистки компонентов при этом не происходит.
Реальный процесс кристаллизации, будучи осуществленным в раз личных условиях, представляет собой сложное сочетание рассмотрен ных случаев. Однако физико-химическое рассмотрение показывает, что для разделения или очистки компонентов может использоваться процесс направленной кристаллизации, в котором для большей эффективности необходимо предусмотреть: 1) охлаждение твердой фазы для предотвращения выравнивания ее химического состава; 2). перемешивание жидкой фазы для обеспечения выравнивания ее химического состава; 3) поддержание интенсивности охлаждения и скорости кристаллизации на таком уровне, который позволял бы использовать первые два условия и не приводил бы к возможности протекания бездиффузионной кристаллизации.
Обратимся к случаю направленной кристаллизации при отсут ствии диффузии в твердой фазе и беспрепятственной диффузии в жидкой фазе.
Рассмотрим двухкомпонентную систему А —В, в которой компо ненты А и В взаимно растворяются в жидком и твердом состояниях. Пусть Сж — концентрация компонента В в жидкой фазе, Ств — в твердой фазе и т — количество кристаллизующейся жидкости.
При |
температуре t + dt масса компонента В в жидкости будет C m . |
|
При |
бесконечно малом понижении температуры до t концентрация |
|
жидкой фазы уменьшается на йСж, т. е. станет Сж — йСж, |
а масса |
|
ее уменьшится до т — dm. Твердая фаза с концентрацией |
Ств при |
этом образуется в количестве dm. Следовательно, масса компонента В
.в жидкости при температуре t составит (Сж — dCx) (т — dm), а в твердой фазе CTBdm. Так как общая масса компонента В остается
постоянной, то |
получаем |
|
(Сж |
— dCx) (т — dm) + CTBdm = Сжт. |
(1.1) |
Отсюда, пренебрегая бесконечно малой второго порядка (<2СЖdm), имеем
__ |
асж___ |
/ т г>\ |
т ~ С ТВ- С ж > |
|
|
а после интегрирования |
|
|
С |
=с„ |
|
1п^> |
dC, |
(1-3) |
|
||
m |
|
|
где т 0 и Сж= С0 отвечают началу кристаллизации и Ств, Сжи т — любой температуре ниже температуры начала кристаллизации.
Полученное выражение известно под названием «уравнение Рэлея».
ю
Конец кристаллизации отвечает условию т = О, вследствие чего
СЖ = С 0
оо.
Ств Сж
Это означает, что интеграл расходится и, следовательно, подын тегральная функция является неограниченной. Поэтому
— L—---->оо и Сж--•> Ств. Другими словами, при отсутствии диффузии
Ств '-'ж
в твердой фазе и беспрепятственной диффузии в жидкой фазе имеет
место разделение компонентов А и В. |
|
|
||||||||||||
При |
исходном содержании |
компонента |
|
|||||||||||
В в жидкой фазе Сж= С0, в твердой |
|
|||||||||||||
фазе получаем непрерывный |
набор кон |
|
||||||||||||
центраций |
компонента |
В |
от Ств |
при |
|
|||||||||
Сж— Со в начале кристаллизации («го |
|
|||||||||||||
ловная фракция») |
до |
СТВ, |
практически |
|
||||||||||
равного |
|
Сж в |
конце |
кристаллизации |
|
|||||||||
(«хвостовая фракция»). |
|
|
|
|
С0) |
|
||||||||
Для определения Ств (при Сж = |
|
|||||||||||||
необходимо |
знать |
функциональную |
|
|||||||||||
связь между Ств и Сж. Если принять, |
|
|
||||||||||||
что Ств = 1г0Сж, где |
к 0 — равновесный |
|
||||||||||||
коэффициент распределения, |
то |
в го |
|
|||||||||||
ловной |
фракции |
при k 0 < |
1 имеем |
Ств |
|
|||||||||
(при Сж = |
С0) < |
С0 |
(происходит |
очи |
Рис. 3. Схема распределения |
при |
||||||||
стка |
по |
компоненту В) |
и |
при |
k 0 > 1 |
месей на границе раздела твердой |
||||||||
имеем Ств |
(при |
Сж = |
С0) > |
|
С0 (проис |
и жидкой фаз при различных усло |
||||||||
|
виях кристаллизации: |
|
||||||||||||
ходит |
очистка |
по |
компоненту |
А). В |
а — равновесные условия (скорость |
|||||||||
хвостовой фракции возможно получе |
кристаллизации пренебрежимо |
ма |
||||||||||||
ла); б — установившиеся • условия |
||||||||||||||
ние чистых компонентов А (при k 0 5> 1) |
(конечная скорость кристаллизации) |
|||||||||||||
или В |
(при k 0 < |
1). |
что чем больше отличается k 0 от единицы, |
тем |
||||||||||
Нетрудно видеть, |
||||||||||||||
больше отличается СТБ (при Сж = С0) |
от С0. |
при |
||||||||||||
Теперь |
обратимся |
к |
случаю |
направленной кристаллизации |
протекании диффузии в жидкой фазе с конечной скоростью.
При недостаточной диффузионной подвижности атомов в жидкой фазе и недостаточном ее перемешивании, когда ■кристаллизация происходит с конечной скоростью, в жидкой фазе вблизи фронта кристаллизации образуется особый так называемый диффузионный слой. Если k 0 < 1, то в результате оттеснения фронтом кристалли зации компонента, которым обеднена твердая фаза, в жидкости возникает градиент концентрации, экспоненциально убывающий в направлении от фронта кристаллизации; при k 0 >-1 градиент концентрации имеет противоположное направление (рис. 3).
В связи с рассмотренным явлением оказывается удобным разли чать эффективный коэффициент распределения k и коэффициент рас пределения на границе раздела k*. Величина k tпредставляет собой
11
отношение концентрации в твердой фазе Ств, кристаллизующейся непосредственно на фронте кристаллизации, к концентрации в объеме жидкой фазы Сж за пределами диффузионного слоя, имеющего тол щину, превышающую б (рис. 3). Практически это расстояние может быть определено как несколько 1 значений D/f, где D — коэффи циент диффузии примеси и / — скорость кристаллизации. Вели чина k* определяется как отношение концентрации в твердой фазе Ств к концентрации в жидкой фазе на поверхности фронта кристаллиза
ции Сж. Величина k , как и величина k, зависит от скорости кри сталлизации, концентрации примеси и ориентации растущего кри сталла. Однако в отличие от величины k значение k* не изменяется
взависимости от интенсивности перемешивания. Величина k ближе
кединице, чем величина k*, но в условиях интенсивного перемеши вания величина k может приближаться к величине k*. Если скорость кристаллизации мала, то k* я=* k 0. Поверхностная активность (ад сорбция) и гетерогенные химические процессы на фронте кристалли зации и другие явления могут приводить к существенному разли чию k* и k Q.
Наиболее удовлетворительную теорию, объясняющую изменение величины k с изменением скорости кристаллизации и рассматри вающую роль диффузионного слоя вблизи фронта кристаллизации, создали Бартон, Прим и Слихтер [33]. Они вывели формулу
k*
^ k* + (1 — £ *)ехр (— А) ’ (1-4)
которая чаще используется в виде |
|
k = _________ ^ ________ |
(1.5) |
ka + (1 — *0)ехр (— А) ’ |
|
где Д = 8//D — безразмерный критерий, называемый |
приведенной |
скоростью кристаллизации. |
|
На рис. 4 зависимость коэффициента распределения от скорости кристаллизации представлена графически. Можно видеть, что с уве личением скорости кристаллизации при данной интенсивности пере мешивания значение k приближается к единице. Применение пере мешивания, уменьшающего толщину диффузионного слоя, приводит
к |
приближению |
величины k |
при данной скорости кристаллизации |
к |
величине k 0. |
|
|
|
Из графика |
зависимости |
от /, который в случае |
подтверждения теории Бартона—Прима—Слихтера должен быть прямолинейным, по тангенсу угла наклона (тангенте) может быть определена величина 8/D, незначительно зависящая от /.
И. А. Мержанов и А. Г. Аникин [34 ] зависимость эффективного коэффициента распределения от скорости кристаллизации предста-
1 Диффузионный и гидродинамический аспект этого вопроса неоднократно обсуждался, например, в работах [8, с. 33 и 22, с. 45, 49, 70, 80].
12
вили в ином виде: |
|
± __L_ t L |
|
|
|
|
||
|
|
|
|
|
( 1.6) |
|||
|
|
k |
k0 |
D ’ |
|
|
|
|
где у — отношение плотностей в твердой и жидкой фазах. |
|
|||||||
Как было показано в работе |
[35] формулы |
(1.5) |
и (1.6) могут |
|||||
быть приведены 1 к виду |
|
|
|
|
|
|
|
|
k = |
k 0 [1 + |
(1 — k 0)f8/D] |
|
|
(1.7) |
|||
и |
|
|
|
|
|
|
|
|
k |
= |
k 0 [1 + |
k 0fy/D] |
|
|
|
(1.8) |
|
соответственно. Отсюда |
видно, |
что при |
k 0 <; 1 |
по |
формуле |
(1.8) |
||
величина k возрастает |
с увеличением |
/ значительнее, чем |
по |
|||||
f, см/с |
|
|
|
f, см/с |
|
|
|
fS/д |
|
|
|
fS/D |
Рис. 4. Зависимость коэффициентов |
распределения k от приведен |
||||
ной |
скорости кристаллизации fb/D |
и скорости кристаллизации f |
|||
при Ь/D = 100 для |
различных значений равновесных коэффициен |
||||
|
тов распределения k 0 (цифры |
на |
кривых): |
||
ч |
а — при &0 < 1; б — при |
kQ > 1 |
|||
|
|
|
|
|
|
формуле (1.7). Экспериментальные |
данные |
подтверждают справед |
|||
ливость формул (1.5) и (1.7). |
|
диффузионных процессов |
|||
Другим |
следствием |
ограниченности |
в жидкости и связанного с этим образования диффузионного слоя является концентрационное переохлаждение (иногда это явление называют диффузионным или конституционным переохлаждением). При повышенных скоростях кристаллизации переохлаждение может достигать значительной степени, и тогда станет возможной бездиффузионная кристаллизация. Рассмотрим оба эти явления.
1 Используются следующие преобразования: дифференцирование k по / и под-
становки d/ ^ |
— 0 = /, k |
k0 |
и k = k0 |
dk |
- f -щ. df. |
13
На рис. 5 схематически показано возникновение концентрацион ного переохлаждения. На рис. 5, а показано положение ликвидуса и солидуса, на рис. 5, б — концентрационный градиент вблизи фронта кристаллизации, а на рис. 5, в — произвольный температурный
Расстояние от фронта кристаллизации
Расстояние от фронта кристаллизации
Рис. 5. Схема, поясняющая возникновение концентрационного переохлаждения при k 0 < I и k 0 > 1:
а — относительное расположение линий ликвидуса и солидуса; б — градиент концентрации в расплаве вблизи фронта кристал лизации; в — относительное расположение линий температур ного градиента в расплаве и начала кристаллизации различ ных сдоев (/, 2, 3, 4, 5) расплава в прифронтовой области
градиент и линия начала кристаллизации различных слоев диффу зионного слоя. Между этими двумя линиями расположена область переохлаждения. Естественно, что вероятность концентрационного переохлаждения возрастает при увеличении скорости кристаллиза-
14
ции, так как это влечет за собой рост градиента концентрации на фронте кристаллизации. Регулирование градиента температур и перемешивание расплава делают возможным предотвращать кон центрационное переохлаждение.
На рис. 6 схематически показано относительное расположение линий изобарного потенциала (свободной энергии по Гиббсу), ана лиз которых позволяет выявить возможность бездиффузионной кри сталлизации. При переохлаждении жидкости любого данного со-
фузионной кристаллизации |
при |
переохлаждениях ниже линий tA t0 |
||||
равных |
изобарных потенциалов |
G (свободных энергий) для |
жидкой |
|||
|
и твердой фаз |
при |
k 0 < 1 (а) |
и k 0 > 1 (б) |
|
|
става до температур |
tAt 0 величина изобарного потенциала твердого |
|||||
раствора G'TB |
будет |
больше |
изобарного |
потенциала |
расплава G* |
(см. состав I), и кристаллизация в этом интервале температур может происходить только диффузионным путем. При переохлаждении жидкости ниже линии tAt 0 встречаемся с обратным соотношением изобарных потенциалов жидкой и твердой фаз (состав II). Пере охлаждение ниже линии равных значений изобарных потенциалов жидкой и твердой фаз делает кинетически более выгодным бездиффузионный процесс кристаллизации, приводящий жидкость в метастабильное G?B, а не стабильное G'c состояние.
Концентрационное переохлаждение, помимо снижения эффекта очистки, при повышенных скоростях направленной кристаллизации из-за бездиффузионной кристаллизации приводит также к искривле
15
нию фронта кристаллизации, дендритному механизму роста, ячеистой структуре («соты»), зарождению избыточных центров кристаллиза ции, механическому захвату примесей и неоднородности их рас пределения.
КОНЦЕНТРАЦИОННАЯ ЗАВИСИМОСТЬ КОЭФФИЦИЕНТОВ РАСПРЕДЕЛЕНИЯ, ВЫРАЖЕННЫХ ЧЕРЕЗ АБСОЛЮТНЫЕ И ОТНОСИТЕЛЬНЫЕ КОНЦЕНТРАЦИИ КОМПОНЕНТОВ
Коэффициент распределения компонентов между соприкасающимися фазами — основная физико-химическая характеристика процессов разделения и очистки ве ществ и, в частности, кристаллизационных процессов.
В общем случае коэффициент распределения является функцией концентрации компонентов в фазах и, следовательно, изменяется в процессе разделения компонен тов. Известны попытки создания 'теории зонной перекристаллизации, учитываю щей эти изменения; так как, однако, это приводит к сложным выражениям, то обычно использ.уется допущение о постоянстве коэффициента распределения в концентра ционном интервале очистки. Коэффициенты распределения могут выражаться через абсолютные концентрации компонента
|
h |
— Схв/Сж |
(1.9) |
и через относительные концентрации |
|
||
ос = |
Ств (1 — Сж)/Сж (1 — Си). |
(1.10) |
|
где Схв и Сж — молярные |
доли |
компонента соответственно |
в твердой и жидкой |
фазах (имеется |
в виду двухкомпонентная |
система). |
Заметим, что физико-химические расчеты могут быть облегчены с помощью таблиц перевода абсолютных концентраций в относительные и обратно.
Высказывалось мнение, что при втором способе выражения коэффициент рас пределения меньше зависит от концентрации, чем при первом способе [22, 36, 37]. При этом подразумевалась аналогия с ректификационными процессами в системах, близких к идеальным. Другими словами, в большом интервале концентраций, на
чиная с нулевой концентрации, когда Схв = 0, справедливо неравенство |
|
||
|
| da/dCTB I < |
| dk0!dCTB | |
( 1.11) |
Известно, между тем |
[38], что при переходе от диаграммы состояния в коорди |
||
натах давление — состав |
к диаграмме в |
координатах температура — состав |
пара |
метр а теряет прежний физический смысл, |
и его постоянство никак не связано с бли |
||
зостью растворов к идеальным. |
|
|
Поэтому представляло интерес исследовать характер концентрационных зави симостей коэффициентов k0 и а без каких-либо допущений о взаимодействии компо
нентов в растворе [39].
Из равенств (1.9) и (1.10) следует, что
|
|
|
|
а — bp |
бхв |
|
( 1. 12) |
||
|
|
|
|
|
1 -— Схв |
|
|
||
и |
|
da |
_ |
k0— 1 |
|
1 |
dk0 |
|
|
|
|
|
(1.13) |
||||||
|
dCTв |
|
(1 — Ств)2 |
1 — Схв dCjв |
|||||
|
|
|
|||||||
Отсюда можно |
получить |
некоторые |
сведения |
о |
связи между |
коэффициентами а |
|||
и А0. При Схв -> 0 имеем |
|
lim |
а — а0 = |
К0, |
(1.14) |
||||
|
|
|
|
||||||
а при Схв 1 |
имеем kn |
1 |
|
СхВ+ ° |
|
|
|
|
|
и |
|
|
|
|
|
||||
|
|
П т |
а = |
a, = |
1 — |
|
|
(1.15) |
схв + 1
16
Следовательно, в области больших концентраций разница между а и ! г 0 тем больше, чем больше в этой области меняется k0 при изменении Схв.
Полагая Схв -> 0 в уравнении (1.13), получим
( d a \ |
= / dkо \ |
\d.CrB)o |
\dCTBy о |
+ Ко — 1• |
(1.16) |
В этом случае неравенство (1.11) выполняется, если
/ dkg \ |
„ |
/ dkn \ |
(1.17) |
\ Щ - Х + |
-1 < |
\dCTBJ о |
|
Ко~ |
|
т. е. если значения (dk0/d.CTB)0 лежат в заштрихованных областях рис. 7, а. Та ким образом, предположение о меньшем изменении коэффициента а с концентра цией по сравнению с изменением коэффициента k0 может не выполняться уже в об-
Рис. |
7. Условия выполнимости соотношения |
da |
\ dk о |
(заштрихован- |
dC„ |
\~dC~ |
|||
ные |
области) при малых (а) и больших |
-'тв | |
I “~тв . |
второго ком |
(б) концентрациях |
понента
ласти самых, малых концентраций. Отметим, что попадание точек в первый и третий квадранты рис. 7, а свидетельствует о немонотонном характере зависимости коэф фициента k 0 от концентрации. Однако может оказаться, что неравенство (1.11)
будет, выполняться вплоть до самых высоких концентраций. Из уравнения (1.13) следует, что при Схв -> 1
Условие выполнения неравенства (1.11) — попадание точек в заштрихованную об ласть рис. 7, б.
Зависимость коэффициента k 0 от концентрации легко было бы получить
если бы существовало общее уравнение, связывающее между собой концентрации линий ликвидуса и солидуса при одинаковых температурах. Такое уравнение пока отсутствует. При ограниченной растворимости компонентов в жидком или твердом состояниях коэффициенты распределения — разрывные функции концентрации. В областях, где коэффициент распределения — непрерывная и монотонная функция концентрации, его можно аппроксимировать выражением
|
|
* 0 = |
< ? в + Ь, |
(1.19) |
где п^> 0. |
Учитывая, |
что |
|
|
|
kg = |
Кд при СхВ = |
0 И kg = 1 ПрИ С?В = L, |
( 1. 20) |
2 В. |
Н. Вигдорович |
|
17 |
получим для систем с неограниченной растворимостью компонентов в твердом и жидком состояниях
* 0 = 0 |
— Ко) Ств + |
Ко- |
|
( 1.21) |
|
Графически эта зависимость представлена на рис. 8 сплошными линиями. |
|||||
Согласно уравнение (1.21), имеем |
|
|
|
|
|
dk, |
|
|
•4—1 |
|
( 1. 22) |
Ж7 = п ( 1 ~ К 0)С:ТВ |
|
||||
^ = п ( |
п - |
1 ) ( 1 |
п—2 |
(1.23) |
|
- к 0) с:тв |
• |
||||
dC: |
|
|
|
|
|
Если Ко^> 1, то k0 монотонно убывает с ростом Схв; |
в этом случае при п^> 1 |
||||
кривая k0 (Схв) обращена выпуклостью в сторону от оси Схв, при |
1 к оси Схв, |
||||
К,.<* |
при |
/С0 <1 1 — соотношения |
обратные. |
||
|
На рис. 8 штриховыми линиями показано |
||||
|
изменение |
коэффициента а |
с концентрацией. |
Ю
|
|
|
|
0 |
0,2 ОМ ОМ 0$Ст8 0' |
0,2 ОМ 0,6 0,8 Сг |
Рис. |
8. |
Концентрационная |
зави |
Рис. 9. |
Концентрационные зависимости коэффициен |
|
симость коэффициентов k 0 (сплош |
тов k Qи а в системе германий—кремний: а, б — для |
|||||
ные |
линии) и а (пунктирные ли |
кремния в германии; в, г — для германия в кремнии; |
||||
нии) |
в |
предположении, что |
= |
а, в — из диаграммы состояния' |
б, г — по уравне |
|
|
|
нию Ван-Лаара; |
/ — k Q\ 2 — а |
= 0 |
- * о ) Ств + К0 и |
» >0 |
|
|
|
(цифры на кривых) |
|
|
|
Особенно сильное возрастание коэффициента наблюдается, |
в соответствии с урав |
|||
нениями (1.14) |
и (1.15), |
при больших значениях п. |
концентрацией коэффи |
|
В |
качестве |
примера |
рассмотрим, как изменяются с |
циенты k0 и а в системах германий—кремний и висмут—сурьма с неограниченной
растворимостью компонентов в жидком и твердом состояниях (рис. 9 и 10). Кроме данных, полученных из диаграмм состояния на рисунках, представлены значения
коэффициентов k0 и а, |
вычисленные по уравнениям Ван-Лаара [31]. |
Из рисунков видно, |
что а является функцией концентрации не только в случае |
регулярных растворов, которыми являются рассмотренные системы, но и в случае идеальных растворов.
18
Здесь умест»-^- остановиться яа одном ошибочном утверждений относительно связи между коэффициентами k0 и а в идеальных растворах. Согласно А. Н. К,ир-
гинцеву [36], в таких растворах одинаковое изменение температуры всегда будет вызывать большее изменение величины k 0, чем величины а.
5
2
1 "
|
О |
0,2 ОА |
0,6 |
0,8 |
1,0 |
О 0,2 ОА |
0,6 |
0,6 |
1,0 |
|||
|
|
|
|
|
^тв |
|
|
|
|
|
^тв |
|
|
Рис. |
10. |
Концентрационные зависимости коэффициентов k 0 и а |
|||||||||
|
в системе висмут—сурьма: |
а, |
б — для |
сурьмы |
в висмуте; в, |
|||||||
|
г — для |
висмута в |
сурьме; |
а, в — из диаграммы |
состояния; |
|||||||
|
|
б, |
г — по |
уравнению Ван-Лаара; |
1 — £0; |
2 — а |
|
|||||
|
Из уравнений Ван-Лаара следует, что |
|
|
|
|
|
||||||
|
|
|
|
/?0 = ехр |
М П - Т ) |
|
|
(1.24) |
||||
|
|
|
|
|
|
|
|
R TJ' |
|
|
|
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
а = ехр ГM 7 V -I) |
М Т 2- Г ) |
-1 |
(1.25) |
||||||
|
|
|
|
|
RT{T |
|
|
RT2Tx |
J ’ |
|||
гДе |
Т — температура, |
значение |
которой лежит в интервале между темпера |
|||||||||
|
турами плавления чистых компонентов Т г и Т 2; |
|
||||||||||
|
^1 и Х2 — скрытые теплоты плавления этих компонентов; |
|
||||||||||
|
R — универсальная газовая |
постоянная. |
|
|
|
|||||||
|
Производные от k0 и а по температуре связаны уравнением |
|
||||||||||
|
|
|
da |
|
|
1 |
^2/^1 |
dk0 |
|
|
(1.26) |
|
|
|
|
W |
|
ехр |
^2 (Tg |
Т) |
W ’ |
|
|
||
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
RT2T |
|
|
|
|
|
из |
которого следует, |
что условие |
|
|
|
|
|
|
|
|||
|
|
|
|
|
I da |
|
|
dk0 |
|
|
|
(1-27) |
|
|
|
|
|
\ d f |
|
|
dT |
|
|
|
2 |
19 |