大型结冰风洞制冷系统蒸发压力预测与降温过程优化研究_张平涛.pdf

2023 年第 1 期总第 251 期低温工程CRYOGENICSNo.1 2023SumNo.251大型结冰风洞制冷系统蒸发压力预测与降温过程优化研究张平涛1王文瑄2郭向东1吕宏宇2陈良2侯予2(1中国空气动力研究与发展中心结冰与防除冰重点实验室绵阳 621000)(2西安交通大学能源与动力工程学院西安 710049)摘要:构建了液氨蒸发器热力学仿真模型,进行了蒸发压力以及风洞降温的预测与验证。采用BP 神经网络算法对试验数据和仿真数据进行训练,建立了氨制冷系统蒸发压力快速预测模型,分析了试验高度和试验风速对风洞降温速度的影响,提出了变风速降温控制方式,对比分析了不同降温控制过程的降温时间及系统能耗。结果表明:对于高空环境试验,先进行风洞减压再降温,可以有效缩短降温时间、减小系统能耗;试验风速 75 m/s 时风洞系统降温速度最快,高于 75 m/s 时可以通过降温过程变风速控制方法降低系统能耗。关键词:结冰风洞 制冷系统 降温过程 蒸发器中图分类号:TB651,TB663文献标识码:A文章编号:1000-6516(2023)01-0072-07收稿日期:2022-08-06;修订日期:2023-02-11基金项目:国家科技重大专项(HT-J2019-III-0010-0054)。作者简介:张平涛,男,38 岁,硕士,工程师。通信作者:陈良,男,38 岁,博士,教授。Study on evaporation pressure prediction and cooling-down method ofrefrigeration system in a large icing wind tunnelZhang Pingtao1Wang Wenxuan2Guo Xiangdong1Lv Honggu2Chen Liang2Hou Yu2(1Key Laboratory of Icing and Anti/De-Icing,China Aerodynamics Research and Development Center,Mianyang 621000,China)(2School of Energy and Power Engineering,Xian Jiaotong University,Xian 710049,China)Abstract:A thermodynamic simulation model of liquid ammonia evaporator is established,and the evaporation pressure and wind tunnel cooling are predicted and validated against experi-mental data.BP neural network algorithm is used to train the test and simulation data,and a fastprediction model of circulation tank pressure of the ammonia refrigeration system is developed.The influence of the test altitude and wind speed on the cooling rate of the wind tunnel is ana-lyzed,and a variable-wind-speed cooling control method is proposed.The cooling time and energyconsumption of different cooling control processes are compared and analyzed.The results showthat the fastest cooling and relatively lower energy consumption can be achieved when the windspeed of 75 m/s is selected as the wind speed condition of the test section during the coolingprocess.Key words:icing wind tunnel;refrigeration system;cooling process;evaporator第 1 期大型结冰风洞制冷系统蒸发压力预测与降温过程优化研究1引言作为飞行器结冰与防除冰试验研究的基础设施平台,大型结冰风洞是中国航空工业发展的重要保障。制冷系统是结冰风洞的重要组成部分,通过热交换器冷却气流为试验提供低温环境。中国已建成的大型结冰风洞用液氨直接蒸发制冷,风洞气流场温度主要取决于液氨蒸发压力和流量1。液氨由多台定频泵从低压循环桶供给至风洞蒸发器,液氨流量通过氨泵启停实现分级调节,液氨蒸发压力的控制是风洞气流场温度及制冷系统负荷调节的关键,与低压循环桶压力相关。为了满足飞行器不同试验工况要求,需要在试验过程中对制冷系统蒸发压力进行快速预测,实现风洞气流场温度快速准确调节和系统能耗降低。目前,大型制冷系统控制参数预测方法已在大型空调系统等应用领域得到了充分研究。Kusiak 等人2通过动态神经网络算法建立了空调系统的能耗预测模型,优化供气温度和压力设定值,与传统的控制策略相比,所提出的模型节省了高达 30%的能耗。王炳文3以夏热冬暖地区某商场中央空调冷源系统为研究对象,以大量冷源系统实际运行数据为驱动,提出了冷源系统运行参数优化的数据挖掘方法。周璇4提出基于支持向量回归机的空调逐时负荷滚动预测算法,解决空调系统节能优化控制的问题。WeiXiupeng5等提出一种数据驱动的方法,用于优 化HVAC 系统的总能耗。Zeng Yaohui6等建立了典型气象年气象参数条件下某办公建筑的逐时冷负荷的支持向量机和 BP 神经网络预测模型,比较了两种模型的预测精度。上述研究表明人工智能算法建立制冷系统运行预测模型是进行制冷系统运行优化的有效手段,现有研究主要针对大型空调系统开展,而大型结冰风洞制冷系统制冷温度跨度更大、负荷变化更加剧烈,需要针对风洞试验工况要求,进一步结合热力学仿真开展主要控制参数快速预测方法研究。本研究建立了制冷系统蒸发压力热力学预测模型,结合BP 神经网络方法实现了不同蒸发器出口气流温度、试验模拟高度、试验风速等试验参数下的低压循环桶蒸发压力快速预测。2制冷系统蒸发压力热力学预测模型2.1蒸发器模型控制方程本研究的大型结冰风洞系统流程图如图 1 所示。针对热交换器部分,建立一维流动换热模型,连续方程、能量方程和压降计算方程分别为7-8:t+(u)x=0(1)(h)t+(uh)x=(2)dpdx+f(u2)2d(3)式中:为换热流体密度,kg/m3;u 为流体平均速度,m/s;为单位换热量,W/m3;h 为焓值,J/kg。换热流体与管壁间的热流密度以及压降使用传热关联式进行计算。图 1大型结冰风洞系统流程图Fig.1Schematic of an icing wind tunnel system换热器低温管路两相流动传热计算包括压降特性与换热特性6,9。在压降计算部分,由于本研究在风洞试验温度范围内氨气液粘度比值均小于 35。根据 Whalley10所提出的 当 L/G 1 000 时,使 用Friedel 关联式所提出的分相模型11,研究中首先计算单相流体的压降,然后通过液相系数 2LO修正因子计算两相压降。Pf=22LOLfLOG2/diL(4)2LO=A1+3.24A2A3Fr0.045We0.035(5)A1=(1-x)2+x2LfGOGfLO()(6)A2=(1-x)0.224x0.78(7)A3=LG()0.91GL()0.191-GL()0.7(8)式中:pf为压降,kPa;L 为流过的管道长度,m;2LO为液相系数;A1、A2和 A3为两相影响因子;x 为干度;为动力粘度,Pas;下标 L 为液相,G 为液相质量流量,(kg/m2s)。37低温工程2023 年采用 Blasius 公式计算换热器管内氨液紊流流动对应的单相摩擦因子。采用叠加法(Chen s mod-el12)进行管内两相流动传热系数计算。在管内流动沸腾的情况下,核态沸腾部分的传热系数(S)被抑制,而对流蒸发的传热系数(E)被加强。E=2.35/(1/Xtt+0.213)0.736(9)S=1/(1+2.53 10-6 Re1.17)(10)Xtt=1-xx()0.087 5LG()0.125GL()0.5(11)h=FhL+ShNB(12)式中:hL由 Dittus-Boelter 公式计算,hNB由 Cop-per 方程式计算。大型换热器的高度落差造成各模块入口氨液过冷度不同,高度落差越大过冷度越大。在这种情况下,需要考虑管内液体单相流动换热的情况,单相流动摩擦因子和换热系数与上述叠加法中的单相部分计算方法一致,即分别由 Blasius 公式和 Dittus-Boelter公式计算。2.2传热器空气侧流动传热计算关联式空气侧为强制对流换热,可按照翅片管式换热器空气侧流动公式进行计算13,其中摩擦压降的计算公式为:P=2LfG2/dc(13)f=0.026 7ReF1dcPtPl()F2FPdc()F3(14)式中:Fp为翅片间距,m;Pt为翅片横向间距,m;Pl为翅片纵向间距,m;dc为空气侧水力直径,m;Redc为空气侧雷诺数。传热系数计算方面还需要考虑翅片对换热的强化作用14。hs=oho=oj (VCp)/Pr2/3(15)j=0.086ReP3dcNP4FPdc()P5FPDh()P6FPPt()-0.93(16)式中:o为强化因子,j 为换热因子,Dh为空气侧水力直径,m;N 为管排数。强化因子计算详见文献14。2.3模型验证图 2a 中针对 10 组试验工况,利用制冷系统蒸发压力热力学模型进行低压循环桶压力预测,使用试验数据进行模型验证。结果表明低压循环桶压力预测值与试验数据均方根偏差为 0.08 105Pa,最大偏差为 0.13 105Pa。将 3 组工况的仿真模型结果与降温试验数据进行对比,模拟降温趋势与试验降温趋势符合,对比结果如图 2b 所示。图 2计算模型验证Fig.2Validation of calculation model3基于 BP 神经网络模型的蒸发压力预测BP(back propagation)神 经 网 络 是 1986 年 由Rumelhart 和 McClelland 为首的科学家提出的概念,是一种按照误差逆向传播算法训练的多层前馈神经网络,是目前广泛应用的神经网络模型。BP 神经网络算法以误差平方为目标函数、采用梯度下降法计算目标函数的误差最小值,其计算流程主要包括了两部分:正向结果计算和反向误差传递。BP 神经网络的隐含层节点数对 BP 神经网络预测精度有较大影响:当节点数太少时,网络不能很好的学习,需要增加训练次数,训练的精度也受影响;当节点数太多时,训练时间增加,网络容易过拟合。最佳隐含层节点数选择参考如式(17),并在范围内取训练最优值:l(m+n)+a(17)式中:l 为隐含层节点数;m 为输出层节点数;n为输入层节点数;a 为 010 之间的常数。47第 1 期大型结冰风洞制冷系统蒸发压力预测与降温过程优化研究通过制冷系统蒸发压力热力学预测模型获得了数据样本,结合试验数据搭建 BP 神经网络算法模型。数据样本来源包括 50 组试验数据和 50 组模型预测数据,以正交设计的 25 组模型预测数据作为测试集。模型训练及预测均在 LabVIEW 平台实现,可在软件中补充试验数据,对模型进行重新训练。本研究选用均方误差(MSE)作为 BP 神经网络模型训练的评价标准。在训练次数相同的条件下,进行不同隐藏层神经元个数对训练结果的影响分析,如图 3a 所示。结果表明训练误差随着神经元个数变化呈现先减小后增大的趋势,存在过拟合情况和最优隐藏层神经元个数。因此,在预测模型中给定训练次数,通过模型预测误差最小化实现隐含层节点数的自动寻优。现有模型的自动寻优训练结果如图 3b 所示。结果表明 80%的预测值的绝对偏差小于 0.1 105Pa,满足结冰风洞制冷系统蒸发压力控制需求。图 3BP 神经网络优化Fig.3Optimization of BP neural network4风洞降温分析风洞降温受试验高度和试验风速的影响,降温时间和降温能耗在不同条件下会产生差异,针对试验高度和试验风速对风洞的影响进行研究,提出降温控制方案。4.1试验高度对风洞降温的影响在试验风速分别为 37 m/s、93 m/s 和 150 m/s,低压循环桶压力为 2.45 105Pa,外界温度为 7 的初始条件下,研究了试验高度分别为 0、2、3、5、7 km条件下的风洞降温情况,忽略试验中抽气减压的时间。仿真计算结果如图 4 所示,可以看出随着高度的升高,风洞回路中的稳态温度呈现降低趋势。图 4不同高度下的蒸发器出口温度降温曲线Fig.4Evaporator outlet temperature at different altitudeduring the cooling-down process and wind speed风机功率-高度-速度关系见图 5,结果表明在相同的风速情况下,试验高度越高,风机功率越低,风洞系统的热负荷更小。以降温到-7 为标准统计各试验高度下的降温时间,结果统计如图 6 所示。图 6表明降温时间与试验高度变化呈负相关,随着高度越高,洞体内部压力下降,系统热负荷下降,降温所需时间变短。高度变化会改变空气压力影响换热,高度升高也会减小风机功率,两者共同影响造成了降温时间变化程度越来越小。考虑风洞降温时的制冷系统能57低温工程2023 年耗和风机能耗,各高度下降温至-7 时的能耗结果统计如图 7,结果表明高度越高,系统降温系统能耗随之减小。当试验风速较高时,导致风机能耗及系统热负荷显著上升,试验高度的影响更加明显。因此,对于高空环境试验,先进行风洞降压再降温,可以有效缩短降温时间、减小系统能耗。图 5主试验段风机在不同试验高度条件下的功率随风速的变化曲线Fig.5Variations of fan power with wind speedunder different tested altitude图 6不同试验风速下的风洞降温时间随试验高度的变化Fig.6Variations of cooling-down time with testedaltitude at different wind speed4.2试验风速对风洞降温的影响在蒸发器试验高度分别为 0 km、3 km 和 5 km,低压循环桶压力为 2.45 105Pa,外界温度为 7 的初始条件下,研究了试验风速分别为 37、56、75、100、130、150 m/s 条件下的风洞降温情况。仿真计算结果如图 8 所示,可以看出随着风速的升高,风洞回路中的稳态温度呈现降温趋势,且降温时间变长。以降温到-7 为标准统计各风速下的降温时图 7不同试验风速下系统能耗随试验高度的变化Fig.7Variations of system power with testedaltitude at different wind speed图 8不同风速下蒸发器出口温度降温规律,试验高度Fig.8Temporal variations of temperature at evaporatoroutlet under different wind speed and tested altitude间。从图 9a 可以看出在低于 100 m/s 的风速下,风洞降温时间变化不明显,降温时间随着风速的增大而降67第 1 期大型结冰风洞制冷系统蒸发压力预测与降温过程优化研究低,而在风速大于 100 m/s 时降温时间变化相反。同时对各风速在对应时间下消耗的能量进行研究,计算结果如图所示,结果表明总消耗能量随着风速的增大而增大,当风速大于 100 m/s 时,能耗快速上升,这是由于高风速下风机功率及其产生的热负荷快速上升所导致的。根据风机功率-高度-速度关系图,可以看出在相同高度下风机功率与风机速度呈现指数级上升趋势,极大地增加了风洞的热负荷。图片同时显示了试验高度对降温时间的影响情况,与前一节中的结论一致。同时统计各风速在对应时间下消耗的能量,结果如图 9b 所示。结果表明功率变化规律与降温时间变化规律存在明显差别,总消耗能量随着风速的增大而增大,当风速大于 100 m/s 时,能耗快速上升,这是由于高风速下风机功率及其产生的热负荷快速上升所导致的。由此可见,选取75 m/s 风速作为降温过程试验段风速条件时,可以实现最速降温,同时维持较低能耗。图 9试验风速优化Fig.9Optimization of experimental wind speed根据上述结果,继续进行变风速工况降温的研究,针对风速大于 75 m/s 的工况(试验高度 5 km),首先在 75 m/s 降温到-7 的情况下,再切换到指定风速进行降温。仿真计算结果如图 10 所示。结果图 10不同试验风速下变风速降温过程Fig.10Temperature variations during the cooling-downprocess with variable-wind speed designfor different tested speed77低温工程2023 年表明较原工况变风速工况在前段降温较快,待降到-7 提高风速后,温度升高,降温稳定后变风速工况与原工况降温趋势逐渐一致。同样对各工况降到-7 (变风速工况调整风速后降温到-7 )下消耗的能量(制冷系统功率和动力系统功率)进行统计研究,计算结果见图 11。结果表明变风速对降温时间影响不大,但是大幅减少了能耗,且随着试验风速的升高,能耗降低幅度增大。上述结果为经济性最优降温策略提供了指导方案。图 11变风速工况与原工况系统能耗随试验风速的变化Fig.11Variations of system power consumption ofvariable-speed and fixed-speed cooling-down processes5结论本研究针对大型结冰风洞复杂试验工况及制冷系统蒸发压力快速准确调节问题,建立了蒸发压力热力学预测模型,使用试验数据完成了模型验证。采用BP 神经网络模型,在 LabVIEW 平台开发了制冷系统低压循环桶压力预测软件,实现了制冷系统调节过程低压循环桶压力快速预测。结果表明:(1)试验风速 75 m/s 时风洞系统降温速度最快。并且 当 试 验 风 速 分 别 为 150 m/s、93 m/s 和37 m/s时,试验高度为 7 km 的降温时间相比于高度为 0 km 时的降温时间能缩短 36%、20%和 11%。(2)对于高空环境试验,先进行风洞降压再降温,可以有效缩短降温时间、减小系统能耗。(3)试验风速选择小于 75 m/s 时风洞系统能耗最低。当试验高度为 5 km 时,试验速度为 75 m/s 的降温能耗相比于速度为 150 m/s 时的降温能耗能降低 83%。(4)针对 100 m/s、112 m/s、130 m/s 和 150 m/s的风速进行变风速工况研究,发现降温能耗分别节省了 27%、36%、50%和 45%。因此对于高风速试验,可以通过降温过程变风速控制方法降低系统能耗。参考文献1 郭向东,张平涛,赵献礼,等.大型结冰风洞热流场符合性验证J.实验流体力学,2020,34(5):79-88.Guo Xiangdong,Zhang Pingtao,Zhao Xianli,et al.The complianceverification of thermodynamic flowfield in the large icing wind tunnelJ.Journal of Experiments in Fluid Mechanics,2020,34(5):79-88.2 Kusiak A,Xu G.Modeling and optimization of HVAC systems using adynamic neural network J.Energy,2012,42(1):241-250.3 王炳文.基于数据挖掘的中央空调冷源系统运行参数优化研究D.广州:华南理工大学,2019.Wang Bingwen.Optimization of operating parameters for chiller plantof central air conditioning system based on data mining D.Guang-zhou:South China University of Technology,2019.4 周璇,杨建成.基于支持向量回归机的空调逐时负荷滚动预测算法 J.中南大学学报(自然科学版),2014,45(3):952-957.Zhou Xuan,Yang Jiancheng.An algorithm for hourly load rollingforecasting of air conditioning system based on SVR J.Journal ofCentral South University(Science and Technology),2014,45(3):952-957.5 Wei X,Kusiak A,Li M,et al.Multi-objective optimization of theHVAC(heating,ventilation,and air conditioning)system perform-ance J.Energy,2015,83:294-306.6 Zeng Y,Zhang Z,Kusiak A.Predictive modeling and optimization ofa multi-zone HVAC system with data mining and firefly algorithmsJ.Energy,2015,86:393-402.7 Anderson J D.Computational Fluid Dynamics:the Basics with Appli-cations M.New York;McGraw-Hill.19958 Mller-Steinhagen H,Heck K.A simple friction pressure drop corre-lation for two-phase flow in pipes J.Chemical Engineering and Pro-cessing:Process Intensification,1986,20(6):297-308.9 Wang H,Fang X.Evaluation analysis of correlations of flow boilingheat transfer coefficients applied to ammonia J.Heat Transfer Engi-neering,2016,37(1):32-44.10Beattie D,Whalley P B.A simple two-phase frictional pressure dropcalculation method J.International Journal of Multiphase Flow,1982,8(1):83-87.11Friedel L.Improved friction pressure drop correlation for horizontaland vertical two-phase pipe flow C.Proc of European Two-PhaseFlow Group Meet,Ispra,Italy,1979,18(2):485-491.12Chen J C.A correlation for boiling heat transfer to saturated fluids inconvective flow J.IndEngChemPrOcDesDev,1966,5(3):322-329.13Horvat A,Leskovar M,Mavko B.Comparison of heat transfer condi-tions in tube bundle cross-flow for different tube shapes J.Interna-tional Journal of Heat and Mass Transfer,2006,49(5):1027-1038.14钱颂文.换热器设计手册 M.北京:化学工业出版社,2002.Qian Songwen.Heat Exchanger Design Manual M.Beijing:Chem-ical Industry Press,2002.87