如图1是刮削法制取冰浆装置的系统模型，其核心组件:压缩机、蒸发器(冰晶器)、冷凝器、热力膨胀阀及干燥过滤器与油分离器。在冰晶器内得到冰浆并通过泵的作用运送至储冰装置中；在冷凝器中，制冷剂将所得到的热量交换给冷凝器外壁的冷凝液体，最后通过热力膨胀阀的节流作用重新回到冰晶器中。

图1. 刮削法制冰系统

根据上述制冷剂的循环过程中绘制出制冰系统的温熵(T-S)图和系统的压焓图如图2和图3所示，过程1-2中制冷剂在冰晶器内等温定压吸热过程(焓增)，过程2-3中制冷剂在压缩机内的压缩升温过程(焓增)，过程3-4中高压制冷剂在冷凝器等温定压放热过程(焓减)，过程4-1中制冷剂在热力膨胀阀不可逆的绝热膨胀降压过程(等焓) [ 10 ]。

工况4处的制冷剂状态为饱和液态工质，工况2处的制冷剂状态为饱和气态工质，整个制冰系统的简化模型如图4所示。

对整体系统建立㶲分析根据上图所示的制冰系统简化得到的热力学模型和系统的温熵图、压焓图可知，图中各状态点的㶲值：

图2. 系统的温熵图

图3. 系统的压焓图

图4. 制冰系统的简化模型

e x . i = h i − h 0 − T 0 ( S i − S 0 ) (1)

式中 T 0 为环境温度，K；㶲值是在相对于环境条件进行相关数据计算的， h i 是图3中压焓图各点处的焓值KJ/kg， h 0 是环境下的焓值，KJ/kg； S i 和 S 0 分别是各点的熵和工质在环境下的熵值，KJ/kg.K。

制冰系统的制冷量：

q L = h 2 − h 1 (2)

q L 为系统的制冷量，KJ/kg。

压缩机耗功：

w c = h 3 − h 2 (3)

w c 为压缩机耗功，KJ/kg。

制冰系统的制冷系数为系统的制冷量和压缩机耗功的比值：

ε = q L w c (4)

ε 为系统的制冷系数。

由热力学第二定律可知，制冰系统的冷量㶲：

e x , q L = q L ( 1 − T 0 T L ) (5)

T L 为低温热源温度，K；系统的输出为冷量，其值为负， e x , q L 为系统的冷量㶲，KJ/kg。

通过系统㶲效率计算公式得到：

η e x = e x , q L w c = e x , q L q L q L w c = e x , q L q L ε (6)

式中 η e x 为系统的㶲效率。

对系统各组分进行㶲平衡分析，系统各组分的㶲如图4所示。

1) 冰晶器的㶲损分析

通过㶲守恒得到冰晶器的㶲损为冰晶器的进口处的㶲值和出口处的㶲值以及冷量㶲的差值：

e x , L , E V = e x 4 − e x 1 − e x , q L (7)

式中 e x , L , E V 为冰晶器的㶲损量，KJ/kg。

冰晶器的㶲效率及㶲损系数：

η e x , E V = e x , q L e x 4 − e x 1 = 1 − e x , L , E V e x 4 − e x 1 (8)

ξ E V = e x , L , E V w c (9)

式中 η e x , E V 为冰晶器的㶲效率， ξ E V 为冰晶器的㶲损系数。

2) 压缩机的㶲损分析

通过对制冰循环系统的㶲简化分析得到压缩机的㶲损：

e x , L , C = w c + e x 1 − e x 2 (10)

式中 e x 2 为压缩机出口处的㶲值。

压缩机的㶲效率及㶲损系数：

η e x , C = e x 2 − e x 1 w c = 1 − e x , L , C w c (11)

ξ c = e x , L , C w c = 1 − e x 2 − e x 1 w c (12)

式中 η e x , C 为系统的㶲效率， ξ c 为压缩机的㶲损系数。

3) 冷凝器的㶲分析

在忽略制冷剂流动阻力时，一般不考虑冷凝器的㶲效率，只考虑广义㶲损系数，冷凝过程中放出的热量㶲：

e x , q H = q H ( 1 − T 0 T H ) = ( H 2 − H 3 ) ( 1 − T 0 T H ) (13)

式中 e x , q H 为冷凝器向外界释放的热量㶲。

q H 为冷凝器向环境放出的热量KJ/kg∙K， T 0 为高温热源温度，K， T H 为外界环境温度，K。

冷凝器的㶲效率：

η e x , E V = e x , q H e x 2 − e x 3 = 1 − e x , L , C O N e x 2 − e x 3 (14)

式中 e x 3 冷凝器出口处的㶲值。

4) 膨胀阀的㶲损分析

制冷剂在膨胀阀中是节流过程，节流前后焓值相同，如图3中的压焓图，即 h 4 = h 1 ，膨胀阀的㶲损：

e x , L , T H = e x 3 − e x 4 = h 3 − h 4 − T 0 ( s 3 − s 4 ) = T 0 ( s 4 − s 3 ) (15)

e x 4 为膨胀阀出口处的㶲值，KJ/kg， s 3 和 s 4 分别为3、4点处的焓值，KJ/kg∙K。

膨胀阀的㶲效率及㶲损系数：

η e x , T H = e x 4 e x 3 (16)

ξ T H = e x , L , T H w c (17)

η e x , T H 为膨胀阀的㶲效率， ξ T H 为膨胀阀的㶲损系数。

蒸发汽化后的制冷剂，在冷凝器内的一定压力的条件下被冷凝成液体时相应的饱和温度被称为冷凝温度，冷凝温度的高低取决于实验过程的环境条件和压缩机允许的排气压力和温度；冷凝温度过高(40℃以上)会使得冷凝效果达不到预期；同时，如果冷凝温度过低(0℃及以下)，参加循环的制冷剂工质大大减少，换热效率低，因此在制冷剂为R410A，环境温度2℃，蒸发温度为−25℃时，设计冷凝温度分别为10℃、20℃、30℃、40℃进行实验，表1是冷凝温度下系统的各个参数。

表1. 冷凝温度下的系统参数

冷凝温度变化时，各部分㶲损及整体的㶲效率如图5所示。

在冷凝温度升高时，系统的制冷量逐渐降低，而压缩机的做功增加，使得制冷系数和㶲效率降低；在冷凝温度有10℃提高到40℃过程中，冰晶器的㶲损系数降低了7%，压缩机耗功从37 kJ/kg增加至65 kJ/kg，而冰晶器的㶲损都在5 kJ/kg左右；在冷凝温度升高时，冷凝器的㶲损系数从20%提高到42%，在冷凝温度达到30℃以上时，冷凝器的㶲损系数保持在42%左右，冷凝温度过高，制冷剂工质在冷凝器内换热不充分，导致了㶲损系数较高，因此降低冷凝温度可以有效的提高系统的效率。

图5. 冷凝温度变化时的㶲效率及各部分㶲损系数

制冰系统的实验环境是影响制冰系统效率的重要因素，郑州市地处中原地区，一年四季分明，将郑州市选择为环境温度的实验条件能够较为科学的进行相关实验，并且能够减少偶然因素引起的误差，初春和秋末的平均气温为22℃，夏季温度较高平均温度为32℃，初冬和冬末的平均温度为12℃和2℃ [ 11 ]；因此在冷凝温度40℃，制冷剂为R410A，蒸发温度为−25℃时，设计环境温度分别为2℃、12℃、22℃、32℃进行实验，环境温度变化时系统的㶲效率及各部分的㶲损系数如图6所示。

图6. 环境温度变化时的㶲效率及各组分㶲损系数

环境温度从2℃升高至32℃时，冰晶器的㶲损系数从7.7%上升到34.8%，而冷凝器的㶲损系数从43%降低至13%；制冷剂工质在冰晶器内蒸发吸热，经过压缩机后在冷凝器内释放热量进行冷却，环境温度主要是通过影响换热过程对㶲损系数产生影响，在2℃至32摄氏度时，提高环境温度可以使㶲效率提升了约2%，环境温度对系统的效率影响较小。

目前制冷剂的种类选择已经成为制冰系统以及制冷系统所面对的重要问题，R22制冷剂是目前使用最广泛地制冷剂，属于HCFC类制冷剂，但是R22制冷剂会对环境造成极大的破坏，因此寻求合理的制冷剂R22的替代物是当前的研究热点之一 [ 12 ]。R410A制冷剂是一种新型环保的制冷剂，是两种烷的混合物(R32和R125)，其不会破坏臭氧层，R410A制冷剂的GWP比R22的GWP高了16%，性能系数COP大约提升8%，因此R410A被视为R22的中长期的制冷剂替代物 [ 13 ]。R407C是一种三元非共沸混合物，其 ODP值为0，GWP值较高，它是直接替代R22的合适工质，可不换压缩机，但R407C的主要缺点是传热性能较差 [ 14 ]。R404A制冷剂可以在一定程度上替代R22，但其使用成本过高，不易大量使用。

R1234ze属于氢氟烯烃制冷剂，ODP为0，R1234ze的GWP值远小于其他几种工质的值，R1234ze制冷剂的不足之处是该制冷剂的单位制冷量较低，而R152A制冷剂具有较高的单位容积的制冷量和能效比，但是制冷剂R152A具有较大的可燃性，可能存在较大的安全隐患 [ 15 ] [ 16 ] [ 17 ] [ 18 ] [ 19 ]。因此目前有学者提出使用混合工质NCUR01 (R1234ze:R152A = 40:60)，可以使得含有40%的R1234ze不燃组分加入到易燃组分R152A中去，能够有效地提升新型制冷剂的安全性能，而且NCUR01制冷剂ODP值也为0，对环境影响较小，是可以替代R22制冷剂的可行性新型制冷剂 [ 20 ]。因此在冷凝温度40℃，环境温度2℃，蒸发温度为−25℃时，设计制冷剂种类分别为R410A、R404A、R407C、R1234ze、R152A、R1234yf进行实验，表2是不同制冷剂时的系统参数。

表2. 不同制冷剂的系统参数

制冷剂工质变化时，系统的㶲效率及各部分的㶲损如图7所示。

图7. 制冷剂变化时的㶲效率及各部分的㶲损系数

当制冷剂为R410A、R407C和混合工质NCUR01时，系统制冷系数在2.2左右，㶲效率均为18%左右，适合作为制冰系统的工质。冰晶器的㶲损系数波动幅度只有2%左右，不同的制冷剂种类对冰晶器的㶲损影响不大。而新型的混合制冷剂工质NCUR01在使用过程中冷凝器的㶲损系数只有20%，相较于大多数制冷剂这一数值较低，但压缩机㶲损系数较高达到35%，这是由于新型制冷剂在实际应用过程中对压缩机的要求更高，因此在使用这一制冷剂时要注意与压缩机匹配。

制冷系统的蒸发温度是影响制冷系统运行效率的关键因素，蒸发温度是工质在工作压力下的沸点。蒸发温度的高低决定了蒸发压力的大小，同时影响着压缩机在实际工作工程中的效率，因此在冷凝温度40℃，制冷剂为R410A，环境温度2℃，设计蒸发温度分别为−25℃、−15℃、−5℃、5℃，表3是蒸发温度变化时系统的参数。

表3. 蒸发温度下的系统参数

蒸发温度变化时，系统的㶲效率及各部分的㶲损如图8所示。

图8. 蒸发温度变化时的㶲效率及各部分的㶲损系数

在蒸发温度为−25℃冷凝器的㶲损系数为43.2%，造成了多数能量的损失；在蒸发温度为5℃时，冰晶器的㶲损系数达到了13.4%，㶲损系数最高，蒸发温度升高时，冰晶器内的蒸发吸热过程收到负面影响，降低其换热效率；同时较高的蒸发温度会使压缩机的排气温度升高，会影响压缩机的正常工作，蒸发温度同时也不能过低，蒸发温度−25℃以下时，系统的工作效率低下，不利于相关实验的开展。

由于单因素实验存在一定的局限性，因此为确保实验的可靠性，在单因素的基础上选取较为合适的实验因素水平进行相关的正交试验设计，在单因素的基础上对制冷剂进行筛选，NCUR01，R410A以及R407C，这三种制冷剂安全环保，且实验过程中，系统的㶲效率和制冷系数较高，因此制冷剂选择分别为NCUR01、R410A、R407C；冷凝温度是将制冷剂工质由饱和气态冷却为液态的温度，冷凝温度过高(40℃或50℃以上)会使得冷凝效果达不到预期；同时，如果冷凝温度过低(0℃及以下)，会有大量制冷剂积存于冷凝器中，实际参加循环的工质大大减小，因此冷凝温度选择为10℃、20℃、30℃；环境温度参考第三节中的郑州市一年四季的平均温度，不考虑冬季极低的(接近0℃及以下)温度水平，因此环境温度的选择为12℃、22℃、32℃；蒸发温度的高低决定着蒸发压力的高低，同时这一参数影响着压缩机的工作状况，蒸发温度一般比冷却水的温度低10℃，因此蒸发温度选择为−25℃、−15℃、−5℃。通过对单因素实验结果以及相关文献的分析，发现冷凝温度在40℃以上以及10摄氏度以下会使制冷系统的工作效率大幅度降低，因此选择冷凝温度为10℃、20℃、30℃；选择制冷剂NCUR01、R410A、R407C是在多种制冷剂在同种工况下进行比较，选择了三种制冷效率较高的制冷剂种；选择环境温度是根据郑州市一年内各季度平均温度进行选择，制冷系统在温度极低的情况下基本不会选择使用，因此选择的环境温度为12℃、22℃、32℃，蒸发温度的高低决定了系统内蒸发器的效率，因此通过单因素分析得到使得系统效率较高的温度为−25℃、−15℃、−5℃，上述各因素及其相应的水平可以使得正交实验比较充分的进行，对上述各因素及水平进行四因素三水平正交实验设计，抽签数据如表4所示。

表4. 抽签数据

通过上述抽签表得到混合水平的正交实验设计表L₉(3⁴)，因实验的评价指标为系统的㶲效率和制冷系数，是双评价指标，因此评价指标采用模糊数学中的隶属度函数对该实验结果进行综合评定，且两种评价指标的权重各占50%，正交试验的极差分析表5所示。

通过正交实验的综合评分数据得到各因素对实验结果的影响主次分别是：A (制冷剂种类)，B (冷凝温度)，D (蒸发温度)，C (环境温度)，通过计算得到系统的最优实验方案是：制冷剂为NCUR01、冷凝温度10℃、蒸发温度−5℃、环境温度12℃，在这个实验条件下，系统的制冷系数ε = 4.97，㶲效率h = 47.1%。由正交试验的极差分析得到，在制冷剂发生改变时，系统的综合评分变化幅度为1.5，会对结果产生较大的影响；提高冷凝温度时，系统的综合评分降低了1.1，且冷凝温度越高，综合评分的结果越小，因此在进行实验过程中，应尽可能的降低冷凝温度来提高系统的制冷系数和㶲效率；而提高蒸发温度时，系统的综合评分提升了0.7，蒸发温度对实验结果的影响比前两者带来的影响较弱；当环境温度从12℃升高至22℃时，系统的综合评分降低了0.3，但是在环境温度从22℃到32℃时，综合评分只有0.01左右的波动，环境温度对实验结果产生的影响最小。

表5. 正交试验的极差分析