Файл: Лодиз, Р. Рост монокристаллов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 14.10.2024

Просмотров: 173

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

398 Р. ПАРКЕР. МЕХАНИЗМЫ РОСТА КРИСТАЛЛОВ

ние кинетических явлений или поверхностной энергии на тем­

пературу движущейся границы раздела фаз, концентрацию на

ней и форму этой границы; влияние гидродинамических

пото­

ков в расплаве, растворе или паре (перемешивание или

естест­

венная конвекция из-за разности плотностей) на условия теп­ лопередачи или диффузии; влияние примесей на температуру

поверхности раздела фаз или на концентрацию на ней

и ее

форму;

совместное действие

потоков

тепла и

вещества

при

росте из

концентрированных

растворов;

наконец,

влияние

сило­

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

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

Здесь

же мы обсудим

эти

эффекты

в духе

задачи Стефана,

т. е. их

влияние на краевые

условия

и

сами

уравнения пере­

носа.

 

 

 

 

 

 

Плоскость, сфера, цилиндр. /. Влияние

примеси. Как уже от­

мечалось, теплоперенос

и диффузия вещества

описываются оди­

наковыми дифференциальными уравнениями (9.1) и (9.36), ре­ шения которых должны удовлетворять сходным граничным ус­ ловиям (9.3), (9.4), (9.37) и (9.38). До сих пор в задачах рассматривался только один из двух процессов переноса. При исследовании направленной кристаллизации считалось, что кри­ сталл растет из чистого расплава. При решении же задачи о кристаллизации цилиндра из пересыщенного раствора тепловые эффекты не учитывались. Но, как показал Франк [54], характер роста может определяться совместным действием обоих процес­ сов переноса. В частности, расплав,из которого растет кристалл, может содержать примесь в таком количестве, что она, накап­ ливаясь на фронте кристаллизации, приведет к снижению на нем температуры плавления. При кристаллизации из раствора температура у фронта роста может из-за выделения теплоты кристаллизации повыситься настолько, что равновесная кон­ центрация там изменится. При одновременном учете обоих про­ цессов значения температуры и концентрации, входящие в гра­ ничные условия, меняются, хотя форма граничных условий ос­ тается прежней. Уравнения переноса также сохраняют свой вид,


I I I . ЗАДАЧА СТЕФАНА

399

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

независимости коэффициентов

D

и к от температуры и состава

среды.

 

Франк [54] исследовал кристаллизацию шара, растущего

из

переохлажденного расплава по закону f/ s в присутствии при­

меси, которая

снижает Г п л

на поверхности

раздела

фаз на ве­

личину

KiCs,

где

Cs есть

концентрация

примеси в расплаве у

фронта

роста, а

Кг — постоянная. Таким

образом,

требуется

одновременно

решить четыре уравнения:

1)

константы роста Аг

и Хс связаны

соотношением, учитывающим,

что они

описывают

перемещение одной и той же сферической

поверхности:

 

 

XTKI'^KCD'1';

 

(10.1)

2)

понижение

температуры

плавления

на поверхности шара

Тп

равно

 

 

 

 

 

 

Та. ш =

Ты — KiCs;

 

(10.2)

3) уравнение

(9.32) и 4) его аналог для

роста

из раствора, где

константы роста выражены через температуру

и концентрацию

на поверхности шара. Как показал Франк, эту систему уравне­ ний можно решить графически с помощью прибора с градуи­ рованными шкалами. Как показывает решение задачи для кон­ кретного числового примера, даже если концентрация примеси в объеме расплава невелика, то примесь все же способна на­ капливаться у фронта кристаллизации в заметном количестве, поскольку коэффициент диффузии мал. В итоге равновесная температура у фронта роста существенно понизится, что при­ ведет к снижению эффективного переохлаждения и в свою оче­ редь замедлит рост кристалла. Такая ситуация складывается, когда кристалл слабо захватывает или вовсе не захватывает примесь. При этом анализе предполагалось, что концентрация примеси или «растворителя» невелика. Подобным же образом можно проанализировать другой предельный случай — рост из раствора. Однако решение такой задачи показывает, что по­ правка, учитывающая изменение температуры вследствие вы­

деления теплоты

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

в общем мала по сравнению

с поправкой, вносимой в решение

задачи о росте из расплава

при учете влияния

примеси.

 

Задачи, подобные рассмотренной Франком, исследовал Иванцов [63, 77, 81]. В частности, он решил задачу о кристал­ лизации бинарного сплава из перегретого расплава для плос­ кого фронта [81]. Им'проведен такой же анализ для кристалли­ зации сплава со сферическим фронтом роста из переохлажден­

ного

расплава [77]. В обоих

случаях фронт перемещался

со

скоростью, пропорциональной

t1'. Решая задачу о направлен­

ной

кристаллизации, подобную

классической задаче Стефана

о


400 Р. ПАРКЕР. МЕХАНИЗМЫ РОСТА КРИСТАЛЛОВ

росте из чистого расплава, Иванцов нашел, по-видимому, не решая количественно самосогласованную задачу, результат, со­

впадающий с данными, полученными Франком для

сферы: из-

за малой подвижности примеси

ее концентрация у фронта роста

возрастает,

так что

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

снижается.

В качестве

примера

Иванцов

взял случай, когда

тугоплавкий

компонент сплава легче переходит в кристаллическое состояние (это эквивалентно ситуации, рассмотренной Франком, когда у

поверхности кристалла

возрастает концентрация легкоплав­

кого компонента, слабо

переходящего в кристалл). В итоге

перед фронтом кристаллизации формируется слой расплава с концентрацией легкоплавкого компонента выше равновесной, который вследствие этого оказывается переохлажденным, хотя на большом удалении от поверхности кристалла расплав пере­ грет, а температура на самой этой поверхности остается равно­ весной. Это явление обнаружили независимо от Иванцова и Франка Раттер и Чалмерс [82], назвавшие его концентрацион­ ным переохлаждением в отличие от обычного переохлаждения, связанного с понижением температуры. Речь о проведенном этими авторами анализе устойчивости формы кристалла еще пойдет ниже.

Любов и Темкин [83] тоже

исследовали задачу

о направлен­

ной

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

сплава

в

формулировке,

предложен­

ной

Иванцовым [81],

однако

они

не учитывали

распределение

температуры и ограничились

анализом распределения примеси,

предположив, что закон перемещения фронта известен пол­ ностью, вплоть до константы роста. Ими было рассчитано рас­ пределение примеси в кристалле при движении фронта кри­

сталлизации

со скоростью, пропорциональной

f/ a

и t.

2.

Учет

влияния

различия

плотностей

твердой и жидкой

фаз.

Гидродинамические потоки

в расплаве

или

пересыщенном

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

говорится

в гл. V I I I . Различие плотностей жидкой и

твердой

фаз легко учесть в рамках задачи Стефана*; его влияние

на

рост

кристалла

рассматривается в настоящей главе подробнее,

чем

в гл. V I I I .

 

 

 

Ясно, что при кристаллизации из расплава, когда для обра­ зования кристалла данного объема обычно требуется на не-


I I I . ЗАДАЧА СТЕФАНА

401

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

Подобный эффект в случае плоского фронта затвердевания исследовали Карслоу и Егер [51], а также Хорвей [84]; последний к тому же проанализировал влияние эффекта на кристаллиза­ цию цилиндра и шара из переохлажденного расплава. Видоиз­

мененное

уравнение

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

(9.1)

в

этом

случае за­

пишется следующим

образом:

 

 

 

 

 

 

 

дТ,

 

дТ,

дгТ.

п

дТ,

 

(10.3)

 

 

dt

1

дг

дг2

г

дг

 

 

 

 

 

Здесь

п =

0 для

плоского фронта, п =

1 для

цилиндра и п =

= 2

для

шара;

и — радиальная

скорость

течения

жидкости.

Движение жидкости, которая считается несжимаемой, описы­ вается двумя дополнительными уравнениями: уравнением непре­ рывности

 

 

 

r~n-§F

(>-"•") = 0

 

 

(Ю.4)

и уравнением Навье — Стокса

 

 

 

 

 

ди

.

ди

' д2и .

I ди/дг

и \1

I

до

, , А г ч

dt

'

дг

 

 

 

 

 

 

где р — давление

жидкости, рь — ее

плотность,

a

v — кинема­

тическая вязкость.

 

 

 

 

 

 

Хорвей

показал, что взятое

в квадратные скобки

выражение

в уравнении (10.5) равно нулю, следовательно,_ поток не за­ висит от вязкости v; кроме того, движение жидкости потенци­ ально. Наряду с условием (9.4) должно выполняться условие баланса масс на фронте кристаллизации:

dR

 

4 г - « д о ] .

< 1 0 - 6 )

Здесь ps — плотность кристалла. Решение

сформулированной

задачи показало, что все три рассмотренные формы растут по

закону

R = 2^(%bt)'h.

Константа роста

|3

в случае

плоского

фронта

выражается

в виде р = (1 -4-е )~'Р';

Р " находится реше­

нием уравнения

cL

лл

— Т^) Ln4*$'e®'

erfc ф'), полностью со­

впадающего по

форме

с уравнением (9.22);

параметр

е, отра­

жающий различие плотностей твердой и жидкой фаз, равен (ps/рь) — 1. Хорвей показал, что р для шара и цилиндра слабо зависит от параметра 8 при малых его значениях.

Из решения задачи, однако, следует, что давление жидкости при направленной кристаллизации очень велико, поскольку на­ чальная скорость жидкости (при t = 0) бесконечно велика,