R513A的饱和液相黏度特性研究

doi:10.11949/0438-1157.20210819

摘要/Abstract

摘要：

工质黏度对于制冷系统的设计优化至关重要。R513A因其良好的环保特性和热力循环性能，在冷水/热泵机组中有望成为R134a的主要替代制冷剂。为了进一步探究R513A的液相黏度特性，基于毛细管法设计搭建了制冷剂液相黏度测试系统，采用R134a作为标准液体对毛细管黏度计进行标定。在温度范围253.15~333.15 K内，开展了R513A的液相黏度实验研究，结果表明R513A的液相黏度略低于R134a。采用R-K多项式方程、硬球模型结合混合规则与实验数据进行关联，R-K方程计算值与实验值的平均绝对偏差（AAD）和最大绝对偏差（MAD）分别为0.71%、1.65%，硬球模型计算值与实验值的平均绝对偏差和最大绝对偏差分别为2.02%和3.39%。上述两种模型能够较好地预测R513A的饱和液相黏度，研究结果可为R513A的替代应用研究提供参考依据。

关键词: 制冷剂, R513A, 黏度, 混合物, 实验验证, 模型预测

Abstract:

The working fluid viscosity is very important for the design and optimization of refrigeration systems. R513A is expected to be the main alternative refrigerant for R134a in chiller/heat pump due to its good environmental protection characteristics and thermal cycle performance. In order to further explore the viscosity characteristics of R513A, a liquid viscosity measurement system for refrigerant was designed and built based on the capillary method, and the capillary viscometer was calibrated using R134a as the standard liquid. The saturated liquid viscosity of R513A was measured in the temperature range of 253.15—333.15 K. The results showed that the saturated liquid viscosity of R513A was slightly lower than that of R134a. The saturated liquid viscosity of R513A were predicted by using R-K polynomial equation, hard sphere model with mixing rule, the average absolute deviation and maximum absolute deviation were 0.71% and 1.65% by R-K equation, the average absolute deviation and maximum absolute deviation were 2.02% and 3.39% by RSH method. It can be seen that the above two models can predict the saturated liquid viscosity of R513A well, and the research results can provide important references for the alternative applications of R513A.

Key words: refrigerant, R513A, viscosity, mixtures, experimental validation, model prediction

中图分类号:

TB 61⁺2

陈裕博, 杨昭, 武晓昆, 吕子建, 张勇. R513A的饱和液相黏度特性研究[J]. 化工学报, 2021, 72(11): 5502-5509.

Yubo CHEN, Zhao YANG, Xiaokun WU, Zijian LYU, Yong ZHANG. Study on the saturated liquid viscosity characteristics of R513A[J]. CIESC Journal, 2021, 72(11): 5502-5509.

图/表 11

参考文献 30

1	Kasaeian A, Hossein S M, Sheikhpour M, et al. Applications of eco-friendly refrigerants and nanorefrigerants: a review[J]. Renewable and Sustainable Energy Reviews, 2018, 96: 91-99.
2	Birmpili T. Montreal Protocol at 30: the governance structure, the evolution, and the Kigali Amendment[J]. Comptes Rendus Geoscience, 2018, 350(7): 425-431.
3	Mota-Babiloni A, Makhnatch P, Khodabandeh R, et al. Experimental assessment of R134a and its lower GWP alternative R513A[J]. International Journal of Refrigeration, 2017, 74: 682-688.
4	张治平, 李宏波, 周宇, 等. 混合工质R513A替代R134a应用于离心式冷水机组的理论分析及试验研究[J]. 制冷与空调, 2018, 18(6): 85-89.
	Zhang Z P, Li H B, Zhou Y, et al. Theoretical analysis and experimental study on mixed refrigerant R513A as a substitute for R134a applied to centrifugal chiller[J]. Refrigeration and Air-Conditioning, 2018, 18(6): 85-89.
5	The Chemours Company FC, LLC. Properties, uses, storage, and handling[EB/OL]. [2017]. .
6	Sun J, Li W H, Cui B R. Energy and exergy analyses of R513A as a R134a drop-in replacement in a vapor compression refrigeration system[J]. International Journal of Refrigeration, 2020, 112: 348-356.
7	Yang M, Zhang H, Meng Z F, et al. Experimental study on R1234yf/R134a mixture (R513A) as R134a replacement in a domestic refrigerator[J]. Applied Thermal Engineering, 2019, 146: 540-547.
8	Velasco F J S, Illán-Gómez F, García-Cascales J R. Energy efficiency evaluation of the use of R513A as a drop-in replacement for R134a in a water chiller with a minichannel condenser for air-conditioning applications[J]. Applied Thermal Engineering, 2021, 182: 115915.
9	王勇, 杜国良, 王发忠, 等. R513A与R1234ze在螺杆冷水机组上替代R134a的试验研究[J]. 制冷与空调, 2019, 19(9): 77-82.
	Wang Y, Du G L, Wang F Z, et al. Experimental research of replacing R134a by R513A and R1234ze on screw chillers[J]. Refrigeration and Air-Conditioning, 2019, 19(9): 77-82.
10	Diani A, Rossetto L. R513A flow boiling heat transfer inside horizontal smooth tube and microfin tube[J]. International Journal of Refrigeration, 2019, 107: 301-314.
11	Diani A, Rossetto L. Characteristics of R513A evaporation heat transfer inside small-diameter smooth and microfin tubes[J]. International Journal of Heat and Mass Transfer, 2020, 162: 120402.
12	Han L Z, Zhu M S, Li X Y, et al. Viscosity of saturated liquid for 1,1,1,2-tetrofluoroethane[J]. Journal of Chemical & Engineering Data, 1995, 40(3): 650-652.
13	Sun L Q, Zhu M S, Han L Z, et al. Viscosity of difluoromethane and pentafluoroethane along the saturation line[J]. Journal of Chemical & Engineering Data, 1996, 41(2): 292-296.
14	刘志刚, 吴江涛, 吕萍. 一种适合于挥发性液体黏度测量的毛细管黏度计[J]. 热科学与技术, 2003, 2(4): 365-369.
	Liu Z G, Wu J T, Lyu P. A new capillary viscometer for volatile liquid[J]. Journal of Thermal Science and Technology, 2003, 2(4): 365-369.
15	吴江涛. 高精度流体热物性测试实验系统的研制及二甲醚热物理性质的研究[D]. 西安: 西安交通大学, 2003.
	Wu J T. Development of the new thermophysical properties measurement system and research of thermophysical properties of dimethyl ether [D]. Xi'an: Xi'an Jiaotong University, 2003.
16	孟现阳, 吴江涛, 刘志刚. 甲醇与蓖麻油混合物运动黏度的实验研究[J]. 工程热物理学报, 2006, 27(S1): 79-81.
	Meng X Y, Wu J T, Liu Z G. The liquid kinematic viscosity of the binary mixtures of methanol and castor oil[J]. Journal of Engineering Thermophysics, 2006, 27(S1): 79-81.
17	潘江, 吴江涛, 刘志刚. 甲基叔丁醚的饱和液相黏度的实验研究[J]. 西安交通大学学报, 2006, 40(9): 1024-1027.
	Pan J, Wu J T, Liu Z G. Viscosity of saturated liquid methyl tert-butyl ether[J]. Journal of Xi'an Jiaotong University, 2006, 40(9): 1024-1027.
18	袁晓蓉. 混合制冷剂黏度特性研究[D]. 杭州: 浙江大学, 2015.
	Yuan X R. Investigation of the viscosity characteristics for refrigerant mixtures[D]. Hangzhou: Zhejiang University, 2015.
19	许晨怡, 叶恭然, 郭豪文, 等. 制冷剂R1336mzz(E)液相黏度理论与实验研究[J]. 化工学报, 2021, 72(6): 3261-3269.
	Xu C Y, Ye G R, Guo H W, et al. Experimental and theoretical study on liquid viscosity of R1336mzz(E)[J]. CIESC Journal, 2021, 72(6): 3261-3269.
20	Meng X Y, Qiu G S, Wu J T, et al. Viscosity measurements for 2,3,3,3-tetrafluoroprop-1-ene (R1234yf) and trans-1,3,3,3-tetrafluoropropene (R1234ze(E))[J]. The Journal of Chemical Thermodynamics, 2013, 63: 24-30.
21	Meng X Y, Wen C Y, Wu J T. Measurement and correlation of the liquid viscosity of trans -1-chloro-3,3,3-trifluoropropene (R1233zd(E))[J]. The Journal of Chemical Thermodynamics, 2018, 123: 140-145.
22	Meng X Y, Zhang J B, Wu J T. Compressed liquid viscosity of 1,1,1,3,3-pentafluoropropane (R245fa) and 1,1,1,3,3,3-hexafluoropropane (R236fa)[J]. Journal of Chemical & Engineering Data, 2011, 56(12): 4956-4964.
23	孟现阳, 孙裕坤, 曹法立, 等. R32和R32/R1234yf混合工质黏度实验测量及模型[J]. 制冷学报, 2018, 39(2): 39-47.
	Meng X Y, Sun Y K, Cao F L, et al. Viscosity measurements and modeling for R32 and binary mixture of R32/R1234yf [J]. Journal of Refrigeration, 2018, 39(2): 39-47.
24	Bair S, Laesecke A. Viscosity measurements of R32 and R410A to 350 MPa[J]. International Journal of Refrigeration, 2017, 83: 157-167.
25	Cui J W, Yan S M, Bi S S, et al. Saturated liquid dynamic viscosity and surface tension of trans -1-chloro-3,3,3-trifluoropropene and dodecafluoro-2-methylpentan-3-one[J]. Journal of Chemical & Engineering Data, 2018, 63(3): 751-756.
26	Cui J W, Bi S S, Meng X Y, et al. Surface tension and liquid viscosity of R32+R1234yf and R32+R1234ze[J]. Journal of Chemical & Engineering Data, 2016, 61(2): 950-957.
27	Dang Y G, Kamiaka T, Dang C B, et al. Liquid viscosity of low-GWP refrigerant mixtures (R32+R1234yf) and (R125+R1234yf)[J]. The Journal of Chemical Thermodynamics, 2015, 89: 183-188.
28	中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 黏度测量方法: [S]. 北京: 中国标准出版社, 2009.
	General Administration of Quality Supervision, Inspection and Quarantine of the People's Republic of China, Standardization Administration of the People's Republic of China. Methods of viscosity measurement: [S]. Beijing: Standards Press of China, 2009.
29	Lemmon E W, Huber M L, Mclinden M O. NIST Reference Fluid Thermodynamic and Transport Properties -REFPROP, Version9.1 [DB]. National Institute of Standards and Technology, 2013.
30	Zhao G J, Bi S S, Fröba A P, et al. Liquid viscosity and surface tension of R1234yf and R1234ze under saturation conditions by surface light scattering[J]. Journal of Chemical & Engineering Data, 2014, 59(4): 1366-1371.

T/K	p/MPa	ρ_L/(kg·m^-3)	ρ_V/(kg·m^-3)	t/s	ν_Ref/(m²·s^-1)	((ν_exp-ν_Ref)/ν_Ref)/%
263.17	0.202	1327.1	10.049	579.18	2.2868×10^-7	-1.00
268.16	0.250	1311.1	12.082	553.03	2.1680×10^-7	-0.49
273.18	0.298	1294.7	14.443	536.44	2.0578×10^-7	0.77
278.16	0.339	1278	17.137	510.65	1.9567×10^-7	0.09
283.13	0.424	1261	20.213	490.06	1.8630×10^-7	0.10
288.12	0.497	1243.5	23.736	476.4	1.7752×10^-7	1.41
293.17	0.586	1225.3	27.797	456	1.6920×10^-7	0.81
298.15	0.676	1206.7	32.35	438.43	1.6150×10^-7	0.46
303.15	0.775	1187.5	37.535	419.87	1.5422×10^-7	-0.53
308.1	0.878	1167.7	43.353	405.35	1.4739×10^-7	0.81
313.13	0.996	1146.8	50.057	392.65	1.4081×10^-7	-0.82

T/K	p/MPa	ρ_L/(kg·m^-3)	ρ_V/(kg·m^-3)	t/s	ν_Ref/(m²·s^-1)	((ν_exp-ν_Ref)/ν_Ref)/%
263.17	0.202	1327.1	10.049	579.18	2.2868×10^-7	-1.00
268.16	0.250	1311.1	12.082	553.03	2.1680×10^-7	-0.49
273.18	0.298	1294.7	14.443	536.44	2.0578×10^-7	0.77
278.16	0.339	1278	17.137	510.65	1.9567×10^-7	0.09
283.13	0.424	1261	20.213	490.06	1.8630×10^-7	0.10
288.12	0.497	1243.5	23.736	476.4	1.7752×10^-7	1.41
293.17	0.586	1225.3	27.797	456	1.6920×10^-7	0.81
298.15	0.676	1206.7	32.35	438.43	1.6150×10^-7	0.46
303.15	0.775	1187.5	37.535	419.87	1.5422×10^-7	-0.53
308.1	0.878	1167.7	43.353	405.35	1.4739×10^-7	0.81
313.13	0.996	1146.8	50.057	392.65	1.4081×10^-7	-0.82

T/K	p/MPa	ρ_L/ (kg·m^-3)	ρ_V/ (kg·m^-3)	k	t/s	ν_exp/(mm²·s^-1)	η_exp/(μPa·s)
253.15	0.156	1288	7.7757	0.9940	571	0.2233	287.59
258.20	0.190	1273	9.4493	0.9926	543.3	0.2111	268.76
263.09	0.225	1258.1	11.329	0.9910	519.28	0.2005	252.23
268.16	0.272	1242.4	13.577	0.9891	499.35	0.1915	237.94
273.19	0.325	1226.4	16.144	0.9868	480.04	0.1828	224.13
278.16	0.380	1210.3	19.047	0.9843	462.13	0.1745	211.22
283.13	0.455	1193.7	22.355	0.9813	445.03	0.1666	198.83
288.21	0.531	1176.3	26.205	0.9777	428.5	0.1588	186.76
293.18	0.621	1158.7	30.484	0.9737	413.09	0.1514	175.40
298.15	0.709	1140.6	35.331	0.9690	399.62	0.1447	165.09
303.13	0.798	1121.8	40.830	0.9636	386.22	0.1380	154.84
308.18	0.892	1102	47.148	0.9572	375.19	0.1322	145.72
313.17	1.018	1081.6	54.231	0.9499	364.62	0.1265	136.84
318.14	1.136	1060.4	62.237	0.9413	354.97	0.1211	128.38
323.12	1.295	1038	71.37	0.9312	345.41	0.1155	119.87
328.19	1.460	1013.9	82.012	0.9191	337.07	0.1102	111.72
333.16	1.640	988.69	94.023	0.9049	329.5	0.1050	103.81

T/K	p/MPa	ρ_L/ (kg·m^-3)	ρ_V/ (kg·m^-3)	k	t/s	ν_exp/(mm²·s^-1)	η_exp/(μPa·s)
253.15	0.156	1288	7.7757	0.9940	571	0.2233	287.59
258.20	0.190	1273	9.4493	0.9926	543.3	0.2111	268.76
263.09	0.225	1258.1	11.329	0.9910	519.28	0.2005	252.23
268.16	0.272	1242.4	13.577	0.9891	499.35	0.1915	237.94
273.19	0.325	1226.4	16.144	0.9868	480.04	0.1828	224.13
278.16	0.380	1210.3	19.047	0.9843	462.13	0.1745	211.22
283.13	0.455	1193.7	22.355	0.9813	445.03	0.1666	198.83
288.21	0.531	1176.3	26.205	0.9777	428.5	0.1588	186.76
293.18	0.621	1158.7	30.484	0.9737	413.09	0.1514	175.40
298.15	0.709	1140.6	35.331	0.9690	399.62	0.1447	165.09
303.13	0.798	1121.8	40.830	0.9636	386.22	0.1380	154.84
308.18	0.892	1102	47.148	0.9572	375.19	0.1322	145.72
313.17	1.018	1081.6	54.231	0.9499	364.62	0.1265	136.84
318.14	1.136	1060.4	62.237	0.9413	354.97	0.1211	128.38
323.12	1.295	1038	71.37	0.9312	345.41	0.1155	119.87
328.19	1.460	1013.9	82.012	0.9191	337.07	0.1102	111.72
333.16	1.640	988.69	94.023	0.9049	329.5	0.1050	103.81

拟合参数	R1234yf^[30]	R134a
k₀	-5.93	2.2484
k₁	0.03342	-0.01585
k₂	-8.626×10^-5	4.0803×10^-5
k₃	1.077×10^-7	-3.7285×10^-8
k₄	4.145	31.7549
n	0.906	7.1102