ГОСТ Р 8.987-2020 Государственная система обеспечения единства измерений. Стандартные справочные данные. Хлорбензол. Теплофизические свойства (плотность, теплоемкость, энтальпия, энтропия, скорость звука, коэффициенты теплопроводности и вязкости) в диапазоне температуры от тройной точки не выше 700 К при давлениях не более 100 МПа

ГОСТ Р 8.987-2020

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

Государственная система обеспечения единства измерений

СТАНДАРТНЫЕ СПРАВОЧНЫЕ ДАННЫЕ. ХЛОРБЕНЗОЛ

Теплофизические свойства (плотность, теплоемкость, энтальпия, энтропия, скорость звука, коэффициенты теплопроводности и вязкости) в диапазоне температуры от тройной точки не выше 700 К при давлениях не более 100 МПа

State system for ensuriug the uniformity of measurements. Standard reference data. Chlorobenzene. Thermophysical properties (density, heat capacity, enthalpy, entropy, speed of sound, thermal conductivity and viscosity coefficients) over a temperature range from the triple point to 700 К with pressures up to 100 MPa

ОКС 07.030

Дата введения 2021-02-01

Предисловие

1 РАЗРАБОТАН Федеральным государственным унитарным предприятием "Всероссийский научно-исследовательский институт метрологической службы" (ФГУП "ВНИИМС")

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

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

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

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

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

Настоящий стандарт распространяется на жидкий и газообразный хлорбензол и устанавливает методы расчетного определения значений стандартных справочных данных плотности , энтальпии h, энтропии s, изобарной теплоемкости , изохорной теплоемкости , скорости распространения звука w, коэффициента динамической вязкости и коэффициента теплопроводности для хлорбензола как в однофазных областях, так и на линии насыщения.

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

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

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

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

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

3 Методические основы разработки стандартных справочных данных

3.1 Основные физико-аналитические модели, принятые для расчетного определения значений термодинамических свойств хлорбензола

В настоящем стандарте приведены основные физико-аналитические модели, принятые для расчетного определения значений термодинамических свойств хлорбензола, разработанные в соответствии с ГОСТ 8.566, ГОСТ Р 8.614 на основе теоретически и практически обоснованного фундаментального уравнения состояния (ФУС), выражающего свободную энергию Гельмгольца в зависимости от температуры Т и плотности .

Безразмерную свободную энергию Гельмгольца представляют в виде суммы идеально-газовой части и избыточной составляющей и вычисляют по формуле

. (1)

Для придания наиболее строгого подхода к ФУС в части учета особенностей термодинамической поверхности хлорбензола и расширения его экстраполяционных возможностей избыточную часть свободной энергии Гельмгольца представляют в виде разложения в ряд по степеням (см. [1]) приведенной температуры и приведенной плотности с оптимизируемыми полиномиальными экспоненциальными членами и вычисляют по формуле

, (2)

где ;

;

, - параметры приведения, в качестве которых принимают значения температуры и плотности фторбензола в критической точке (в [1] приняты кмоль/м, К).

Для определения значений параметров ФУС по формуле (2) и расширения его функциональных возможностей при нахождении значений коэффициентов ФУС, учитывают разнородные экспериментальные данные о термодинамических свойствах хлорбензола:

- о р, v, T-данные*;

________________

* Текст документа соответствует оригиналу. - .

- упругости насыщенных паров ,

- плотности насыщенной жидкой и газовой фаз;

- теплоемкости насыщенной конденсированной фазы ;

- изобарной теплоемкости;

- энтальпии h;

- скорости распространения звука w.

Корректность в описании термодинамической поверхности хлорбензола при обработке экспериментальных данных достигается путем ввода системы ограничений, накладываемых в виде неравенств на термодинамическую поверхность. В число основных видов вводимых ограничений включают (см. [1]):

- критические условия;

- правило Максвелла;

- контроль кривизны идеальных кривых;

- положительность значений теплоемкостей;

- правило прямолинейного диаметра;

- контролирование знаков производных для различных термодинамических величин.

Эти ограничения обеспечивают "физическую" форму поверхности состояния и улучшают экстраполяционные возможности уравнения.

Определение коэффициентов ФУС выполняют с применением алгоритма, представленного в [1], реализующего метод случайного поиска с возможностью возврата в начало процедуры поиска при неудачном шаге. При этом алгоритм модифицируют введением элементов детерминированного поиска на шаге корректировки величины шага поиска и выбора направления поиска.

В алгоритме [1] применяют аддитивный критерий оптимальности - минимизируемый функционал, представленный в соотношении (3), который образуют путем сложения выходных параметров, преобразованных к безразмерным слагаемым. Это осуществляют с помощью введения нормирующих множителей - весовых коэффициентов. Нормирование вводят для объединения нескольких выходных параметров - термодинамических свойств, имеющих в общем случае различную физическую размерность. Минимизируемый функционал содержит слагаемые, ответственные за точность аппроксимации результатов измерений разнородных данных о термодинамических свойствах, а также ограничения, накладываемые в виде неравенств на термодинамическую поверхность. Алгоритм представляют следующим соотношением:

, (3)

где W - весовой коэффициент для каждой опытной точки;

F - функция, используемая для минимизации отклонений.

Например, для изохорной теплоемкости данных, функцию вычисляют по соотношению

. (4)

Квадратичные функции для других термодинамических свойств имеют аналогичный вид. - функция, учитывающая различные ограничения на область изменения переменных.

Весовой коэффициент W для каждой выбранной экспериментальной точки назначают индивидуально с учетом типа данных, области состояний и требуемой точности. Типичное значение W для данных р, , T и давления насыщенных паров составляет 1, для теплоемкости - 0,5, для скорости звука - 1.

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

В минимизируемый функционал включают несколько слагаемых, каждое из которых ответственно за определенную категорию обрабатываемых термодинамических характеристик (см. [1]).

Для расчетного определения значений термодинамических свойств используют известные дифференциальные соотношения термодинамики (5)-(10).

Коэффициенты и показатели степени при температуре и плотности оптимизированной формуле (2) представлены в [1]. Процедура построения ФУС более подробно описана в [1].

Термодинамические свойства хлорбензола рассчитывают по следующим соотношениям:

- плотность

; (5)

- энтальпия

; (6)

- энтропия

; (7)

- изохорная теплоемкость

; (8)

- изобарная теплоемкость

; (9)

- скорость звука

, (10)

где нижний индекс при коэффициенте показывает частную производную по соответствующей переменной.

За термодинамическое начало отсчета при составлении таблиц термодинамических свойств хлорбензола принимают состояние насыщенной жидкой фазы при температуре К. Значения энтальпии и энтропии определены в [1] как кДж·кг, кДж·кг·К.

3.2 Коэффициенты переноса

3.2.1 Коэффициент вязкости

Табличные значения коэффициентов переноса определяют по эмпирическим уравнениям, разработанным на основе наиболее надежных экспериментальных данных и апробированным на практике.

Для расчетов значений коэффициента динамической вязкости используют уравнение, предложенное [1].

, (11)

где - вязкость разреженного газа при нулевой плотности;

- избыточная вязкость.

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

; (12)

, (13)

где вязкость , мкПа·с;

М=112,557 - масса киломоля, кг/кмоль;

Т - температура, K;

- линейный масштабный параметр потенциала Леннарда - Джонса [1], нм;

- энергетический масштабный параметр;

К; - приведенный эффективный интеграл столкновений, аппроксимированный уравнением (13);

Т*- приведенная температура .

Избыточную вязкость аппроксимируют уравнением

, (14)

где ;

.

Плотность рассчитывают по фундаментальному уравнению состояния (2).

В таблице 1 приведены параметры уравнений для расчета вязкости разреженного газа.

Таблица 1 - Параметры уравнений (12) и (13) для расчета вязкости разреженного газа

Поиск коэффициентов и показателей степени при температуре и плотности уравнения (14), а также коэффициентов приведенного эффективного интеграла столкновений (13) осуществляют методом случайного поиска с возвратом при неудачном шаге [1]. При поиске коэффициентов вводятся ограничения на форму поверхности состояния, обеспечивающие "правильные" знаки производных и тем самым улучшают экстраполяционные возможности уравнения. Более подробно метод описан в [1].

В таблице 2 приведены коэффициенты и показатели степени уравнения (14).

Таблица 2 - Коэффициенты и показатели степени уравнения (14)

Средняя неопределенность описания коэффициента динамической вязкости по уравнениям (11)-(14) составляет 1,28%. Уравнения (11)-(14) не учитывают критическую аномалию вязкости, которая не исследована для хлорбензола. Результаты сравнения с имеющимися экспериментальными данными представлены в [1].

3.2.2 Коэффициент теплопроводности

Расчетное определение значений коэффициента теплопроводности хлорбензола проводят согласно [1] на основе применения эмпирического уравнения вида:

, (15)

где - теплопроводность разреженного газа при нулевой плотности;

- избыточная теплопроводность;

;

;

, - опорные значения плотности и температуры (принимают критические значения =632,35 К; =3,24 кмоль/м).

Теплопроводность разреженного газа определяют по уравнению

, (16)

где - вязкость разреженного газа при нулевой плотности, мкПа·с.

Избыточную теплопроводность аппроксимируют уравнением

. (17)

Поиск коэффициентов и показателей степени при температуре и плотности уравнений (16) и (17) осуществляют методом случайного поиска с возвратом при неудачном шаге, который подробно описан [1]. Также вводятся ограничения, обеспечивающие "правильный" знак производных. Коэффициенты и показатели степени уравнений (16) и (17) представлены в таблице 3.

Таблица 3 - Коэффициенты и показатели степени уравнений (16) и (17)

Средняя неопределенность описания экспериментальных данных в жидкой фазе уравнениями (15)-(17) составляет соответственно 3,2%, в газовой фазе - 2,31%. Уравнения (15)-(17) не учитывают критическую аномалию теплопроводности, которая экспериментально не исследована для хлорбензола. Рассчитанные значения коэффициентов динамической вязкости и теплопроводности в однофазной области и на линии насыщения представлены в [1] в таблицах Б.3, Б.4.

4 Анализ и отбор экспериментальных данных

4.1 Данные о термодинамических свойствах хлорбензола

Для оценки точности расчетных значений и анализа разработки расчетных уравнений проводят анализ и отбор разнородных экспериментальных данных о термодинамических свойствах хлорбензола, полученными различными авторами. Затем проводят сравнительный анализ экспериментальных данных с данными, полученными по ФУС. Результаты сравнения данных о термодинамических свойствах хлорбензола со значениями, рассчитанными по ФУС, представлены в [1], в таблице Б.1 и на рисунках Б.1-Б.3, на которых показан характер отклонений экспериментальных данных от расчетных значений.

4.2 Данные о коэффициентах переноса хлорбензола

Для анализа разработки расчетных уравнений и оценки точности расчетных данных о динамической вязкости хлорбензола используют массив экспериментальных данных с чистотой образца 99,92%, представленный в [1]. Динамическую вязкость хлорбензола исследуют в паровой, жидкой, критической и сверхкритической областях параметров состояния, в интервале давлений от 0,1 до 40 МПа, в диапазоне температур от нормальной температуры кипения до 623,15 К. Измерения проводят методом вынесенного капилляра.

Для обработки коэффициента теплопроводности в переменных температура - плотность, используют уравнение состояния (1) и массив экспериментальных данных, представленный [1]. Средняя неопределенность описания экспериментальных данных в жидкой фазе (см. [1]) уравнением составляет соответственно 3,2%, в газовой фазе (см. [1]) - 2,31%.

5 Оценка достоверности расчетных значений свойств хлорбензола

5.1 Результаты оценки достоверности расчетных значений термодинамических свойств хлорбензола

Величины неопределенности расчетных значений термодинамических свойств хлорбензола оценивают по результатам сравнения с наиболее надежными экспериментальными данными. Оценки, представленные в [1], даны: для жидкой фазы: , , для газовой фазы: , , для сверхкритического флюида: , включая критическую область: , .

В таблице 3 (см. [1]) дана оценка неопределенностей расчетных значений термодинамических свойств хлорбензола, а поля неопределенностей приведены в таблицах В.5-В.7. В настоящем стандарте оценка неопределенностей расчетных значений термодинамических свойств и поля неопределенностей хлорбензола представлены в таблицах 4-5, 6-7.

На рисунках Б.1-Б.3 (см. [1]) продемонстрированы отклонения экспериментальных данных о плотности жидкой фазы, о давлении насыщенных паров, о скорости звука хлорбензола от рассчитанных данных по ФУС. На рисунке Б.4 (см. [1]) показан ход идеальных кривых хлорбензола. На диаграммах Б.5-Б.7 (см. [1]) продемонстрированы поверхности состояния основных термодинамических свойств, построенные по ФУС. Вид этих поверхностей свидетельствует о хороших интерполяционных и экстраполяционных свойствах хлорбензола, разработанного ФУС.

Таблица 4 - Оценки неопределенности расчетных значений термодинамических свойств хлорбензола

Таблица 5 - Поля неопределенности расчета плотности

Таблица 6 - Поля неопределенности расчета изобарной теплоемкости

Таблица 7 - Поля неопределенности расчета скорости распространения звука w

5.2 Результаты оценки достоверности расчетных данных о коэффициентах переноса

Поля неопределенностей расчета коэффициентов переноса представлены в таблицах В.8, В.9 (см. [1]). В настоящем стандарте - в таблицах 8, 9. На диаграммах Б.8, Б.9 (см. [1]), продемонстрированы поверхности состояний интерполяционных и экстраполяционных возможностей уравнений по вязкости и теплопроводности. Анализ результатов отклонений экспериментальных данных показывает, что вероятная неопределенность расчета коэффициента теплопроводности хлорбензола составляет 3,0-5,0%. Неопределенность экспериментального определения вязкости составляет 1,2%.

Таблица 8 - Поля неопределенности расчета коэффициента теплопроводности

Таблица 9 - Поля неопределенности расчета коэффициента динамической вязкости

Оценки достоверности расчетных значений теплофизических свойств хлорбензола на линии равновесия "жидкость - газ" представлены в таблице В.10 (см. [1]). В настоящем стандарте поля неопределенности расчета теплофизических свойств хлорбензола на линии равновесия "жидкость - газ" приведены в таблице 10.

Таблица 10 - Поля неопределенности расчета теплофизических свойств хлорбензола на линии равновесия "жидкость - газ"

Итоговые значения термодинамических свойств и коэффициентов переноса хлорбензола в однофазной области и на линии насыщения представлены в [1], в таблицах стандартных справочных данных Б.3, Б.4.

Библиография

Электронный текст документа

и сверен по:

, 2020