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

книги из ГПНТБ / Бокштейн Б.С. Термодинамика и кинетика диффузии в твердых телах

.pdf
Скачиваний:
19
Добавлен:
27.10.2023
Размер:
18.93 Mб
Скачать

Заметим

что

общее число

узлов

постоянно,

следовательно,

CA + CV = const

и

V С„ =

- V

СА. Тогда вместо второго

уравне-

'

'•'V

------

 

 

~

 

/

-W C o ,

откуда

(с уче-

ния (316)

получим

 

/0=

kTLvv у Со

 

Са У

 

 

том

С0 С

СА),

 

 

 

 

 

 

 

 

(318)

 

■kT

Lvv

 

 

 

 

 

 

 

 

 

Cv

 

 

 

 

 

 

 

 

 

Поскольку

Laa =

Lvo,

находим связь между. коэффициентом

самодиффузии

(DAA)

и

коэффициентом

диффузии

вакансии

(Uvv).

Г)

_ Cv В

vv-

 

 

 

 

 

 

 

 

(^^)

и АА qJ

u

 

 

 

 

 

 

 

 

 

Коэффициент Dm можно найти экспериментально, например из опытов по кинетике изменения электросопротивления избыточных вакансий после закалки (гл. III), а затем рассчитать DAA. Следует только иметь в виду, что рассчитанное значение будет отличаться от истинного, поскольку в условиях пересыщения решетки вакансиями

Cv =f= Cpv. Однако результат (319) справедлив и при Си <=>Ср, когда

CJCA я» Npv {Nl — атомная доля равновесных вакансий; СА я» С — общему числу узлов), так что действительно коэффициент самодиффу­

зии

D

= NpD

(320)

U АА

V W

 

Диффузия в бинарных системах

Система уравнений (304) для потоков атомов растворителя (Л) и

растворенного вещества (В) имеет вид:

 

 

iA =

— LAA\ / (цл —

^ а в Д (Н-в — Вв). 1

(321)

/ в =

L A b V

( И л — щ )

^ b b V ( Ц в

В в ) -

1

 

Все потоки определены по отношению к средней скорости переме­

щения частиц:

 

 

 

 

 

 

CdVi

V

 

 

 

 

 

2j Ck^k

 

 

 

(322)

 

к

 

 

 

 

V ■

 

 

 

 

 

 

 

где

 

 

Cv — число вакантных

узлов в единице объема;

 

 

 

Ck — число узлов,

занятых -атомами

&-того сорта;

 

 

vv и

vk — соответствующие

скорости в

неподвижной

 

С = С0 + ^

 

(лабораторной) системе отсчета;

 

Ck — общее число узлов в единице объема, прини-

 

 

k

 

маемое постоянным.

 

 

 

 

 

 

 

Для линейной задачи вектора скоростей совпадают или противо­

положны.

 

 

 

 

 

 

 

Уравнения (321) переходят в первый закон Фика при выполнении

следующих условий:

а) потоки атомов обоих сортов

не взаимодей-

160

ствуют,

L ab

= 0; б) вакансий находятся в равновесии,

Cv — С„,

так

чт0

п

== 0

(достаточно

менее

жесткого требования V

=

0 ■

концентрация вакансий

постоянна); в) раствор —- идеальный,

так

что для обоих компонентов рг =

р; +

kT In Сг и V

= kT V In Сг.

Действительно, если все эти условия выполняются, то оба уравне­

ния

(321) можно записать в виде /г =

Dt V Сг, где

 

 

 

D. =

feT^==feT«,..

 

 

 

 

 

(323)

 

W

 

 

 

 

 

 

 

Мы нашли, таким образом, условия применимости первого закона

Фика. „ Откажемся от первого ограничения. Пусть потоки атомов взаимо­

действуют: Lab =h 0. Остальные условия сохраним. Проанализируем выражения для потоков и jB, полагая, что градиенты химических потенциалов атомов связаны уравнениями Гиббса Дюгема:

(324)

СдУЦл -Ь CBV pa = 0.

Из выражений (321) и (323) следует, что

/-.— « ■ ( ' б т - т г ) ч С л ’ ]

(325)

 

Переход от выражения (321) к (325), как и в случае самодиффузии, эквивалентен учету корреляции. В уравнениях (325) учтено, что ско­ рость обмена атомов сорта А с вакансиями зависит также от скорости обмена атомов сорта В с вакансиями, поскольку в соответствии с вы-

ражением (303) Ьдв = Asa

А40

или ^ ав ~~ Ааа ц ву»

а кинетические коэффициенты LAb и LBv

отражают взаимодействие

потоков атомов обоих сортов и вакансий.

Отказ от третьего ограничения эквивалентен введению термо­

динамической активности. В реальном растворе рг-

= рГ +

kT In щ.

Подстановка этого выражения в уравнения

для

потоков

приведет

в конечном счете к результату:

 

 

 

 

 

 

 

(326)

где у- =

а-1С■—- коэффициент активности.

предположение о том,

Как

было подробно показано в гл. III,

что вакансии находятся в равновесии, во многих случаях оказывается неправильным, причем неравновесная концентрация (избыток или недостаток) вакансий может существовать достаточно долго.

Авторами монографии [2 ] была введена характеристика мощности источников или стоков (ав) как разница между числом вакансии, соз­ даваемых и исчезающих в единице объема в единицу времени. По су-1

11 Заказ № 737

.

ществу о0 — это скорость образования или исчезновения вакансий; соответственно можно написать

С — Ср„

(327)

Т у

где Cv — число вакансий в единице объема;

Ср — равновесное значение концентрации вакансий; т0 — среднее время жизни неравновесных вакансий.

Авторы работы 12] определили а0 следующим образом:

(У,

 

(328)

где

р0 — химический потенциал вакансий, равный kT In (Cv/Cp).

При

С

__qp

сравнительно небольших пересыщениях In (C jC p)r - v

v

так

что

Cpv

 

ov — Ak cv- c p

(329)

 

cp

 

т. e. cr0 с точностью до множителя Ak совпадает с введенной в гл. III величиной пересыщения (недосыщения) решетки вакансиями. Сравне­

ние выражений (328) и (329) показывает, что А = —Cp/krv, т. е. коэффициент А зависит от плотности источников и стоков, подвиж­ ности, энергии образования вакансий и от температуры. Если в кри­ сталле устанавливается равновесная концентрация вакансий, то их химический потенциал (ра) и мощность источников (ап) обращаются в нуль.

Для атомов выполняется закон сохранения, так что

 

— div (Cfttift) = — d iv /° .

(330)

Индекс «0» при jk означает, что речь идет о потоке атомов й-того сорта в неподвижной системе отсчета. Однако вакансии могут не подчи­ няться закону сохранения. Для них

~ = — div / “+ div (Си),

(331)

где v — средняя скорость перемещения атомов (скорость смещения решетки) в неподвижной системе отсчета, определяемая по формуле (322), а С = Cv + £ Ck.

k

Очевидно, член div (Cv) в выражении (331) описывает возникнове­ ние или исчезновение вакансий. Авторы монографии [2] показали, что

сг„ = div (Си).

(332)

162

В диффузионной зоне, т. е. в области, где возникают большие

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

Эффект Киркендалла

В 1947 году Киркендалл после пяти лет экспериментальной ра­ боты надежно показал, что метки из молибденовой проволоки, поме­ щенные на границе раздела между медью и латунью, содержащей около 30% Zn, сдвигаются в сторону латуни. Величина сдвига ока­ залась пропорциональной корню квадратному из времени и доста­ точно большой, так что ее нельзя было объяснить изменением объема 1

Киркендалл пришел к выводу, что поток атомов цинка из латуни превосходит встречный поток атомов меди в латунь, — именно это приводит к сдвигу меток. Таким образом, компоненты бинарного сплава диффундируют с разной скоростью; иначе говоря, их пар­ циальные коэффициенты диффузии различны. Как справедливо отме­ чено в работе [28], в 1947 г. эта концепция была совершенно новой и не очевидной.

Эффект Киркендалла наблюдали впоследствии в разных системах и с разными метками. Он имеет весьма общий характер. Некоторые результаты будут приведены ниже.

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

зии компонентов, если известен коэффициент взаимной диффузии (D) и скорость сдвига меток (v), совпадающая со скоростью течения ре­ шетки, как целого.

Согласно Даркену, если объем образца остается постоянным, то суммарные потоки обоих компонентов 2* (с учетом течения) равны

друг другу, т. е. выполняется условие

 

 

- D

^ + v C ^ - D ^ + vC*

(333)

Обозначим общее число мест в единице объема через С : С ^

Сг +

С2

(при точной записи следует учитывать концентрацию вакансий

Св),

тогда C-JC = Nt — молярной доле компонента и

 

 

v =

(D1- D 2) d^ .

(334)

Коэффициент взаимной диффузии можно определить в соответ­ ствии с первым законом Фика как коэффициент пропорциональности

t При добавлении меди к латуни меняется плотность.

2 Отметим, что суммарный поток определен по отношению к неподвижной (лабо­ раторной) системе координат. Скорость течения решетки в этой системе (о) совпа­

дает со средней атомной скоростью (322).

 

11*

163

между суммарным потоком любого компонента и соответствующим градиентом концентрации, т. е.

D —

i ? +oC'

- D* i i + va

(335)

dCj

dC2

 

dx

dx

 

Подставляя в выражение (335) значение v из (334) и полагая число

узлов (мест) в единице объема постоянным, вследствие чего

=

dC,

ПОЛуЧИМ

 

 

 

= ----:

 

 

 

D = D 1N 2 ~\~ D%Nх.

 

 

(336)

Уравнения (334)

и (336)

позволяют рассчитать D 1 u D 2 n o D n v

в плоскости с концентрациями Сг и Сг.

 

Табл. 27 составлена на основании данных измерений, относящихся

к металлам с г. ц. к.

и о. ц.

к. решетками. В последних эффект Кир-

кендалла также наблюдается; раньше считали, что это не так. Ско­ рость сдвига меток рассчитывали по величине сдвига (|) и времени (():

если |

\'1, то

v = |/2 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а 27

Парциальные коэффициенты диффузии

в некоторых системах

 

 

Система

t. °C

Ni

T, r

 

Парциальные коэффи­

d j d 2

CM

 

циенты, см2/с

Ag—Au

940

0,635

100

0,029

7,3 • 10" 9;

1,7 -10-9

4,3

Znp—Ag

400

0,492

18

0,0037

1,2• 10~7;

3.10“ 8

4,0

Mge—А1

425

0,51

194

0,00835

3

-10~s;

1,7 ■10"9

1,8

Au—Pt

1020

0,967

100

0,0178

5,6-10~9; 7,3• 10-11

77

BeNi

1200

0,99

8

0,0221

7,8 • 10“ 7; 5 • 10"9

157

Cu—Ni

1054

0,899

313

0,0140

3,8 • 10"10; 2,9-10-10

1,3

Zna —Cu

785

0,225

24

0,0015

5,1 • 10"9; 2,2 • 10" 9

2,3

Zna —Cu

890

0,28

36

0,012

6,3

■10-8;

1,2 • 10" 8

5,3

Zrio—Cu

395

0,84

213

0,12

4,7

• 10-8;

1,0

-10-3

47

U -T i

950

0,82

36

0,0118

4,7

-10“ 9;

1,2

• 10“ 9

3,9

U—Mo

1050

0,90

100

0,0126

3,4

-10-8;

1,3

-10“ ®

2,6

U—Nb

800

0,995

96

0,0166

4

• 10"10; 2 • 10"11

20

U—Zr

950

0,875

36

0,0159

7,7

• 10-a; 6,5 10~10

11,8

U—Zr

1075

0,95

36

0,0300

3,7

10“ 8;

2,8

10-

12,8

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

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

164

непосредственно не фигурируют. Обсуждение связй между эффектом Киркендалла и процессами возникновения — исчезновения вакан­

сий было дано в работе [3].

На рис. 47 изображено положение инертных меток в плоскости х = 0 до начала диффузионного отжига (а) и после некоторого вре­ мени t, когда метки сдвинулись и х фО (б). Назовем плоскость, в кото­ рой расположены метки, — плоскостью Киркендалла (К)\ она сме­ щается в лабораторной системе координат.

Рис. 47. Положение меток в плоскости х — О

Рис. ^48.

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

(вертикальная полоска) при i = 0 (а) :'и сдвиг

(jv

= jA

jB)

в эффекте Киркендалла

меток в результате эффекта Киркендалла (б)

 

 

 

 

Потоки через плоскость К можно

записать

без учета течения:

jA = — Da WCa и jB = —DB4 C B. ПосколькуDA ф DB, то и jA ф jB,

если еще учесть, что V СА = —V Св .

Если общее число узлов постоянно в любом элементе объема (это предположение, которое было сделано и в выводе Даркена, исклю­ чает, например, возможность порообразования), то в соответствии с теорией Онзагера в бинарной системе с вакансиями поток последних

h = - U л + /в) = Ф а - D b ) V C a .

(337)

Если Da ~>Db, то атомов А уходит слева направо больше, чем справа приходит атомов В, и число вакансий слева от первоначаль­ ной плоскости раздела (в фазе А) возрастает. Это эквивалентно по­ току вакансий, направленному справа налево (рис. 48). В результате в фазе А возникает избыток вакансий по сравнению с термодинами­ чески равновесной их концентрацией. Лишние вакансии должны ис­ чезнуть. Это возможно, например, в результате переползания дисло­ каций, так как при этом сохраняется число узлов в единице объема. Но это приводит к исчезновению атомных плоскостей и соответственно плоскость К приближается к левому краю образца, изображенного на рис. 47. Справа от первоначальной границы раздела (в фазе В) атомные плоскости появляются, чтобы восполнить недостаток вакан­ сий. Результат — тот же: плоскость К удаляется от правого края кристалла.

165

В расчете на единицу площади плоскость К смещается в единицу времени на объем /0Й (й — атомный объем; как и у Даркена, прини­ мается, что все частицы имеют одинаковый атомный объем), т. е. скорость ее перемещения

Vk = dM = j 0Q = (DA- D B)Q v C A

 

или

 

 

vk = (DA- D B) 4 N A,

(338)

где

— атомная доля компонента А в сплаве.

 

 

Полученный результат совпадает с уравнением (334) для скорости

сдвига меток в теории Даркена. Таким образом, оба подхода (сдвиг меток определяется течением решетки или потоком вакансий) физи­

чески идентичны.

Бардин и Херринг [3 ] обсудили также механизм возникновения или исчезновения точечных дефектов, который приводит к замещению одних плоскостей решетки другими и смещению границы раздела. Поскольку это смещение одинаково вблизи поверхности, около гра­ ниц зерен и внутри образца, они сделали предположение, что решаю­ щую роль в нем играют дислокационные источники и стоки вакан­ сий. Предложенный ими механизм возникновения вакансий анало­ гичен механизму Франка—Рида и был рассмотрен в гл. III. Как видно из рис. 17, дислокация, закрепленная в точках А и В, восходит при поступлении к ней вакансий (плоскость скольжения перпендику­ лярна плоскости рисунка). Если есть избыточные вакансии, они анни­ гилируют по мере восхождения дислокации, которая все сильнее изгибается (положение С2, С8 . . .), пока (в положении С5) образуется петля, ограничивающая область, в которой исчезает атомная пло­ скость. Оставшийся сегмент приводит к повторению процесса. Так работают источники Бардина—Херринга.

Отклонения от расчетных формул чаще всего связаны с образова­ нием пор, которые также являются внутренними стоками вакансий. Соответственно поры образуются там, где скорость исчезновения ва­ кансий (o0) максимальна — эта область всегда находится вблизи пло­ скости К со стороны быстрого компонента.

Процесс образования пор в диффузионной зоне обсуждался в гл. III. Этот процесс был весьма детально изучен в работах Гегузина [127], который разделил эффекты Френкеля (порообразование) и Киркендалла. Последний, таким образом, сводится к поглощению избыточных вакансий дислокациями, что приводит к диффузионному восхождению дислокаций и сопровождается смещением исходной границы раздела между компонентами диффузионной пары.

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

166

приложением внешнего давления всестороннего сжатия 1 можно по­ давить эффект Френкеля.

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

(и-о ** 0).

Как было показано выше, пренебрежение взаимодействием потоков означает фактически неучет фактора корреляции. Маннинг [115] показал, что введение этого фактора, который он назвал «вакансионным ветром», увеличивает поток быстрого компонента (Л) и его пар­ циальный коэффициент диффузии (DA). Соответственно предсказы­ ваемое смещение плоскости К увеличится примерно на 28% в г. ц. к. (/~4 = 1,28) и еще в большей степени в других (в кристалле с решет­ кой алмаза — вдвое). Влияние этого фактора на коэффициент взаим­ ной диффузии меньше: если DA и DB различаются в три раза или меньше, то поправка при диффузии в г. ц. к. решетке меньше 7% (в кристалле с решеткой алмаза меньше 25%). Величина поправки возрастает с приближением NA к NB и с увеличением DAIDB. Макси­ мальное ее значение равно /~ 4, где f — фактор корреляции.

Химический потенциал вакансий не равен точно нулю, иначе равнялась бы нулю скорость их образования (или исчезновения).

Согласно выражению (328), сг0= —А величина может быть

велика даже при очень малых р0, так как А — большая величина.

Действительно, А = Cpvlkxv^ 1024/10~16-10~2 = 1042 ~-а~ г |р '~с~Р~

(ПРИ

Т ^ Т пл).

что пересыщение вакансиями составляет 1%,

тогда

Примем,

\ i j T = k\n (C jC l) ^ 1СГ~16.10—2=

10“ 18 эрг7град

и, следовательно,

| а01«=* 1024

вакансий/м3-с.

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

вакансий физиче­

Предположение о равновесной

ски означает, что все 1024 вакансий каждую секунду успевают исче­ зать на дислокационных стоках. Формально же это значит, что в урав­

нении (331)

—div jv -)- div(Cu) член dCJdt в левой

части мал

и в диффузионной зоне выполняется соотношение

 

Oj, = div (Си) =

div / S-

(339)

Однако если по каким-либо причинам компенсация (т. е. сток избы­ точных вакансий) не обеспечивается, то в расчеты следует вводить поправки.

Если раствор, образующийся при взаимной диффузии, не является идеальным, то парциальные коэффициенты (D£) в уравнениях Даркена содержат термодинамический множитель Dt — utkT (1 -f- + д In yjd In C£), так что если смешивание компонентов приводит

' Н г у е н Ч о н г Б а о . Исследование влияния малых давлений всестороннего сжатия на процессы в диффузионной зоне. Автореф. канд. дис. Харьков, 1969.

* 1 эрг = 10"7 Дж.

167

к уменьшению энтальпии системы (теплота смешения отрицательна),^ то усиливается тенденция ко взаимной диффузии и парциальные коэф­

фициенты растут.

Таким образом,

для подробного анализа процесса взаимной диф­

фузии необходимо, кроме D и v, знать еще зависимость коэффициента активности от концентрации. Анализ [133] показал, что результаты расчета по уравнениям Даркена достаточно хорошо согласуются

сопытом во всяком случае для твердых растворов с г. ц. к. решеткой

ив пределах погрешности существующих методов определения D для металлических систем (20—30%). Согласие хуже для твердых

растворов с о. ц. к. решеткой.

Особый интерес представляет анализ ситуации, когда скорость

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

2,

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

компонентов отличаются не слишком сильно {N1

 

N 2). Именно это

наблюдалось в работах [134, 7; с. 56]: в первой работе

— для взаим­

ной диффузии в [S-латуни (г. ц. к. решетка), а во второй

для систем

типа Nb—V, Mo—Ti и т. д. (о. ц. к. решетка). Совместное решение

уравнений Даркена для D и v привело в обоих случаях к отрицатель­ ным значениям одного из парциальных коэффициентов (D2).

Действительно, рассмотрим случай, когда v/W N l ^>DIN2 (именно так было в указанных работах). Из формулы (336) следует, что

Учтем теперь выражение (334) и получим

 

D2= D -

v N 2

(340)

VЛД

 

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

парциальных коэффициентов.

этой монографии (совместно с

В этой связи одним из авторов

с Л. С. Швиндлерманом) ситуация,

складывающаяся при эффекте

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

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

168

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

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

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

потенциалов

компонентов (V цг),

однако в них

учитывают только

энтропийный

член

(feTVlnC;)

с

поправкой

на

неидеальность

(1 + д In y jd

In С;), что и приводит к уравнениям Даркена. В общем

случае V ц,-

содержит также член,

пропорциональный градиенту

давления — vt V р,

где vL— парциальный объем (если парциальные

объемы всех

компонентов одинаковы, то vt =

Q).

Этот член, как

правило, не учитывают, считая р = const. При наличии осмотиче­

ского давления это неверно.

в выражение для

Рассмотрим процесс взаимной диффузии, введя

движущей силы (302) член, учитывающий градиент давления:

k T ^ X n C i - v ^ p .

(341)

Напишем выражения для потоков в неподвижной системе коорди­ нат, как и Даркен, пренебрегая перекрестными членами (L12 = = L21 = 0) и полагая, что V Сг = —V С2:

ii =

- L

п ( ^ v C x +

QVp),

 

 

 

 

 

1

 

-

 

 

й =

^22

( — £ “ V C

2 4 - n v p j .

 

 

Из условия постоянства объема j\

/2

= 0 следует, что

V p :

UT

t-22/^2

Д г Д 1 VCi.

 

 

 

 

 

р (^-11+ t,22)

 

 

Так как согласно выражению (323)

Ln

= D fii/kT, то

Vp —

kT

P 2

D\

X7 c

 

 

Q D A + B A V

 

 

(342)

(343)

(344)

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

169

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