Научная статья на тему 'Уравнение линии упругости хладагента R236еа'

Уравнение линии упругости хладагента R236еа Текст научной статьи по специальности «Физика»

CC BY
114
27
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ELASTIC LINE / УРАВНЕНИЕ СОСТОЯНИЯ / EQUATION OF STATE / ЛИНИЯ ФАЗОВОГО РАВНОВЕСИЯ / THE LINE OF PHASE EQUILIBRIUM SCALING THEORY / ХЛАДАГЕНТ R236EA / R236EA / HEXAFLUOROPROPANE / ЛИНИЯ УПРУГОСТИ / МАСШТАБНАЯ ТЕОРИЯ / ГЕКСАФТОРПРОПАН

Аннотация научной статьи по физике, автор научной работы — Полторацкий М. И., Рыков С. В., Свердлов А. В.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

The equation of elasticity R236EA

The paper proposed the equation of elasticity hexafluoropropane, taking into account the requirements of scale hypotheses and behaviors vapor pressure near the critical and triple points. A comparison with the experimental data on the pressure in the line of elasticity. Revealed good agreement with experiment in the vicinity of the critical point, in contrast to the equations of other authors. Information about the relative deviation is presented in graphical form On the basis of the proposed equation calculated detailed tables, depending on the temperature of the pressure on the line of phase equilibrium. Also in the pre-table shows calculated values of the first and second derivatives of pressure on the elasticity line. The article provides a detailed algorithm for finding the coefficient of elasticity line. The possibility of using the proposed equation in the construction of scale and wide-range equations of state.

Текст научной работы на тему «Уравнение линии упругости хладагента R236еа»

УДК 536.71

Уравнение линии упругости хладагента R236EA

Полторацкий М.И.,

Канд. техн. наук Рыков С.В. [email protected] Университет ИТМО 191002, Россия, Санкт-Петербург, ул. Ломоносова, 9 Свердлов А.В. [email protected] Air Comfort & Fire Safety Europe

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

Ключевые слова: линия упругости, уравнение состояния, линия фазового равновесия, масштабная теория, хладагент R236ea, гексафторпропан.

The equation of elasticity R236EA

Poltoratskiy M.I., Ph. D. Rykov S.V. [email protected] ITMO University 191002, Russia, St. Petersburg, Lomonosov str., 9 Sverdlov A.V. [email protected] Air Comfort & Fire Safety Europe

The paper proposed the equation of elasticity hexafluoropropane, taking into account the requirements of scale hypotheses and behaviors vapor pressure near the critical and triple points. A comparison with the experimental data on the pressure in the line of elasticity. Revealed good agreement with experiment in the vicinity of the critical point, in contrast to the equations of other authors. Information about the relative deviation is presented in graphical form On the basis of the proposed equation calculated detailed tables, depending on the temperature of the pressure on the line of phase equilibrium. Also in the pre-table shows calculated values of the first and second derivatives of pressure on the elasticity line. The article provides a detailed algorithm for finding the coefficient of elasticity line. The possibility of using the proposed equation in the construction of scale and wide-range equations of state.

Key words: elastic line, the equation of state, the line of phase equilibrium scaling theory, R236ea, hexaf-luoropropane.

1,1,1,2,3,3-гексафторпропан (^Y-236ea, номер CAS 431-63-0) был предложен в качестве альтернативы для хлорфторуглеродов (ХФУ) и гидрохлорфторуглеродов (ГХФУ). HFC236ea имеет нулевой

64

потенциал разрушения озона (ODP) (Ozon Depletion Potential) и низкий потенциал глобального потепления (GWP) (Global warming potential). По этой причине, HFC236ea используется в качестве замены для R-114 для применения в высокотемпературных тепловых насосах. С недавнего времени, ГФУ-236ea используется в качестве рабочей жидкости в микроканальных теплообменниках и низкотемпературных органических системах, работающих по циклу Ренкина, а также для утилизации низкопотенциальной теплоты. Все это вызывает повышенный интерес к данному холодильному агенту и требует уточнения расчетных зависимостей для нахождения равновесных свойств R236ea. При этом надо иметь в виду, что согласно современным представлениям [1], при расчете термодинамических свойств жидкостей и газов надо учитывать особенности термодинамической поверхности в окрестности критической точки. Именно этим и обусловлено использование при расчете термодинамических таблиц ССД таких важных в техническом плане веществ как аммиак, холодильные агенты R134a, R218, R23 уравнений, удовлетворяющих масштабной теории критических явлений. До настоящего времени давление на линии упругости R236ea рассчитывалось без учета особенностей критической точки. В данной работе в качестве исходного уравнения для расчета давления на линии упругости используется зависимость [2, 3]:

( „ 7 Л

1.2л I |2-a I |2-а+Д v^ п i

Ps = Рс еХР ~а01 tT 1 + а1Т + а2 |Т| + а3 |Т| + Zj aiX

2 I I ™3|

V /'=4

(1)

где - постоянные коэффициенты; рс - критическое давление; СИ - критический индекс изохорной теплоемкости; А - «неасимптотический» критический индекс; п I — массив из натуральных чисел;

т = Т/Тс-\^ = Т/Тс.

Зависимость (1) обеспечивает асимптотически правильное поведение линии упругости в окрестности тройной точки:

1п Р*1Рс =-у (2)

и удовлетворяет требованиям масштабной теории при Т —> Тс :

Р$ Т л | |2-а | |2-а+А

-= \ + а1т + а2\т\ К ••• (3)

Рс

Уравнение линии упругости в форме (1) выгодно отличается от уравнений типа уравнения Вагнера [4-6]:

Рс Т

( т

1- —

Т

V с J

(4)

Действительно, в структуре (4) есть коэффициенты, которые играют ведущую роль как вблизи тройной точки, так и в критической области. Из зависимостей (2) и (3) непосредственно следует, что уравнение (1) построено таким образом, что характер поведения линии упругости при описании области у тройной точки определяет коэффициент а0 , а в окрестности критической точки - коэффициенты щ,

a2 и a.

Важным обстоятельством является и то, что выбор выражения линии упругости в форме (1), как показано в [7, 8], дает возможность количественно и качественно верно передать поведение первой:

(1Т

а0т

= Рв

{ гр Л

—х-2

Т

V

) Рс -- + —ехр

т т:

X

{ п Л

V * у 7

X

(5)

^ I |1-а ^ * I |1-а+Д V-1 • п \ -1

- 2 - а а2 т - 2 - а + Д а2 т + ¿^п 1 а^х

/=4

Л

и второй производных функции р Т по температуре от тройной точки до критической. При этом

давление насыщенного пара в окрестности тройной точки уточняется по уравнению Клапейрона-Клаузиуса:

1 г t

и уравнению идеального газа

р Т ? Дй ' Рз=КТ Р~-

(6)

Здесь р - плотность насыщенного пара, а функция Г X - «кажущейся» теплота парообразова-

ния.

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

Значения критических индексов выбраны в соответствии с рекомендациями [9]: ОС =0,11; А = 0,51. Критические параметры приняты равными значениям р =34,2 бар и Тс =412,44 К, которые

согласуются как с уравнением Клапейрона-Клаузиуса (5), так и с фундаментальным уравнением состояния [10, 11]:

Р Р,Т = X <м Ё|ДР|

8+1+А,г/(3

ап х +

п=о

"з Л

(7)

Т + ДГ1пса + ДГсаХ X С,,х/ Ар

/=1 у=о

которое использовано при построении термического уравнения состояния хладагента R236ea.

Здесь ап X - масштабные функции свободной энергии Гельмгольца Р р,7^ ; % СО,^ - кроссо-

верная функция; X - масштабная переменная; 6 - критический индекс критической изотермы.

Поиск коэффициентов уравнения (1) осуществлялся на массиве экспериментальной и расчетной информации [12-14] путем поиска минимума функционала F :

Л'! -1

/ = 1

Г) пРасч _ „эксп

V,. , У В. / Р,

(8)

В результате получены следующие значения искомых параметров уравнения (1): а0 = 13,7; а = 8,587824476; а2 = 172,2216673; а3 = 45,56289106; а4 = - 202,4047127; а5 = - 43,53179291; а6 = -80,807200; а7 = - 41,50773797; п 4 =2; п 5 =4; п 6 = 5; п 7 = 6.

2

Относительные отклонения между расчетом по уравнению (1) и экспериментальными и табличными данными представлены на рис. 1, 2. На рис. 3-5 построены зависимости давления на линии упругости, первой и второй производных. Алгоритм нахождения коэффициентов уравнения (1) представлен на рис. 6, 7. Данный алгоритм реализован в пакете MathCAD для использования в учебных целях [21, 22].

Предлагаемое уравнение линии упругости (1) точнее передает экспериментальные данные [12, 13], чем уравнение [14], во всем температурном диапазоне. Уравнение линии упругости, предложенное в работе, с более высокой точностью описывает критическую область, чем уравнения работ [13, 14] (рис. 2). Это позволяет использовать уравнение (1) при построении как масштабных [ 15-18], так и широкодиапазонных уравнений состояния [10, 11, 19, 20].

ТК= 412.44 РК= 34.2

oPSj% 0,8 0,5 0,3 0,0 -0,3 -0,5 -0,8 -1,0

♦ ♦

4

J AJA А J 1 1 ,,-J 1 * * * 1 » * *» ♦ ♦

kiAJ А А ja * * ♦♦ * * ♦

А ■ Л А А

А А

i i

♦ 1 □ 2

3

190 200 210 220 230 240 250 260 270 280 290 300 310 320 330 340 350 360 370 380 390 400 410 Т,К

Рис. 1. Отклонение давления ps, рассчитанного по уравнению от данных:

(1) - [12]; (2) - [13]; (3) - [14]

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

5.% 2,0 1,5 1,0 0,5 0,0 -0,5 -1,0 -1,5 7.1

■ ♦ 1 ■ 2

■ ■ ■

■ ■ ■

■ ■ ■

—- /

———^

Ф

.5 250 275 300 325 350 375 400 Т,К

Рис. 2. Отклонение давления рц, рассчитанного по уравнению от данных, рассчитанных по уравнениям (1) - [13]; (2) - [14]

Р^бар 30 25 20 15 10 5 0 и

К

Ю 2( ю 2; Ю 11 Ю 260 280 300 320 340 360 380 400 Т,

Рис. 3. Линия упругости ps = /(Т), рассчитанная по уравнению (1)

рЧ 0,7 0,6 0,5 0,4 0,3 0,2 ОД 0 • 2( ёРэск

К) 250 300 350 400 Т,К

Рис. 4. Первая производная р'я

Р = 0,04 0,035 0,03 0,025 0,02 0,015 0,01 0,005 0 2С

)0 250 300 350 400 Т,К

Рис. 5. Вторая производная р"

Рис. 6. Блок-схема поиска коэффициентов уравнения (1) (начало)

Рис. 7. Блок-схема поиска коэффициентов уравнения (1) (окончание)

Таблица

Термодинамические свойства, рассчитанные на основе уравнения линии упругости хладагента R236ea

Т,К Ps,6ap P's P"s т,к Ps,6ap P's P"s т,к Ps,6ap P's P"s

190 0.00266 0.00029 2,8Е-05 264 0,518526 0.02423 0.00091 340 7,222 0,1869 0,0037

192 0.00331 0.00035 3,2Е-05 266 0,568832 0.0261 0.00096 342 7,6032 0,1943 0,0038

194 0.00408 0.00042 3,8Е-05 268 0,622985 0.02807 0.00101 344 7,9994 0,2019 0,0039

196 0,005 0,0005 4,4Е-05 270 0,681194 0.03015 0.00107 346 8,411 0,2097 0,004

198 0,0061 0,0006 0,00005 272 0,74367 0.03234 0.00112 348 8,8385 0,2178 0,0041

200 0.00739 0,0007 5,7Е-05 274 0,810632 0.03464 0.00118 350 9,2822 0,226 0,0042

202 0.00892 0.00083 6,5Е-05 276 0,882301 0.03705 0.00123 352 9,7425 0,2344 0,0043

204 0.01071 0.00097 7,4Е-05 278 0,958904 0.03957 0.00129 354 10,22 0,243 0,0044

206 0,0128 0.00112 8,4Е-05 280 1,040674 0.04222 0.00135 356 10.715 0.2518 0,0045

208 0.01522 0,0013 9,5Е-05 282 1,127845 0.04498 0.00141 358 11,227 0.2608 0,0046

210 0.01802 0,0015 0,00011 284 1,220658 0.04786 0.00147 360 11,758 0,27 0,0047

212 0.02125 0.00173 0,00012 286 1,319357 0.05086 0.00153 362 12,308 0,2795 0,0048

214 0.02495 0.00198 0,00013 288 1.42419 0.05399 0.0016 364 12.876 0.2891 0,0049

216 0.02917 0.00226 0,00015 290 1.53541 0.05725 0.00166 366 13.464 0.2989 0,005

218 0.03399 0.00257 0.00016 292 1,653273 0.06064 0.00173 368 14,072 0.3089 0.0051

220 0.03946 0.00291 0.00018 294 1.77804 0.06415 0.00179 370 14,7 0.3192 0.0052

222 0.04564 0.00328 0.0002 296 1,909975 0.0678 0.00186 372 15.349 0.3297 0.0053

224 0.05262 0.0037 0.00022 298 2.049347 0.07159 0.00193 374 16.019 0.3403 0.0054

226 0.06046 0.00415 0.00024 300 2,196429 0.07552 0.002 376 16,71 0.3513 0.0055

228 0.06925 0.00465 0.00026 302 2,351498 0.07958 0.00207 378 17.424 0.3624 0.0056

230 0.07908 0.00519 0.00028 304 2,514838 0.08379 0.00214 380 18,16 0.3738 0.0058

232 0.09004 0.00578 0.00031 306 2,686733 0.08814 0.00221 382 18,92 0.3855 0.0059

234 0.10223 0.00642 0.00033 308 2,867477 0.09263 0.00229 384 19.703 0.3975 0.0061

236 0.11575 0.00711 0.00036 310 3,057366 0.09728 0.00236 386 20.51 0,4099 0.0063

238 0.13072 0.00787 0.00039 312 3,256702 0.10208 0.00244 388 21.342 0,4226 0.0065

240 0.14725 0.00868 0.00042 314 3,465793 0.10704 0.00252 390 22.201 0,4357 0.0067

242 0.16547 0.00955 0.00045 316 3,684951 0.11215 0.0026 392 23.086 0,4494 0.007

244 0,1855 0.01049 0.00049 318 3,914496 0.11742 0.00268 394 23.999 0,4637 0.0073

246 0.20749 0.0115 0.00052 320 4,154752 0.12286 0.00276 396 24.941 0,4787 0.0077

248 0.23156 0.01258 0.00056 322 4,406049 0.12847 0.00284 398 25.914 0,4945 0.0082

250 0.25788 0.01374 0.0006 324 4,668724 0.13424 0.00293 400 26,92 0.5116 0.0088

252 0.28658 0.01498 0.00064 326 4.94312 0.14019 0.00302 402 27.961 0.5301 0.0097

254 0.31784 0.0163 0.00068 328 5,229585 0.14631 0.00311 404 29.041 0.5505 0.0108

256 0.35183 0.0177 0.00072 330 5,528472 0.15261 0.0032 406 30.165 0.5737 0.0125

258 0.38871 0.01919 0.00077 332 5,840143 0.15909 0.00329 408 31,339 0.6011 0.0151

260 0,42866 0.02078 0,00082 334 6,164961 0.16576 0.00338 410 32,574 0.6357 0,0201

262 0,47187 0,02245 0,00086 336 6,503298 0,17261 0,00347 412 33,893 0,6886 0,0389

Список литературы (References)

1. Лысенков В.Ф., Попов П.В., Рыков В.А. Параметрические масштабные уравнения состояния для асимптотической окрестности критической точки // Обзоры по теплофизическим свойствам ве-ществ.-М.:Изд-во ИВТАН. -1992. № 1(93). - 78 с.

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

3. Кудрявцев Д.А., Рыков В.А., Устюжанин Е.Е. Расчет линии упругости перфторпропана в пакете MathCad // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2014. № 3. С. 78-89.

4. Ustyuzhanin E.E., Shishakov V.V., Abdulagatov I.M., Popov P.V., Rykov V.A., Frenkel ML. Scaling Models of Thermodynamic Properties on the Coexistence Curve: Problems and Some Solutions// Russian Journal of Physical Chemistry B, 2012, Vol. 6, No. 8, P. 912-931.

5. Устюжанин Е.Е., Шишаков В.В., Абдулагатов И.М., Рыков В.А., Попов П.В. Давление насыщения технически важных веществ: модели и расчеты для критической области // Вестник Московского энергетического института. 2012. № 2. С. 34-43.

6. Устюжанин Е.Е., Шишаков В.В., Абдулагатов И.М., Попов П.В., Рыков В.А., Френкель М.Л. Скейлинговые модели для описания термодинамических свойств на линии насыщения: проблемы и некоторые решения // Сверхкритические флюиды: Теория и практика. 2012. Т. 7. № 3. С. 30-55.

7. Рыков С.В., Рябова Т.В. Расчет линии фазового равновесия аммиака в пакете MathCad // Научный журнал НИУ ИТМО. Серия: Холодильная техника и кондиционирование. 2013. № 2. С. 8.

8. Рыков С.В., Самолетов В.А., Рыков В.А. Линия насыщения аммиака // Вестник Международной академии холода. 2008. № 4. С. 20-21.

9. Киселев С.Б. Масштабное уравнение состояния индивидуальных веществ и бинарных растворов в широкой окрестности критических точек // Обзоры по теплофизическим свойствам веществ.- М.: Изд-во ИВТАН. 1989. №2(76). 150 с.

10. Рыков С.В. Фундаментальное уравнение состояния, учитывающее асимметрию жидкости // Научно-технический вестник Поволжья. 2014. № 1. С. 33-36.

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

12. V.A. Gruzdev, R.A. Khairulin, S.G. Komarov, S.V. Stankus, Int. J. Thermophys. 29 (2008) 546556.

13. D R. Defibaugh, K.A. Gillis, MR. Moldover, J.W. Schmidt, L A. Weber, Fluid Phase Equilib. 122 (1996) 131-155.

14. Lemmon, E.W., Huber, M.L., McLinden, M.O. NIST Standard Reference Database 23: Reference Fluid Thermodynamic and Transport Properties-REFPROP, Version 9.0, National Institute of Standards and Technology, Standard Reference Data Program, Gaithersburg, 2010

15. Рыков С.В., Кудрявцева И.В. Выбор параметров масштабного уравнения в физических переменных // Научно-технический вестник Поволжья. 2014. № 3. С. 34-37.

16. Рыков С.В., Рыков В.А. Обобщенная модель масштабного уравнения, основанная на феноменологической теории критических явлений // Фундаментальные исследования. 2014. № 11-9. С. 19161920.

17. Рыков С.В., Кудрявцева И.В. Непараметрическое масштабное уравнение и феноменологическая теория критических явлений // Фундаментальные исследования. 2014. № 9 (8). С. 1687-1692.

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

19. Рыков С.В. Фундаментальное уравнение состояния, учитывающее асимметрию жидкости // Научно-технический вестник Поволжья. 2014. № 1. С. 33-36.

20. КудрявцеваИ.В., Рыков В.А., Рыков С.В. Ассиметричное единое уравнение состояния R134a // Вестник Международной академии холода. 2008. № 2. С. 36-39.

21. Практикум по работе в математическом пакете MathCAD: Учеб. пособие / С.В. Рыков, И.В. Кудрявцева, С.А. Рыков, В.А. Рыков. - СПб.: Университет ИТМО, ИХиБТ, 2015. - 87 с.

22. Методы оптимизации в примерах в пакете MathCAD 15. Ч. I: Учебн. пособие / И.В. Кудрявцева, С.А. Рыков, С.В. Рыков, Е.Д. Скобов - СПб.: НИУ ИТМО, ИХиБТ, 2014. - 164 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.