ГОСТ Р 8.989-2020. Национальный стандарт Российской Федерации. Государственная система обеспечения единства измерений. Стандартные справочные данные. Бензол жидкий и газообразный. Термодинамические свойства, коэффициенты динамической вязкости и теплопроводности при температурах от 280 К до 725 К и давлениях до 100 МПа

(утв. и введен в действие Приказом Росстандарта от 29.04.2020 № 182-ст)

Утвержден и введен в действие Приказом Федерального агентства по техническому регулированию и метрологии от 29 апреля 2020 г. № 182-ст

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

ГОСУДАРСТВЕННАЯ СИСТЕМА ОБЕСПЕЧЕНИЯ ЕДИНСТВА ИЗМЕРЕНИЙ

СТАНДАРТНЫЕ СПРАВОЧНЫЕ ДАННЫЕ. БЕНЗОЛ ЖИДКИЙ И ГАЗООБРАЗНЫЙ

ТЕРМОДИНАМИЧЕСКИЕ СВОЙСТВА, КОЭФФИЦИЕНТЫ ДИНАМИЧЕСКОЙ ВЯЗКОСТИ И ТЕПЛОПРОВОДНОСТИ ПРИ ТЕМПЕРАТУРАХ ОТ 280 К ДО 725 К И ДАВЛЕНИЯХ ДО 100 МПа

State system for ensuring the uniformity of measurements. Standard reference data. Liquid and gaseous benzene. Thermodynamic properties, dynamic viscosity and thermal conductivity at temperatures from 280 K to 725 K and pressures up to 100 MPa

ГОСТ Р 8.989-2020

ОКС 07.030

Дата введения 1 февраля 2021 года

Предисловие

1 РАЗРАБОТАН Главным научным метрологическим центром «Стандартные справочные данные о физических константах и свойствах веществ и материалов» (ГНМЦ «ССД»)

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 180 «Стандартные справочные данные о физических константах и свойствах веществ и материалов»

3 УТВЕРЖДЕН И ВВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 29 апреля 2020 г. № 182-ст

4 ВВЕДЕН ВПЕРВЫЕ

Правила применения настоящего стандарта установлены в статье 26 Федерального закона от 29 июня 2015 г. № 162-ФЗ «О стандартизации в Российской Федерации». Информация об изменениях к настоящему стандарту публикуется в ежегодном (по состоянию на 1 января текущего года) информационном указателе «Национальные стандарты», а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет (www.gost.ru)

1 Область применения

Настоящий стандарт распространяется на жидкий и газообразный бензол и устанавливает методы расчетного определения значений стандартных справочных данных по плотности \rho , энтальпии h, энтропии s, изобарной теплоемкости cp, изохорной теплоемкости cv, скорости звука w, коэффициента динамической вязкости \mu и коэффициента теплопроводности \lambda как в однофазных областях (газ, жидкость и флюид), так и на линии фазового перехода газ — жидкость (линии насыщения), а также значений давления на линии насыщения ps.

2 Нормативные ссылки

В настоящем стандарте использованы нормативные ссылки на следующие стандарты:

ГОСТ 8.566 Государственная система обеспечения единства измерений. Межгосударственная система данных о физических константах и свойствах веществ и материалов. Основные положения

ГОСТ Р 8.614 Государственная система обеспечения единства измерений. Государственная служба стандартных справочных данных. Основные положения

Примечание - При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет или по ежегодному информационному указателю «Национальные стандарты», который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя «Национальные стандарты» за текущий год. Если заменен ссылочный стандарт, на который дана недатированная ссылка, то рекомендуется использовать действующую версию этого стандарта с учетом всех внесенных в данную версию изменений. Если заменен ссылочный стандарт, на который дана датированная ссылка, то рекомендуется использовать версию этого стандарта с указанным выше годом утверждения (принятия). Если после утверждения настоящего стандарта в ссылочный стандарт, на который дана датированная ссылка, внесено изменение, затрагивающее положение, на которое дана ссылка, то это положение рекомендуется применять без учета данного изменения. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, рекомендуется применять в части, не затрагивающей эту ссылку.

3 Общие положения

Стандартные справочные значения (ГОСТ Р 8.614, ГОСТ 8.566) \rho , h, s, cp, cv, w и ps рассчитаны по единому для жидкой и газовой фаз фундаментальному уравнению состояния (ФУС) — зависимости свободной энергии (функции Гельмгольца) F от плотности \rho и температуры T

\frac{F(\rho,T)}{RT} = f(\omega,\tau) = f_{0}(\omega,\tau) + f_{r}(\omega,\tau) , (1)

где f, f0, fr — безразмерные полная свободная энергия, идеально-газовая и неидеальная составляющие свободной энергии, соответственно;

\omega — относительная плотность, \omega = \rho\rho_{\text{кр}} ;

\tau — относительная температура, \tau = TT_{\text{кр}} .

Значения плотности \rho_{\text{кр}} и температуры Tкр бензола в критической точке приведены в таблице А.1 приложения А.

Уравнение для идеально-газовой составляющей свободной энергии имеет следующий вид

f_{0} = ln(\omega) + \alpha_{1} + \alpha_{2}\tau^{- 1} + \alpha_{3}ln(\tau^{- 1}) + \sum_{i = 4}^{6} \alpha_{i}ln[1 - exp(- \delta_{i}\tau^{- 1})] . (2)

Коэффициенты {_{\alpha_{i}}} и параметры {_{\delta_{i}}} уравнения (2) приведены в таблице А.3 приложения А.

Уравнение для неидеальной составляющей свободной энергии имеет следующий вид

f_{r} = \sum_{j = 1}^{14} b_{j}\phi_{j} , (3)

где

(4)

В формулах (3) и (4) bj — коэффициенты уравнения состояния, значения которых вместе с показателями степеней rj, tj, lj и параметрами gj, \alpha_{j} , \beta_{j} , \epsilon_{j} , приведены в таблице А.2 приложения А.

Плотность \omega в однофазных областях при заданных значениях давления p и температуры T определяют из решения следующего уравнения

\pi = \omega\tau(1 + A_{0})z_{\text{кр}} , (5)

где \pi = pp_{\text{кр}} ;

z_{\text{кр}} = 10^{3}p_{\text{кр}}(\rho_{\text{кр}}RT_{\text{кр}}) .

Значения давления pкр и фактора сжимаемости zкр в критической точке, а также газовой постоянной R бензола приведены в таблице А.1 приложения А.

Плотности газовой \omega\prime\prime и жидкой \omega\prime фаз на линии насыщения при заданной температуре T определяют из условий фазового равновесия в результате решения следующей системы уравнений:

, (6)

где \phi_{r}(\tau,\omega) — безразмерная неидеальная составляющая изобарно-изотермического потенциала (потенциала Гиббса)

\phi_{r} = f_{r} + A_{0} + ln(\omega) . (7)

Давление на линии насыщения ps определяют по формуле (5) для \omega\prime\prime .

Энтальпия, энтропия, изобарная и изохорная теплоемкости и скорость звука как в однофазных областях (для T и \omega ), так и на линии насыщения (для T, \omega\prime или T, \omega\prime\prime ) вычисляют по следующим формулам:

h = h₀ + A₃RT; (8)

s = s₀ + RA₄; (9)

cp = cv + R(1 + A2)2/(1 + A1); (10)

cv = cv0 + A5R; (11)

w = [103RTcp(1 + A1)/cv]0,5, (12)

где h0, s0, cv0 — энтальпия, энтропия и изохорная теплоемкость в идеально-газовом состоянии.

Термодинамические свойства в идеально-газовом состоянии определяют по формулам, полученным из f_{0}(\tau,\omega) с привлечением табличных данных (см. [1]).

c_{v0} = R[\alpha_{3} + \sum_{i = 4}^{6} \alpha_{i}E_{i}D_{i}^{2}] ; (13)

h_{0} = RT[1 + \alpha_{3} + \alpha_{2}\Theta + \sum_{i = 4}^{6} \alpha_{i}E_{i}D_{i} + \frac{\Delta h_{0}}{RT}] ; (14)

s_{0} = R\{\alpha_{3}(1 - ln\Theta) - \alpha_{1} + \sum_{i = 4}^{6} \alpha_{i}[E_{i}D_{i} - ln(1 - E_{i})] + \Delta s_{0}R - ln\omega\} , (15)

где \Theta = \tau^{- 1} ;

Ei и Di — функции от \Theta , имеющие следующий вид:

E_{i} = exp( - \delta_{i}\Theta),D_{i} = \delta_{i}\Theta(1 - E_{i}) . (16)

Коэффициенты {_{\alpha_{i}}} в формулах (13) - (15) и параметры {_{\delta_{i}}} в формуле (16), а также значения энтальпии \Delta h_{0} и энтропии \Delta s_{0} приведены в таблице А.3 приложения А. Значения \Delta h_{0} и \Delta s_{0} добавлены в выражения для h0 и s0 для удобства сравнения с табличными данными (см. [1]).

Комплексы A0 — A5 в формулах (5) - (12) определяют по следующим соотношениям, полученным из уравнения (3) для fr с использованием известных дифференциальных уравнений термодинамики

A_{0} = \sum_{j = 1}^{14} b_{j}\phi_{j}X_{j} ; (17)

A_{1} = \sum_{j = 1}^{14} b_{j}\phi_{j}[X_{j}(X_{j} + 1) + U_{j}] ; (18)

A_{2} = \sum_{j = 1}^{14} b_{j}\phi_{j}[X_{j}(Y_{j} + 1)] ; (19)

A_{3} = \sum_{j = 1}^{14} b_{j}\phi_{j}[X_{j} - Y_{j}] ; (20)

A_{4} = - \sum_{j = 1}^{14} b_{j}\phi_{j}[Y_{j} + 1] ; (21)

A_{5} = - \sum_{j = 1}^{14} b_{j}\phi_{j}[Y_{j}(Y_{j} + 1) + Q_{j}] , (22)

где

; (23)

; (24)

; (25)

(26)

Стандартные справочные значения коэффициента динамической вязкости бензола рассчитывают по следующему уравнению (см. [1]):

\mu = \mu_{0} \cdot (1 + B_{\mu} \cdot \tilde{\rho}_{\mu}) + \Delta\mu(\omega,\tau) , (27),

где \mu_{0} , B_{\mu} и \Delta\mu — соответственно, коэффициент динамической вязкости в состоянии разреженного газа, второй вязкостный вириальный коэффициент и составляющая коэффициента динамической вязкости, отражающая поведение вязкости при высоких плотностях; \tilde{\rho}_{\mu} = \rho M .

Составляющие \mu в (27) определяют по следующим уравнениям:

\mu_{0} = \frac{0,021357 \cdot \sqrt{MT}}{\sigma^{2} \cdot \Omega(T_{r})} ; (28)

\Omega(T_{r}) = exp[\sum_{i = 0}^{3} a_{i}(lnT_{r})^{i}] . (29)

В уравнениях (28, 29) относительная температура T_{r} = T(\epsilon k) ; параметр приведен в таблице А.1 приложения А, а коэффициенты {ai} уравнения (29) приведены в таблице А.4 приложения А.

B_{\mu} = 0,6022137\sigma^{3} \cdot [\sum_{i = 0}^{6} b_{i}T_{r}^{- 0,25i} + b_{7}T_{r}^{- 2,5} + b_{8}T_{r}^{- 5,5}] . (30)

Параметр \sigma приведен в таблице А.1 приложения А; коэффициенты {bi} уравнения (30) приведены в таблице А.5 приложения А.

\Delta\mu(\omega,\tau) = (\omega^{23}\tau^{12}) \cdot \{c_{0}\omega^{2} + \frac{c_{1}\omega}{c_{2} + c_{3}\tau + c_{4}\omega} + \frac{c_{5}\omega + c_{6}\omega^{2}}{c_{7} + c_{8}\omega^{2}} \} . (31)

Коэффициенты {ci} уравнения (31) представлены в таблице А.6.

Стандартные справочные значения коэффициента теплопроводности бензола рассчитывают по следующему уравнению (см. [1]):

\lambda = \lambda_{0} + \Delta\lambda + \Delta\lambda_{\text{кр}} , (32)

где \lambda_{0} , \Delta\lambda и \Delta\lambda_{\text{кр}} — соответственно, коэффициент теплопроводности в состоянии разреженного газа, избыточная по отношению к \lambda_{0} составляющая коэффициента теплопроводности и аномальная составляющая коэффициента теплопроводности в околокритической области.

\lambda_{0} = \frac{\sum_{i = 1}^{3} a_{i}\tau^{i - 1}}{\sum_{i = 4}^{6} a_{i}\tau^{i - 4}} ; (33)

\Delta\lambda = \sum_{i = 1}^{5} (b_{1i} + b_{2i}\tau) \cdot \omega^{i} ; (34)

(35)

\Delta\chi = [\frac{\chi(\tau,\omega) - \chi(\tau_{ref},\omega)T_{ref}T^{- 1}}{\text{Г}}] ; (36)

\xi = \xi_{0}\Delta\chi^{\nu\gamma} ; (37)

\tilde{\Omega} = \frac{2}{\pi} [(1 - \frac{c_{v}}{c_{p}})arctg(\xi q_{D}) + \frac{c_{v}}{c_{p}} (\xi q_{D})] ; (38)

\tilde{\Omega}_{0} = \frac{2}{\pi} \{1 - exp[\frac{- 1}{(\xi q_{D})^{- 1} + \frac{1}{3} (\xi q_{D}^{- 1}\omega^{- 1})^{2}}]\} ; (39)

\chi(\tau,\omega) = \frac{\omega z_{c}}{\tau[1 + A_{1}(\tau,\omega)]} . (40)

Значения коэффициентов {ai} и {b1i}, {b2i} уравнений (33, 34) приведены в таблице А.7 и А.8 приложения А. Универсальные теоретически обоснованные постоянные kB = 1,380658·10-2; R0 = 1,02; v = 0,63; \gamma = 1,239 . Значения подгоночных параметров для бензола \xi_{0} , Г, qD, Tref приведены в таблице А.1 приложения А. В выражениях (35), (38) - (40) теплоемкости cp, cv и расчетный комплекс A1 определяются по формулам (10) - (11) и (18), соответственно; коэффициент динамической вязкости \mu(\tau,\omega) определяется по формулам (28) - (31).

Рассчитанные стандартные справочные данные контрольных значений термодинамических и переносных (\mu,\lambda) свойств бензола приведены в таблицах Б.2 (линия насыщения) и В.1 (однофазные области).

4 Расширенные неопределенности расчетных значений стандартных справочных данных по свойствам бензола

Расширенные неопределенности с доверительной вероятностью 95% расчетных значений термодинамических свойств: плотности \delta\rho = \Delta\rho\rho , скорости звука \delta w = \Delta ww , изохорной \delta c_{v} = \Delta c_{v}c_{v} и изобарной \delta c_{p} = \Delta c_{p}c_{p} теплоемкостей, а также давления насыщения \delta p_{s} = \Delta p_{s}p_{s} определяются в соответствии с оценками, приведенными в [1].

Для околокритической области авторы работы (см. [1]) вместо \delta\rho приводят значения \delta p = \Delta pp = 0,2% , поэтому значения \delta\rho определяют в соответствии с теорией переноса ошибок по формуле

\delta\rho = (\frac{1 + A_{0}}{1 + A_{1}}) \cdot 0,2 . (41)

Расширенные неопределенности расчетных значений энтальпии \Delta h и энтропии \delta s = \Delta ss определяют в соответствии с теорией переноса ошибок через значение \delta\rho по следующим выражениям:

\Delta h = 0,1 + RT|(\frac{\partial A_{3}}{\partial\omega})_{\tau}|\omega\delta\rho/100 , кДж/кг; (42)

\delta s = \{0,01 \cdot s_{0}(\tau) + R|\omega(\frac{\partial A_{4}}{\partial\omega})_{\tau} - 1|\delta\rho\}s^{- 1} , %. (43)

В формулах (42) - (43) A3 и A4 - расчетные комплексы (20) - (21); s_{0}(\tau) рассчитывают по формуле (15), но без учета ln(\omega) .

Расширенные неопределенности с доверительной вероятностью 95% расчетных значений коэффициента динамической вязкости \delta\mu = \Delta\mu\mu и коэффициента теплопроводности \delta\lambda = \Delta\lambda\lambda определяют в соответствии с оценками, приведенными в [1].

Расширенные неопределенности расчетных значений стандартных справочных данных бензола представлены в таблицах 10 и 11 (см. [1]), где для всех теплофизических свойств, кроме энтальпии, приведены относительные величины неопределенностей \delta A = 100 \cdot \Delta AA , %. Для энтальпии приведена абсолютная величина \Delta h , кДж/кг.

Приложение А (обязательное)

ОСНОВНЫЕ ФИЗИЧЕСКИЕ ПАРАМЕТРЫ И КОЭФФИЦИЕНТЫ УРАВНЕНИЙ ДЛЯ ОПРЕДЕЛЕНИЯ ЗНАЧЕНИЙ СТАНДАРТНЫХ СПРАВОЧНЫХ ДАННЫХ ПО СВОЙСТВАМ БЕНЗОЛА

Таблица А.1

Основные физические параметры бензола

Молярная масса M, кг/кмоль	78,1118
Газовая постоянная R, кДж/(кг·К)	0,1064432
Параметры в тройной точке:
давление p_{\tau} , МПа	0,004785
температура T_{\tau} , К	278,674
Параметры в критической точке:
давление pкр, МПа	4,894
температура Tкр, К	562,02
плотность \rho_{\text{кр}} , кг/м3	304,792
фактор сжимаемости zкр	0,268405066
Параметры для расчета \mu :
\sigma , нм	0,540
, К	412,0
Параметры для расчета \Delta\lambda_{\text{кр}} :
\xi_{0} , нм	0,216
Г	0,0569
qD, нм	0,620
Tref, К	843,0

Таблица А.2

Коэффициенты, показатели степеней и параметры уравнения для неидеальной составляющей ФУС бензола [см. уравнения (3) и (4)]

Таблица А.3

Коэффициенты уравнений (2) и (13) - (16) для термодинамических свойств бензола в идеально-газовом состоянии, энтальпия \Delta h_{0} и энтропия \Delta s_{0}

Таблица А.4

Коэффициенты \alpha_{i} уравнения (29) для \mu_{0} бензола

Таблица А.5

Коэффициенты bi уравнения (30) для B_{\mu} бензола

Таблица А.6

Коэффициенты {ci} уравнения (31) для \Delta\mu бензола

Таблица А.7

Коэффициенты ai уравнения (33) для \lambda_{0} бензола

Таблица А.8

Коэффициенты bji уравнения (34) для \Delta\lambda бензола

Приложение Б (обязательное)

ТАБЛИЦЫ КОНТРОЛЬНЫХ СТАНДАРТНЫХ ЗНАЧЕНИЙ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ БЕНЗОЛА НА КРИВОЙ НАСЫЩЕНИЯ

Таблица Б.1

Обозначения и размерности теплофизических свойств и их неопределенностей, представленных в таблицах Б.2 и В.1 приложений Б, В

Таблица Б.2

Контрольные стандартные значения теплофизических свойств бензола на кривой насыщения

Окончание таблицы Б.2

Приложение В (обязательное)

ТАБЛИЦЫ КОНТРОЛЬНЫХ СТАНДАРТНЫХ ЗНАЧЕНИЙ ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ БЕНЗОЛА В ОДНОФАЗНОЙ ОБЛАСТИ

Таблица В.1

Контрольные стандартные значения теплофизических свойств бензола в однофазной области

БИБЛИОГРАФИЯ