ГОСТ 34898—2022
Окончание таблицы 4
р
о
м
-•i41
eo4"
м
01
О
LL
О
СО
II
-ti41
Коэффициент сжимаемости Z(f2, р0)
У
Компонент
t2 =
20 °С
54
Кислород
0,99325
0,99434
0,99438
0,99466
55
Диоксид углерода
0,97509
0,98024
0,98040
0,98161
Для достижения согласованности этот пакет реализует, в частности, уравнение GERG—2004
(см. [60]), с помощью которого вычислены значения для метана, этана, пропана, н-бутана, 2-метил-
пропана (изобутана), водорода, гелия, аргона, азота, кислорода и диоксида углерода. Другие источни
ки, используемые RefProp®: для этена (см. [44]), пропена (см. [45]), сероводорода, оксида углерода и
диоксида серы (см. [41]), аммиака (см. [64) и неона (см. [62]). Все значения приведены для
стандарт ного давления р0, равного 101,325 кПа. Следует учитывать, что перечисленные значения
коэффици ента сжимаемости частично взаимоконвертируемы с коэффициентами суммирования в
ГОСТ 31369 (см. также [1]), и значения для температуры 0 °С не идентичны наиболее
соответствующим значениям Z в таблице 3. Последнюю цифру для значений, приведенных в таблице
4, следует рассматривать как помощь в интерполяции.
8.4 Неопределенности чистых компонентов
8.4.1 Неопределенность вторых вириальных коэффициентов
Используя формулу (31), можно математически вывести стандартную неопределенность u(Sj) в
коэффициенте суммирования Sjдля чистого компонента, заданную выражением
u(sj)/sj =u(Bjj)/2-Bjj,(40)
где u(Bjj) — стандартная неопределенность второго вириального коэффициента В;у.
Однако это тождество может вводить в заблуждение, поскольку, как объяснено в 7.4, формула (31)
основана на двух приближениях, при которых формула (20) одновременно усекается на втором члене
и линеаризуется таким образом, что преобразуется в формулу (25).
Поэтому для получения истинных оценок неопределенностей для коэффициентов суммирования
необходимо учитывать вклады из трех источников, только первый из которых учтен в формуле (40).
Для данного вклада значения и(В^)для большинства компонентов, рассматриваемых вГОСТ31369
(см. также [1]), — (см. [28]). Для тех немногих случаев, когда данные по Вуу и u(Bjj) отсутствовали
(см. [28]), первая составляющая величины
u(Sj)
установлена равной 0,01.
8.4.2 Ошибка усечения
Несмотря на их разное концептуальное происхождение, влияние усечения и линеаризации на
неопределенность можно объединить и рассматривать совокупно. Это может быть достигнуто путем
оценки Zy(n, следовательно, sy), которую вычисляют:
a) исходя из двухчленного вириального разложения по давлению р [формула (25)], и
b
) от более строгого трехчленного (или более) вириального разложения по обратному молярному
объему MV или молярной плотности D [формула (20)], и устанавливают полученное значение как ошиб ку
усечения (точнее, смещение усечения) в t/(Zy).
Для целей данного сравнения второй вириальный коэффициент Вууможно оценить из корреляции
(см. [30], [формула (37)], а третий вириальный коэффициент С^ из корреляции (см. [65]) следующим
образом:
(
r
C
)
j
=(
0
,
01407
-
0
,
02676 .m)+ (0’02432%°801770
“ ) +^^- ^^
(R ■Тс)ТгТгТг
- (
0
’
00313
\ °
f
228
Ю>. (
41
)
Тг
Эта процедура соответствует рекомендованной процедуре (см. [28]), где она объясняется более
подробно. Критические значения температуры Тс, давления рс и ацентрического фактора со, необходи
мые для решения по формуле (41), взяты из таких источников, как и для формулы (37) (см. [32], [33]).
Решение по формуле (20) для Zyс этой целью требует итерационной процедуры. Второй вклад в u(sj)
затем оценивают окончательно по разности [см. формулы (20) (25)], где второй и третий вириальные
коэффициенты оценивают по формулам (37) и (41) соответственно.
21