Расчет коэффициентов активности компонентов бинарного раствора с помощью уравнения Гиббса – Дюгема

– Задача. С помощью уравнения Гиббса – Дюгема рассчитать коэффициент активности железа в расплавах Fе–Аl при 1873 К по экспериментальным значениям коэффициента активности алюминия.

– Исходные данные. Зависимость коэффициента активности алюминия (g_Al) от концентрации в расплавах Fе–Аl при 1873 К [6] приведена в таблице.

Результаты расчета коэффициентов активности железа
в расплавах Fе–А1 при 1873 К с помощью уравнения Гиббса – Дюгема

x Al	x Fe		-ln gAl
			2,38		2,38
0,1	0,9	0,111	2,20	0,01	2,72	0,005
0,2	0,8	0,250	1,95	0,05	3,05	0,047
0,3	0,7	0,428	1,60	0,18	3,26	0,165
0,4	0,6	0,667	1,25	0,37	3,47	0,374
0,5	0,5	1,00	0,82	0,64	3,30	0,695
0,6	0,4	1,50	0,50	1,02	3,12	1,01
0,7	0,3	2,33	0,24	1,44	2,67	1,58
0,8	0,2	4,00	0,10	1,90	2,50	2,00
0,9	0,1	9,00	0,02	(2,33)	2,0	2,48
1,0		¥		(2,80)		2,84
________________ * Уравнение (1.31). ** Уравнение (1.32).

– Теория. Парциальные термодинамические характеристики компонентов бинарного раствора связаны уравнением Гиббса – Дюгема. Для парциальных энергий Гиббса имеем: . Учитывая, что и dx ₁ = - dx ₂, это уравнение можно записать в виде

(1 – х ₂)(¶ ln a ₁/¶ х ₂) _{T , р} = - x ₂(¶ ln a ₂/¶ х ₂) _{T , р} .

– Отсюда

(1 – х ₂)(¶ lng₁/¶ х ₂) _{T , р} = - x ₂(¶ lng₂/¶ х ₂) _{T , р} . (1.30)

– Таким образом, если известен коэффициент активности одного из компонентов, например g₂, то с помощью уравнения (1.30) можно вычислить g₁. Проинтегрировав уравнение (1.30), получим

– или

. (1.31)

– Решение. Определение коэффициента активности железа в расплавах Fe–Аl по заданной концентрационной зависимости коэффициента активности алюминия сводится к вычислению интеграла (1.31). Воспользуемся графическим методом интегрирования. С этой целью строим график зависимости отношения x _Аl/ x _Fe от lng_Al (см. рисунок). Заштрихованная на рисунке площадь равна интегралу

x Al/ x Fe

и значению lng_Fe при x _Al = x _Fe = 0,5. Рассчитанная таким образом концентрационная зависимость коэффициента активности железа приведена в таблице вмсте с исходными данными.

– Примечания. 1. В области низких концентраций железа графическое интегрирование уравнения (1.31) затруднено по двум причинам. Во-первых, в этой области невысока точность определения величины g_Al, во-вторых, отношение x _Al/ x _Fe ® ¥ при x _Fe ® 0, поэтому при низких значениях х _Fe (x _Fe < 0,2; х _Al > 0,8) коэффициенты активности железа получены экстраполяцией зависимости lng_Fe к точке x _Fe = 0.

– 2. Более удобное для расчета интегральное уравнение Гиббса – Дюгема может быть получено с помощью вспомогательной функции Даркена [7]:

a₂ = lng₂/(1 - х ₂)².

– Тогда

. (1.32)

– В данном случае подынтегральное выражение ограничено во всем концентрационном интервале 0 £ x ₂ £ 1. Функция a_Al и результаты расчета концентрационной зависимости lng_Fe по уравнению (1.32) приведены в таблице.

– 3. Более точно интегралы в уравнениях (1.31) и (1.32) можно вычислить численными методами (трапеций, Симпсона и др.) с помощью ЭВМ.

– 4. Для ряда систем коэффициент активности g₂ известен лишь в ограниченной области концентраций от до . Например, для системы Fе–С при 1873 К величина g_Fe известна в интервале 0,79 £ x _Fe £ 1. В этом случае величину g_C можно рассчитать по уравнению

.

– Для расчета концентрационной зависимости g_C по этому уравнению необходимо знать величину при x _Fe = 0,79 ().