管壳式换热器的分段设计方法.docxVIP

管壳式换热器的分段设计方法 1 管壳式换热器的设计 设计方法粗糙，设计精度低是传统管壳式换热器的制造成本和运营成本之一。国外对管壳式换热器设计的研究已经取得了显著成果,其中影响最大的是HTRI和HTFS两家公司,但出于商业机密原因,设计方法一直没有公开。在管壳式换热器发展的几十年中,设计的重点和难点在于壳程流体流动与传热具有复杂性和不确定性,结果难以准确预测,所以准确计算壳程传热膜系数和压降成为管壳式换热器设计的一个重要研究对象。目前公开的设计方法结果都不很理想,其中Bell法是应用最为广泛,也是计算相对最准确的设计方法,这种方法把换热器作为一个整体来考虑,但未能涉及到换热器内部的温度场和随温度场变化的一系列参数,尤其是温度变化较大时,流体物性参数的影响会很大,这也是换热器设计不确定性的主要原因,基于上述问题,本文采用分段设计方法,将换热器按照管壳程结构划分成有限个单元,分别进行设计计算,有效减小换热器设计误差。 随着计算机的发展,工业设计大都采用计算机程序进行计算,VB6.0作为Windows系统下的软件开发工具,具有开发效率高、设计灵活、界面友好、交互性强等优点,本文运用VB6.0编程语言,开发了管壳式换热器设计软件,不仅节省计算时间,而且提高了设计精度,还能为后续软件升级提供保证。 2 设计原则 2.1 本问题的解决方案 在换热器设计中,各设计目标之间存在一定的矛盾性,这使得构造目标函数最优点集的隶属函数非常困难,流体物性参数随温度的变化是产生不确定性的主要因素,而传统的换热器设计方法未能从根本上解决这些问题。本文采用分段计算方法进行设计[8~10],将换热器按照管壳程的几何结构和流动形式分为有限个连续的换热单元,流体物性取每个单元平均温度对应的物性参数,有效降低了物性变化带来的影响。如果壳程数为Ns,管程数为n,折流板数为Nb,那么总换热单元数为(Nb+1)Nsn个;如图1所示,单壳程双管程换热器,折流板数为2,则共有6个换热单元。分别求出每个单元流体的进出口温度,得到各单元的传热系数和压降,进而求出整个换热器的总传热系数和压降,完成工艺计算,可以看出,分段设计计算最重要的步骤是要得到每个单元的温度分布。 2.2 冷凝器温度分布计算 采用有限差分法计算单元温度分布,不仅简单快捷,而且精度很高。雷俊杰等对换热器无相变传热温度分布计算做了详细研究,但对于冷凝温度分布,需要借助专业物性软件生成的冷凝曲线来计算,利用专业物性软件实际上只是得到热流体几个温度点和气相分率及对应的焓值,对于特定介质,其冷凝曲线是确定的,因此可以拟合出焓值-温度和气相分率-温度间的关系式,进而求解温度分布。这里以1-2型冷凝器设计为例(管程冷凝),介绍冷凝温度分布计算,如图2所示,换热器折流板数为Nb,管程数为n,模型做以下假设: (1)冷却介质质量流速及定压比热沿整个换热器是常量。 (2)壳程流体在横截面上均匀混合,温度相等。 对于无相变设计,当折流板的数目足够多时,子单元内的平均温度差近似等于两流体靠近前端管箱一侧的温度差,但对于冷凝计算,若直接用靠近前端管箱一侧的温差代替对数平均温差,误差较大,需引入温差修正系数Ft,壳程流体传热速率方程: 已知管程进出口温度和管箱侧的壳程温度,则换热器设计型计算和校核型计算时的所有工艺温度都可以得到。根据物性软件拟合焓值-温度和气相分率-温度间的关系式: 从左向右依次进行设计计算,边界条件:S=0时,tⅠ=to,tⅡ=ti,Tx=To 以壳程第x个单元为研究对象,管程上方子单元内,壳程流体速率方程为: 管程流体能量方程为: 联立式(1)和(5)可得Ⅰ管程第x单元温度t1;同理,可得下方Ⅱ管程第x单元温度。 开始计算时需初选换热系数U及温差修正系数Ft,如果最右端单元温度,或者,则重新选择换热系数,返回迭代计算;如果换热系数超出经验值范围,则调整Ft值,返回迭代计算,直至满足条件。计算得到单元温度分布后,代入式(3),可得到单元气相分率。 2.3 单元流型图判断 对于无相变传热,根据每个单元的雷诺数Re选择对应的传热与压降公式进行计算。对于有相变传热(冷凝传热或沸腾传热),需首先计算流型参数Cg和均相流体积分率Rlh,根据流型图判断单元流体流动状态,并选择公式计算;对于有相变传热压降计算,需根据气相分率判断流动模型,y≥0.7时为均相流,y0.7时为分离流,依次完成单元设计计算。 管程平均传热系数: 式中αti——管内单元传热膜系数,W/(m2·K) Ati——管内单元换热面积,m2 式中αsj——管外单元传热膜系数,W/(m2·K) Asj——管外单元换热面积,m2总传热系数: 管程压降由直管压降、管箱内回弯压降、流体进出管箱压降3部分组成: 壳程压降由折流板间压降与进出口局部压降两部分组成: 3 软件开发