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

3159

.pdf
Скачиваний:
3
Добавлен:
15.11.2022
Размер:
3.17 Mб
Скачать

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

В случае, если коэффициент диффузии D есть функция времени, а не координат, то второе уравнение Фика можно записать как

c

2c

 

 

D(t)

 

.

 

x2

t

 

Это значит, что можно использовать все решения, применявшиеся для постоянного значения коэффициента диффузии, при этом заменив D(t) его средним значением

t

D(t) D(t)dt .

(1.68)

0

 

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

Рассмотрим одномерную задачу, в которой коэффициент диффузии есть функция координаты, т.е. D D(x). Второе уравнение Фика в этом случае имеет вид:

 

 

 

c

 

c

 

D c

 

2c

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

D

 

.

(1.69)

 

 

 

 

 

 

x x

x2

 

 

 

t

x

x

 

 

 

 

Член

D

 

делает

решение

этого

 

уравнения

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

неоднородным, и получить его в замкнутом виде трудно, либо невозможно.

Предположим, что причиной диффузии является градиент концентрации c x. В этом случае зависимость от

41

координаты D D(x) однозначно связана с зависимостью от

концентрации. Рассмотрим решение

для D D(c).

Можно начальное условие

записать через одну

переменную

x

при этом c есть функция только . В

t1/ 2

 

 

 

этом случае уравнение (1.69) перейдёт в обычное однородное дифференциальное. Используя переменный параметр , получим:

 

 

 

 

c

 

dc

 

 

 

 

 

1

 

 

x

 

 

 

dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 t3/ 2 d

 

 

 

 

и

 

t

 

 

d t

 

 

 

 

 

 

 

 

 

 

 

 

c

 

 

 

 

 

 

dc

 

 

 

 

 

1

 

 

 

 

dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

d x

 

 

 

 

 

 

t1/ 2 d

 

 

 

 

Подставляя последние выражения в уравнение (1.69),

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x dc

 

D dc

 

 

1 d

 

 

 

 

dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

3/2

 

 

 

 

 

 

1/2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2t

 

d

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

d

t d

d

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dc

 

 

 

 

 

 

 

 

d

 

dc

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D

 

 

 

 

 

 

 

.

 

 

(1.70)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 d

 

 

d

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Переход от уравнения (1.69) к (1.70) сделал Больцман.

Для

определения

 

 

 

 

 

D D(c)

 

 

 

 

 

 

 

 

Матано

рассмотрел

полубесконечную диффузионную пару (см. раздел 1.4.2.3.), которая описывается следующими начальными условиями:

c c0 при x 0,t 0; c 0 при x 0,t 0;

Так как точка x 0 исключена из рассмотрения при t 0 и исходная концентрация не зависит от расстояния, то,

42

отбросить разрыв при x 0, начальное условие можно записать только через :

c c0 при ; c 0 при ;

Уравнение (1.70) содержит только полные дифференциалы, поэтому можно его сократить на 1 слева и

d

справа и проинтегрировать от c 0 до c c , где c - любая концентрация 0 c c0

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

 

1 c

dc c c

dc

 

 

 

xdc Dt

 

 

 

Dt

 

.

(1.71)

 

 

 

2

0

dx c 0

dx c c

 

Последнее равенство в выражении (1.71) следует из того

в бесконечной системе, dc 0 как при c 0, так и при c c0 . dx

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

c0

 

xdc 0.

(1.72)

0

 

Это условие определяет плоскость, на которой

x 0.

Имеется в виду плоскость, разделяющая вещества A и B (см.

рис.1.1.) при t 0.

При таком условии, зависимость коэффициента диффузии от концентрации можно получить графическим интегрированием и дифференцированием c(x) с помощью

соотношения

 

 

 

 

 

 

 

1 dx

c

 

D(c )

 

 

 

 

xdc.

(1.73)

 

 

 

2t dc c 0

 

Экспериментальные данные по диффузии из одного вещества ( A) в другое (B) можно представить графически,

43

как показано на рисунке на рисунке 1.4. Точка x 0 (штриховая линия на рисунке) определяется в соответствии с уравнением (1.72) и называется поверхностью Матано. Эта линия, которая делит заштрихованную область на две равные по площади половины.

Величину D(c ) для выбранного

c c0 рассчитывают по

измеренной площади

 

участка,

заштрихованного

в

клетку

(интеграл

c

), и обратной величине наклона касательной к

xdc

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

кривой (производной

 

 

 

) в этой точке. Этот метод получил

 

 

 

 

dc c

 

 

 

 

 

 

 

название метода Больцмана – Матано.

 

c c0 1 и

Ошибки

в вычислениях

D(c)

велики

при

c c0 0,

так

как в окрестности

этих точек

площадь

очень

мала, а

производная

 

dc

 

очень

велика. Для уменьшения

 

dx

 

 

 

 

 

 

 

 

 

 

 

ошибки кривую зависимости концентрации от расстояния обычно строят по точкам методом наименьших квадратов – в результате она должна быть похожа на кривую, представленную на рис. 1.4.

Рис. 1.4. Поверхность Матано расположена так, чтобы заштрихованные участки слева и справа от нее были равны. Площадь, заштрихованная в клетку, и касательная – величины,

используемые для расчета D при c 0,2c0

44

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

При расчете коэффициента диффузии в методе Больцмана – Матано по сути определяется коэффициент

~

взаимной диффузии D, который описывает процесс гомогенизации сплава. Это ставит новую задачу определения

~

связи между D и коэффициентами самодиффузии компонентов A и B при том же составе сплава.

~

На первый взгляд, кажется, что D является средним значением коэффициентов диффузии двух компонентов. Однако, если оба компонента диффундируют с различными скоростями, то необходимо ввести дополнительный параметр, определяющий эту разницу. Впервые этот параметр экспериментально определил Киркендалл. Он намотал на прямоугольный пруток латуни (сплав меди с цинком) молибденовую тонкую проволоку и затем нанёс на этот образец слой чистой меди. (рис. 1.5). Молибден практически нерастворим ни в меди, ни в латуни и поэтому никакого участия в диффузии не принимает. Таким образом, был изготовлен образец с неравновесной концентрацией цинка и начальными условиями, соответствующими системе пары полубесконечных тел.

Рис. 1.5. Сечение диффузионной пары в эксперименте Киркендалла

45

Образец был подвергнут серии последовательных отжигов, в результате которых цинк диффундировал из латуни в медь, а медь из внешней зоны диффундировала в латунь. После каждого отжига от образца отрезался кусок и измеряли расстояние между витками молибдена (d ). Было установлено, что d монотонно убывает со временем. Поскольку общий объём образца не менялся, то такая кинетика возможна в случае, если поток цинка превосходит поток меди через одну и ту же плоскость. Сдвиг молибденовых меток в латунь направлен в сторону, противоположную диффузионному потоку цинка.

Каков механизм движение меток? Даркен показывает это на следующем простом примере. Всем известна диффузия чернил в воде. Под диффузионным движением понимается движение молекул или атомов. Представим ту же ситуацию когда чернила вылиты в реку, у которой имеется собственное движение – течение. Тогда для наблюдателя на берегу скорость движения чернил является результатом сложения скорости потока воды и скорости диффузии чернил. Чтобы измерить скорость диффузионного потока необходимо выбрать систему отсчёта, связанную с водой реки и тем самым исключить роль течения. Даркен предлагает связать эту новую систему отсчёта с щепкой, плывущей по течению. В случае нашей задачей роль щепки играют молибденовые метки, которые дрейфуют (текут) вместе с кристаллической решёткой - аналогом реки. Их положение определяет истинную поверхность раздела системы медь - латунь

Диффузионные потоки компонентов A и B в системе отсчёта, связанной с метками, не равны (JA JB ). Но в неподвижной системе отсчёта, связанной, например, с внешней поверхностью образца (берегом реки) потоки компонентов A и B ( JA и JB ) выравниваются за счёт течения решётки (Vk ) .

46

JA

JA VkcA JB (JB VkcB),

(1.74)

где cA

и

cB

-

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

компонентов,

рассчитываемые как число атомов на единицу объёма.

JA

 

В такой

 

постановке

 

задачи

потокам

и

JB соответствуют собственные коэффициенты диффузии DA

и

 

 

 

 

 

 

 

 

 

 

 

~

 

DB , а потокам JA

и JB коэффициент взаимной диффузии D.

 

Воспользовавшись первым уравнением Фика запись

(1.74) можно представить как

 

 

 

 

 

 

 

 

 

 

D

cA

V c

 

D

cB

V c

 

(1.75)

 

 

 

 

x

 

 

 

A

x

k

A

B

k

B

 

 

При постановке этой задачи было оговорено, что объём системы в процессе отжига не меняется. Следовательно, общее

число атомов в единице объёма

c cA cB

есть

величина

постоянная,

а cA

x cB x.

Заменив

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

молярными

долями

компонентов:

NA cA c

и

NB cB c,

получим из уравнения (1.75) выражение для скорости течения решётки.

 

V (D

 

D )

NA

 

(1.76)

 

 

x

 

k

A

 

B

 

 

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

скорость

течения

решётки

и

соответственно дрейфа молибденовых меток определяется разницей собственных коэффициентов диффузии.

~

Для установления связи между D и собственными коэффициентами диффузии подставим уравнение (1.75) в уравнение (1.74)

~

J

A

V c

A

 

J

B

V c

B

 

D

 

 

k

 

 

 

k

(1.77)

cA

x

 

cB

x

 

 

 

 

 

 

Учитывая вывод (1.76) выражение для коэффициента взаимной диффузии приобретает вид

47

~

DBNA

(1.78)

D DANB

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

~

заданного состава путем измерения Vk и D. В разбавленном

~

растворе NB 0 и D DB , т.е. коэффициент взаимной диффузии определяется скоростью диффузии атомов примеси. Если NA NB , но DA DB коэффициент взаимной диффузии

~

D DANB , т.е. скорость перемешивания определяется диффузией более быстрого компонента.

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

 

x

t

(1.79)

Тогда скорость дрейфа меток

 

V

 

dx

 

x

(1.80)

 

2t

k

 

dt

 

 

При рассмотрении метода Больцмана – Матано (раздел 1.7) было введено понятие поверхности Матано (уравнение 1.72), как границы раздела диффузионной пары, изменяющей своё пространственное положение по мере увеличения времени отжига. Выражениями (1.79) и (1.80) по сути и определяется её положение и скорость перемещения, т.е. метки всегда находятся на поверхности Матано. Основным свойством этой поверхности является постоянство концентрации диффундирующего вещества на неё для всех t 0. В эксперименте Смигельскаса и Киркендалла

48

концентрация цинка на этой поверхности соответствовала содержанию 22,5% Zn. Они нашли, что за 56 дней при температуре равной 785оС метки сдвинулись на 0,0125 см. Для этого состава Даркен получил следующие коэффициенты собственной диффузии:

DCu 2,2 10 9 см2/сек;

DZn 5,1 10 9 см2/сек;

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

~

Nb V ) и опытные значения D становятся такими, что для их объяснения одному из собственных коэффициентов диффузии приписывается отрицательное значение, что не имеет физического смысла.

Бардин и Херринг объясняют сдвиг молибденовых меток возникновением в зоне взаимной диффузии вакансионного

потока (Jv ), компенсирующего разность потоков компонент:

Jv (JA JB )

Этот поток направлен в строну более быстрого компонента. В нашем случае в сторону потока Zn, т.е. в латунь.

Со стороны меди возникает избыток вакансий по сравнению с их термодинамически равновесной концентрацией, что приводит к образованию вакансионных пор. Лишние вакансии должны исчезнуть в диффузионной зоне, например, поглотиться дислокациями за счёт их переползания. При таком движении дислокаций число узлов кристаллической решётки в единице объёма сохраняется, однако, переползание приводит к исчезновению атомных экстраплоскостей. Если до начала диффузии (t 0) метки

49

находились в плоскости раздела диффузионной пары, созданной экспериментатором (т.н. плоскость Киркендалла), то в процессе отжига они сместятся в сторону цинка, оставаясь теперь всё время на поверхности Матано (x 0 ).

Таким образом, теория Бардина – Херринга в место течения решётки создаёт механизм вакансионного дрейфа меток. Их сдвиг есть следствие стремления системы установить равновесную концентрацию вакансий.

Проведенные ими расчёты показали, что скорость смещения плоскости Киркендалла (VK )можно выразить как

V

(D

D )

NA

,

(1.81)

 

K

A

B

x

 

что совпадает с уравнением (1.76) для скорости течения решётки. Поэтому сдвиг меток за счёт течения решётки или потока вакансий физически идентичны.

50

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]