袁曉蓉 高贊軍 王學(xué)會(huì) 徐象國 陳光明 韓曉紅

(浙江大學(xué)制冷與低溫研究所 杭州 310027)

一種適用于低沸點(diǎn)混合制冷劑粘度測量裝置的研制

袁曉蓉 高贊軍 王學(xué)會(huì) 徐象國 陳光明 韓曉紅

(浙江大學(xué)制冷與低溫研究所 杭州 310027)

針對(duì)混合制冷劑等揮發(fā)性混合物粘度測量的需要，研制了一種新的適用于測量低沸點(diǎn)混合物粘度的旋轉(zhuǎn)式毛細(xì)管粘度計(jì)。該新型粘度測量裝置在壓力容器內(nèi)嵌入旋轉(zhuǎn)式毛細(xì)管粘度計(jì)，將旋轉(zhuǎn)法升液和壓力容器承壓結(jié)合起來，避免了傳統(tǒng)密封型毛細(xì)管粘度計(jì)由于抽放氣的升液方式而導(dǎo)致混合物成分的變化，可以實(shí)現(xiàn)在較高壓力下循環(huán)測量揮發(fā)性混合物溶液的粘度。采用R22和R290對(duì)旋轉(zhuǎn)式毛細(xì)管粘度計(jì)進(jìn)行了標(biāo)定，并用R410A對(duì)粘度計(jì)的測量精度和性能進(jìn)行了評(píng)價(jià)，粘度測量值與文獻(xiàn)值最大相對(duì)偏差為0.81%。

混合制冷劑；粘度；旋轉(zhuǎn)式毛細(xì)管粘度計(jì)；低沸點(diǎn)液體；高壓

我國是制冷劑生產(chǎn)和消耗大國，近年來制冷空調(diào)行業(yè)發(fā)展迅速，制冷劑帶來的環(huán)境問題也引起強(qiáng)烈關(guān)注。人們開始尋找環(huán)境友好型的制冷劑用以替代傳統(tǒng)制冷劑。然而，制冷劑的替代是個(gè)嚴(yán)格而慎重的過程，需要對(duì)制冷劑相關(guān)的各個(gè)方面進(jìn)行一系列研究測試，獲得大量而全面的數(shù)據(jù)，才能被很好地評(píng)估和論證，并且滿足工業(yè)應(yīng)用的需要，最終完成替代過程。

粘度是流體的基本屬性之一，在擬定和測試液體行為的各種理論中具有相當(dāng)?shù)膬r(jià)值[1]。粘度是研究流動(dòng)與換熱的重要基礎(chǔ)數(shù)據(jù)，當(dāng)有傳熱發(fā)生時(shí)，流體熱物性與傳熱是相互影響的，在流體的熱物性中，粘度對(duì)溫度的變化最為敏感，它對(duì)傳熱過程的影響遠(yuǎn)大于其他物性的影響，這必將對(duì)流動(dòng)中的傳熱和壓降特性產(chǎn)生重要影響[2]。制冷劑的液相粘度測量一直是國際上熱物性研究關(guān)注的內(nèi)容之一，目前國內(nèi)外研究學(xué)者已提出了多種不同的方法對(duì)不同的制冷劑粘度進(jìn)行了測量。表1給出了國內(nèi)外對(duì)制冷劑粘度的實(shí)驗(yàn)研究情況。從表1中可以看出這些研究主要集中在對(duì)純質(zhì)的粘度研究上，而對(duì)混合制冷劑的粘度，特別是液相粘度的研究非常少。

劉志剛等[11]研制了一種帶壓力容器的玻璃毛細(xì)管粘度計(jì)。該粘度計(jì)利用旁通管升液，但由于采用抽放氣的方式，測量揮發(fā)性混合物溶液粘度存在一定的局限性，這是因?yàn)槌榉艢獾纳悍绞綍?huì)改變混合物組分。美國NIST熱物理部的Ripple[3]提出了一個(gè)類似Ubbelhode粘度計(jì)的高壓毛細(xì)管粘度計(jì)。粘度計(jì)能承受超過1MPa的蒸汽壓力，整體精度約為3%，測量精度可與玻璃毛細(xì)管相媲美。吳江濤等[12]在此基礎(chǔ)上，結(jié)合玻璃毛細(xì)管粘度計(jì)計(jì)時(shí)泡的特點(diǎn)，研制了一種利用翻轉(zhuǎn)方式來升液的金屬毛細(xì)管粘度計(jì)，避免了放氣對(duì)混合物成分的影響，可以實(shí)現(xiàn)揮發(fā)性混合物或純質(zhì)溶液體系粘度的測量，但受目前國內(nèi)加工水平的限制，不銹鋼計(jì)時(shí)泡和毛細(xì)管內(nèi)表面粗糙度、光潔度和線性度難以達(dá)到要求，使得表面張力對(duì)流體流動(dòng)的影響加強(qiáng)，實(shí)驗(yàn)結(jié)果的不確定度加強(qiáng)。

在前人的工作基礎(chǔ)上，針對(duì)混合制冷劑等揮發(fā)性混合物粘度測量的需要，將旋轉(zhuǎn)法升液和壓力容器承壓結(jié)合起來，研制了一種新的適用于低沸點(diǎn)混合物粘度測量的旋轉(zhuǎn)式毛細(xì)管粘度計(jì)。

1 新型粘度測量裝置的提出

1.1 毛細(xì)管粘度計(jì)的測量原理與方法

毛細(xì)管粘度計(jì)測量的基本原理是 Hagen-Poiseuille定律。19世紀(jì)法國科學(xué)家泊肅葉對(duì)牛頓流體在毛細(xì)管中的流動(dòng)作了理論推導(dǎo)及實(shí)驗(yàn)研究。他假設(shè):1)流體是不可壓縮的；2)流體是牛頓流體；3)管子足夠長、直線狀、內(nèi)徑均勻一致；4)在管壁處無滑動(dòng)；5)流動(dòng)為穩(wěn)定流；6)流動(dòng)為層流。

根據(jù)描述粘性流體在管子中流動(dòng)的哈根-泊肅葉定律有:

式中:η為流體的動(dòng)力粘度；R為毛細(xì)管的半徑；p為毛細(xì)管兩端的壓力差，對(duì)于重力型毛細(xì)管粘度計(jì)而言，即為重力和浮力作用之差；Q為流經(jīng)毛細(xì)管的流量；L為毛細(xì)管長度；t為流體的流動(dòng)時(shí)間；V為在t時(shí)間內(nèi)流經(jīng)毛細(xì)管的流體體積。

泊氏公式是建立在外力完全用于克服內(nèi)摩擦力的假設(shè)上的，但實(shí)際上，需考慮動(dòng)能修正和末端修正。經(jīng)動(dòng)能修正和末端修正后的哈根－泊肅葉公式成為:

式中:ρL為待測流體液相密度；ρV為待測流體氣相密度；h為毛細(xì)管內(nèi)液柱的高度；g為重力加速度；m為動(dòng)能修正系數(shù)；n為末端修正系數(shù)。η，t，V，L，R代表的含義與式(1)相同。

毛細(xì)管法可以分為絕對(duì)測量和相對(duì)測量。對(duì)于絕對(duì)測量，按照式(2)，不僅要求對(duì)毛細(xì)管的尺寸、流量、壓力進(jìn)行精確的測量，還需要用實(shí)驗(yàn)方法獲得動(dòng)能修正系數(shù)m及末端修正系數(shù)n。由于毛細(xì)管粘度計(jì)的尺寸在加工中難以控制，制成后也很難測量，因此絕對(duì)測量法目前僅用于粘度基準(zhǔn)的建立和一些特殊場合。相對(duì)測量法則不需要測量這些參數(shù)，而是通過與已知粘度的標(biāo)準(zhǔn)液進(jìn)行比較來測量。即先用已知粘度的流體準(zhǔn)確標(biāo)定毛細(xì)管粘度計(jì)常數(shù)A和B，在實(shí)驗(yàn)中只需準(zhǔn)確測量待測流體的流動(dòng)時(shí)間即可由式(3)得到運(yùn)動(dòng)粘度ν，再由液體的密度計(jì)算出動(dòng)力粘度η。

將式(2)轉(zhuǎn)換為:

對(duì)同一粘度計(jì)而言，其幾何尺寸是固定的，因此，V、R、L＋nR、h為常數(shù)，g、n、m可視為常數(shù)，于是毛細(xì)管粘度計(jì)常數(shù)A和B的值也是固定的。

1.2 新型粘度測量裝置的研制

新型粘度測量系統(tǒng)如圖1所示，主要包括:毛細(xì)管粘度計(jì)、壓力容器、真空系統(tǒng)、恒溫系統(tǒng)、測量系統(tǒng)及輔助管路。

由于混合制冷劑是一種高揮發(fā)性的低沸點(diǎn)工質(zhì)，因此必須在一個(gè)能耐高壓的密閉容器中進(jìn)行。壓力容器的設(shè)計(jì)壓力為10 MPa，各部件之間通過不銹鋼螺栓連接，考慮到實(shí)驗(yàn)觀測的需要，釜體的左右兩側(cè)設(shè)計(jì)有直徑為40 mm的觀察窗，如圖2所示。

流體在毛細(xì)管中的流動(dòng)狀態(tài)是與毛細(xì)管的結(jié)構(gòu)參數(shù)直接相關(guān)的，對(duì)于毛細(xì)管長度而言，希望盡可能長，以減少起始段的影響，但是受恒溫槽高度的限制，所使用的毛細(xì)管長度L為125 mm。對(duì)于毛細(xì)管內(nèi)徑的確定，主要依據(jù)是保證Hagen-Poiseuille定律的適用，即毛細(xì)管內(nèi)的流動(dòng)屬于層流:

根據(jù)式(5)可以得到V的計(jì)算公式:

為了便于計(jì)算，設(shè)計(jì)中不考慮動(dòng)能損失和末端效應(yīng)，則:把式(11)代入式(9)后可以得到:

k值取為1.0(在離臨界點(diǎn)較遠(yuǎn)時(shí)是可行的)，于是運(yùn)動(dòng)粘度可以用下式計(jì)算:

則式(12)可化簡為:

上式中L和g的值分別為125 mm和9800 mm/ s2，h的值與L的值則相同，ν的選取則要依據(jù)待測流體的粘度而定，取最小值。考慮到目前制冷工質(zhì)等低沸點(diǎn)流體的運(yùn)動(dòng)粘度一般大于0.1 mm2/s，于是可以得到毛細(xì)管的半徑:

由于所需毛細(xì)管內(nèi)徑太小，廠家難以精確加工，實(shí)際使用的毛細(xì)管半徑據(jù)廠家描述為0.175 mm，毛細(xì)管半徑經(jīng)初步標(biāo)定為0.2 mm，此時(shí)按照式(14)計(jì)算出Re為1960＜2000，可認(rèn)為毛細(xì)管內(nèi)流動(dòng)處于層流狀態(tài)。

對(duì)于計(jì)時(shí)泡的體積，按運(yùn)動(dòng)粘度最小時(shí)，流體經(jīng)過毛細(xì)管的最短時(shí)間不能小于300 s來確定計(jì)時(shí)泡的體積。根據(jù)式(11)和式(13)可以得到計(jì)時(shí)泡體積V與測量時(shí)間t之間的關(guān)系:

由于體積的精確加工有相當(dāng)大的困難，粘度計(jì)的計(jì)時(shí)泡設(shè)計(jì)體積為9000 mm3，盡管減少了測量的時(shí)間，但滿足實(shí)驗(yàn)要求。

為保證旋轉(zhuǎn)升液時(shí)，計(jì)時(shí)球能被被測樣本填滿，上貯液器的體積應(yīng)大于計(jì)時(shí)球體積，設(shè)計(jì)的上貯液器體積為20400 mm3，下貯液器體積與上貯液器一樣。旋轉(zhuǎn)式毛細(xì)管粘度計(jì)已獲得國家發(fā)明專利授權(quán)[13]，其實(shí)物圖如圖3所示。

表1 國內(nèi)外制冷劑粘度研究情況Tab.1 Domestic and international research on refrigerant viscosity

圖1 實(shí)驗(yàn)裝置示意圖Fig.1 The schematic layout of the experimental system

圖2 壓力容器與粘度計(jì)組合剖面圖Fig.2 Sectional view of pressure vessel and viscometer

圖3 旋轉(zhuǎn)式粘度計(jì)實(shí)物圖Fig.3 Photo of rotatable capillary viscometer

2 實(shí)驗(yàn)裝置的可靠性驗(yàn)證

采用R22和R290對(duì)毛細(xì)管粘度計(jì)進(jìn)行了標(biāo)定，并用R410A對(duì)粘度計(jì)的測量精度進(jìn)行了評(píng)價(jià)。

2.1 樣品來源

實(shí)驗(yàn)用的R22、R290和R410A如表2所示。所有試劑在實(shí)驗(yàn)前沒有進(jìn)行進(jìn)一步的提純。

2.2 毛細(xì)管粘度計(jì)常數(shù)A和B的標(biāo)定

由毛細(xì)管粘度計(jì)相對(duì)測量法的原理和式(3)可知，在測量工質(zhì)粘度之前，必須要利用標(biāo)準(zhǔn)液體對(duì)毛細(xì)管粘度計(jì)的常數(shù)A和B進(jìn)行標(biāo)定。

表2 樣品參數(shù)Tab.2 Parameters of the sample

實(shí)驗(yàn)選用R22和R290作為標(biāo)準(zhǔn)液，通過測量R22和R290液相流經(jīng)粘度計(jì)計(jì)時(shí)泡的時(shí)間以及對(duì)應(yīng)的飽和溫度值，根據(jù)飽和溫度查詢得到粘度文獻(xiàn)值(NIST REF-PROP 9.0)，回歸得到毛細(xì)管粘度計(jì)常數(shù)A、B。

經(jīng)過粗大誤差判斷，并將同一溫度點(diǎn)的實(shí)驗(yàn)數(shù)據(jù)取平均值后得到最終用于標(biāo)定毛細(xì)管粘度計(jì)常數(shù)的實(shí)驗(yàn)數(shù)據(jù)，初步標(biāo)定的粘度計(jì)常數(shù):

A＝4.523293e－10m2/s2

B＝1.971138e－6m2

選擇兩種標(biāo)準(zhǔn)液和多個(gè)溫度點(diǎn)的數(shù)據(jù)來擬合A值的目的是為了減少測量中的偶然誤差。

表3給出了R22和R290的液相粘度測量結(jié)果，表中T和t是測量值。ρl、ρg、νref為NIST數(shù)據(jù)，νcal為式(3)計(jì)算值。從表3中可以看出計(jì)算值νcal與文獻(xiàn)值νref(NIST REF-PROP 9.0)差值的平均值為0.000001429 mm2/s，最大值為0.00108 mm2/s，平均相對(duì)偏差為0.0014%，最大相對(duì)偏差為－0.65%。

2.3 R410A液相粘度的測量

為了評(píng)價(jià)粘度計(jì)的測量精度和性能，對(duì)溫度為278～288 K的R410A粘度進(jìn)行了測量，并分別與文獻(xiàn)值(NIST REF-PROP 9.0)進(jìn)行了比較，結(jié)果如表4所示。從表4中可以看出，R410A的粘度測量值νcal與文獻(xiàn)值νref(NIST REF-PROP 9.0)差值的平均值為0.000696667 mm2/s，最大值為0.00107 mm2/s，平均相對(duì)偏差為0.5367%，最大相對(duì)偏差為0.81%。

2.4 實(shí)驗(yàn)不確定度分析

關(guān)于實(shí)驗(yàn)不確定度分析，實(shí)驗(yàn)中被測量的物理量是溫度T和時(shí)間t。溫度測量的不確定度包括溫度測量系統(tǒng)的不確定度和恒溫槽的溫度波動(dòng)度。鉑電阻溫度計(jì)的測溫不確定度為±1 mK。對(duì)于時(shí)間的測量，計(jì)時(shí)器精度是0.1 s，實(shí)驗(yàn)時(shí)間偏差最大約±3 s，最短測量時(shí)間t＝295.9 s，因此△t/t＝±1.0%。根據(jù)式(3)，由于修正項(xiàng)B/t2在運(yùn)動(dòng)粘度中所占比例很小，可將其忽略，因而運(yùn)動(dòng)粘度的測量不確定度約等于時(shí)間測量的不確定度，即運(yùn)動(dòng)粘度的測量不確定度為±1.0%。又由動(dòng)力粘度的計(jì)算公式:

液相密度的不確定度一般小于±1.0%，因此實(shí)驗(yàn)動(dòng)力粘度的不確定度為±2.0%。

表3 R22和R290粘度實(shí)驗(yàn)值與標(biāo)準(zhǔn)值的比較Tab.3 Comparison of R22 and R290 viscosity experimental value and the standard value

表4 R410A粘度實(shí)驗(yàn)值與標(biāo)準(zhǔn)值的比較Tab.4 Comparison of R410A viscosity experimental value and the standard value

3 小結(jié)

1)在毛細(xì)管粘度測量原理的基礎(chǔ)上，提出了一種低沸點(diǎn)混合制冷劑粘度測量裝置，該裝置的主要特點(diǎn)是在壓力容器內(nèi)嵌入旋轉(zhuǎn)式毛細(xì)管粘度計(jì)，將旋轉(zhuǎn)法升液和壓力容器承壓相結(jié)合，這樣不僅可以解決由于放氣導(dǎo)致混合物成分的變化和待測流體溫度的變化等問題，而且可以承受高壓，循環(huán)測量。

2)采用R22和R290對(duì)毛細(xì)管粘度計(jì)進(jìn)行標(biāo)定，并對(duì)R410A的粘度進(jìn)行了測量，粘度測量值與文獻(xiàn)值最大相對(duì)偏差為0.81%。測量結(jié)果表明，該新型粘度計(jì)具有良好的可操作性與較高的精度。

3)在隨后的研究中將給出更大溫區(qū)范圍內(nèi)的混合制冷劑的液相粘度。

[1] Phillips T W，Murphy K P.Liquid Viscosity of Halocarbons[J].Journal of Chemical and Engineering Data，1970，15(2):304-307.

[2] 師晉生.管內(nèi)定壁溫定型流傳熱有效能損失分析[C]//中國工程熱物理學(xué)會(huì)第十一屆年會(huì)論文集工程熱力學(xué)與能源利用(上冊(cè)):13-16.

[3] Ripple D.A compact high-pressure capillary viscometer [J].Rev.Sci.Instrum.，1992，63(5):3153-3155.

[4] Oliveira C M B P，Wakeham W A.The Viscosity of Liquid R134a[J].International Journal of Thermophysics，1993，14(1):33-44.

[5] Geller V Z，Paulaitis M E，Bivens D B，et al.Viscosities of HFC-32 and HFC-32/Lubricant Mixtures[J].International Journal of Thermophysics，1996，17(1):75-83.

[6] 孫立群，朱明善，韓禮鐘，等.HFC-32飽和液粘度的試驗(yàn)研究[J].工程熱物理學(xué)報(bào)，1996，17(4):393-396.(Sun Liqun，Zhu Mingshan，Han Lizhong，et al.Experimental Research of the viscosity of saturated liquid HFC-32 [J].Journal of Engineering Thermophysics，1996，17 (4):393-396.)

[7] Fr?ba A P，Leipertz A.Thermophysical properties of the refrigerant mixtures R410A and R407C from dynamic light scattering(DLS)[J].International Journal of Thermophysics，2003，24(5):1185-1206.

[8] 吳江濤.高精度流體熱物性測試實(shí)驗(yàn)系統(tǒng)的研制及二甲醚熱物理性質(zhì)的研究[D].西安:西安交通大學(xué)能源與動(dòng)力工程學(xué)院，2003.

[9] 何茂剛，張穎，王小飛，等.傾斜管式液體粘度測量方法理論及實(shí)驗(yàn)研究[J].工程熱物理學(xué)報(bào)，2007，28(3): 382-384.(He Maogang，Zhang Ying，Wang Xiaofei，et al.A method for liquid viscosity coefficient measurement using inclined tube[J].Journal of Engineering Thermophysics，2007，28(3):382-384.)

[10]Zhelezny V P，Sechenyha V V，Semenyuka Yu V，et al. An experimental investigation and modelling of the viscosity refrigerant/oil solutions[J].International Journal of Refrigeration，2009，32:1389-1395.

[11]劉志剛，吳江濤，呂萍.一種適合于揮發(fā)性液體黏度測量的毛細(xì)管黏度計(jì)[J].熱科學(xué)與技術(shù)，2003，2(4):365-369.(Liu Zhigang，Wu Jiangtao，Lü Ping.A new capillary viscometer for volatile liquid[J].Journal of Thermal Science and Technology，2003，2(4):365-369.)

[12]吳江濤，畢勝山，劉志剛.一種適用于低沸點(diǎn)液體粘度測量的金屬毛細(xì)管粘度計(jì)[J].工程熱物理學(xué)報(bào)，2004，25 (6):921-924.(Wu Jiangtao，Bi Shengshan，Liu Zhigang.A new metal capillary viscometer for the volatile liquid[J].Journal of Engineering Thermophysics，2004，25 (6):921-924.)

[13] 浙江大學(xué).一種旋轉(zhuǎn)式毛細(xì)管粘度計(jì):中國，CN201110394660.9[P].2012-3-28.

[14]Lemmon E W，McLinden M O，Huber M L.NIST REFPROP，version 9.0[CP/DK].USA:National Institute of Standards and Technology，2010.

Development of a Viscosity Measuring Device for Low Boiling Point Refrigerant Mixtures

Yuan Xiaorong Gao Zanjun Wang Xuehui Xu Xiangguo Chen Guangming Han Xiaohong

(Institute of Refrigeration and Cryogenics，Zhejiang University，Hangzhou，310027，China)

According to the need for viscosity measurement of refrigerant mixtures，a new rotatable capillary viscometer is developed and designed.It is adequate for the viscosity investigation of the low boiling point liquid mixtures.A rotatable capillary viscometer is embedded in a pressure vessel.It combines the rotation method and the pressure vessel to realize the continuous experiment and eliminates the error brought from the deflation，which can repeat the viscosity measurement of the low boiling point liquid mixtures in the high pressure. The rotatable capillary viscometer was calibrated by R22 and R290 and checked by R410A.According to the results，the relative deviation of the viscometer has been controlled to be less than 0.81%.

refrigerant mixtures；viscosity；rotatable capillary viscometer；low boiling point liquid；high pressure

TB61＋2；TB65

A

0253－4339(2014)05－0001－06

10.3969/j.issn.0253－4339.2014.05.001

韓曉紅，女(1973－)，副教授，浙江大學(xué)制冷與低溫研究所，(0571)87953944，E-mail:hanxh66＠zju.edu.cn。研究方向:制冷劑的替代及關(guān)鍵技術(shù)研究(主要包括:混合制冷劑的汽液相平衡、制冷劑的循環(huán)性能、制冷劑與潤滑油的相容性研究、制冷劑/潤滑油與材料的相容性、制冷劑的密度、制冷劑的粘度，制冷劑在管內(nèi)的傳熱物性研究等)，制冷用換熱器的數(shù)值模擬與軟件開發(fā)，換熱器新型材料及換熱器的結(jié)構(gòu)設(shè)計(jì)；液化天然氣的物性，如天然氣的汽液相平衡，天然氣的密度等；熱驅(qū)動(dòng)制冷系統(tǒng)的開發(fā)；高熱流密度的換熱設(shè)備及系統(tǒng)的研究與開發(fā)。

國家自然科學(xué)基金(51176166)資助項(xiàng)目。(The project was supported by the National Natural Science Foundation of China(No. 51176166).)

2013年12月2日

About the corresponding author

Han Xiaohong(1973－)，female，Ph.D./associate professor，Institute of Refrigeration and Cryogenics，Zhejiang University，(0571)87953944，E-mail:hanxh66＠zju.edu.cn.Research fields:alternative refrigerant research(vapor-liquid phase equilibrium，solubility and miscibility between refrigerant and lubricating oil，compatibility between refrigerant and material，density，cycle performance，and viscosity，heat transfer characteristics of new refrigerants，etc)；structural optimization and material of the heat exchanger；thermophysics of natural gas；new low-grade，heat-driven refrigeration systems；heat sink and air-conditioning systems for the high heat-flux space，etc.