На правах рукописи

РЫКОВ АНДРЕЙ ВЛАДИМИРОВИЧ

МЕТОД ПОСТРОЕНИЯ ФУНДАМЕНТАЛЬНOГО УРАВНЕНИЯ

СОСТОЯНИЯ И ТЕРМОДИНАМИЧЕСКИЕ ТАБЛИЦЫ ПЕРФТОРПРОПАНА

Специальность 01.04.14 – Теплофизика и теоретическая теплотехника

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

Санкт-Петербург 2 2013

Работа выполнена в Институте холода и биотехнологий, ФГБОУ ВПО «Санкт-Петербургский национальный исследовательский университет информационных технологий, механики и оптики»

Научный руководитель: доктор технических наук, профессор Самолетов Владимир Александрович

Официальные оппоненты: Коваленко Анатолий Николаевич доктор технических наук, ведущий научный сотрудник ФТИ им. А.Ф.Иоффе Прошкин Станислав Станиславович кандидат технических наук, доцент, доцент кафедры физики СПб МНИУ ГУ

Ведущая организация: ООО «ЛЕННИИХИММАШ»

Защита диссертации состоится 18 декабря 2013 г. в 14 часов на заседании диссертационного Совета Д 212.227.08 при Санкт-Петербургском национальном исследовательском университете информационных технологий, механики и оптики по адресу: 191002, Санкт-Петербург, ул. Ломоносова, д. 9, тел./факс: (812) 315-30-15.

С диссертацией можно ознакомиться в библиотеке университета.

Автореферат разослан 15 ноября 2013 г.

Ученый секретарь Владимир Алексеевич диссертационного Совета Рыков

Общая характеристика работы

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

Цель работы Разработка метода построения масштабного уравнения в физических переменных в рамках феноменологической теории Мигдала и обоснование перехода от фундаментального масштабного уравнения в физических переменных, содержащего интегралы от дифференциальных биномов, к фундаментальному уравнению состояния скейлингового типа, на основе которого все термодинамические функции описываются алгебраическими выражениями, не содержащими интегралов. Разработка метода построения масштабных и единых уравнений, удовлетворяющих МТ, в которых нелинейные параметры являются универсальными величинами с точностью до универсальности критических индексов и коэффициента пропорциональности x0, связывающего линию насыщения и обобщенную масштабную переменную. Апробация масштабных уравнений состояния на примере описания термодинамической поверхности трифторметана; фундаментального асимметричного единого уравнения на примере описания равновесных свойств аргона и фундаментального единого уравнения состояния на примере описания равновесных свойств перфторпропана.

Разработка на основе фундаментального единого уравнения состояния перфторпропана точных термодинамических таблиц на линии насыщения, в широкой окрестности критической точки, а также в однофазной области в диапазоне параметров состояния по температуре от 130 K до 500 K и по давлению от 0,0001 МПа до 70 МПа.

Научная новизна Впервые разработано масштабное уравнение состояния в физических переменных на основе гипотезы Бенедека, которому дано строгое обоснование в рамках феноменологической теории Мигдала. Отработана методика расчета параметров масштабных функций через критические индексы и значение масштабной переменной на линии насыщения.

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

Предложенные масштабные функции с универсальными нелинейными параметрами апробированы на примере построения фундаментального асимметричного единого УС аргона. Показано, что включение в структуру УС таких масштабных функций позволяет улучшить расчетные характеристики УС как в области низких температур, так и в широкой окрестности КТ.

Разработана методика построения фундаментального единого уравнения состояния перфторпропана, в котором все нелинейные параметры универсальны в той же мере, в какой универсальны критические индексы. Показано, что полученное уравнения состояния позволяют количественно верно рассчитать термические и калорические данные, как в широкой окрестности критической точки, так и в регулярной части термодинамической поверхности. Причем при описании регулярной части термодинамической поверхности предложенное фундаментальное единое уравнение состояния не уступает аналитическим УС по точности описания как термических, так и калорических свойств R218.

Разработаны новые таблицы термодинамических свойств перфторпропана в однофазной области и на линии насыщения, в том числе в окрестности критической точки.

Автор защищает:

- новое масштабное уравнение состояния в физических переменных, обоснованное в рамках гипотезы Бенедека и феноменологической теории Мигдала;

- метод выбора структуры масштабной функции свободной энергии в переменных плотность-температура и алгоритм расчета ее нелинейных параметров;

- систему взаимосогласованных уравнений, описывающих линию фазового равновесия, и выведенное на ее основе новое, физически обоснованное, уравнение для «кажущейся» теплоты парообразования;

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

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

- фундаментальное асимметричное единое уравнение состояния аргона, структурно включающее асимптотическую масштабную функцию свободной энергии с универсальными нелинейными параметрами и имеющее рабочую область: по температуре – от тройной точки и до 1200 K, по плотности – от 0 до 3,33 c ;

- фундаментальное единое уравнение состояния перфторпропана, которое удовлетворяет масштабной гипотезе и не имеет ни одного индивидуального нелинейного параметра в сингулярной составляющей УС;

- таблицы термодинамических свойств перфторпропана, рассчитанные в диапазоне параметров состояния: по температуре от 130 K до 500 K и по давлению от 0,0001 МПа до 70 МПа;

Практическая ценность работы Разработанные методы расчета термодинамических свойств рабочих веществ в однофазной области, на линии фазового равновесия и в области метастабильных состояний оформлены в виде пакета прикладных программ на языке Фортран и в пакете MathCad и могут быть использованы для разработки масштабных и единых неаналитических УС в физических переменных.

На основе полученного единого неаналитического УС разработаны подробные термодинамические таблицы R218 как в регулярной части термодинамической поверхности, так и в широкой окрестности критической точки.

Предложенные методы расчета равновесных свойств технически важных веществ, включая криогенные газы и жидкости и хладагенты используются в учебном процессе при обучении бакалавров, специалистов и магистров по направлениям: 140700 – «ядерная энергетика и теплофизика», 141200 – «холодильная, криогенная техника и системы жизнеобеспечения», 190600 – «эксплуатация транспортно-технологических машин и комплексов».

Основные результаты работы докладывались и обсуждались на IV международной научно-технической конференции «Низкотемпературные и пищевые технологии в XXI веке» (Санкт-Петербург, 2009); научно-технической конференции с международным участием «Холодильные агенты на все времена.

Евроожидания и Российский опыт» (Санкт-Петербург, 2010); международной конференции с элементами научной школы для молодежи «Инновационные разработки в области техники и физики низких температур» (Москва, 2010);

научно-технической конференции с международным участием «Холод – 2011.

Проэкология и энергосбережение» (Санкт-Петербург, 2011); V международной научно-технической конференции «Низкотемпературные и пищевые технологии в XXI веке» (Санкт-Петербург, 2011); международной конференции с элементами научной школы для молодежи «Инновационные разработки в области техники и физики низких температур» (Москва, 2011); научно-технической конференции с международным участием «Киотский протокол за чертой года – экологические доминанты и императивы будущего индустрии холода»

(Санкт-Петербург, 2012); II международной научно-технической конференции «Современные методы исследований теплофизических свойств веществ»

(Санкт-Петербург, 2012); научно-технической конференции с международным участием «25-летие Монреальского протокола по озоноразрушающим хладагентам в контексте экологической бивалентности и доминирующей реальности» (Санкт-Петербург, 2013); XIX International Conference of Chemical Thermodynamics in Russia (Москва, 2013); IX International scientific conference «Modern problem of refrigeration equipment and technology» (Одеса, 2013); VI международной научно-технической конференции «Низкотемпературные и пищевые технологии в XXI веке» (Санкт-Петербург, 2013).

Публикации Основные результаты диссертации опубликованы в 23 печатных работах, в том числе 4 работы в изданиях, рекомендованных ВАК РФ.

Структура и объем работы.

Диссертация состоит из введения, четырех глав, выводов, списка литературы (145 наименований) и приложения. Содержание работы изложено на страницах машинописного текста, содержит 85 рисунков и 11 таблиц.

СОДЕРЖАНИЕ РАБОТЫ

Для описании широкой области параметров состояния, включая включая область критических и метастабильных состояний в основном используются два подхода. Первый связан с тем, что для описания окрестности критической точки применяется масштабное уравнение в той или иной форме. Так для описания асимптотической окрестности критической точки обычно используют уравнение, предложенное в 1965 г. Вайдомом:

где масштабная функция химического потенциала h ( x ), являющаяся однородной функцией масштабной переменной x = /, выбирается таким образом, чтобы передавать характерные особенности термодинамической поверхности в области сильно развитых флуктуаций плотности.

Для описания регулярной части термодинамической поверхности используют при этом УС или системы согласованных УС, не приводящие к возникновению особенностей в критической точке, например, вириального вида.

Альтернативные подход связан с использованием уравнений состояния, широкодиапазонного или единого, неаналитического вида, которое передает основные закономерности в поведении и регулярной и нерегулярной части термодинамической поверхности.

Развитию как первого, так и второго подхода посвящено большое число работ, среди которых можно отметить работы Абдулагатова И.М., Алтунина В.В., Анисимова М.А., Геллера В.З., Вассермана А.А., Клецкого А.В., Рабиновича В.А., Спиридонова Г.А., Сычева В.В., Лысенкова В.Ф., Платунова Е.С., Рыкова В.А., Кудрявцевой И.В., Яковлевой М.В., Мигдала А.А., Мартынеца В.Г., Паташинского А.З., Матизена Э.В. и многих других исследователей. Этой же области посвящены труды многих иностранных исследователей: Скофилда, Литстера, Хо, Сенжерса, Леммона и др.

Несмотря на достигнутые успехи в области построения уравнений состояния, «работающих» в окрестности критической точки в соответствии со степенными законами, вытекающими из современной теории критических явлений, многие задачи еще требуют своего решения.

Одна из них связана с выбором структуры масштабных функций, например h ( x ), входящей в (1), или масштабной функции свободной энергии Гельмгольца, с которой h ( x ) связанна зависимостью:

Эта задача особенно актуальна для масштабных уравнений состояния в физических переменных, которые известны также как непараметрические уравнения состояния скейлингового вида. Причем предпринято несколько попыток построить такое уравнение на основе гипотезы Бенедека. В работе рассмотрены два пути решения этой проблемы. Показано, что она имеет решение в рамках феноменологической теории Мигдала:

где ( m ) = m + 3m +... ; µ = ( c pc ) ( µ (, T ) µ 0 ( T ) ) ; c – критическая плотность; pc – критическое давление; – плотность; T – абсолютная температура;

µ 0 ( T ) – регулярная функция температуры; и – критические индексы; m – переменная, определяемая на основе равенства: 3 – постоянный коэффициент.

Изотермическая сжимаемость KT выбрана в соответствии с гипотезой Бенедека:

Масштабное уравнение в виде:

Уравнение линии насыщения найдем из равенства µ = 0 :

Так как µ = ( F ) T, то выражение для свободной энергии F, рассчитанное на основе (5), имеет вид:

Показано, что на линии x = x1 изотермическая сжимаемость обращается в нуль в каждой точке, за исключением критической точки, т.е. x = x1 определяет геометрическое место точек, удовлетворяющих равенствам, доказанным в рамках МТ Рыковым: ( T s ) = 0 ( p ) T = 0, а следовательно, уравнение x = x1 задает на термодинамической поверхности линию расходимости Cv и качественно верно, в соответствии с ( p ) T = 0, передает спинодаль.

Установлена связь уравнения состояния (7) с известными масштабными уравнениями состояния. Показано, что вычисление интегралов, входящих в (7), путем замены переменной и интегрирования по частям приводит в первом приближении к следующему выражению для свободной энергии:

Масштабная функция в (8) принадлежит к классу функций, исследованных Лысенковым и Рыковым:

В этом классе рассчитана, получившая наибольшее распространение, функция В отличие от масштабных непараметрических уравнений состояния, разработанных на основе феноменологической теории Мигдала, термодинамические функции, рассчитанные на основе (8)(10) имеют простую структуру и не содержат интегралов от дифференциальных биномов.

Линейные коэффициенты A1, A2, B1 устанавливались по известной методике, а для поиска нелинейных параметров x1, x2 и x3 были исследованы несколько алгоритмов. Показано, что оптимальным является алгоритм, построенный на системе равенств:

где h, h', f и hl, hl', fl – масштабные функции, рассчитанные на основе (10) и линейной модели (ЛМ) Скофилда-Литстера-Хо, соответственно.

Результаты сравнения масштабных функций h, h', f с рассчитанными по предложенной методике параметрами 1 = x1 / x0, 2 = x2 / x0, 3 = x3 / x0, и hl, hl', fl показали, что расхождения между ними не превосходят 0,9% во всем рабочем диапазоне линейной модели.

Масштабная функция (10) с приведенными выше параметрами i { 1, 2, 3} использована при разработке фундаментальных УС хладагента R23, аргона и перфторпропана. При апробации ниже приведенных фундаментальных УС учтено, что при их построении нельзя использовать правило Максвелла (в силу расходимости производной ( p ) T на линии сингулярности Cv ), поэтому в функционал, на основе которого устанавливаются коэффициенты этих УС, должна быть включена составляющая, обеспечивающая согласованность +,, ps, Ts -данных на линии фазового равновесия. Кроме того, линия фазового равновесия играет важную роль в формировании уравнений в физических переменных, удовлетворяющих МТ, в структуру которых входит обобщенная масштабная переменная x = / s, где функция s определяется через линию насыщения Ts по формуле s = x0 ( Ts / Tc 1). С целью оптимального выбора структуры системы уравнений для расчета линии фазового равновесия (ЛФЗ), на основе которой рассчитывается x были исследованы наиболее известные модели поведения линии насыщения в широкой окрестности критической точки:

В результате получены, во-первых, две новые системы уравнений для ЛФЗ, удовлетворяющие соответственно а) или б), во-вторых, новая физически обоснованная форма для «кажущейся» теплоты парообразования.

Новые уравнения для ЛФЗ были апробированы при описании равновесных свойств +,, ps, Ts и r * аргона от тройной точки до критической.

Предложено масштабное уравнение состояния:

где a2 ( x ) = B2 ( x + x1 ) ( x + x2 ) + C2. Введение функции a2 ( x ), обусловлено тем, что в разложении функции a1 ( x ) на критической изохоре ( x ), помимо членов, пропорциональных x 2 и x, присутствует член, пропорциональный x 1, и его надо компенсировать. В результате границы рабочей области параметров состояния масштабного уравнения (12) составили: по плотности 0,55 с 1,45 с ; по температуре Ts T 1,15Tc. Апробация уравнения (12) проводилась на примере описания равновесных свойств R23.

С целью расширить рабочую область масштабных уравнений состояния, предложено уравнение:

Рассмотрено фундаментальное асимметричное уравнение состояния Кудрявцевой:

которое апробировано на примере описания равновесных свойств аргона в диапазоне параметров состояния от разряженного газа до границы жидкость-твердое тело, от тройной точки до 1200 К, включая окрестность критической точки и метастабильную область. Во всем указанном диапазоне основной массив разнородной экспериментальной информации передан в пределах экспериментальной погрешности. Показано, что в околокритической области погрешность описания уравнением (14) опытных данных Воронеля и Анисимова с соавторами в десятки раз меньше, чем кроссоверным уравнением Rizi-Abbaci (2012 г.).

Масштабная функция с универсальными нелинейными параметрами (10) включена в выражение для свободной энергии Гельмгольца:

где f (, t ) – новая кроссоверная функция f (, t ) = exp( a 2 / 5 ) / T, обеспечивающая качественно верное описание вириальных коэффициентов.

Количественный анализ фундаментального УС (15) осуществлялся на основе надежной экспериментальной информации о термических и калорических данных перфторпропана, относящихся к однофазной области и линии фазового равновесия.

Схема поиска коэффициентов единого УС (15) существенно упрощена по сравнению со схемой, используемой при поиске параметров фундаментального асимметричного УС (14). Это обусловлено тем, что в отличие от уравнения состояния (14), УС (15) не содержит неизвестных нелинейных параметров.

Поэтому сначала на основе имеющейся экспериментальной информации о равновесных свойствах перфторпропана составляется матрица плана и на ее основе составляется функционал из условия минимума которого и определяются коэффициенты фундаментального единого УС (15). В функционале (16) Фp, ФCv, ФB обеспечивают передачу однофазной части термодинамической поверхности, а Фps и Фµ ± согласованность свойств µ ± Ts ± и +,, ps, Ts -данных на линии фазового равновесия.

В массив экспериментальной информации, формирующей матрицу плана, включены p T -данные Брауна (345,15 К 360,15 К, 1968 г.) Баpышева (133,15 К 340,15 К, 1980 г.), Клофмара (173,15 К 340,15 К, 2012 г.). Средние квадратические отклонения рассчитанных по фундаментальному неаналитическому УС (14) равновесных свойств и экспериментальных данных Баpышева В.П., Клофмара, Беляевой О.В. и др., Пономаpёвой О.П. и др., Рыкова В.А., Рябушевой Т.И. и др., Brown G., Pase E. и др., Growder G. и др. и т.д. составили в однофазной области: плотности – = 0,34%; давления – p =0,9%; изохорной теплоемкости – cv =1,2%; и на линии фазового равновесия: плотности на паровой ветви кривой сосуществования – =0,74%; плотности на жидкостной ветви кривой сосуществования - + =0,14%; «кажущейся» теплоты парообразования – r* =0,41 %; скорости звука (Hallewell, 2012 г.) – w =0,23%. Среднее квадратическое отклонение от опытных данных Барышева в области плотной жидкости составляет по плотности 0,052%. В области газа погрешность несколько выше, но не превосходит значения, рассчитанные по аналитическим уравнениям состояния Барышева и Рябушевой. Отклонения рассчитанных по уравнению состояния (15) равновесных свойств перфторпропана от экспериментальных данных представлены на рис. 1–4. Из рис. 1–4 видно, что фундаментальное неаналитическое УС (15) и качественно, и количественно воспроизводит поведение равновесные свойства перфторпропана с меньшей погрешностью, чем уравнение Леммона-Спана (2006 г.), и в окрестности критической точки, и в области плотной жидкости, и особенно, в окрестности тройной точки.

Рис. 1. Относительные отклонения плотности перфторпропана, рассчитанные по УС (15), от экспериментальных данных Brown на изохорах: 618,71 кг/м3; 564,6 кг/м3; 479,9 кг/м3; 375, кг/м3; 187,8 кг/м3. 1 – расчет по уравнению (15); 2 – расчет по уравнению Lemmon-Span.

Рис. 2. Относительные отклонения давления перфторпропана, рассчитанные по УС (15), от экспериментальных данных Барышева. Изотермы: 133,15 К; 153,15 К; 173,15 К; 193,15 К;

213,15 К; 233,15 К; 253,15 К; 273,15 К; 293,15 К; 313,15 К. 1 – расчет по уравнению (15);

Рис. 3. Зависимость Cv перфторпропана от температуры на изохоре 503,8 кг/м3. 1 – экспериментальные точки Рыкова; 2 – расчет по уравнению (15); 3 – расчет по уравнению Lemmon-Span; 4, 5 – пересечение с линией фазового равновесия.

Рис. 4. Относительные отклонения скорости звука перфторпропана в однофазной области, рассчитанные по уравнению (15), от экспериментальных данных Hallewell на изотермах:

1 –248,2 К; 2 – 259,8 К; 3 – 270,2 К; 4 – 272,6 К; 5 – 283,2 К; 6 – 296,1 К; 7 – 303,9 К.

1. В рамках феноменологической теории Мигдала и гипотезы Бенедека выведено новое непараметрическое уравнение скейлингового вида в переменных плотность и температура, которое воспроизводим все степенные законы, следующие из современной теории критических явлений. Оно удовлетворяет условиям (p/v)T=0 и (T/s)P=0, и качественно верно описывает область метастабильных состояний вещества, в частности, спинодаль.

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

3. Разработан метод выбора структуры масштабной функции свободной энергии в переменных плотность-температура и алгоритм расчета ее нелинейных параметров.

4. Разработана система взаимосогласованных уравнений линии фазового равновесия, удовлетворяющая условиям правила криволинейного диаметра в окрестности критической точки и передающая линию фазового равновесия от тройной точки до критической с экспериментальной погрешностью. Последнее показано на примере описания +,, ps и r *, r*-данных Ar, хладагентов R и R23.

5. На основе предложенной системы взаимосогласованных уравнений линии фазового равновесия выведено новое уравнение для «кажущейся теплоты»

парообразования. Обоснована структура этого уравнения в широкой окрестности критической точки и рассчитаны его параметры.

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

Полученное фундаментальное масштабное уравнение состояния апробировано на примере описания равновесных свойств хладагента R23. Показано, что его рабочая область равна: по плотности 0,55 / c 1,45, по температуре Ts T 1,15 Tc.

7. Разработано новое асимметричное масштабное уравнение состояния трифторметана, строго удовлетворяющее преобразованию равенства химических потенциалов на линии фазового равновесия, с рабочей областью: по плотности 0,4 c 1,7 c, по температуре Ts T 1,17Tc.

8. На базе рассчитанной масштабной функции свободной энергии в переменных и T разработан метод построения фундаментального единого неаналитического уравнения состояния, выполняющего все требования, которым удовлетворяют уравнения состояния вириального типа. Разработана новая кроссоверная функция, обеспечивающая описание термических и калорических свойств в более широкой области параметров состояния без потери точности описания экспериментальной информации.

9. Масштабная функция с новыми универсальными значениями параметров апробирована на примере построения фундаментального асимметричного единого уравнения состояния аргона. Показано, что предложенное УС ни по одному показателю не уступает известным уравнениям состояния Ar, а по некоторым превосходит. Например, оно лучше, чем фундаментальное асимметричное УС Кудрявцевой, описывает поведение второго вириального коэффициента в области низких температур. Это позволяет сделать вывод о том, что установлены оптимальные значения параметров асимптотической масштабной функции свободной энергии, и она может быть рекомендована для разработки сингулярных составляющих не только масштабных, но и широкодиапазонных уравнений состояния.

10. Метод построения фундаментального единого неаналитического УС, удовлетворяющего МТ, применен к описанию равновесных свойств перфторпропана. Установлено, что предложенное УС позволяет описать с малой погрешностью калорические и p T -данные и в околокритической области, и в регулярной части термодинамической поверхности. Важным обстоятельством является то, что все нелинейные параметры предложенного уравнения состояния являются универсальными с точностью до универсальности критических индексов и значения масштабной переменной (как классической, так и обобщенной) на линии насыщения.

11. На основе фундаментального неаналитического УС данной работы рассчитаны таблицы равновесных свойств R218 (в диапазоне 130 K T 500 K и 0,0001 МПа p 70 МПа) как в регулярной части термической поверхности, так и в окрестности критической точки. Таблицы включают, h, s, Cv, C p, w, 12. На основе полученных результатов можно сделать вывод, что разработанный метод построения фундаментальных масштабных и единых УС, удовлетворяющих МТ, в переменных плотность-температура обоснован в рамках феноменологической теории Мигдала, и фундаментальные единые уравнения состояния, с разработанными на его основе масштабными функциями, передают регулярную и нерегулярную части термодинамической поверхности с погрешность, аналогичной аналитическим УС и масштабным уравнениям, соответственно.

1. Рыков А.В. Кудрявцева И.В., Рыков В.А. Ассиметричное масштабное уравнение состояния R23 // Вестник Международной академии холода. 2012. № 4.

С. 26–28.

2. Рыков С.В., Рыков А.В. Метод построения уравнения состояния, учитывающего особенности критической области // Сборник трудов II Международной научно-технической конференции «Современные методы исследований теплофизических свойств веществ». СПб. 2012. – С. 253–258.

3. Кудрявцева И.В., Рыков А.В., Рыков В.А. Модифицированное уравнение линии насыщения, удовлетворяющее требованиям масштабной теории [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. – № 2.

4. Кудрявцева И.В., Рыков А.В. Uniform nonanalytic equation of state and thermodynamic tables R218 // XIX International Conference on Chemical Thermodynamics in Russia (RCCT-2013). – М.: MITHT Publisher, 2013. – 468 р.

5. Кудрявцева И.В., Рыков А.В. Линия фазового равновесия и новое уравнение для «кажущейся» теплоты парообразования // Сборник тезисов IX International scientific conference «Modern Problems of Refrigeration Equipment and Technology». Одеса. 2013. С. 229–231.

6. Кудрявцева И.В., Рыков А.В., Рыков В.А. Метод расчета равновесных свойств жидкости и газа в широкой области параметров состояния, включая критическую // Сборник тезисов IX International scientific conference «Modern Problems of Refrigeration Equipment and Technology». М. – 2013. – С. 223–225.

7. Кудрявцева И.В., Рыков А.В., Рыков В.А. Метод расчета равновесных свойств сверхкритических флюидов, используемых в СКФ-технологиях [Электронный ресурс]: Электронный научный журнал «Процессы и аппараты пищевых производств». НИУ ИТМО, 2013. № 2.

8. Кудрявцева И.В., Рыков А.В., Рыков В.А. Модифицированное уравнение линии насыщения, удовлетворяющее требованиям масштабной теории [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

9. Кудрявцева И.В., Рыков А.В., Рыков В.А. Непараметрическое уравнение состояния скейлингового вида и метод псевдокритических точек [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

10. Кудрявцева И.В., Рыков А.В., Рыков В.А. Новое уравнение для «кажущейся» теплоты парообразования [Электронный ресурс]: Электронный научный журнал «Процессы и аппараты пищевых производств». НИУ ИТМО, 2013. № 2.

С. 30.

11. Кудрявцева И.В., Рыков А.В., Рыков В.А., Рыков С.В. Единое неаналитическое уравнение состояния перфторпропана, удовлетворяющее масштабной теории критических явлений // Вестник международной академии холода. 2013. № 3. С. 22–26.

12. Кудрявцева И.В., Рыков А.В., Рыков С.В. Анализ структуры непараметрического уравнения состояния скейлингового вида [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

13. Кудрявцева И.В., Рыков А.В., Рыков С.В. Непараметрическое уравнение состояния скейлингового вида и метод псевдокритических точек // Сборник тезисов IX International scientific conference «Modern Problems of Refrigeration Equipment and Technology». Одесса. 2013. С. 225–226.

14. Кудрявцева И.В., Рыков А.В., Рыков С.В. Описание метастабильной области непараметрическими уравнениями состояния скейлингового вида [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

15. Рыков А.В., Кудрявцева И.В., Рыков С.В. Непараметрическое масштабное уравнение состояния, не содержащее дифференциальных биномов [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

16. Рыков А.В., Кудрявцева И.В., Рыков С.В. Уравнение линии насыщения, удовлетворяющее модифицированному правилу криволинейного диаметра [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

17. Рыков А.В., Кудрявцев Д.А., Рыков В. А. Метод выбора масштабных функций в переменных плотность и температура // Научно-технический вестник Поволжья. 2013. № 5. С. 111–116.

18. Рыков А.В., Кудрявцева И.В., Рыков В.А. К вопросу описания термодинамической поверхности, включая критическую область, уравнениями состояния в физических переменных [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 1.

19. Рыков А.В., Кудрявцева И.В., Рыков С.В. Анализ экспериментальной информации о равновесных свойствах R218 на основе неаналитического уравнения состояния [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 1.

20. Рыков А.В., Кудрявцева И.В., Рыков С.В. Непараметрическое масштабное уравнение состояния, не содержащее дифференциальных биномов [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

21. Рыков А.В., Кудрявцева И.В., Рыков С.В. Уравнение линии насыщения, удовлетворяющее модифицированному правилу криволинейного диаметра [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 2.

22. Рыков А.В., Кудрявцева И.В., Рыков С.В. Уравнения линии насыщения и упругости хладона R218 // Вестник Международной академии холода. 2013. № 4. С. 54–57.

23. Рыков С.В., Кудрявцева И.В., Рыков А.В., Курова Л.В. Метод построения фундаментального уравнения состояния, учитывающего особенности критической области [Электронный ресурс]: Электронный научный журнал «Холодильная техника и кондиционирование». НИУ ИТМО, 2013. № 1.

НИУ ИТМО. 197001, Санкт-Петербург, Кронверкский пр., 49.

ИИК ИХиБТ. 191002, Санкт-Петербург, ул. Ломоносава, 9.