Вычисление свойств тройной смеси n2–Ar–o2

Для вычисления свойств тройной смеси предварительно определяем теплофизические свойства чистых компонентов в зависимости от давления p и состава смеси x₁, x₂, x₃, используя в расчетах коэффициенты, приведенные в табл. 2.

Таблица 2

Значения коэффициентов в формулах для вычисления свойств чистых компонентов тройной смеси n2–Ar–o2

Коэффициент	N2	Ar	O2	Коэффициент	N2	Ar	O2
i	1	2	3	i	1	2	3
FAi	0,289	0,208	0,257	FNi	5,46	1,00	14,80
FBi	0,00460	0,00423	0,00628	FTi	256,7	283,3	228,5
FCi	0,00943	0,00350	0,00540	FIi	302,8255	345,1872	364,7802
FDi	0,507	0,413	0,384	FHi	5,92	5,98	6,05
FEi	25,160	34,095	36,360	FPi	24867,26	26767,32	31222,48
FKi	0,211	0,247	0,255	FFi	6811,50	7364,81	8177,97
FLi	0,006725	0,009410	0,008580	FSi	4,030	3,885	4,186
FMi	114	142	131

При этом температуры кипения чистых компонентов (К)

, (13)

а температура жидкой смеси на линии насыщения (К)

, (14)

где

Свойства чистых компонентов:

плотность насыщенного пара (кг/м³)

; (15)

плотность насыщенной жидкости (кг/м³)

; (16)

теплота испарения (Дж/кг)

; (17)

поверхностное натяжение (Н/м)

; (18)

динамическая вязкость насыщенного пара (Пас)

; (19)

динамическая вязкость насыщенной жидкости (Пас)

. (20)

Свойства смеси:

молярная масса смеси в паровой фазе (кг/кмоль)

, (21)

где M_i – молярные массы чистых компонентов; M₁ = 28 кг/кмоль, M₂ = 39,94 кг/кмоль, M₃ = 32 кг/кмоль;

молярная масса смеси в жидкой фазе (кг/кмоль)

; (22)

поверхностное натяжение (Н/м)

; (23)

плотность смеси в паровой фазе (кг/м³)

; (24)

плотность смеси в жидкой фазе (кг/м³)

; (25)

теплота фазового перехода пар–жидкость (Дж/кг)

; (26)

теплота испарения смеси (Дж/кг)

, (27)

где

.

По уравнениям (13) – (27) в подпрограмме ssmes(p,x1,x2) вычисляем параметры смеси, коэффициенты из табл. 2 задаем в подпрограмме smes.