1. 引言
在可持续发展战略和能源环境问题的世界背景下,绿色能源的开发相对于一般能源的利用能极大降低环境污染。而吸附式制冷技术作为一种低品位能源驱动的绿色制冷技术给冰箱、空调、热泵等行业带来新的机遇。
因此,吸附式制冷设备逐渐成为研究热点,其研究领域主要包括制冷工质对、物理与化学吸附的热动力学特性、复合吸附剂的吸附机理及热动力学特性、吸附床技术以及吸附式制冷系统的热力学循环等 [1]。由于吸附床的传热传质能极大影响吸附式制冷设备的性能,绝大多数制冷设备的热力学分析通常建立吸附床的理论计算或数值模拟的基础上 [2] - [7]。此外,吸附式制冷设备的运行参数也是一个影响性能的关键因素 [8],在此基础上,控制水汽从吸附床导向冷凝端的通道结构特征或许也能一定程度地影响设备的运行效率。
在本文中,研究重点主要在于不同的水汽通道结构中,设备运行状态下的热力学特性。其中,通过观察通道结构的温度分布以及冷凝分量和冷凝端内换热工质的温度分布来对不同水汽通道结构中的吸附式制冷设备的热力学特性进行分析。
2. 制冷设备结构及工作原理
本文中研究的吸附式制冷设备结构如图1所示,主要的研究对象为解吸后的气态工质在设备的工作部件中发生的热交换。为了降低具备多样的热循环制冷系统的复杂程度,设备由机箱外壳构建接近真空的工作环境,设备内部的主要工作部件为吸附床,冷凝床。
Figure 1. Structure of an adsorption cooling system
图1. 吸附式制冷设备结构
其中,基于该吸附式制冷设备的工作原理,如图2所示,吸附床主要发挥吸附气态工质(水汽)和高温下解吸水汽的作用,而冷凝床则包含低温下冷凝水汽及高温下蒸发冷凝水达到制冷效果的两个主要功能。具体来说,吸附床及冷凝床之间利用气压差驱动水汽循环,并通过适当密度的管式翅片来达到热量的交换的目的。然而除了冷凝床中冷凝水汽发生的换热过程之外,用于隔离冷凝床与吸附床的水汽通道也会与在其表面冷凝后的水汽产生热交换,以上两个过程主要发生在吸附剂的等量温变过程,在该过程中,吸附剂所释放的热量主要为:
(1)
其中,
为吸附床释放的水汽质量,
为水汽的定压比热容,
为通道结构上冷凝的水汽质量分量,
为冷凝端换热结构上冷凝的水汽质量分量。
为吸附剂在通道结构上等量降温过程前后的温差,
为吸附剂在冷凝端等量降温过程前后的温差。
此外,整个工作环境中处于真空,热传导为主要换热方式,因此可以近似得到水汽通道表面的换热量为:
(2)
其中,
为水汽通道表面的导热系数,
为通道表面的有效换热面积,
为导热方向的温度梯度。
在管式翅片的换热结构中,管径与翅片的几何参数,如翅片高度、间距、厚度均会显著影响制冷设备的换热性能 [7],因此制冷功率(SCP)与性能系数(COP)同样会在不同热循环周期中产生浮动。与之相比,板式翅片换热结构 [9] 内吸附质与换热器的接触面积更大,可以更明显地提高设备换热性能,但是由于结构紧凑易于造成传质堵塞,影响设备的可靠性,因而管式翅片结构更普遍用于吸附式制冷设备。
3. 仿真方案
制冷设备的结构热设计与分析由于其复杂的热循环过程,通常建立在理论计算或数值模拟的基础上 [10] [11] [12],具体主要包括吸附工质对的热力学特性、设备制冷的基本热力学循环以及回质/类回质的热力学循环 [4]。本文主要围绕设备制冷的基本热力学循环开展热仿真分析,其中,吸附床与冷凝床间的通道结构为主要分析对象。
Figure 2. Working principle of the adsorption cooling
图2. 吸附式制冷原理
该热仿真分析主要基于一项通用的工程热设计软件FloEFD,其属于SIEMENS旗下并广泛应用于航天航空、机械船舶、制冷空调、消费电子产品以及汽车等工业领域 [13]。与其它CFD软件相比,具备网格自适应、流态自动识别求解以及MCAD系统兼容的特点,能够极大简化工程师对于CFD分析的流程,同时便于在原始CAD模型下增加结构复杂度。
由于该吸附式制冷设备模型结构较为复杂,并存在算力限制,在尽可能保证仿真结果准确性的前提下,对该模型进行简化处理:
1) 忽略设备内部密闭空间的辐射换热;
2) 忽略冷凝床的流道结构内换热工质的沿程损失与局部损失,同时流道内为充分发展流动的对流换热;
3) 降低模型几何复杂度,对复杂的几何特征如倒角,零部件配合等进行简化处理。其中,冷凝床内管式翅片换热结构简化为几字型矩形截面流道,吸附床内则等效为近似体积大小的恒温热源,解吸后的水汽由热源近通道表面释放,如图3所示。
此外,对应的热分析模型的边界条件根据吸附床与冷凝床不同的热循环要求设置如表1所示,其中重力为Z轴方向。同时,设备外壳材料默认为铝金属,密度为2688.9 kg/m3,导热系数为240 W/(m∙k);换热管道材料为铜,密度为8940 kg/m3,导热系数为387.7 W/(m∙k)。
网格质量与边界条件同样决定着热仿真的准确性,因此对于常规的热仿真方案来说,划分网格并判断网格的收敛性是一项重要的任务。然而FloEFD作为一款面向工程的仿真软件,其自带的网格划分功能可以满足工程案例中网格的精密要求 [12]。在该仿真方案中,网格的初始级别设置为2,并选用高质量的六面体网格进行网格自动划分,其中原始结构总网格数约为131,046,流体网格数约为81,918,固体网格数约为49,128,接触固体的流体网格约为53,983。
Figure 3. Simplified structure model of the cooling system for thermal analysis
图3. 简化的制冷设备热分析模型结构
Table 1. Parameters of the boundary condition
表1. 边界条件的参数设置
4. 计算结果与分析
该热设计仿真的目的为分析特定吸附式制冷设备中的水汽通道结构对其基本热力学循环的影响。因此,根据仿真方案,设计了如图3所示的三组通道结构,并对其开展热仿真分析。
对于吸附剂来说,在设备工作的循环过程中主要有等压解吸,等量降温,等压吸附以及等量升温四个主要循环过程 [1]。在该仿真结果的分析中,研究的换热形式主要发生在吸附剂等量温变的过程。不同的水汽通道结构会在不同程度上影响水汽进入冷凝器的效率,同时其壁面温度会随着在通道壁面冷凝的水汽含量变化而变化,从而可能影响设备整体的制冷效率,具体如图4及图5所示。
可以看到,当采用矩形孔板结构时,温度分布基本一致,高温集中在中心区域,且温度向两侧呈梯度分布。由于受重力影响,水汽相变释放的潜热更集中于通道结构的Z轴正方向,而最大冷凝的物质分量的位置却由于结构的不同而发生改变。由于该冷凝物质分量仅仅以通道结构为参考,因此无法在不同结构间进行冷凝质量大小的比较。特别地,当孔板结构变为单一的矩形通道时,温度分布发生变化,如图4(c)所示,且平均温度降低了约16%,而发生冷凝的水汽逐渐发散且主要集中于结构的X轴正方向,如图5(c)所示,这可能是由于水汽直接与冷凝管接触而减少了在水汽通道上相变潜热的释放。
Figure 4. Temperature profiles of 3 types of tunnel structure
图4. 三种通道结构温度分布
Figure 5. The amount of vapor condensed on the walls of 3 types of tunnel structure
图5. 三种通道结构壁面冷凝水汽分量
吸附剂在该热力学过程中处于等量降温的过程,释放的热量Q1被吸附设备内的换热结构吸收,如公式(1)所示。由于实际的真空工作环境,在分析过程中仅考虑热传导的换热形式。因此根据热传导的基本表达式,水汽通道结构的栅栏面积减少,有效换热面积减少,因此在通道结构上的吸收的相变潜热减少,如公式(2)所示。
结合公式(1)和公式(2),可以知道,在减少通道结构的有效换热面积后,其表面冷凝的水汽质量分量也会减少,因此更多热量导向冷凝端,管内的换热工质会吸收更多的热量导致温度上升,如图6所示。从仿真结果来看,尽管不采用孔板结构时或许更有利于水汽在冷凝管表面的相变过程,但会导致热量集中在通道结构的局部。由于仿真模型的简化,导致冷凝端的有效换热面积大幅度减少,吸附剂参与工作循环的质量因而被极大消减,导致较大部分的热量通过机箱外壳产生了损失。该部分热损失造成冷凝管进出口水温的变化不敏感,因此热量在通道结构的局部集中和冷凝管进出口水温变化在该结果中无法反映对吸附设备的制冷效率的影响。
由于设备的内部工作环境为低压真空,吸附床不断产生水汽使得设备内部气压升高,该过程为等量升温的过程。在大量水汽发生冷凝后,冷凝端气压降低,因此在等量降温的过程中,吸附床解吸的水汽会不断流向冷凝端。在流动过程中,水汽通道在此起到控制流量的作用。在该仿真过程中,通道结构均为较薄的矩形孔板,由于通道较短,对于水汽的分压和分流效果并不明显,因此结构上的温度变化也不显著,同时由于在低压环境中主要流动流体的为水汽,所以发生汽蚀现象 [14] 的概率也较低。值得注意的是,当采用矩形孔板结构时,冷凝管内进出口的工质温差无明显变化,这是由于大量的热损失导致的,因此无法通过进出口水温直接计算制冷设备的效率和性能。但对于低品位热源驱动的吸附式制冷设备,单位制冷量为优化设计中的主要参数,而由于该仿真控制管内换热工质流速,不考虑吸附与脱附时间,可以知道SCP (Specific Cooling Power)由制冷量与吸附器中吸附剂的质量决定的 [15]。在仿真结果中,制冷量可以从吸附剂在冷凝端有效冷凝的质量体现,因此比起矩形多孔的通道结构,单一矩形孔板的通道可能提供更大的制冷量,从而提高设备的单位制冷量。
Figure 6. Temperature profiles of heat transfer working fluid in the condensers of 3 types of tunnel structure
图6. 三种通道结构下冷凝管内换热工质温度分布
总的来说,尽管得到了初步的仿真结果,但并未平衡吸附质的质量流量,这容易导致最终结果与实际失真。同时,由于简化的热分析模型不包含吸附床及冷凝端的具体参数特征,如换热器结构、肋片间距、吸附剂和吸附质的质量等,导致冷凝端的换热面积不足并小于整体设备外壳的有效换热面积,从而对水汽的冷凝热吸收不充分,导致冷凝流道内的温差变化不明显,无法验证特定功率下的吸附式制冷设备的制冷效率受水汽通道结构的具体影响程度。因此还原具体的设备结构参数特征以及对热力学循环的瞬态分析将有利于探究水汽通道结构对设备性能的影响。
5. 结论
该报告针对特定的吸附式制冷设备设计了三组热仿真分析方案,主要研究对象为冷凝段与吸附床间的水汽通道结构。通过简化的热分析模型对不同通道结构的工作环境下进行了热力学仿真,计算了结构的表明温度变化以及冷凝物质分量,对比了对应通道结构对冷凝管内工质换热的影响。热分析模型可在增加复杂度的情况下,结合实验对其进行验证,验证过后的模型可用于对吸附式制冷设备的理论设计和分析。
基金项目
广州市基础与应用基础项目(No.202102021306 & No.202102080648);广州开发区国际科技合作项目(No.2020GH08)。
NOTES
*通讯作者。