换热系统优化中换热与阻力的协同

中国工程热物理学会传热传质学学术会议论文编号：153001换热系统优化中换热与阻力的协同王怡飞陈群1清华大学工程力学系,热科学与动力工程教育部重点实验室,北京，100084（Tel:[phone]E-mail:[e-mail]）摘要：本文在对换热系统流动阻力和传热过程分析的基础上，建立了系统的阻力约束和传热约束，提出了一种考虑流动阻力的换热系统优化方法。通过实验测量换热系统实验台的结构参数，以变频泵的总能耗最小为目标，对换热系统运行参数进行优化，得到了不同管网阻力特性条件下各变频泵的最优运行频率。优化结果表明，在给定条件下，某一支路的阻力系数增大，其最优分配的流量则会相应减小，各变频泵转速不同程度地增大。最后，通过比较最优工况和可行工况下各换热器和整个换热系统的熵产及变频泵的总功耗，明确了将流动阻力通过熵产引入优化准则的方法并不总能够获得不同应用场合下的最优结果，而建立系统的换热和阻力的广义约束能够很好地实现换热系统优化中换热与阻力的协同。关键词：换热系统，节能，协同优化，热阻，流动阻力0.前言换热器结构或布局的优化可以增强换热效果，但因为强化换热的同时增加了流动阻力[1,2]，导致强化换热未必能达到节能的效果。因此，减小流动阻力在换热系统节能优化中具有重要意义[3]。然而，很多关于换热系统的研究中并没有考虑流动阻力的影响[4,5]，使其优化结果并非最优[6]。也有很多学者尝试从换热过程、换热部件和换热系统三个层次平衡减小流动阻力和强化换热这两个相互矛盾的优化目标。在换热过程层次，研究人员大多致力于实现提高传热系数和减小流动阻力的平衡[7]。例如，Paoletti等人采用Bejan数来衡量传热过程中流动和传热所引起的不可逆性[8]，Mahmud则推导得到了不同对流换热过程中的Bejan数[9]。Cheng等人则采用与流动阻力和换热量相关的目标函数来优化毫米量级的空气槽道的尺寸参数[10]。这些方法本质上是将传热过程中与流动和传热的相关参数作为评价准则进行优化。然而，对于不同应用场合的传热过程，某个固定的评价参数未必能总对应最优解。因此，除了将流动阻力引入优化准则外，研究人员也尝试在给定流动阻力的前提下提高换热系数。例如，Polley等人在给定的流动压差的条件下优化了换热过程中的传热系数[11]。在换热部件层次中，单个换热器的优化注重于平衡强化换热以减少初投资和减小流动阻力以减少运行成本。同样，大量评价参数被作为准则函数引入到换热器的性能优化中。例如，由压差引起的熵产和由传热引起的熵产之和被用来优化换热器的结构及翅片的尺寸[12,13]，以达到运行成本和初始投资均为较低的目的[14]。除此之外，很多方法将流基金项目：国家自然科学基金(资助号：51422603)和全国优秀博士学位论文作者专项基金(资助号：201150)动阻力作为前提，优化换热器的热导以求得最大换热量[15]。然而，绝大多数的换热器都是热力系统的部件[16]，其流动阻力不仅仅只来自于换热器内部，而且来自于换热系统的输配管网[17]，因此在不完整的阻力分析的基础上得到的结果并非最优。由于过程和部件层次上的优化存在流动阻力考虑不完全和局部优化结果并非系统最优解的局限，因而很多学者提出了换热系统的流动阻力模型，并将流动阻力作为准则函数用于换热系统的优化。例如，Zhu等人将管壳式换热器网络的流动阻力作为准则函数优化其运行成本[18]，Kim等人建立了换热器网络压差与体积流量的关系并将其作为准则函数用于优化冷冻水网络的优化[19]。这些研究为分析换热系统阻力对系统性能优化的影响提供了参考，但是这些方法实际上仍然是将流动阻力引入准则函数，而系统的约束仅仅来自于传热过程分析。另外，也有学者将流动阻力作为系统优化的约束。例如，Frausto-Hernandez等人在分别给定了冷、热流体的最大可行压差的前提下对换热系统进行了优化[20]。然而，事先给定最大可行压差并没有从物理上提供有效的约束，并且其数值的给定往往依赖经验。因此，换热系统优化中仍然缺乏有效的流动阻力约束。除此之外，人们通常采用能量守恒和传热方程构建换热系统中的传热约束。该方法不可避免需要引入中间温度[21]，增加了系统优化的复杂度，甚至影响了最终的优化结果。近年来，Chen等人提出了基于火积耗散热阻的换热系统优化方法，消除了中间温度，建立了换热系统中结构参数、需求参数和运行参数间的直接联系[22-25]。本文在对换热系统流动阻力和传热过程分析的基础上，建立了系统的流动阻力约束和传热约束，提出了一种考虑流动阻力和变频泵运行特性的换热系统优化方法。通过实验测量换热系统实验台的结构参数的基础上，以变频泵总能耗最小为目标，得到不同管网阻力系数的条件下各变频泵的最优运行频率。优化结果表明，增大某一支路的阻力系数，其最优分配的流量则会相应地减小，各变频泵转速则会不同程度地增大。最后，本文通过分析最优工况和可行工况中各换热器和系统的熵产及变频泵总功耗的大小，明确了将流动阻力引入优化准则的方法未必能够得到不同应用场合下的最优结果，而将流动阻力引入系统优化的约束条件能够很好地实现换热系统中换热与阻力的协同。1.换热器网络实验台及测量仪器图1是换热器网络实验台的结构示意图[25]，实验台的组成包括两个逆流换热器1和2，3个变频泵及对应管网上3个开度可调节的调节阀，1个恒温热水水浴箱和1台制冷机，换热器及管道外表面均用保温材料包裹。在变频泵的驱动下，热量从恒温水浴箱经过换热器1和2，最终传递到制冷机的蒸发器。图2是换热器网络实验台的照片[25]。每个支路上均安装量程为20L/min精度为±0.5%的涡轮流量计以测量各支路的体积流量，每个变频泵进出口均接有量程为350kPa精度为±0.2%的压差传感器以测量其压差。在图1所示的每个温度节点均安装了2个精度为±0.2oC的T型热电偶（Omega生产）以测量各节点的温度。支路1和3中的工质为水，支路2中的工质为R142b。本文取水得密度为1000kg/m3，定压比热容为4196J/(kgK)。支路2处装有量程为1.8Mpa的绝对压力表，测量其支路的绝对压力为0.47MPa，在该压力下，R142b为饱和液态，查得在300K±10K的温度范围内，其密度为1100kg/m3，定压比热容为1304J/(kgK)（NIST-Refprop）。HeatExchanger1HeatExchanger2EvaporatorTankTw1,iTw1,oTr1,oTr1,iTr2,iTr2,oTw2,oTw2,iTwe,iTwe,oTeTeVSP1VSP2VSP3DPT1DPT2DPT3VFM2VFM1VFM3PGMV1MV2MV3DACSChiller图1换热器网络实验台示意图图2换热器网络实验台照片2.系统分析及其中特征参数的实验测量2.1换热器网络传热过程分析基于火积耗散热阻，图1所示的换热器网络的传热过程可分为4个部分：（1）换热器1中热水被冷却；（2）换热器2中冷冻水被加热；（3）被冷却的热水与热水箱中热水混合；（4）被加热的冷冻水在蒸发器中被冷却。这四个部分中的火积耗散热阻为[22]:11,1,1,1,11111,12,211exp111,22exp1wiworirohxppkATTTTRQmcmckA，(1)22,2,2,2,22222,23,322exp111,22exp1rirowiwohxppkATTTTRQmcmckA，(2)1,1,1,1122wiwompTTRQmc，(3),,3,3exp121,22exp1eweiweoeeeeepeekATTTRQmckA，(4)其中，R火积耗散热阻，T是工质温度，ξ是换热器的流型系数，k是换热器的传热系数并且假定为常数，A是换热器的传热面积，m是质量流量，cp是水定压比热容。下标hx1/hx2表示换热器1/2，数字1/2/3表示1/2/3循环支路，m表示混合过程（3），e代表蒸发器，w/r代表工质水/R142b，i/o代表换热器的进口/出口。假设不考虑漏热，则Tr1,o等于Tr2,i，Tr1,i等于Tr2,o,Tw2,o等于Twe,i，且Tw2,i等于Twe,o，基于公式(1)~(4)，可以得到换热器网络中关于热水箱水温Th、蒸发温度Te、传热量Q和不同过程的火积耗散热阻的直接关系如下：12hemhxhxeTTQRRRR。(5)2.2换热器网络的流动阻力分析当变频泵在转速ω下运行时，变频泵的压头H和质量流量m有以下关系[26,27]220122mmHaaa，(6)其中，a0/a1/a2是泵的特性参数，可以通过实验测得。不考虑机械损失，变频泵消耗的电功率可表示为PmgH。(7)换热器网络中，管网的压头可以表示为[28]2222iisdsiiiifLmHHHHKdgS。(8)其中，Hs是静压头，Hd是动压头，g是重力加速度，ρ是工质的密度，S是管道截面积，d是管道直径，f是沿程损失系数，L是管道长度，K是具体结构（比如90°弯头、突扩、突缩和阀门等的局部损失系数，下标i代表第i段管道。对于一个给定的管网，管道结构、管道长度、管道截面积和流体的密度均为已知，假设沿程损失系数为常数，那么式(8)可以变为：2sHHbm，(9)2212ifLbKgSd。(10)联立公式(6)和(9)，则可以得到通过管网流动阻力分析而建立的约束2220122smmaaaHHbm。(11)2.3部件特征参数的实验测量参考文献[25]中对换热器网络实验台的测量方法，分别测得了给定转速下换热器1、2和蒸发器的热导值为182.6、143.6和56.0W/K，表1给出了各变频泵的特征参数，表2给出了各管网的特征参数。表1各变频泵的特征参数VSPa0a1a2δmaxVSP10.01202.63-1304.740.6%VSP20.00811.50-1266.213.7%VSP30.00932.07-717.081.3%表2各管网的特征参数LoopHs(m)bδmax100%80%60%40%20%Loop10.572623.63213.04050.45135.76468.96.9%Loop20.101205.21349.01553.51818.52144.23.7%Loop30.36413.0548.7741.3990.71297.16.8%3.换热器网络优化方法对于上述换热器网络实验台，在给定换热量的情况下，如何确定每个变频泵的运行参数以使泵的总能耗最低是一个典型的优化问题。基于对换热网络中的传热过程分析和管网阻力分析，约束式(5)(11)可以通过拉格朗日乘子（Lagrangemultipliers）法构建方程为：212,220,1,2,21,2,3ihemhxhxeiisiiiiiiiiiiiiiiiPTTQRRRRHHbmmmHaaai，(12)其中，α和β都是拉格朗日乘子。式(12)分别对mi和ωi求变分可依次得到1121,1111,112,111111121111111221,111exp11exp122=0expexp1ppkAmckAaamgHdmQkAmkAmckA，(13)1122,2111111222,2111,222,222222222222222,222222222,222exp112exp1expexp122exp112exp1expexpppppkAmckAkAkAmckAaamgH