CC BY
115
УДК 621.512/.513:536.7.001.57 ББК 31.762:22.317.24в631
Е. Н. Туголуков, Е. С. Егоров
МЕТОДИКА МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ ТЕРМОДИНАМИЧЕСКИХ ПРОЦЕССОВ ПОРШНЕВОГО КОМПРЕССОРА
E. N. Tugolukov, E. S. Egorov
METHODS OF MATHEMATICAL MODELING OF THERMODYNAMIC PROCESSES OF PISTON COMPRESSOR
Предложен подход к математическому моделированию термодинамических процессов сжатия газа или пара в цилиндре поршневого компрессора с учетом параметров реальных газов. Отличительной особенностью предлагаемого подхода является отказ от уравнений состояния, которые используются в традиционных моделях. Эти уравнения заменяются значениями, полученными интерполяцией экспериментальных данных для отдельных веществ, из которых состоит рабочее тело. Предложена математическая модель для определения термодинамических параметров рабочего тела в цилиндре поршневого компрессора с учетом теплообмена с конструкционными элементами.
Ключевые слова: математическое моделирование, параметры реальных газов, поршневой компрессор, фазовые переходы.
An approach to the mathematical modeling of thermodynamic processes of compression of gas or vapor in the cylinder of piston compressor with use of parameters of real gases is presented.
A distinctive feature of this approach is an abandonment of equations of state, which are used in traditional models. Theses equations are replaced with the values obtained by interpolation of the experimental data for the individual substances of which the operational body is made. The mathematical model to determine the thermodynamic parameters of the operational body in the cylinder of piston compressor based on heat exchange with the structural members is developed.
Key words: mathematical modeling, parameters of real gases, piston compressor, phase transitions.
Введение
Данное исследование посвящено разработке методики моделирования термодинамического состояния рабочего тела в цилиндре поршневого компрессора. Процесс сжатия рабочего тела в поршневом компрессоре в реальных условиях сопровождается целым комплексом взаимосвязанных явлений. К ним относятся: изменение объема, давления и температуры рабочего тела; изменение массы рабочего тела вследствие протечек через уплотнения; изменение фазового состава и массы газовой фазы при сжатии влажных паров; влияние тепловых эффектов фазовых переходов; теплообмен рабочего тела с конструкционными элементами.
Данная проблема имеет большое практическое значение и неоднократно рассматривалась разными авторами. В зависимости от учитываемых параметров предложено множество разнообразных моделей, основанных на использовании подходящих для каждого конкретного случая уравнений состояния. Это уравнения состояния для идеальных газов или вириальные уравнения, которые, как правило, дают хорошие результаты в узкой области изменения параметров рабочего тела и практически непригодны для расчетов в областях с фазовыми переходами. При описании поведения газа в областях с фазовыми переходами или в окрестностях критических точек использование таких моделей приводит к систематическому накоплению ошибок вычисления и сводит на нет их применение для реальных расчетов. Именно поэтому решение данной проблемы остается одной из самых трудных и актуальных задач термодинамики.
Отличительной особенностью предлагаемого подхода является отказ от уравнений состояния, использующихся в традиционных моделях, и их замена на значения, полученные интерполяцией экспериментальных данных для отдельных веществ, из которых состоит рабочее тело. Таким образом, при моделировании предлагается использовать только основные законы термодинамики, закон сохранения энергии и экспериментальные данные.
В таком подходе есть свои плюсы и минусы. Минус в том, что мы отказываемся от хорошо известных уравнений состояния, описывающих поведение газов и жидкостей для большого количества веществ, и тем самым делаем шаг назад. Очевидно, что это приведет к излишней работе по поиску и обработке данных по термодинамическим свойствам веществ в различных агрегатных состояниях.
Но, с другой стороны, при таком подходе: во-первых, увеличивается точность вычислений (т. к. вычислительные процедуры носят, как правило, итерационный характер, что приводит к накоплению ошибок вычислений при использовании классических уравнений состояния), а во-вторых, порой это единственный путь получить адекватное описание процессов с фазовыми переходами первого или второго рода с использованием параметров реальных газов.
К настоящему времени накоплен огромный экспериментальный материал по свойствам веществ, представленный во многих справочниках и их электронных аналогах. Это, например, справочник Н. Б. Варгафтика, выдержавший множество изданий [1, 2]; электронные базы данных по теплофизическим свойствам веществ [3-5]. Современный уровень развития вычислительной техники позволяет получать прямой доступ к ним посредством использования облачных серверов [6]. Это позволяет достаточно легко преодолеть указанные выше недостатки.
Методика моделирования термодинамических процессов с использованием параметров реальных газов
Предлагаемая методика включает следующие этапы:
1. Поиск экспериментальных данных.
Исходя из поставленных условий задачи, для каждого отдельного вещества, которое участвует в рассматриваемом процессе, производим поиск табличных зависимостей параметров состояния.
2. Обработка табличных зависимостей.
Далее эти таблицы аппроксимируются аналитическими функциями. Если изначально удается оценить диапазон значений для температуры и давления, в котором будет решаться поставленная задача, то табличные данные берутся из этого диапазона с некоторым запасом. Это сильно упрощает таблицы, т. к. в большинстве справочников диапазон измеренных интервалов большой: например, давление принимает значения от 0,01 до 1000 атм, а температура - от 100 до 3000 К.
Как известно, при фазовых переходах некоторые теплофизические характеристики, такие как плотность, вязкость, теплоемкость, претерпевают значительные скачкообразные изменения. В связи с этим возникает сложность аппроксимации данных характеристик и адекватного использования значений в математической модели, и поэтому для каждой задачи выбирается наиболее подходящий алгоритм аппроксимации (линейная, кубическая, сплайн-интерполяция), позволяющий без каких-либо существенных искажений использовать значения в данных областях.
3. Составление математической модели.
Математическая модель включает полученные аппроксимационные зависимости, основные уравнения термодинамики и законы сохранения энергии. Таким образом, уравнения состояния заменяются аппроксимационными зависимостями табличных термодинамических параметров. В зависимости от конкретных условий в математической модели могут присутствовать дополнительные уравнения, учитывающие такие факторы, как теплообмен с конструкционными элементами, утечки рабочего тела через уплотнения, фазовые переходы, химические реакции и т. д.
4. Решение уравнений модели.
Решение уравнений разработанной математической модели осуществляется итерационным методом, в результате чего определяются искомые параметры рабочего тела и энергетические характеристики текущих процессов.
Рассмотрим этапы 1-4 более подробно на примере задачи определения термодинамических параметров водяного пара в цилиндре поршневого компрессора.
Этап 1. Поиск экспериментальных данных
На этом этапе (в соответствии с предложенной выше схемой) в качестве примера выбираем таблицы для теплоемкости воды и водяного пара [1] при различных значениях температуры и давления. Пример такой таблицы приведен на рис. 1.
р, бар 1, °с 0,1 1 10 20 40 60 80 100
О 4,218 4,217 4,212 4,207 4,196 4,186 4,176 4,165
50 1,929 4,181 4,179 4,176 4,172 4,167 4,163 4,158
ЮО 1,910 2,038 4,214 4,211 4,207 4,202 4,198 4,194
120 1,913 2,007 4,243 4,240 4,235 4,230 4,226 4,221
140 1,918 1,984 4,283 4,280 4,275 4,269 4,263 4,258
160 1,926 1,977 4,337 4,334 4,327 4,320 4,313 4,307
180 1,933 1,974 2,613 4,403 4,395 4,386 4,378 4,370
200 1,944 1,975 2,433 4,494 4,483 4,472 4,461 4,450
220 1,954 1,979 2,316 2,939 4,601 4,586 4,571 4,557
240 1,964 1,985 2,242 2,674 4,763 4,741 4,720 4,700
260 1,976 1,993 2,194 2,505 3,582 4,964 4,932 4,902
280 1,987 2,001 2,163 2,395 3,116 4,514 5,25 5,20
300 1,999 2,010 2,141 2,321 2,834 3,679 5,31 5,70
320 2,011 2,021 2,126 2,268 2,649 3,217 4,118 5,79
340 2,024 2,032 2,122 2,239 2,536 2,943 3,526 4,412
350 2,030 2,038 2,125 2,235 2,504 2,861 3,350 4,043
360 2,037 2,044 2,127 2,231 2,478 2,793 3,216 3,769
365 2,040 2,048 2,128 2,227 2,462 2,759 3,134 3,655
370 2,043 2,050 2,128 2,222 2,446 2,725 3,072 3,546
375 ■ 2,046 2,053 2,127 2,218 2,428 2,690 3,018 3,446
380 2,049 2,056 2,127 2,212 2,412 2,657 2,964 3,356
385 2,052 2,059 2,126 2,207 2,396 2,627 2,913 3,274
390 2,056 2,061 2,125 2,202 2,381 2,600 2,867 3,201
395 2,059 2,065 2,125 2,200 2,369 2,575 2,826 3,137
400 2,062 2,068 2,126 2,197 2,358 2,553 2,789 3,078
405 2,066 2,071 2,127 2,195 2,349 2,534 2,756 3,025
410 2,069 2,074 2,128 2,193 2,340 2,517 2,727 ' 2,979
415 2,072 2,077 2,129 2,192 2,334 2,501 2,700 2,936
Рис. 1. Пример таблицы для теплоемкости воды и водяного пара (до обработки)
Так как данные получены экспериментально, то значения представляют собой дискретный набор значений теплоемкости в зависимости от термодинамических параметров состояния.
Этап 2. Обработка табличных зависимостей
На этом этапе обрабатываем исходную таблицу и решаем вопрос, каким способом лучше всего произвести вычисление для промежуточных значений. Спектр возможностей достаточно велик - это разнообразные подходы по интерполяции и аппроксимации.
Сложность в обработке данных таблиц состоит в том, что при фазовом переходе происходит существенный скачок значений теплофизических характеристик: так, например, значения теплоемкости для воды и пара (рис. 2) в непосредственной близости от линии фазового перехода отличаются в несколько раз.
Р, бар Т С 0,1 1 10 20 40 60 80 100
0 ' '■ 11 Д'1 7 12 1 2 1,20 4 4 66 4 , ! 10
50 00 ’Г. т ! 7 7-. 470 2 4 72 2 4 67 ;.1ё?
100 "! .91 2,038 | 4 2 14 И ,1 1 П‘7 2 2:02 13 ;,1 ул 1
120 1,0 ■! 3 ,00 7 4,24 4 73
140 1,9; а ; МИЛ 4,2 я 4 7? 6:3
160 1,023 , -д77 :■ 33 2 :22
100 1,933 2,613 4,403 ! .о у 4..386
200 1 ,044 18 4 .4 7 2
220 47154 ; /Я 7 4 7222; .30 15 8 Й 4.57 1 2 ,6 61
240 1 .444 22424 '6 4.-4 5 = .. 4 7
260 ■: ,7 7 0 4ЙЗ 2,1 04 2,60 2 7: ; 2 4,9 64 33 4 ,302
280 “,777? ■:_»_] : ; ! р. 2 7:4 4 3,11 25 3.2
300 1 4! 4 9 4 4!! 7 4 ~П 7 4 г ,63 4 22 7 : 7: 2
320 И 2 л=4; 2, \ 24 7. 'Л 7 2.74 0 27 2 ) 7 4,11В ■ :
340 2.024 2.032 2.122 2.239 2,536 2,943 3,526 4,412
350 2.03 2.038 2.125 2,235 2,504 2,861 3,35 4,043
360 2,037 2.044 2,127 2,231 2,478 2,793 3,216 3,769
365 2,04 2.048 2.128 2,227 2,462 2,759 3,134 3,655
370 2.043 2,05 2.128 2,222 2,446 2,725 3,072 3,546
375 2.046 2,053 2.127 2,218 2,428 2,69 3,018 3,446
380 2,049 2,056 2,127 2,212 2,412 2,657 2,964 3,356
385 2.052 2,059 2,126 2,207 2,396 2,627 2,913 3,274
390 2,056 2,061 2.125 2,202 2,381 2,6 2,867 3,201
395 2.059 2,065 2,125 2,2 2,369 2,575 2,826 3,137
400 2,062 2,068 2.126 2,197 2,358 2,553 2,789 3,078
405 2.066 2,071 2,127 2,195 2,349 2,534 2,756 3,025
410 2.069 2,074 2.128 2,193 2,34 2,517 2,727 2,979
415 2.072 2,077 2.129 2,192 2,334 2,501 2,7 2,936
Рис. 2. Пример таблицы для теплоемкости воды и водяного пара (после обработки)
Разобьем исходную таблицу на 2 части по линии фазового перехода. Принцип разбиения виден на рис. 2. Так, область, выделенная светлым серым, включает фазовый переход, а не выделенная область содержит параметры вещества только в одном фазовом состоянии. Обозначим эти зоны как I (светлый серый) и II. Для упрощения вычислений целесообразно для зоны II использовать метод простой линейной интерполяции. Изменение величин происходит равномерно, и данный метод дает высокую точность при вычислениях. Для зоны I линейная интерполяция неприемлема в связи с наличием разрыва первого рода. В результате исследований было установлено, что в данном случае целесообразно применение сплайн-интерполяции. На рис. 3 приведены примеры построения поверхности для зоны I и зоны II.
Для
данных
области
V VI
Рис. 3. Пример обработки экспериментальных данных
Этап 3. Составление математической модели
На этапе 3 составляем математическую модель для определения термодинамических параметров рабочего тела в цилиндре поршневого компрессора. Рассмотрим этот этап более детально.
При движении поршня изменяется объем рабочего пространства, масса рабочего тела в рабочем и дополнительном пространствах, фазовый состав рабочего тела, температура газовой и жидкой фракций, давление.
Кроме того, при наличии фазового перехода учитывается теплота фазового перехода и изменение количества газообразной фракции. В связи с большим количеством факторов, подлежащих учету, данная задача была разбита на подзадачи, что позволило, последовательно усложняя систему уравнений, составить наиболее полную модель рассматриваемой задачи.
Для рабочего пространства в области насыщенного пара запишем следующую систему уравнений.
Работа сжатия газа:
Лр = PДV. (1)
Изменение теплосодержания газа:
Qg = С (Т + АТ,Р + АР) (Т + АТ)(ш& -ш/ -ши)-С (Т,Р) Тш& . (2)
Изменение теплосодержания жидкой фракции:
Ql = С, (Т + ЛТ,Р + ЛР) (Т+ЛТ)(ш1 + шу -шии)-С, (Т,Р) Тш,. (3)
Теплота фазового перехода:
(4)
Уравнение баланса:
Ар - Qf + Qg + й = 0 •
(5)
Условие насыщения:
Т + АТ =Тнас (Р+АР) .
(6)
При решении данной задачи только с учетом того, что пар насыщенный, сложности не возникает, и нет необходимости в использовании уравнения состояния для идеального газа. Зависимости давления от температуры берутся непосредственно из экспериментальных данных. Если начальное состояние газа не соответствует состоянию насыщения, то необходимо добавление дополнительной зависимости. Это либо добавление уравнения Клайперона - Менделеева или его аналога, или его замена на полученные экспериментальным путем табличные зависимости удельного объема от давления и температуры.
Дополнительное условие (уравнение Клайперона - Менделеева)
Здесь Р, V, Т - соответственно текущие давление, объем и температура; А Р, А Т - соответственно изменения давления и температуры, вызванные изменением объема А V; шг, ш^ ши, ш, ши - соответственно текущие массы газовой фракции; рабочего тела, совершившего фазовый переход; газообразного рабочего тела, вытесненного в дополнительное пространство; жидкой фракции; жидкой фракции, вытесненной в дополнительное пространство; Су (Т, Р), С1 (Т, Р) - соответственно теплоемкость газовой и жидкой фракций при температуре Т
и давлении Р; г (Т) - удельная теплота фазового перехода как функция температуры; г - сжимаемость; Rc = R / ц ; ц - молекулярная масса рабочего тела.
Размерность всех величин дана в системе СИ.
Следующий этап составления системы уравнений - это описание рабочего пространства вне области насыщения. Для этого используются те же соотношения (1)-(7), в которых величины Qfi шш,, ш,и, отвечающие за фазовые переходы, равны нулю.
Работа нагнетания:
(7)
(8)
Изменение теплосодержания газа:
(9)
Изменение теплосодержания жидкой фракции:
(10)
Теплота фазового перехода:
(11)
Теплота с нагнетаемым рабочим телом:
&п = С, (Т,Р+АР)Т ши + С, (Т,Р+АР)Т шы
(12)
Уравнение баланса:
Ара - Qfa ~ + &іа = 0
(13)
Условие насыщения:
Т +АТ + Т (Р + АР).
а а нас V '
(14)
Дополнительное условие:
(Р + АР )У&а = 2аШаЯс (Та +АТа )•
(15)
Здесь индекс а соответствует дополнительному объему; Тг - температура рабочего тела на выходе из дополнительного объема.
Следующий этап позволяет учитывать изменения температуры конструкционных элементов в области насыщенного пара.
Для учета влияния теплообмена с конструкционными элементами запишем уравнения для процессов сжатия и расширения рабочего тела в области насыщенного пара с использованием энтальпии процесса.
Процесс расширения (процесс сжатия моделируется аналогично, меняются лишь знаки перед соответствующими компонентами уравнений).
Работа расширения газа:
Таким образом, предложенная модель (1)-(22) учитывает все основные факторы данного процесса.
Этап 4. Решение уравнений модели
На последнем этапе - этапе 4 переходим к решению уравнений составленной модели. Фрагмент блок-схемы решения уравнений математической модели, относящийся к процессу сжатия в цилиндре поршневого компрессора, представлен на рис. 4. Здесь, в зависимости от начальных условий рассматриваемого процесса, задаются начальные данные (блок 2), зависящие от формы рабочего пространства, рабочего тела, характеристик процесса. В первую очередь это такие величины, как давление, температура, начальный объем, занимаемый рабочим телом, состав фаз (количество рабочего тела, находящегося в жидкой и газообразной форме), величина перемещения поршня А х.
В зависимости от используемого рабочего тела, для заданных веществ ищутся следующие таблицы (блок 2): таблица зависимости давления от температуры и объема, зависимость удельной теплоемкости рабочего тела. В случае, когда присутствует фазовое равновесие, использует-
(16)
Изменение теплосодержания газовой фракции:
(17)
Изменение теплосодержания жидкой фракции:
= 11(Т-АТ) (т, + ту)-11 (Т) т,.
(18)
Теплота, отданная конструкционными элементами за время А т:
& = аК (*п (т)-^ )Ат
(19)
(20)
Условие насыщения:
Т -АТ =Тнас (Р-АР).
(21)
Дополнительное условие:
V + А^ )нас = ^(Т -АТ )•
(22)
ся таблица зависимости температуры от давления для насыщенного пара. В зависимости от особенностей процесса опционально используются также таблицы значений вязкости, плотности, коэффициента сжимаемости рабочего тела как функции основных термодинамических величин.
1
Началоработы программы I
О
2
Ввод данных
Ро, То, Го, Ак
3
Аппроксимация
табличных
данных
___ 4
Расчет теплоемкости Со = К(Ро, То)
5
Цикл по і =1 ... N.
N = И Ак
V, = Го + і Ак S
Расчет Ті и Рі
Расчет
теплоемкости
С, = К(Рі, Ті)
Давление, температура, начальный объем, занимаемый рабочим телом, шаг поршня
Т аблицы экспериментальных данных
9
Вывод результатов
Решение
системы
уравнений (23)-(27)
10 ' г
I Конец работы I программы
6
7
8
Рис. 4. Блок-схема расчета параметров системы на этапе сжатия рабочего тела в цилиндре поршневого компрессора
Далее (блок 3) эти табличные зависимости аппроксимируются аналитическими функциями. Если изначально удается оценить диапазон изменения температуры и давления в рассматриваемом процессе, то табличные данные берутся из этого диапазона с некоторым запасом. Это сильно упрощает таблицы, т. к. в большинстве справочников диапазон измеренных интервалов очень большой, например давление принимает значения от 0,01 до 1000 атм, а температура - от 100 до 3000.
На этапе ввода исходных данных определяется значение теплоемкости рабочего тела С0 как функции начальных значений температуры Т0 и давления Р0 (блок 4).
Далее в цикле (блоки 5-8) происходит вычисление текущих значений объема рабочего тела при смещении поршня на величину А h.
VI = V) + I А h S ^ - площадь поршня), где I = и А h, L - перемещение поршня; I - количество шагов.
Для текущего значения объема рабочего тела решается следующая система уравнений.
Работа сжатия газа:
Лр = РА V. (23)
Изменение теплосодержания газа:
Qg = С (Т+АТ,Р+АР) (Т+АТ)(ш!? -ш} -ши)-С (Т,Р) Тш?. (24)
Изменение теплосодержания жидкой фракции:
Ql = С, (Т + АТ,Р + АР) (Т + АТ)(ш, + шу -шии)-С, (Т,Р) Тш,. (25)
Теплота фазового перехода:
Q'=г+А^] ш>. (26)
Уравнение баланса:
-Лр -Qf + Qg + & = °. (27)
Используя аппроксимационную зависимость температуры и давления рабочего тела от удельного объема, вычисляем новые значения давления и температуры (блок 7):
Т = Т +АТ и Р = Р0 +АР .
Далее, на основе пересчитанных значений для температуры и давления, пересчитываем значения для теплоемкости в соответствии с табличными значениями (блок 8). Эти вычисления повторяются до окончания цикла по 1.
В результате определяем термодинамические параметры рабочего тела поршневого компрессора при любом положении поршня.
Сравнение результатов решения
Результаты математического моделирования процесса сжатия реального газа в поршневом компрессоре с использованием данного подхода и с использованием уравнений для идеального газа приведены в [7]. Анализ полученных результатов показывает, что в областях при низком давлении и без фазовых переходов расчетные кривые совпадают, но расходятся в областях с фазовыми переходами. Это объясняется тем, что уравнения идеального газа в данных областях менее точные, чем в предложенном нами подходе.
Заключение
Таким образом, разработана методика математического моделирования процесса сжатия реального газа или пара в цилиндре поршневого компрессора с учетом теплообмена с конструкционными элементами. Особенностью разработанной методики является замена уравнения состояния для рабочего тела на аппроксимационные зависимости соответствующих табличных данных, приводимых в справочной литературе.
СПИСОК ЛИТЕРАТУРЫ
1. Варгафтик Н. Б. Справочник по теплофизическим свойствам газов и жидкостей / Н. Б. Варгаф-тик. М.: Наука, 1972. 720 с.
2. Vargaftik N. B. Handbook of physical properties of liquids and gases. Pure substances and mixture / N. B. Vargaftik, Yu. K.Vinogradov, V. S.Yargin. New York: Begel House, 1996. 1358 p.
3. Lemmon E. W. REFPROP - Reference Fluid Thermodynamic and Transport Properties / E. W. Lemmon, M. L. Huber, M. O. McLinden // CD ROM of the Seventeenth Symp. on Thermophys. Prop. Boulder, Colorado, USA, 2009. 599 p.
4. Александров А. А. Web-версия справочника «Теплофизические свойства рабочих веществ теплоэнергетики» / А. А. Александров, К. А. Орлов, В. Ф. Очков [Электронный ресурс]: http://twt.mpei.ac.ru/rbtpp/index.html.
5. NIST Reference Fluid Thermodynamic and Transport Properties Database (REFPROP): Version 9.1 [Электронный ресурс]: http://www.nist.gov/srd/nist23.cfm.
6. Очков В. Ф. Облачный сервис по свойствам рабочих веществ для теплотехнических расчетов / В. Ф. Очков, К. А. Орлов, Л. М. Френкель, А. В. Очков, В. Е. Знаменский // Теплоэнергетика. 2012. № 7. С. 79-86.
7. Туголуков Е. Н. Подход к моделированию термодинамических процессов с использованием параметров реальных газов / Е. Н. Туголуков, Е. С. Егоров // Математические методы в технике и технологиях «ММТТ-25»: c6. тр. Междунар. науч. конф. Саратов, 2012. Т. 8. С. 106-108.
REFERENCES
1. Vargaftik N. B. Spravochnik po teplofizicheskim svoistvam gazov i zhidkostei [Reference on thermophysical parameters of gases and liquids]. Moscow, Nauka Publ., 1972. 720 p.
2. Vargaftik N. B., Vinogradov Yu. K., Yargin V. S. Handbook ofphysical properties of liquids and gases. Pure substances and mixture. New York, Begel House, 1996. 1358 p.
3. Lemmon E. W., Huber M. L., McLinden M. O. REFPROP - Reference Fluid Thermodynamic and Transport Properties. CD ROM of the Seventeenth Symp. on Thermophys. Prop. Boulder, Colorado, USA, 2009. 599 p.
4. Aleksandrov A. A., Orlov K. A., Ochkov V. F. Web-versiia spravochnika «Teplofizicheskie svoistva rabochikh veshchestv teploenergetiki» [Web-version of the reference "Thermo-physical properties of the working substances of power engineering"]. Available at: http://twt.mpei.ac.ru/rbtpp/index.html.
5. NIST Reference Fluid Thermodynamic and Transport Properties Database (REFPROP): Version 9.1. Available at: http://www.nist.gov/srd/nist23.cfm.
6. Ochkov V. F., Orlov K. A., Frenkel' L. M., Ochkov A. V., Znamenskii V. E. Oblachnyi servis po svoist-vam rabochikh veshchestv dlia teplotekhnicheskikh raschetov [Cloud service on the properties of the working substances for thermodynamic calculations]. Teploenergetika, 2012, no. 7, pp. 79-86.
7. Tugolukov E. N., Egorov E. S. Podkhod k modelirovaniiu termodinamicheskikh protsessov s is-pol'zovaniem parametrov real'nykh gazov [Approach to modeling of thermodynamic processes using parameters of real gases]. Matematicheskie metody v tekhnike i tekhnologiiakh «MMTT-25». Cbornik trudov Mezhdunarod-noi nauchnoi konferentsii. Saratov, 2012, vol. 8, pp. 106-108.
Статья поступила в редакцию 22.11.2013, в окончательном варианте - 6.12.2013
ИНФОРМАЦИЯ ОБ АВТОРАХ
Туголуков Евгений Николаевич — Тамбовский государственный технический университет; д-р техн. наук, профессор, профессор кафедры «Техника и технологии производства нанопродуктов»; [email protected].
Tugolukov Evgeniy Nickolaevich - Tambov State Technical University; Doctor of Technical Sciences, Professor; Professor of the Department "Techniques and Nanoproduction Technologies"; [email protected].
Егоров Евгений Сергеевич — Тамбовский государственный технический университет; аспирант кафедры «Автоматизированное проектирование технологического оборудования»; [email protected].
Egorov Evgeniy Sergeevich - Tambov State Technical University; Postgraduate Student of the Department "Automated Designing of Technological Equipment"; [email protected].
