же состава он может меняться в зависимости от времени, если меняется температура. Коэффициент диффузии может также зависеть от состава, а, следовательно, при наличии градиента концентрации, меняться от точки к точке по объёму образца.
В случае, если коэффициент диффузии 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