王晓坡，杜家浩，陈建林，何茂刚

(西安交通大学能源与动力工程学院，热流科学与工程教育部重点实验室，陕西西安，710049)

吸收式制冷技术能够有效利用余热资源，一直是各界研究的热点。然而，传统的吸收式制冷工质对如H2O/LiBr 或NH3/H2O 等仍存在一些问题[1−3]，如NH3有毒且需要在发生器出口安装精馏装置，增加了系统的运行成本；LiBr溶液易结晶，对金属有腐蚀性且对装置密封性要求比较高。因此，寻找新型吸收式制冷工质对一直是吸收式制冷领域一个重要的研究方向。

由于离子液体具有较高热稳定性和极低饱和蒸气压，以离子液体为吸收剂组成的新型吸收式制冷工质对引起了广泛关注。MARTÍN等[4]比较了CO2与多种离子液体组成的工质对在吸收式制冷循环中的应用潜力，发现CO2/[bmpyrr]Tf2N的性能系数最高可达到0.55；KIM 等[5−6]分析比较了不同HFCs 与[bmim]PF6组成的工质对的性能，发现R134a/[bmim]PF6性能最优；WU 等[7−8]分析了不同HFOs 与[hmim]Tf2N 组成的工质对在吸收式制冷循环中的性能，发现工质对R1234ze(E)/[hmim]Tf2N的性能系数可达0.498；SUN 等[9−10]比较了R1234yf与[emim]BF4，[hmim]BF4，[omim]BF4和[hmim]Tf2N等组成的工质对的性能，发现R1234yf/[hmim]Tf2N的性能最好；张垚等[11]研究了R1234ze(E)与3 种离子液体组成的工质对的性能，发现R1234ze(E)/[omim]PF6的性能系数最大，可达到0.21。

上述研究表明，研究人员研究了以离子液体为吸收剂的工质对，获得了一些有价值的结论。然而，目前的研究主要集中在对某一种或几种工质对的性能分析上，而离子液体的种类繁多，进一步对比分析多种离子液体工质对的性能，对于开发新型吸收式制冷工质对具有重要意义。基于此，本文对烷烃类、烯烃类以及NH3等制冷剂与不同离子液体组成的19 种工质对在吸收式制冷循环中的性能进行分析，探究发生器出口温度、吸收器出口温度、蒸发温度和冷凝温度等参数对系统性能系数及㶲效率的影响。

1 热力学模型的建立

1.1 工质的名称及相关信息

本文研究的制冷剂和离子液体的名称及相关信息如表1所示。

表1 制冷剂和离子液体的名称及相关信息Table 1 Names and related information of refrigerants and ionic liquids

1.2 溶解度计算模型

由于离子液体的饱和蒸气压非常低，通常可以忽略，因此，当制冷剂与离子液体组成的混合物处于相平衡状态时，可认为气相中全部为制冷剂气体，即气相中制冷剂的摩尔分数y1=1。此时，相平衡方程可写为

式中：x1为液相中制冷剂的摩尔分数；γ1为制冷剂的液相活度系数；为制冷剂的饱和蒸气压；p为平衡时系统的总压；φ1为制冷剂的气相逸度系数，

式中：R为气体常数；B1和v1分别为制冷剂的第二维里系数和饱和液相摩尔体积。

制冷剂的液相活度系数γ1由非随机双液模型(NRTL)[12]计算得到：

式中：x1和x2分别为液相中制冷剂和离子液体的摩尔分数；τ12和τ21为二元交互参数。二元交互参数τ与温度T的关系式如下：

式(4)和(5)中，，，和通过文献[13−21]中的气液相平衡数据拟合得到，其结果如表2所示。

表2 制冷剂/离子液体工质对的NRTL模型系数回归结果Table 2 Regressed results of coefficients in NRLT model for refrigerant/ionic liquid working pairs

为了验证NRTL模型对溶解度实验数据的复现性，以R152a/[P(14)666]TMPP为例，图1给出了由NRTL方程计算的制冷剂摩尔分数与实验值之间的偏差(δARD)分布。由图1可见：NRTL方程计算值与实验测量值具有良好一致性，偏差均在5%以内，故可以利用所得到的NRTL方程分析循环的性能。

图1 NRTL模型预测值与实验值对比Fig.1 Comparison of predictive values from NRTL model and experimental data

1.3 循环模型

单效吸收式制冷循环的流程简图如图2所示。在分析时，假设系统处于稳态运行，忽略工质在管道和换热器中的热损失和压降，且出发生器和吸收器的蒸气和稀溶液均处于气液相平衡状态。

图2 单效吸收式制冷循环的流程图Fig.2 Flow chart of single-effect absorption refrigeration cycle

在进行循环的性能分析时，需要用到各状态点的焓。制冷剂/离子液体混合物的焓H可由下式计算：

式中：Hid为理想溶液的焓；HE为二元溶液的过量焓。

理想溶液的焓Hid定义为

式中：H1和H2分别为制冷剂和离子液体的焓。离子液体的焓H2用下式计算得到：

式中：h0依据国际制冷学会的标准选取；当T0=273.15 K 时(p0=psat(T0))，h0=200 J/g；cp,IL和vIL分别为离子液体的比定压热容和比体积。离子液体定压比热容根据GARDAS 等[22]提出的基团贡献法计算：

式中：A，B和D为离子液体中各基团贡献值加和求得的特征参数，具体数值可参考文献[22]。离子液体的比体积vIL由JACQUEMIN 等[23]提出的基团贡献法得到：

式中：vcation和vanion分别为阳离子和阴离子的比体积；Di，Ei和Fi分别为离子液体中各基团贡献值加和求得的特征参数，具体数值可以参考文献[23]。

二元溶液的过量焓HE通过以下式计算：

根据式(6)～(13)可以计算得到各点的焓。在此基础上，根据质量守恒和能量守恒，可以建立单效吸收式制冷循环中各部件的质能方程，得到各部件的热负荷Q。

对于发生器，

对于吸收器，

对于冷凝器，

对于蒸发器，

对于溶液换热器，

式中：Qg，Qa，Qc，Qe和Qshx分别为发生器、吸收器、冷凝器、蒸发器和溶液换热器的热负荷。对于溶液换热器，本文取其能效系数为0.8[24]。h为比焓，各下标分别对应于图2中的各状态点；mr为制冷剂蒸汽的质量流量；mw为进入发生器的稀溶液质量流量，

xw为制冷剂在稀溶液中的质量分数；xs为制冷剂在浓溶液中的质量分数。

单效吸收式制冷循环的性能系数用η表示，其计算公式为

当忽略溶液泵功时，循环㶲效率ECOP可以用下式计算：

式中：EQe和EQg分别表示收益㶲和代价㶲；Tref为环境温度，298.15 K；Tg和Te分别为发生温度和蒸发温度。

2 吸收器出口温度(ta)对溶液质量流量的影响

分析单效吸收式制冷循环的性能时，其计算工况如下：蒸发温度和冷凝温度分别为25 ℃和50 ℃，蒸发器出口工质过热5 ℃，冷凝器出口工质过冷5 ℃，循环制冷量为100 W。

图3所示为循环中进入发生器的稀溶液质量流量随着吸收器出口温度的变化趋势。从图3可见：溶液质量流量随吸收器出口温度增加而增加。吸收器出口温度增加，制冷剂在稀溶液中的质量分数降低，从而导致溶液质量流量随之增大。其中，工质对R143a/[P(14)666]TMPP，R143a/[hmim]Tf2N和R290/[hmim]FEP 的溶液质量流量较大，当吸收器出口温度由35 ℃增加到50 ℃时，上述3种工质对的溶液质量流量增加幅度分别为98.68%，91.53%和91.22%。

图3 溶液质量流量随着吸收器出口温度的变化趋势Fig.3 Trends of solution mass flow with absorber outlet temperature

3 吸收器出口温度(ta)对系统性能的影响

结合工业实际的需要，设定发生器出口温度为80 ℃。图4所示为系统的性能系数随吸收器出口温度的变化趋势。由图4可见：吸收器出口温度增加，引起制冷剂在浓溶液中的溶解度减小，浓溶液和稀溶液中的制冷剂质量分数的差相应减小，系统出现循环倍率单调增大，系统性能系数呈现单调减小的趋势。

此外，从图4还可见：在相同工况下，NH3/[omim]BF4和NH3/[hmim]BF4具有较高的性能系数，而含有[P(14)666]TMPP 的工质对的性能系数普遍偏低。此外，结合溶解度实验数据以及图4(a)和(b)可知，对相同的离子液体，碳链较短的氟化烷烃的溶解度更大，这有利于提升系统性能。同时，降低吸收器出口温度可以显著提升系统性能系数。

图4 系统性能系数随着吸收器出口温度的变化趋势Fig.4 Trends of coefficient of performance of system with absorber outlet temperature

4 发生温度(tg)对系统性能的影响

图5所示为系统性能系数随发生器出口温度的变化趋势。由图5可见：当发生器出口温度较低时(tg<80 ℃)，制冷剂在稀溶液中的质量分数不变，制冷剂在浓溶液中的质量分数减小，溶液的循环倍率减小，导致发生器需要提供的加热量减小，此时，性能系数随发生器出口温度增加而急剧增加。当发生器出口温度达到一定程度后(tg>80 ℃)，溶液质量流量随发生器出口温度的变化率逐渐减小，但此时发生器出口焓却逐渐增大，两者共同作用导致热负荷几乎不变，系统性能系数逐渐保持平稳。因此，在所研究的工质对中，NH3/[omim]BF4具有最高的性能系数，可达到0.62，R290/[hmim]FEP具有最低的性能系数。

图5 系统性能系数随着发生器出口温度的变化趋势Fig.5 Trends of coefficient of performance of system with generator outlet temperature

5 系统㶲效率的影响因素

本文以R161/[P(14)666]TMPP工质对为例，分别研究蒸发温度和冷凝温度对循环㶲效率的影响，结果如图6所示。不同蒸发温度时对应的机组运行所需最低热源温度(t8min)、峰值温度(tpeak)以及系统㶲效率的最大值(ECOPmax)如表3所示。

从表3可见：蒸发温度平均每升高2 ℃可导致系统㶲效率最大时对应的峰值温度(tpeak)降低6 ℃左右，系统㶲效率的最大值(ECOPmax)可以提高1.39倍。从图6(a)可见：随蒸发温度增加，㶲效率先急剧上升，㶲效率达到最大值时的“峰值温度”(图中箭头所指方向)后逐渐降低。

图6 蒸发温度和冷凝温度对系统㶲效率ECOP的影响Fig.6 Influence of evaporation temperature and condensation temperature on system exergy efficiency ECOP

表3 同蒸发温度对应的t8min，tpeak和ECOPmax数值Table 3 Values of t8min,tpeak and ECOPmax at different evaporator temperatures

从图6(b)可见：增加冷凝温度导致系统㶲效率ECOP沿着图中箭头所指方向逐渐降低，一方面，增加冷凝温度会增大溶液质量流量，导致发生器和溶液换热器等由于摩擦和压降等增加不可逆损失；另一方面，冷凝器本身由于存在温差也会增加传热㶲损失。增加冷凝温度会提高发生器对于热源温度的要求，如当冷凝温度由46 ℃升高到54 ℃时，发生器所需的最低热源温度将会由60 ℃升高到70 ℃左右。

6 结论

1)NH3较大的蒸发潜热导致NH3/[omim]BF4具有最高的性能系数，最高可达0.62。而含有离子液体[P(14)666]TMPP 的工质对的循环的性能系数普遍偏低。

2)增加吸收器出口温度导致系统性能系数减小，而提高发生器出口温度，性能系数先增加后逐渐平稳。

3)㶲效率随蒸发温度增加先急剧上升，㶲效率达到最大值时的“峰值温度”逐渐降低，而增加冷凝温度导致系统㶲效率逐渐降低。

4)增加蒸发温度或者减少冷凝温度都会降低发生器对于最低热源温度的要求。