ВУЗ:
Составители:
где x=T*10
-4
.
Дифференцированием (42) по температуре легко
получить выражения для удельной энтальпии, энтропии и
теплоёмкости вещества:
S=a
1
+a
2
(1+lnx)-a
3
/x
2
+4a
7
x
3
+3a
6
x
2
+2a
5
x
H=10
4
(a
2
x-2a
3
/x-a
4
+a
5
x
2
+2a
6
x
3
+3a
7
x
4
)
C
p
=a
2
+2a
3
/x
2
+2a
5
x+6a
6
x
2
+12a
7
x
3
Полиномы (42) используются в большинстве программ физико
-
химического назначения, написанными в СССР и России.
В [5] используется полином (полином JANAF) для аппроксимации
теплоёмкости:
Cp
R
=a
1
+a
2
T+a
3
T
2
+a
4
T
3
+a
5
T
4
(43)
откуда несложно получи
т
ь выражения для других термодинамических
функций:
H
RT
=a
1
+
a
2
2
T+
a
3
3
T
2
+
a
4
4
T
3
+
a
5
5
T
4
+
a
6
T
S
R
=a
1
lnT+a
2
T+
a
3
2
T
2
+
a
4
3
T
3
+
a
5
4
T
4
+a
7
Φ
R
=a
1
(lnT-1)+
a
2
2
T+
a
3
6
T
2
+
a
4
12
T
3
+
a
5
20
T
4
-
a
6
T
+a
7
Полиномы JANAF используются в большинстве программ,
написанных в США. С конца 90-
х
годов формат таблиц JANAF был
изменён и сделан аналогичным (42), тем не менее, большинство
западных программ по-прежнему ориентированы на формат (43).
В обеих случаях приведённые выражения дают относительные
значения H, S, Φ. Для получения полного значения энтальпии к
этому значению следует прибавить энтальпию образования H
f
,
которая приводится отдельно.
Д
ля полиномов [3] вычисление полных
функций (H, S и Φ
)
будет производиться автоматически, если
изменить значение коэффициента a
4
: a
*
4
=a
4
-H
f
*10
-4
.
10.Схема программы расчёта равновесия.
Как следует из вышеизложенного, программа расчёта равнове
с
ия
гетерогенной системы оказывается достаточно сложной даже в
приближении идеальный газ - идеальные растворы - несжимаемые
конденсированные фазы. Она должна содержать несколько
вспомогательных процедур, а в ходе основного цикла (ньютоновских
итераций) выполнять ряд проверок. В настоящее время существует
ряд отечественных и зарубежных программных пакетов,
предназначенных для расчёта термодинамического равновесия.
Перечень наиболее распространённых программ такого рода можно
найти в [2]. Рассмотрим вариант программы:
где x=T*10-4. Дифференцированием (42) по температуре легко
получить выражения для удельной энтальпии, энтропии и
теплоёмкости вещества:
S=a1+a2(1+lnx)-a3/x2+4a7x3+3a6x2+2a5x
H=104(a2x-2a3/x-a4+a5x2+2a6x3+3a7x4)
Cp=a2+2a3/x2+2a5x+6a6x2+12a7x3
Полиномы (42) используются в большинстве программ физико-
химического назначения, написанными в СССР и России.
В [5] используется полином (полином JANAF) для аппроксимации
теплоёмкости:
Cp 2 3 4
R =a1+a2T+a3T +a4T +a5T (43)
откуда несложно получить выражения для других термодинамических
функций:
H a2 a3 2 a4 3 a5 4 a6
RT =a 1+ 2 T+ 3 T + 4 T + 5 T + T
S a3 2 a4 3 a5 4
R=a1lnT+a2T+ 2 T + 3 T + 4 T +a7
Φ a2 a3 2 a4 3 a5 4 a6
R =a1(lnT-1)+ 2 T+ 6 T +12T +20T - T +a7
Полиномы JANAF используются в большинстве программ,
написанных в США. С конца 90-х годов формат таблиц JANAF был
изменён и сделан аналогичным (42), тем не менее, большинство
западных программ по-прежнему ориентированы на формат (43).
В обеих случаях приведённые выражения дают относительные
значения H, S, Φ. Для получения полного значения энтальпии к
этому значению следует прибавить энтальпию образования Hf,
которая приводится отдельно. Для полиномов [3] вычисление полных
функций (H, S и Φ) будет производиться автоматически, если
*
изменить значение коэффициента a4: a4=a4-Hf*10-4.
10.Схема программы расчёта равновесия.
Как следует из вышеизложенного, программа расчёта равновесия
гетерогенной системы оказывается достаточно сложной даже в
приближении идеальный газ - идеальные растворы - несжимаемые
конденсированные фазы. Она должна содержать несколько
вспомогательных процедур, а в ходе основного цикла (ньютоновских
итераций) выполнять ряд проверок. В настоящее время существует
ряд отечественных и зарубежных программных пакетов,
предназначенных для расчёта термодинамического равновесия.
Перечень наиболее распространённых программ такого рода можно
найти в [2]. Рассмотрим вариант программы:
