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

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

Повний опис

Збережено в:
Бібліографічні деталі
Дата:2015
Автори: Левтеров, А.М., Умеренкова, К.Р.
Формат: Стаття
Мова:Russian
Опубліковано: Інститут технічної теплофізики НАН України 2015
Назва видання:Промышленная теплотехника
Теми:
Онлайн доступ:http://dspace.nbuv.gov.ua/handle/123456789/142166
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии / А.М. Левтеров, К.Р. Умеренкова // Промышленная теплотехника. — 2015. — Т. 37, № 2. — С. 32-38. — Бібліогр.: 10 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-142166
record_format dspace
spelling irk-123456789-1421662018-09-30T01:23:04Z Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии Левтеров, А.М. Умеренкова, К.Р. Тепло- и массообменные аппараты В статье приведен оригинальный метод и результаты расчета теплофизических свойств азота, используемого в качестве рабочего тела для транспортных поршневых установок. У статті викладено оригінальний метод та результати розрахунку теплофізичних властивостей азоту, що використовується як робоче тіло для транспортних поршневих двигунів. In article the original method and results of calculation thermo physical properties of the nitrogen used as a working body for transport piston installations is resulted. 2015 Article Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии / А.М. Левтеров, К.Р. Умеренкова // Промышленная теплотехника. — 2015. — Т. 37, № 2. — С. 32-38. — Бібліогр.: 10 назв. — рос. 0204-3602 DOI: https://doi.org/10.31472/ihe.2.2015.04 http://dspace.nbuv.gov.ua/handle/123456789/142166 532+541.11 ru Промышленная теплотехника Інститут технічної теплофізики НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language Russian
topic Тепло- и массообменные аппараты
Тепло- и массообменные аппараты
spellingShingle Тепло- и массообменные аппараты
Тепло- и массообменные аппараты
Левтеров, А.М.
Умеренкова, К.Р.
Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
Промышленная теплотехника
description В статье приведен оригинальный метод и результаты расчета теплофизических свойств азота, используемого в качестве рабочего тела для транспортных поршневых установок.
format Article
author Левтеров, А.М.
Умеренкова, К.Р.
author_facet Левтеров, А.М.
Умеренкова, К.Р.
author_sort Левтеров, А.М.
title Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
title_short Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
title_full Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
title_fullStr Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
title_full_unstemmed Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии
title_sort расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. часть iii. вычисление энтальпии и энтропии
publisher Інститут технічної теплофізики НАН України
publishDate 2015
topic_facet Тепло- и массообменные аппараты
url http://dspace.nbuv.gov.ua/handle/123456789/142166
citation_txt Расчетная оценка теплофизических свойств азота, как рабочего тела поршневого криодвигателя. Часть III. Вычисление энтальпии и энтропии / А.М. Левтеров, К.Р. Умеренкова // Промышленная теплотехника. — 2015. — Т. 37, № 2. — С. 32-38. — Бібліогр.: 10 назв. — рос.
series Промышленная теплотехника
work_keys_str_mv AT levterovam rasčetnaâocenkateplofizičeskihsvojstvazotakakrabočegotelaporšnevogokriodvigatelâčastʹiiivyčislenieéntalʹpiiiéntropii
AT umerenkovakr rasčetnaâocenkateplofizičeskihsvojstvazotakakrabočegotelaporšnevogokriodvigatelâčastʹiiivyčislenieéntalʹpiiiéntropii
first_indexed 2025-07-10T14:19:43Z
last_indexed 2025-07-10T14:19:43Z
_version_ 1837269980178022400
fulltext ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №232 УДК 532+541.11 РАСЧЕТНАЯ ОЦЕНКА ТЕПЛОФИЗИЧЕСКИХ СВОЙСТВ АЗОТА, КАК РАБОЧЕГО ТЕЛА ПОРШНЕВОГО КРИОДВИГАТЕЛЯ. ЧАСТЬ III. ВЫЧИСЛЕНИЕ ЭНТАЛЬПИИ И ЭНТРОПИИ Левтеров А.М., Умеренкова К.Р. Институт проблем машиностроения им. А.Н. Подгорного НАН Украины, ул. Пожарского, 2/10, г. Харьков, 61046 В статье приведен оригиналь- ный метод и результаты расчета теплофизических свойств азота, используемого в качестве рабочего тела для транспортных поршневых установок. У статті викладено оригіналь- ний метод та результати розрахун- ку теплофізичних властивостей азо- ту, що використовується як робоче тіло для транспортних поршневих двигунів. In article the original method and results of calculation thermo physical properties of the nitrogen used as a working body for transport piston installations is resulted. ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ Библ. 10, табл. 2, рис. 2. Ключевые слова: поршневой криодвигатель, азот, теплофизические характеристики, рабочее тело, энтальпия, энтропия, математическая модель. Состояние проблемы. Пневматические си- ловые установки для автомобилей, в которых используется газ высокого давления, стали объ- ектом интенсивных исследований в последние несколько лет и рассматриваются как один из перспективных вариантов для автомобилей ма- лой и средней грузоподъемности. Сюда относят- ся и криогенные силовые установки, в которых используются низкотемпературные пневматиче- ские двигатели, а рабочим телом для них явля- ется жидкий азот или воздух, хранящийся в кри- остатах на борту автомобиля и подогреваемый до температуры окружающей среды в теплооб- менниках [1–3]. Отмечено, что сжиженный газ обеспечивает большую работу на единицу массы рабочего тела по сравнению с аналогичным по- D – массовая плотность; D – параметр приведения для массовой плотности; E – параметр приведения для температуры; F, f − свободная энергия; h –удельная энтропия; I – групповой интеграл; k – постоянная Больцмана; M – молярная масса; N – число частиц; NA – число Авогадро; p – давление; Р – параметр приведения для давления; R = R0/M – газовая постоянная; R0 – универсальная газовая постоянная; u – потенциал взаимодействия; T – температура; V – объем; z – фактор сжимаемости; β = 1/(kT); ε, σ – параметры потенциала межмолекулярного взаимодействия; μ – химический потенциал; η – параметр упаковки; ρ = N/V – плотность числа частиц. Индексы верхние: * – приведенные величины. Индексы нижние: 0 – нулевое приближение; 1 – первый порядок теории возмущений; 2 – второй порядок теории возмущений; st – стандартное; L – жидкость; S – насыщенный пар; V – пар; ДВС – двигатель внутреннего сгорания; РТ – рабочее тело; ХНАДУ – Харьковский национальный автомо- бильно-дорожный университет. ~ ~ ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №2 33 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ тавалась нерешенной. Существовали различного рода прогности- ческие алгоритмы или упрощенные модельные методы, обзор которых приведен в работе [7]. Данные методы расчетов обеспечивали прием- лемые результаты для отдельных термодинами- ческих характеристик в ограниченных областях состояний, однако в целом оказывались неудов- летворительными, поскольку для набора свойств погрешности их описания существенно превы- шали экспериментальные ошибки. Постановка задачи. Такое состояние про- блемы обусловило необходимость развития ста- тистико-механических методов описания свойств молекулярных веществ, в рамках которых ис- пользуется минимум исходных данных и пара- метров. Одним из важных условий успешного использования азота в качестве РТ является соз- дание современных инженерных методов рас- чета его термодинамических свойств. Эти ме- тоды расчета должны удовлетворять одновре- менно таким требованиям, как точность и воз- можность применения в широком диапазоне дав- лений и температур, а также доступность исполь- зования результатов в режиме текущего времени, т.е. в процессе управления режимами работы дви- гателя. Проводимые исследования посвящены применению оригинальной модифицированной схемы термодинамической теории возмущений [7] для описания энтальпии и энтропии азота. Метод и результаты расчета. Cвободная энергия F = N·f системы N частиц, взаимодействующих посредством исходного по- тенциала u(r), при T и ρ имеет вид βf = βf0 + ρ* (I1 + I2/T*)/T* + ..., (1) где ρ* = ρσ3. Величина f0 представляет собой удельную (на одну частицу) свободную энергию системы ТС. I1 и I2 – групповые интегралы пер- вого и второго порядков теории возмущений [8]. Термодинамические свойства азота получе- ны на базе рассчитанных согласно (1) поверх- ностей свободной энергии f(T, V) или F(T, ρ) с использованием стандартных термодинамичес- ких соотношений. Начальным этапом расчетов является определение плотности. Плотности жидкой DL и паровой DV фаз вещества на линии казателем сжатого газа. Жидкий азот получают на воздухораздели- тельных установках из атмосферного воздуха, где его запасы практически не ограничены. Кро- ме того, азот является экологически безопасным и химически инертным газом, что позволяет ис- пользовать криогенные силовые установки даже в закрытых помещениях повышенной пожаро- опасности, в том числе при работе с огнеопас- ными веществами. Первые экспериментальные автомобили с криодвигателями, использующимися в авиаци- онной промышленности уже созданы в США и по мнению специалистов имеют прекрасную перспективу в будущем. Представляют также интерес разработки малолитражных автомоби- лей французской компании Motor Development International. В Украине инженерно-конструктор- ские работы по разработке высокоэффективного автомобильного пневмодвигателя для криоген- ной силовой установки ведутся в ХНАДУ [3]. К преимуществам использования криогенных силовых установок в автомобилях относятся: - пожаробезопасность; - экологическая чистота; - возможность рекуперации энергии; - доступная технология производства; - доступность и относительно низкая стоимость рабочего тела (сжиженный или сжатый азот, воз- дух); - экономическая целесообразность (стоимость 1 км пробега меньше, чем для ДВС) [3, 4]; - возможность массового применения в городах с высокой плотностью населения и автотранс- порта [5]. В настоящее время предполагается разраба- тывать автомобили с пневматической (криоген- ной) силовой установкой для использования в густонаселенных районах городов в качестве маршрутных такси, небольших грузовых авто- мобилей, машин-холодильников, а также специ- альных машин для обслуживания аэропортов [6], шахт и взрывоопасных производств. Анализ исследований и публикаций. Про- блема создания эффективных методов расчета термодинамических свойств веществ, которые одновременно применимы для газообразного и жидкого состояний, до последнего времени ос- ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №234 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ насыщения определяются из условий равенства давлений и химических потенциалов фаз: (2) а плотность вещества D(p,T) [кг/м3] в однофаз- ной области в термодинамическом состоянии p [МПа] и T [K] – из уравнения: p(D,T) – p = 0. (3) Приведем выражения для основных свойств. Давление p = PT*ρ*z. Фактор сжимаемости z = 103p/RTD (4) Свободная энергия F = RT(βf) (5) Внутренняя энергия U = RT(βu) (6) Энтропия S = Rs (s = βu – βf) (7) Энтальпия H = RT(βh) (βh = βu + z) (8) Энергия Гиббса G = RT(βμ) (μ = f + z/β) (9) Здесь η = 0,4177ρ*; ρ* = D/D, T* = T/E, p* = p/P, D = M/(NAσ3), E = ε/k, P = RED. Величины вы- ражены в следующих единицах: R – кДж/(кг·K); p, P – МПа; D, D – кг/м3; T, E – K; F, G – кДж/кг. Для азота: Е = 97,55 K, σ = 3,5996·10-10 м [7], D = = 997,38 кг/м3, P = 8,877 МПа, R = 0,2968 кДж/ (кг·K). Формы χ1[n] и χ2[n] определяются согласно выражениям, приведенным в работе [7]. B выражениях (1) – (9) zN = 103pst/RTD – функ- ция нормировки на стандартное давление pst = = 0,101325 МПа, а величины h0, cp0, s0 – удельные (приходящиеся на одну молекулу) энтальпия, энтропия и изобарная теплоемкость в идеально- газовом состоянии. Для определения h0, cp0, s0 использованы вы- ражения Пассата-Даннера [9]. Энтальпия и энтропия азота рассчитаны в однофазной области и вдоль линий кипения и конденсации. В табл. 1 приведены расчетные значения плотности, энтальпии и энтропии для различных значений давлений и температур. При сравнении с экспериментальными данными [10] к расчет- ным значениям Н добавлена теплота сублимации азота h0 = 247,6 кДж/кг при Т = 0 K, поскольку для идеально-газовой энтальпии, используемой в [10], принята ненулевая константа отсчета. Отклонения расчетных значений Н и S вдоль линий кипения и конденсации от эксперимен- тальных данных [10] (δ = (Нрас/Нэксп) – 1; δ = = (Sрас/Sэксп) – 1) и средние значения модулей этих отклонений показывают, что при- веденные результаты расчетов удовлетворитель- но согласуются с опытными данными. Выводы Анализ полученных результатов показыва- ет возможность успешного применения строгих статистико-механических методов к описанию свойств молекулярных систем в различных фа- зовых состояниях. При использовании предло- женного метода вполне достаточно минималь- ного набора исходных данных (давление, темпе- ратура). Рассчитанные термодинамические ха- рактеристики жидкого и газообразного азота могут быть использованы для разработки пер- спективных видов экологически чистых и пожа- робезопасных транспортных средств с криоген- ными силовыми установками. ~ ~ ~ ~ ~ ~ ~ ~ 0 ~ ~ ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №2 35 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ Табл. 1. Сравнение расчетных и экспериментальных значений плотности, энтальпии и энтропии Т, К D, кг/м3 D [10], кг/м3 H, кДж/кг H [10], кДж/кг S, кДж/(кг·K) S [10], кДж/(кг K) p = 0,1 МПа 70 838,45 837,86 126,36 110,5 2,757 2,628 80 4,353 4,37 328,76 328,3 5,463 5,457 90 3,833 3,85 339,56 338,8 5,59 5,58 100 3,429 3,44 350,19 349,5 5,703 5,693 110 3,104 3,11 360,76 360,2 5,803 5,795 120 2,837 2,84 371,28 370,8 5.895 5,887 200 1,6882 1,69 454,87 454,5 6,429 6,422 250 1.34885 1,35 506,84 506,7 6,661 6,655 300 1,1234 1,12 558,77 558,8 6,85 6,844 350 0,9626 0,96 610,77 610,9 7,011 7,005 p = 5 МПа 70 849,657 847,64 130,28 114,4 2,73 2,60 80 807,857 806,29 148,1 134,8 2,968 2,873 90 761,098 760,87 167,86 155,1 3,2 3,112 100 710,637 711,62 188,69 175,7 3,42 3,329 110 655,410 654,93 210,28 197,6 3,625 3,538 120 591,739 582,07 233,14 222,6 3,824 3,755 200 93,677 93,42 432,42 429,8 5,184 5,168 250 69,877 69,58 492,15 491,1 5,451 5,442 300 56,587 56,35 548,53 548,2 5,657 5,65 350 47,8325 47,67 603,47 603,6 5,827 5,821 p = 10 МПа 70 860,083 856,49 134,36 118,5 2,705 2,575 80 820,748 818,06 154,63 138,4 2,935 2,841 90 777,58 776,24 170,84 158,2 3,161 3,074 100 732,34 732,53 190,9 178,0 3,372 3,282 110 685,119 685,75 211,25 198,3 3,566 3,476 120 635,141 633,03 231,8 219,8 3,745 3,663 200 207,29 199,54 408,64 405,8 4,886 4,869 250 142,62 140,59 478,13 477,1 5,198 5,188 300 112,91 111,69 539,13 539,1 5,421 5,414 350 94,48 93,7 596,92 597,6 5,599 5,596 ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №236 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ Табл. 2. Средние значения модулей отклонений δ энтальпии и энтропии азота Линия кипения Линия конденсации H S H S δ 0,078 0,028 0,031 0,028 Рис. 2. Энтальпия азота на линии кипения и конденсации. Экспериментальные данные [10]. 2 3 4 5 6 60 80 100 120 140 T, K S , кДж/(кг∙К) энтропия на линии конденсации - расчет энтропия на линии конденсации - эксперимент энтропия на линии кипения - расчет энтропия на линии кипения - эксперимент 100 150 200 250 300 350 60 80 100 120 140T, K Н, кДж/кг энтальпия на линии конденсации - расчет энтальпия на линии конденсации - эксперимент энтальпия на линии кипения - расчет энтальпия на линии кипения - эксперимент Рис. 1. Энтропия азота на линии кипения и конденсации. Экспериментальные данные [10]. ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №2 37 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ ЛИТЕРАТУРА 1. Plummer M.C. Cryogenic heat engine experiment / M.C. Plummer, C.P. Koehler, D.R. Flanders // Proc. of 1997 Cryogenic Eng. Conf., Portland, July 1997, USA. – P. 7. 2. Williams J. Frost-Free Cryogenic Heat Exchanger for Automotive Propulsion / J. Williams, C. Knovlen, A.T. Mattick, A. Hertzberg // Proc. of 33rd AIAA/ASEE Joint Propulsion Conf. @ Exhibit, July 6-9, 1997, Seattle, USA. – P. 18. 3. Туренко А.Н. Экологически чистый крио- генный транспорт: современное состояние про- блемы /А.Н. Туренко, А.И. Пятак, И.Н. Кудрявцев // Вестн. Харьк. автомоб.-дорожн. техн. ун-та. – 2000. – Вып. 12-13. – С. 42–47. 4. Сухов А.П. Выхлоп чище воздуха // За рулем. – 2001.– № 2. – С. 40–42. 5. Бондаренко С.И. Криоавтомобиль: будем ез- дить на азоте? / С.И. Бондаренко, И.Н. Кудрявцев // Автоцентр. – 2000.– № 40. – С. 32–33. 6. Бондаренко С.И. Новый вид пожаробезо- пасного и экологически чистого транспортного средства для аэропортов / С.И. Бондаренко, И.Н. Кудрявцев, А.И. Пятак, И.И. Тимченко, А.П. Ку- дряш // Пробл. машиностроения. – 2002. – Т. 5. – № 7. – С. 92–95. 7. Маринин В.С. Теплофизика альтернативных энергоносителей. – Харьков: Форт. – 1999. – 212 с. 8. Левтеров А.М. Расчетная оценка теплофи- зических свойств азота, как рабочего тела порш- невого криодвигателя. Часть I. Математическая модель фазовых равновесий / А.М. Левтеров, К.Р. Умеренкова // Пром. теплотехника. – 2013. – Т. 35, № 4. – С. 90 – 95. 9. Левтеров А.М. Расчетная оценка теплофи- зических свойств азота, как рабочего тела порш- невого криодвигателя. Часть II. Вычисление те- плоемкости / А.М. Левтеров, К.Р. Умеренкова // Пром. теплотехника. – 2014. – Т. 36, № 2. – С. 93 – 100. 10. Термодинамические свойства азота / В.В. Сычев, А.А. Вассерман, А.Д. Козлов. – М.: Издат. стандартов. – 1977. – 720 с. CALCULATION ASSESSMENT OF THERMAL PROPERTIES OF NITROGEN AS THE BODY WORK PISTON CRYOGENIC ENGINE. PART III. THE CALCULATION OF ENTHALPY AND ENTROPY A. Levterov, K. Umerenkova The A. N. Podgorny Institute for Mechanical Engineering Problems of the National Academy of Sciences of Ukraine, Kharkiv, Ukraine, ul. Pozharsky, 2/10, Kharkov, 61046, Ukraine In article theoretical model for description of the thermodynamic properties of liquid and gaseous nitrogen which used as a working body for transport piston power plants is proposed. Mathematical model for description of heat capacity is based on the unified statistical mechanical approach within the thermodynamic perturbation theory without invoking of any empirical parameters. The numerical calculations which accomplished for enthalpy and entropy of the molecular nitrogen agree with the available experimental data. Elaborated calculation procedures are used at creation of new kind of ecological and fireproof vehicle. References 10, tables 2, figure 2. Key words: reciprocating cryogenic engine, nitrogen, thermal characteristics, working body, enthalpy, entropy, mathematical model. 1. Plummer M.C. Cryogenic heat engine expe- riment / M.C. Plummer, C.P. Koehler, D.R. Flanders ISSN 0204-3602. Пром. теплотехника, 2015, т. 37, №238 ТЕПЛО- И МАССООБМЕННЫЕ АППАРАТЫ // Proc. of 1997 Cryogenic Eng. Conf., Portland, July 1997, USA. – P. 7. 2. Williams J. Frost-Free Cryogenic Heat Exchan- ger for Automotive Propulsion / J. Williams, C. Knovlen, A.T. Mattick, A. Hertzberg // Proc. of 33rd AIAA/ASEE Joint Propulsion Conf. @ Exhibit, July 6-9, 1997, Seattle, USA. – P.18. 3. Turenko A.N. Eco-friendly cryogenic trans- port: the current state of the problem / A.N. Tu- renko, A.I. Pjatak, I.N. Kudryavtsev // Vestnik Khar- kov Avtomobilno-doroznogo tekhnicheskogo uni- versiteta – 2000. – Issue. 12–13. – P. 42–47. (Rus.). 4. Suhov A.P. Vyhlop is cleaner than air // Za rulom. – 2001. – № 2. – P. 40–42. (Rus.). 5. Bondarenko S.I. Krioavtomobil: we shall go by nitrogen? / S.I. Bondarenko, I.N. Kudryavtsev // Avto centr. – 2000. – № 40. – P. 32–33. (Rus.). 6. Bondarenko S.I. New a kind of a fireproof and non-polluting vehicle for the airports / S.I. Bonda- renko, I.N. Kudryavtsev, A.I. Pjatak, I.I. Timchenko, A.P. Kudrjash // Problem mashinostroeniya. – 2002. – V. 5. – № 7. – P. 92–95. (Rus.). 7. Marinin V.S. Thermal physics of alternative energy sources. – Kharkov: Fort. – 1999. – 212 p. (Rus.). 8. Levterov A.M. Calculation assessment of thermal properties of nitrogen as the body work piston cryogenic engine. Part I. Mathematical model of phase equilibrium / A.M. Levterov, K.R. Umerenkova // Promyshlennaya teplotekhnika. – 2013. –V. 35, № 4. – P. 90–95. (Rus.). 9. Levterov A.M. Astimation of thermophysical properties of nitrogen as the working body for piston cryo-engine. Part II. The calculation of heat capacity / A.M. Levterov, K.R. Umerenkova // Prom. teplotekhnika. – 2014. –V. 36, № 2. – P. 93–100. (Rus.). 10. Thermodynamic Properties of Nitrogen / V.V. Sychev, A.A. Wasserman, A.D. Kozlov etc. – M: Izdatelstvo Standartov. – 1977. – 350 p. (Rus.) Получено 22.12.2014 Received 22.12.2014