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

книги из ГПНТБ / Процессы взаимной диффузии в сплавах

..pdf
Скачиваний:
22
Добавлен:
25.10.2023
Размер:
17.2 Mб
Скачать

gii

= бі/ +

д In С ;

 

 

(2.81)

 

 

 

 

 

(2.82)

 

 

 

a in

»1

 

Если пренебречь, как

обычно,

величиной

д In c;-

TO

 

 

 

 

 

(2.83)

Некомпенсированный поток вакансий приводит к эффекту Киркендалла, причем, как обычно, скорость смещения кристал­ лических слоев как целых (т. е.— U) определяется формулой

В результате наблюдаемый поток J І 1-ГО компонента дается выражением

Ji=Ji

±U =

Ji-ci'2lJj.

(2.84)

 

J

/ - і

 

Воспользовавшись соотношением (2.78), получим

j=i

I=i л = і

 

8 - і

S

 

V 1-і

( 1 — СІ)Dij - С І V> D у

«ЗА: (2.85)

i=lL

(=1

 

В случае бинарной системы ранее мы получили

h =

ET^Vc i.

т. е. наблюдаемый поток /-го компонента определялся одним па­ раметром— коэффициентом взаимной диффузии D, равным D = C\D2-\-c2Di. В s-компонентной системе наблюдаемый поток определяется аналогичными по физическому смыслу параметра­ ми, которые поэтому также можно назвать коэффициентами вза­ имной диффузии, но число их уже равно не единице, а 51. Дру­ гими словами, можно написать

1 5 - 1

~

 

(2.86)

J i = - — 2

6«Vc />

 

і—і

 

 

 

где

 

 

 

Di, = (1 - СІ) Da - d Z

Du.

(2.87)

Этот результат уже был получен в феноменологической тео­ рии (см^ формулу (1.127)). Новое здесь заключается в том, что теперь Dij можно представить в явном виде через микроскопи­ ческие параметры системы.

§8. Вывод результатов Маннинга

вформализме метода дырочного газа

Легко показать [103], что в общем виде результаты Маннин­ га, изложенные в § 5 гл. I , можно получить и в формализме ме­ тода дырочного газа. Рассмотрим снова бинарную систему Л—В с объемноцентрированной кубической структурой. Тогда плот­ ность потока дырок в узле 0 запишется опять в виде (деформа­ цией решетки мы пренебрегаем)

J(0) = 4"il

{c(0)u0 / +

c(/)u/ 0 }

 

 

 

 

7=1

 

 

 

 

 

1

w m ^ r I і

1

1

д с { 0 ) )

Величины

1 /тоj выражают вероятность

перескока

дырок в еди­

ницу времени между соответствующими узлами, т. е. среднюю частоту таких перескоков. Следуя Маннингу, учтем теперь, что частота перескоков по направлению нескомпенсирозанного пото­ ка вакансий и против него различна.

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

J(0)

= - i - c ( 0 ) | / o v a

 

''аО

 

L0a

 

 

1

/

dc(0)\ 1

 

w

 

аО

+

fop

 

 

 

L 03

Lpo

 

п' +

W X

' LxB0

где индексы а и |3 указывают любой из узлов /, расположенных правее или левее узла 0 соответственно, а 1/т0°а и т. д. означают частоту перескоков в отсутствии некомпенсированного потока вакансий, т. е. в прежних обозначениях

- L = f ( a ) = f(0) + r O a M

Oa

И Т . Д.

Тогда, поскольку

Гоа =

—Гор =

а/2,

 

7 = / А + / Б ,

находим

j ( 0 ) = ^ C ( o ) ; f(0(l+W)+~

 

 

^

(

1

 

+

ш)

 

 

 

 

 

 

f(0)(l-w)

_ і _ _ | _ а ^ )

/ ( 0 ) ( 1 _ ш ) _ / ( 0 )

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

_ | _ е т ( 1 _ ш ) + / ( 0 ) ( 1 +

й ;)

 

 

 

 

 

 

 

 

 

 

=

£ В { 2 а д а _ 2

а % ^ _ ^ /

(

0

)

+

 

8 й / ( 0 ) ш |

 

 

 

=

і - 2а"-с(0)/л(0)

[a in / л (0)

 

5 in с

(0)

,

4

 

+

 

 

 

 

 

 

 

 

+

 

W

 

4 - 2а 2 с(0)/ в (0)

d l n / B ( 0 )

ainc(0 )

 

,

 

w

 

 

 

 

 

 

дх

 

дх

 

 

і

а

 

 

 

 

 

= - 1 { 2 а 2 С ( 0 ) / л ( 0 )

 

 

 

 

 

dlnc(O)

\

дАс

 

j _ _4_

W

СІ

д In с.

 

 

 

д In с .

 

 

 

 

 

 

 

 

1 ^ l n / B ( 0 )

 

a In с (0)

^

.

_4_

 

2 с(0)/ в (0) [ с в \

 

д\псЕ

 

 

 

д In сг.

 

дх '

" w

Учтем теперь, что 2a"cfA = cAD*A,

 

2 а 2 с / в

 

cBDB,

 

 

 

 

d l n f A

a In с _

 

 

 

d l n / B

 

a In с

 

gBB

 

 

d In с л

5 In cA ~

&AA'

 

 

д\псв

 

 

д In

cB

 

 

и £ В В = £ А Л . Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 Г *

дсА

 

* 4

 

 

 

 

*

^с в

 

 

* 4

J = — [ А 4 £лл - д р + сл DA — ЬУ + D B gBB

 

 

+ с в D B w »

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.88)

 

J=±(DA-D]

 

 

дсл

 

 

^{cADA-\-cBD*B)w}.

 

 

(2.89)

 

 

 

дх

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вместе с тем, так как / определяет скорость U в эффекте Киркендалла, исходной причиной которого является неодинако­ вая диффузионная подвижность атома Л и В в отсутствие на­ правленного потока вакансий, то из общих физических сообра­ жений, очевидно, должно быть справедливым соотношение

J = X-Jr[DA- D B ) g A A дх^ ,

причем і не может явно зависеть от температуры и состава. Сравнивая последнее соотношение с формулой (2.89), сразу

получаем

w = Р

•D,

дс.

(2.90)

4 cADA + cBDB ёАА'Э*

'

где р = х — 1 .

.112

£>л=-=

DAgAA(l+CAK),\

(2.91)

DB

= DB gBB

(1 СвК), J

(2.92)

Анализ структуры

правой

части формулы (2.S8) с учетом

явного выражения для w (2.90)

дает

 

где

Таким образом, мы получили результаты Маннинга, изложен­ ные в § 5 гл. I (см. формулы (1.102), (1.103)).

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

§ 9. Взаимная диффузия при наличии градиента температуры и внешнего силового поля

Процесс диффузии при наличии градиента температуры в случае чистых металлов рассматривался рядом авторов [86, 104—108]. В основу этих работ положена методика, разрабо­ танная в [109]. Терімодиффузия в бинарных системах рассмот­ рена в работах [22, 110, 111] на основе описанного выше фор­ мализма модели дырочного газа.

Рассмотрим наиболее общий случай непостоянства по объе­ му состава и температуры в бинарном твердом растворе А— В. Этот случай соответствует взаимной диффузии, усложненной

эффектом

термодиффузии.

 

Вообще

говоря, в силу соотношения сА-\-св-\-с=\

И зависи­

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

бых температурах С-ССА И С<ССВ И, следовательно, приближен­

но можно считать, что состав системы задается независимо от

температуры, т. е. параметры

системы Е\, QA,

QB, т\[А>,

г\(А), г\[вК

•ч^» V A

, v s

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

зависят

только

от

одного параметра

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

сА

 

или

концентра­

ции св),

 

характеризующего

состав системы,

и

не

зависят от

температуры. Это, в частности, означает, что

в

данной

задаче

можно

положить

 

 

 

 

 

 

 

 

с . г + С п = 1,

усв — — \'СА .

 

 

 

 

 

8 Под ред. К. П. Гурова

113

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

сов. В таком диапазоне

изменение указанных

выше

параметров

системы

должно

быть

пренебрежимо

малым.

Для чистых ме­

таллов

в этом

можно

непосредственно убедиться,

например,

для Q, используя связь

между

Q и дебаевской характеристи­

ческой температурой металла в

(см. Приложение I I ) .

Весь формализм вывода выражения для плотности потока

дырок в

узле 0

остается прежним с

тем только

изменением,

что теперь надо учитывать зависимость температуры от гидро­

динамической

координаты.

 

 

 

 

 

 

 

 

 

 

Из исходной формулы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J =

V

{p(0)uo/-|-p(/)u/ol,

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J =

2 r o / { p ( 0 ) / ( / ) - p ( / ) / ( 0 ) J ,

 

 

 

 

 

где

 

 

 

7=1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

/ (/) = ІА (І) +

їв (j),

U (/) = : r

- - Y

°л 0"> ^

(і) схр { -

Щу),

 

 

 

 

 

ф л ( / ) =

^ ) ^ ( / ) '

 

 

 

 

 

 

снова получаем в случае одномерного процесса

(вдоль

оси х)

для объемноцентрированной кубической решетки

 

 

 

 

/ =

2 М0)р(0) ( —

 

 

а Г ^ - ] +

 

 

 

 

 

 

 

 

 

 

 

 

+

2а-/в

(0) р (0) [ — Т х

 

 

 

 

Заметим, что теперь

 

 

 

 

 

 

 

 

 

 

 

 

d\njA

 

і [дсА

( dlnq>AdcA

 

\ dQA дсл\

t

QA

дТ

 

 

 

 

дх

А~\~дх

*dlncA

дх

кТд\псАдх

)~>

kT2dx'

 

 

din

р _

dine

_

gin и

1_ a In с дсА

a In с дТ _

 

 

 

 

 

 

дх

~~

дх

 

дх

сАд

In сА

дх

1

дТ дх

a In У дсА

 

 

 

 

 

 

 

 

 

 

 

 

1

діп

у дТ

 

 

 

 

 

 

 

 

 

сА

din

сА

дх

 

дТ

дх '

Как и выше, пренебрегаем эффектом деформации решетки при изменении состава. Можно также пренебречь тепловым расширением решетки при термодиффузии.

Тогда

 

 

 

 

 

a In р

 

1

 

дЕ1

дсА

 

Ел

дТ

 

 

 

 

 

 

 

 

дх

~~

 

cAkT

 

д In сА

дх

'

kT2

дх "

 

В результате

получим

 

 

 

 

 

 

 

 

 

 

 

 

г

=

1

п* (Л .

д *

g » ^ дсл

 

г

E i ~ Q a

1

д Т \ .

 

J

 

~

° А ( 1

1 +

ТиГ^

)

 

- с ^

- к ?

 

ТЩ

 

 

 

 

 

 

1

г»* (Л

а 1

, , 1 р в \

а с в

 

£ 1 - < Эв і

ат

 

 

 

 

7 D B

 

' " ' " Т Ь Г ^

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ч'л =

фд ехр I

kT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда для

потоков

компонентов

(например,

компонента А)

имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где D A — обычный парциальный коэффициент диффузии этого компонента, a D A — коэффициент его термодиффузии, опреде­ ляемый как

DTA= E I ~ Q A D A .

(2.95)

Здесь следует сделать одно замечание. Если бы мы в разло­ жении в ряды Тейлора сохраняли квадратичные члены, то, как и раньше, все они взаимно сократились бы. Поэтому и в случае термодиффузии поправка к линейному приближению на два по­ рядка выше по параметру малости.

Легко дать физическую интерпретацию последнему соотно­ шению для коэффициента термодиффузии компонента А. Под­ вижность в условиях термодиффузии характеризуется двумя факторами: подвижностью при самодиффузии (т. е. коэффици­ ентом D A ) И результатом конкуренции двух процессов, выра­ жаемым множителем —г ~kT^A • Возникновение таких конкуриру­ ющих процессов объясняется следующим. Пусть в двух сосед­ них узлах температуры различны; скажем, в правом температура выше. Тогда в правом узле концентрация дырок больше, что приводит к направленному потоку атомов слева направо; но в этом узле и «активированных» атомов больше, что приводит к направленному потоку атомов в обратном направлении. Таким образом, возникают два конкурирующих процесса.

Поскольку. D A и D B В общем случае различны, даже в пространственно однородной бинарной системе термодиффузия

8*

115

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

цессом взаимной диффузии.

 

 

 

 

 

 

 

 

 

 

Наличие некомпенсированного

потока

дырок

означает, что

в системе должен возникнуть эффект Киркендалла,

причем

характерная для него скорость определяется

 

по формуле

 

U — — (DA

~DB)d-~

+ {cADTA + cBDTB)

 

 

~r|Г.

 

(2.96)

Наблюдаемый

поток компонента

 

А теперь имеет вид

 

J ' A = J А - Ц- U = - ~ (cBDA

+ c A D B ) ^ +

 

 

 

 

 

 

1 г rJn]

T

i_ п'г\ 1 д Т

-

-

1

пд°А

і

1

~

п г 1

д Т

/ о с т

~-CACB[DA-TDB)YW-

 

 

~

D - ^ -

V

 

D

y g p

(2.97)

где формально введен «эффективный коэффициент термодиф­ фузии»

D T = сАсв {DTA + DTB).

(2.98)

Таким образом, в рассматриваемом процессе перераспреде­ ление состава характеризуется не парциальными коэффициента­ ми диффузии, а едиными для системы в целом диффузионными

параметрами D и DT, т. е. существует далеко идущая

аналогия

со случаем одной лишь взаимной диффузии.

 

Из формулы

(2.97) сразу же видно, что в чистом

металле

в = О,

oj

истинный процесс термодиффузии полностью

компенсируется

возникающим эффектом Киркендалла

и, зна­

чит, результирующий (наблюдаемый) поток массы отсутствует. Рассмотрим температурную зависимость коэффициента тер­

модиффузии DTA.

Имеем

 

 

 

 

 

 

D A = const• exp(<-

^ j j ^ - )

E l w Q A -

 

(2-99)

Запишем

производную

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

от этого

выражения

dDTA

E,+QA

т

і

 

і

+

 

 

D T -

K T ,

D^-YDA

=

Y D A ~ ^ T ^ ~ 4 -

( 2 Л 0 0 )

 

 

 

 

 

 

 

1 DT

 

Заметим,

что всегда

E1-\-QA>kT.

Поэтому

-~-^>0,

т. е.

 

 

 

 

 

 

 

DA

 

с ростом температуры термодиффузионная подвижность атомов, определяемая величиной DTA, растет, но относительная скорость этого роста

1

dDTA

і

(E1

+ QA

DT

dT ~

Т

\

kT

с повышением температуры убывает.

В заключение рассмотрим вопрос о величине термодиффу­ зионного вклада в эффект Киркендалла. Для этого в соотноше­ нии (2.96) следует оценить величину

(cADTA+cBD

T]J_dT\

Т dx

 

В наиболее типичных условиях опыта температура равна по

порядку

величины

1000° К, на

расстояниях 100 мк она меняется

на 10° К,

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

компонентов — на единицу,

т. е.

 

 

 

 

Т

dxj

дх :

101—2

 

 

 

 

 

 

 

 

Тогда оцениваемая

величина приближенно сводится к виду

 

 

 

 

с

\ E I -

Q A \

•10-

 

 

 

 

 

 

 

kT

 

 

 

 

где С: 0,5 -4-1

(здесь для грубых

оценок по порядку

величины

 

C

A D A

+св£>В

0 , 5 - 1 ) .

 

 

 

положено (D*A -

D*B) gAA

 

 

 

 

 

 

 

 

Разность

E\—QA

обычно

имеет тот же порядок, что и они

сами. Наиболее

типичным значением является 50000

кал/моль,

так что

 

 

 

:25,

С

E I - Q A

• 1 0 — S « 0,1 .

 

 

kT

 

kT

 

 

Другими

словами, термодиффузионный вклад в эффект Кир­

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

Рассмотрим теперь влияние внешнего силового поля на вза­

имную диффузию в бинарной системе.

 

Пусть внешнее

силовое поле, действующее на атом сорта А,

равно FA; примем,

что это силовое поле — потенциальное,

т. е.

 

 

F A = - V W A .

 

Величины

F A и

WA будем рассматривать как функции гидро­

динамических

координат. Это значит, что не учитывается

влия­

ние внешнего поля на условия колебаний атомов, т. е. сохра­

няются

гармоническое приближение, спектр

колебаний и все

параметры,

характеризующие

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

энергии

QA,

QB, E I ) .

применимости

такой модели. Во-

Легко установить критерий

первых, полная энергия частицы равна сумме потенциальной энергии в узле «замороженной решетки», колебательной энергии и потенциальной энергии W частицы во внешнем силовом поле. Чтобы сохранялись все параметры колебаний, величина W по абсолютной величине должна быть много меньше средней коле­ бательной энергии частицы (на одну степень свободы), т. е.

\W\4ZkT.

Во-вторых, как

мы

знаем, колебательное

движение

частицы

в гармоническом приближении описывается

уравнением

 

 

 

 

тг =

— аг.

 

 

 

При наличии

внешнего поля

его нужно

заменить уравнением

 

 

 

т г = — аг •-

'.

 

 

 

Следовательно, для сохранения гармонического приближения

должно выполняться

условие

 

 

 

 

 

 

 

 

WW\«

« Й ;

 

 

 

здесь |г|

усредняется

по периоду

колебаний

и по всем

колеба­

тельным

состоянием. На стр. 83 мы уже получили j г | =

(/гГ/а)

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

 

 

 

 

 

 

\yW\^(akT)4*

 

 

(2.101)

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

 

 

a SFI

 

 

 

F дх < 1 .

(2.102)

Как

правило, все перечисленные выше условия выполняются.

Тогда сохраняется весь формализм метода дырочного газа,

развитый ранее. По-прежнему

 

J

= 2 (Р (0) но/ +

Р (/) и/о} = 4Іі го/ (0) / (/) -

с (/) f (0))

И f =

fA+fB.

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

сомножитель,

Величина / л ( / ) ,

который имеет физический смысл вероятности нахождения атома А в двухузельном состоянии при условии, что первоначально он

был в одноузельном состоянии

в узле

/

и по соседству находи­

лась вакансия (дырка). В отсутствие

внешнего поля эта вероят­

ность равна

 

 

 

 

^

(

ОХ

 

tjU) е х

Р \

kT

 

 

Однако при наличии внешнего

поля

следует учитывать, что

в гармоническом приближении средние положения атома А в одноузельном (в узле /) и двухузельном состояниях не совпа­ дают и в последнем случае это положение смещено на '/г і>.

Следовательно, при переходе атома в двухузельное состояние из одноузельното в узле / необходимо совершить работу против

силового поля, равную

- 1 / 3 r / 0 ( v ^ ) = 1 /a r3 -oF.

Вероятность перехода в двухузельное состояние теперь будет иметь вид

 

 

 

 

 

w—)

=

^ Є

Х Р

 

AT

 

В результате

получим

 

 

 

 

 

 

 

 

J =

-

2 1 Г0 ,-

С (0) 1 / , С 4 (/)

ф л (/) ЄХ Р

_

Q-1 (/) -

1 /"-Г ;0Р Л (/)

 

 

 

 

 

 

 

 

 

 

 

 

 

kT

 

 

 

 

- c(s)V 2 Q (0)cp . 4 (0)

exp

-

QA

(0) - V s

r 0 / F A (0)

 

 

 

 

 

 

 

 

 

 

 

 

AT

 

 

 

+

~

S

%

{ с (0 ) •*,в

(/) ф в (/) ехр [-

QB ( 0 ) - V s r 0 / F B

(0)

 

 

AT

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- c ( / ) V 2 c B ( 0 ) 9 s ( 0 ) e x p

QB ( ( Ц - У а У в (0)

 

 

 

 

 

 

AT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Учитывая условие (2.101), сюда вместо F(/) можно подста­

вить

F(0). Кроме того, поскольку r i 0

= — r0 j,

имеем

 

J -

і

|

% {с (0) Ч,сл <fl

*

И ,х Р [ -

< b « + W ' < ° > ]

_

 

 

'

-

о (;) Ч,сл (0) , Р л

(0) « р

[ -

Я д М - ^ ^ М ] }

+

-(-аналогичное выражение с заменой индекса Л индексом 5.

Примем по-прежнему, что мы имеем дело с кубической, объемно-

центрированной решеткой ( Z = 8 )

и что F направлено вдоль по­

ложительного

направления

оси

х.

Тогда

для

четырех

узлов

г о / = + а / 2

и

для

других

четырех

узлов

rxQj

= —al2,

и мы

получим

 

 

 

 

 

 

 

 

 

 

 

 

т

1 Г

,гл

/ v

- ч

/

Q л

(а) — 1 4

a f , (0) \

-

 

 

J

= ~[ас(0) с

л ( « ) Фл ( а ) е х р ( -

 

 

 

4 1 1 )

 

 

 

 

 

 

 

 

 

 

AT

 

 

 

 

 

- а с ( 0 ) сл(р)флф) ехр

 

0л (Р) - У< «fл (0)

 

 

 

 

 

 

AT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-

„сW

с д

(0)г л (0) ехр ( -

О д М + ^ ' ^ О ^

+

 

 

+ а С ( В с л ( о ) Ф , ( о ) е х р ( - ^ ' 0 ' ; ; ; ^ ( 0 ' ) }

+

 

 

 

 

 

 

 

 

 

+

аналогичное

выражение,

где а и р —символы узлов, расположенных

правее и левее узла

0.

Разлагая

все

величины

типа

с (а)

и

с(р)

в ряд Тейлора

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