ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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 (произведение скорости на
I I I . ЗАДАЧА СТЕФАНА |
397 |
в котором через Та. ш обозначена считающаяся |
неизвестной |
температура поверхности раздела фаз на вершине иглы. Иван-
цов [63] |
исследовал |
кристаллизацию иглы |
в форме |
карандаша |
и иглы |
квадратного |
сечения при росте с |
постоянной |
скоростью |
из переохлажденного расплава и определил считавшуюся неиз вестной температуру вершины в функции скорости роста. Ре зультат, полученный Иванцовым, хорошо согласуется с расчетом Фишера для цилиндра с полусферической вершиной.
Хорвей [52] в своем весьма подробном анализе модифициро ванной задачи Стефана (форма и скорость роста заданы, а температура поверхности кристалла считается неизвестной)
исследовал |
кристаллизацию |
ци |
|
|
|
||||||||||
линдра |
со |
сфероидальной |
|
вер |
|
|
|
||||||||
шиной. |
Включив |
рассмотренную |
|
|
|
||||||||||
Фишером |
форму |
|
как |
|
частный |
|
|
|
|||||||
случай |
и |
решив |
задачу |
более |
|
-Кристалл |
|||||||||
точно, он нашел, что коэффи |
|
|
|
||||||||||||
циент |
|
перед |
|
vpt/2-кь |
в выражении |
|
|
|
|||||||
(9.52) |
|
равен |
0,94, |
а |
не |
1,0. |
|
|
|
|
|
||||
4. |
Прочие |
|
формы |
роста. |
|
Уил- |
|
- |
Тигель |
||||||
кокс |
[48] |
исследовал |
так |
назы |
|
||||||||||
ваемую |
«форму |
Чохральского», |
|
|
|
||||||||||
т. е. |
форму |
поверхности |
раздела |
|
Расплав |
||||||||||
фаз при выращивании |
кристалла |
|
|||||||||||||
|
|
|
|||||||||||||
методом |
Чохральского |
|
(фиг. |
9). |
|
|
|
||||||||
Предполагалось, |
что |
кристалли |
|
|
|
||||||||||
зация |
|
идет |
с |
постоянной |
скоро |
Ф и г . 9. |
Схема выращивания |
||||||||
стью, |
расплав |
перегрет, |
а |
темпе |
|||||||||||
кристаллов |
вытягиванием |
из ра |
|||||||||||||
ратура |
фронта |
кристаллизации |
сплава |
по Чохральскому. |
|||||||||||
постоянна |
и |
|
равна |
температуре |
|
|
|
плавления. Численный расчет показал, что поверхность кристал
ла вогнута, форма ее близка к параболоиду, |
а |
кривизна зави |
сит от условий теплообмена кристалла |
с |
окружающей |
средой.
Хорвей [52] провел детальный численный расчет кристалли зации пластины с эллиптической вершиной (двумерный слу чай) .
10.УСЛОЖНЕННАЯ ЗАДАЧА СТЕФАНА;
РА З Л И Ч Н Ы Е ФОРМЫ
Имеется несколько типов усложняющих факторов, которые можно ввести в уже рассмотренную «чистую» задачу Стефана. Разумеется, усложнения вводятся для того, чтобы привести выбранную модель в более точное соответствие с действитель ностью; среди таких факторов надо назвать следующие: влия-