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

книги из ГПНТБ / Вигдорович, В. Н. Совершенствование зонной перекристаллизации

.pdf
Скачиваний:
28
Добавлен:
21.10.2023
Размер:
9.49 Mб
Скачать

печивает беспрепятственного протекания процесса по схеме (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

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