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

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

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

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

Добавлен: 14.10.2024

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

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

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

V I I I

ГИДРОДИНАМИЧЕСКИЕ ЭФФЕКТЫ

32. В В Е Д Е Н И Е

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

сталлов из текучей среды, т.

е. из раствора, пара и

расплава.

В свою очередь жидкую среду

приводят в движение

различные

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

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

сти кристалла имеют более высокие значения, чем в

отсут­

ствие гидродинамического течения. Возрастают и скорости

роста,

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


V I I I . ГИДРОДИНАМИЧЕСКИЕ ЭФФЕКТЫ

fill

способно гораздо сильнее снизить концентрацию примеси у фронта кристаллизации, поскольку диффузионный слой при та­ ком перемешивании становится гораздо тоньше, так как его толщина определяется толщиной гидродинамического погранич­ ного слоя (см., например, [251, 252]).

Хотя гидродинамические эффекты во многих случаях кристал­ лизации играют не менее важную роль, чем теплопроводность и диффузия, все же количественных исследований в этой области проведено, по-видимому, несравненно меньше, чем по диффу­ зионной задаче Стефана и другим аналогичным задачам. [Впро­ чем, к исследованию влияния перемешивания на некоторые кри­ сталлизационные процессы часто прибегают в химической про­ мышленности (см., например,[280, 281]).] Такую недооценку мож­ но объяснить несколькими причинами. Во-первых, если в опытах жидкая фаза специально не перемешивается, то о возможности естественной конвекции часто забывают. Во-вторых, если суще­ ствование потоков жидкости при росте кристалла и учитывалось,, то громоздкость дифференциальных уравнений удерживала ис­ следователей от попыток определить аналитически или хотя бы полуколичественно распределение таких потоков. В-третьих, в опытах по росту кристаллов часто невозможно наблюдать за распределением потоков жидкости, особенно если последняя непрозрачна.

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

ной проблеме, на группы, приняв за

основу классификации

форму или геометрию системы; такая

классификация часто

употребляется в литературе по гидродинамике, теплопроводно­ сти и диффузии (основная литература: Ландау и Лифшиц [212], Гольдштейн [282], Шлихтинг [283], Эккерт и Дрейк [284]; см. также журнал [285]).

33. ОБТЕКАНИЕ ПЛОСКОЙ ПЛАСТИНЫ

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

подобного

рода относятся работы Карлсона [286], Беннемы [182]

и Брайса

[287].

Главным средством исследования распределения потоков

жидкости

служат в принципе уравнения Навье — Стокса (см.


512

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

Ландау и Лифшиц [212], Шлихтинг [283]). Это сложные нелиней­ ные уравнения движения, точные решения которых известны только для особых случаев (потоки в трубе и около вращаю­ щегося диска, очень медленное обтекание шара). Как показал Прандтль, когда поток жидкости обтекает твердые тела, в слое жидкости, прилегающем к их поверхности и называемом погра­ ничным слоем, возникают большие градиенты скорости течения. Рассматривая движение жидкости в этом слое, следует учиты­ вать трение (вязкость); вне пограничного слоя трением можно пренебречь. Такой приближенный анализ позволяет упростить уравнения движения жидкости в пограничном слое, которые все

«со

Ф и г . 40. Гидродинамический пограничный слой при обтекании плоской пластины [283].

еще остаются нелинейными. В случае стационарного двумерного потока несжимаемой жидкости с постоянными характеристиками эти уравнения имеют следующий вид:

 

 

« & + » • & = - < 3 3 1 »

 

 

 

 

& + £ - а

 

 

 

(зз.2)

Здесь

и и

и — составляющие скорости

потока по

осям

х

и у;

v — кинематическая вязкость, v — ц1рь\

Ц — вязкость; pL

— плот­

ность; считается, что давление р известно из решения для по­

тенциального течения невязкой жидкости. При обтекании пла­

стины

(фиг. 40)

dp/dx

= 0 и граничные

условия: и = v

= 0

при

у = 0,

и =

Uoo при у

= оо. Точное решение

этой задачи

все

еще

сопряжено со значительными трудностями; первым ее решил

методом разложения

в ряд Блазиус (см. [283]). Приведем здесь

лишь результат, полученный в [284] более простым методом,

основанным

на

использовании уравнения

количества

движе­

ния пограничного слоя. Распределение

скоростей

описывается г-


V I I I . ГИДРОДИНАМИЧЕСКИЕ ЭФФЕКТЫ

513

в этом приближении формулой

и _ 3

у

 

 

(33.3)

2

6

2

16

 

в которой толщина пограничного слоя б определяется следую­ щим образом:

 

 

± = -Щг

 

 

 

(33-4)

(число Рейнольдса Rex =

UocXJv).

В этом приближении б — рас­

стояние от пластины,

при

котором

и =

и,*,; в действительности

и приближается к их

асимптотически.

В

точном выражении

(33.4) для б, при котором

и = 0,99их,

стоит

коэффициент

~5,0

вместо 4,64. По формуле

(33.4)

толщина

пограничного

слоя

равна нулю у края пластины, встречающего поток, и возрастает пропорционально х"к с удалением от этого края.

Для определения потока тепла к пластине или от нее сна­ чала выводится уравнение энергетического баланса для элемен­ тарного объема жидкости. Это уравнение в принципе следует решать совместно с уравнениями Навье — Стокса. Более простой его вариант в приближении пограничного слоя записывается в следующем виде:

I дТ

дТ\

К, д2Т

/ди\2

 

 

C ^ { u - ^ + V W ! = ~ f W + V \ d y - )

'

( 3 3 - 5 )

Здесь K L — теплопроводность

жидкости, a

CLP — ее

удельная

теплоемкость при постоянном давлении. Кроме граничных усло­

вий, связанных

с уравнениями (33.1)

и (33.2),

должны

выпол­

няться еще и

граничные

условия,

связанные

с уравнением

(33.5) , а

именно

Т = Та.ш

при у = 0 и Т =

при у =

со. Если

свойства

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

скоростей

не связано с полем температур, так что можно

сначала

найти

решение

уравнений (33.1)

и (33.2), а

затем

решить

линейное

уравнение (33.5) относительно Т(х, у). Польхаузен (см. [283]) получил точное решение уравнения (33.5), использовав решения уравнений (33.1) и (33.2). Если тепловой пограничный слой б* тоньше гидродинамического пограничного слоя б, то приближен­ ное решение имеет вид

 

т-тп.ш

_ з у

 

1 / ^ 4 3

 

где

Т „ . ш ~ 2 6,

2 \ 6 , / '

( 3 3 ' 6 )

 

 

 

 

 

 

 

=

 

гг-

 

(3 3 -7)

 

6

1,026

 

( Р г ) / з

 

 

Здесь число

Прандтля

Pr =

V / X L ;

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

Х-L = KJPLCLP-

П О Т О К

тепла

 

через

поверхность

нагретой

17 Зак. 718

 

 

 

 

 

 


514

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

пластины

единичной площади выражается

в виде

 

Q = 4 4 f O . - T V J ;

(33.8)

Воспользовавшись выражениями (33.7) и (33.4), можно записать

Q = 0,332^ ( Г м - Г„.ш ) (Рг)'/ з ( ^ ) ' / 2 ,

(33.9)

т. е. поток Q обратно пропорционален корню квадратному из расстояния от данной точки до края, встречающего поток. Заме­ тим, что, согласно (33.7), в вязкой жидкости, т. е. при Рг 1000, толщина теплового пограничного слоя составляет меньше '/ю толщины гидродинамического пограничного слоя. Для воды при комнатной температуре Рг « 7; для многих газов критерий Рг близок к единице. Для жидких металлов Рг <С 1 (например, для ртути при комнатной температуре Рг « 0,02), но тогда формула (33.9) уже неприменима.

Если происходит перенос вещества, например при росте из разбавленного раствора или при растворении, то [282] поток ве­ щества <3вещ выражается формулой [286]

Фвещ - 0.332D (С„ - Св) (Sc)7 ' ( - ^ ) ' / ! ,

(33.10)

где Sc = v/D — число Шмидта; D — коэффициент диффузии рас­ творенного вещества; С3 — концентрация на поверхности; — концентрация на бесконечности. В случае диффузии простых молекул в воде критерий Sc есть число порядка единицы. Тогда толщина диффузионного пограничного слоя бс запишется в соот­ ветствии с формулами (33.7) и (33.4) в виде

 

^ = w ( S c r ' / s f e ) " ' / 2 -

( З З Л 1 >

или

бс = 4,5D'l3v'uuZ'l2x'h.

Брайс [287],

Карлсон

[286] и Беннема

[182]

пользовались этими

формулами

для толщины диффузион­

ного пограничного слоя около плоской пластины при истолкова­ нии результатов измерения скоростей роста кристаллов в пере­ мешиваемом растворе с известной скоростью протекания иж.

Брайс [287] проанализировал роль диффузионного погранич­ ного слоя при росте кристаллов из раствора и, в частности, во­ прос о том, каким образом следует учитывать существование этого слоя при определении зависимостей скорости роста кри­ сталла от пересыщения из экспериментальных данных. Зависи­ мость скорости роста от пересыщения можно записать, по край­ ней мере при линейной и квадратичной кинетике, в виде

Q =

A(CS-C0)n,

(33.12)