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

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

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

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

Добавлен: 14.10.2024

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

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

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

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

393

Здесь

gi >

а >

| 2

>

Ь >

£3 >

0;

а

\\ b — постоянные.

Следуя

Морсу

и Фешбаху

[72], преобразуем

уравнение

(9.36)

к виду

 

^

y

f

i

- ^

=

^

_

V ± »

^ .

(9.45)

 

 

2 ^

\ К

dl,

h

 

 

 

 

 

Здесь hi (угловые масштабы) представляют собой известные

функции gj. Координата \\ принимает на

поверхности

эллип-

соида

с главными

полуосями длиной (£0

а )

t , (| 0

— о J

г

и

|о^'/ 2

постоянное

значение £о-

Решение,

зависящее

только

от

| i ,

удовлетворяющее уравнению

(9.45)

и

принимающее

значе­

ние Со на такой эллипсоидальной поверхности

и значение

С»

на бесконечности, выражается

формулой

 

 

 

 

 

 

 

 

С а,) = С0 + (С„ -

 

Со) { 1 -

 

 

} ,

 

(9.46)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

f (g,) =

J е-< 2 /*° [(г2

-

а2 ) (*2 -

б 2 ) ] - 7 2

dt.

 

 

 

s.

Чтобы определить £0, обращаются к условию постоянства пото­ ков на границе фаз:

 

\

т - С0 ) У (Ы = -

D т е - С0 ) (^F )

(9-47)

Здесь

Ст,

как и прежде, есть плотность диффундирующего ве­

щества

в

выделяющейся фазе.

Как показал проведенный Хэ-

мом [70] анализ, эллипсоидальный фронт кристаллизации пере­ мещается пропорционально t'!* (то же самое верно и для сфе­ рического, цилиндрического и плоского фронтов, координатные системы для которых представляют собой частные случаи

эллиптической системы координат). Надо

заметить еще и то,

что решением задачи служит эллипсоид с

любым отношением

полуосей а и Ь. Таким образом, если рост

лимитируется диф­

фузией, то при одной и той же концентрации исходного ра­

створа сохраняется

форма

любых

эллипсоидальных частиц.

К этому

выводу

мы

еще

вернемся

при изучении дендритного

роста. Частный

случай кристаллизации эллипсоида — рост эл­

липсоида

вращения (сфероида) рассмотрели Хорвей и Кан [73],

а также

Иванцов [74].

 

 

 

Параболоидальные дендриты. Под «дендритом» понимают кристалл в форме дерева, или с признаками древовидное™, прежде всего ветвления. Параболоидальными называют денд­ риты, для которых форму вершин стволов и ветвей с некото­ рой степенью приближения можно считать параболоидальной.


394

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

 

 

/.

Кристаллизация

параболоида

вращения

из

переохлаж­

денного

расплава.

В исследовании,

ставшем

 

классическим,

Папапетру [75] предположил,

что

вершина дендрита

имеет

форму

параболоида

вращения.

Иванцов

[61, 63]

первым

решил

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

Метод Иванцова, подробно рассмотренный Хорвеем и Каном [73], состоит в том, чтобы найти решение нелинейного диффе­ ренциального уравнения сохранения теплоты или вещества на движущейся границе [см. уравнение (9.4)], удовлетворяющее также уравнению диффузии (9.1) и граничным условиям (9.3), (9.6) и (9.7). Обще решение уравнения (9.4) содержит пять произвольных постоянных и одну произвольную функцию, кото­ рые определяются из граничных условий, из уравнения (9.1) и из соотношений, связывающих постоянные. Для дендрита, ра­ стущего в направлении оси г, распределение температур в рас­ плаве выражается формулой

TL(B)-T„ E l [ - 0 p , g / 2 * J

Здесь

Гоо температура

расплава

при

бесконечно

большом г,

pt — радиус

кривизны

вершины

параболоида,

v — постоянная

скорость

роста, xL — температуропроводность

расплава; функ­

ция

Ei

протабулирована

[76]

[при

малых

х функция

Ei(—x)

«0,577+lrw], а

 

 

 

 

 

 

 

 

 

 

 

 

VI

.

J_\2

+/Z-Vt\2

 

 

 

 

 

 

 

г

 

 

 

 

 

 

Величина

vgtl2%i. определяется

из

уравнения

 

 

 

 

ч

( г п л - и

=

Zlh.

_

 

\ . п ; (-

щ

 

 

 

 

L

 

2 х ,

 

 

 

 

 

 

Как и в случае эллипсоида, решение задачи о росте парабо­ лоида определяет величину произведения vpt, а не значения каждой из этих величин в отдельности. Следовательно, при одинаковых температурных условиях осуществляется одна из двух возможностей: либо медленно растет тупая игла, либо быстро растет острая игла (температура на поверхности и внутри обеих игл равна температуре плавления). Выражение


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

395

(9.49) Иванцов привел также в виде, отвечающем росту из рас­

твора

[63,

77].

 

 

 

 

2.

Кристаллизация

эллиптического

параболоида

из

переох­

лажденного

расплава.

Хорвей и Каи

[73] обобщили

метод

Иван-

цова [61], применив его к расчету роста параболоида с эллипти­ ческим сечением. На фиг. 8 приведена связь безразмерной скорости роста 2Q = vpt/xL с формой (т. е. с отношением полу­ осей) эллиптического сечения; параметр А2 зависит от безраз-

 

 

W'5

Ю~*

Ю'3

Ю'г

0,1

 

1,0

Ф и г .

8.

Зависимость приведенной

скорости роста

2Q

от

приведенного

переохлаждения

Uf для ряда значений эллиптичности

А2

при росте эллип­

 

тических параболоидов из переохлажденного расплава [73].

мерной степени переохлаждения Uf—cL(TUn—T^/L.

 

 

Оказывает­

ся,

что

для

параболического

цилиндра 2

0)

скорость Q

меняется при малых переохлаждениях приблизительно пропор­

ционально

U2f

(см. ниже),

а для

иглы с

круговым сечением

приблизительно

пропорционально

U{.

Как

и прежде,

извест­

но лишь произведение upt.

 

 

 

 

 

Другие формы роста. 1.

Цилиндр

с

параболическим

сечени­

ем. Хорвей

и Кан [73], а также Иванцов

[74] исследовали

задачу

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


396

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

КРИСТАЛЛОВ

радиус кривизны

вершины); как и в трехмерном

случае,

 

Х

пл -

T J = (яр')'/ 2 ер'

erfc

(9.50)

где p'=vpt/2xL

(так называемое число

Пекле Q).

 

2. Многогранники

и

многоугольники.

Кристаллизацию пра­

вильных многоугольников, куба и октаэдра из пересыщенного раствора исследовал Зегер [78], а из переохлажденного рас­ плава — Иванцов [63]. Зегер, в частности, ограничился рассмот­ рением кристаллизации с постоянной скоростью при очень ма­

лых концентрациях и пользовался стационарным

уравнением

диффузии

 

V2C = 0.

 

(9.51)

 

 

 

Поскольку

предполагается,

что скорость

роста грани постоянна

как во времени, так и вдоль поверхности грани,

правая часть

условия непрерывности потока (9.4) или (9.38)

остается по­

стоянной,

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

напротив,

меняется вдоль поверх­

ности кристалла; этим рост гранной формы отличается от роста других уже рассмотренных форм. Разность между максималь­ ной концентрацией на углу квадрата и минимальной концент­ рацией в середине стороны квадрата составляет примерно3,5% от значения интеграла нормальной составляющей градиента

концентрации1 )- Грубая оценка

показывает, что концентрация

у вершины отличается от концентрации в центре грани

куба

примерно на 2,5—3%. Иванцов

[63], также исследовавший

кри­

сталлизацию с очень малой линейной скоростью, установил, что при росте из расплава центр грани — самый горячий участок, а вершина — самый холодный. Этот результат качественно со­ гласуется с выводами Зегера. Однако Иванцов учел зависи­ мость температуры поверхности кристалла от времени, что за­ трудняет установление количественного соответствия его рас­

четов с данными

Зегера.

 

3. Иглы иной

(не параболоидальной)

формы. Фишер [79] (см.

также [80] или [52]) исследовал кристаллизацию цилиндра со сферической вершиной радиусом pt, растущего вдоль своей оси со скоростью v из переохлажденного расплава. Используя при­ ближенное выражение для распределения температуры и пола­ гая, что выделение скрытой теплоты примерно равно теплоотводу, автор получил (см., например, работу Хорвея [52]) вы­ ражение

T-Wn.m-TJ'-^.

(9-52)

'J То есть произведения нормальной составляющей градиента на длину грани. — Прим. ред.


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

397

в котором через Та. ш обозначена считающаяся

неизвестной

температура поверхности раздела фаз на вершине иглы. Иван-

цов [63]

исследовал

кристаллизацию иглы

в форме

карандаша

и иглы

квадратного

сечения при росте с

постоянной

скоростью

из переохлажденного расплава и определил считавшуюся неиз­ вестной температуру вершины в функции скорости роста. Ре­ зультат, полученный Иванцовым, хорошо согласуется с расчетом Фишера для цилиндра с полусферической вершиной.

Хорвей [52] в своем весьма подробном анализе модифициро­ ванной задачи Стефана (форма и скорость роста заданы, а температура поверхности кристалла считается неизвестной)

исследовал

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

ци­

 

 

 

линдра

со

сфероидальной

 

вер­

 

 

 

шиной.

Включив

рассмотренную

 

 

 

Фишером

форму

 

как

 

частный

 

 

 

случай

и

решив

задачу

более

 

-Кристалл

точно, он нашел, что коэффи­

 

 

 

циент

 

перед

 

vpt/2-кь

в выражении

 

 

 

(9.52)

 

равен

0,94,

а

не

1,0.

 

 

 

 

 

4.

Прочие

 

формы

роста.

 

Уил-

 

-

Тигель

кокс

[48]

исследовал

так

назы­

 

ваемую

«форму

Чохральского»,

 

 

 

т. е.

форму

поверхности

раздела

 

Расплав

фаз при выращивании

кристалла

 

 

 

 

методом

Чохральского

 

(фиг.

9).

 

 

 

Предполагалось,

что

кристалли­

 

 

 

зация

 

идет

с

постоянной

скоро­

Ф и г . 9.

Схема выращивания

стью,

расплав

перегрет,

а

темпе­

кристаллов

вытягиванием

из ра­

ратура

фронта

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

сплава

по Чохральскому.

постоянна

и

 

равна

температуре

 

 

 

плавления. Численный расчет показал, что поверхность кристал­

ла вогнута, форма ее близка к параболоиду,

а

кривизна зави­

сит от условий теплообмена кристалла

с

окружающей

средой.

Хорвей [52] провел детальный численный расчет кристалли­ зации пластины с эллиптической вершиной (двумерный слу­ чай) .

10.УСЛОЖНЕННАЯ ЗАДАЧА СТЕФАНА;

РА З Л И Ч Н Ы Е ФОРМЫ

Имеется несколько типов усложняющих факторов, которые можно ввести в уже рассмотренную «чистую» задачу Стефана. Разумеется, усложнения вводятся для того, чтобы привести выбранную модель в более точное соответствие с действитель­ ностью; среди таких факторов надо назвать следующие: влия-