弦月型通道热虹吸提升管流动特性数值计算.docVIP

弦月型通道热虹吸提升管流动特性数值计算 　　摘　要：应用两相流和溶液热力学的基本理论对热虹吸提升管的流动特性进行了分析研究，采用两相流分相模型对其进行了数值计算，得出了对应不同截面、不同浸没高度下加热负荷与农溶液的提升量、蒸汽发生量、浓溶液最大提升高度及稀溶液最小浸没高度的关系曲线，分析结果表明：在小加热负荷下(7.0kw左右)，直径比为19/32的组合套管的提升高度优于其他套管，最大提升高度高于直径比为16/32的组合套管。计算结果与实验结果符合良好，为该类型热虹吸泵的设计和运行提供了依据。 　　关键词：两相流；热虹吸；分相模型；数值计算 　　中图分类号：TK124　文献标志码：A 文章编号：0253-987X(2008)01-0017-06 　　 　　太阳能无泵LiBr吸收式制冷机最主要的特点之一是不使用机械泵，而是采用热虹吸泵将溶液从较低部位的发生器提升到较高部位的汽液分离器，从而实现制冷机中的溶液循环，系统装置如图l所示，随着发生器热源加热负荷的不同，热虹吸管内的两相流流型可能是泡状流、弹状流或环状流，有研究表明，若热虹吸管中可实现弹状流的话，将会获得较好的热虹吸泵工作特性，所以，应尽可能使热虹吸管内处于弹状流动，本文采用分相流动模型对热虹吸管内的流体动力特性进行计算。 　　 　　1　数学模型的建立及求解 　　 　　热虹吸管内的流动相当复杂，为方便分析，做如下假设： 　　 　　(1)垂直热虹吸管内的两相流是稳定的一维两相流，且流型为弹状流； 　　(2)采用分相流动模型进行热虹吸管内的流体动力特性计算，即把汽液两相看作是人为分离的两股流体的流动，它们的速度是常量但不一定相等，两相间处于热力平衡状态； 　　(3)热虹吸管微元控制体中蒸汽和溶液各自占有的容积总和等于整个微元的容积，各微元控制体内有一速度平均值和密度平均值； 　　(4)进入发生器的稀溶液浓度与发生器中的浓溶液浓度保持不变； (5)应用经验公式或简化的概念来表示两相摩擦因子及截面含汽率与独立变量之间的关系； 　　在以上假设的基础上，可求解下面的数学模型。 　　(1)由发生器质量和能量守恒方程，可得蒸汽质量流量Mg和浓溶液质量流量M1与发生器热负荷QG的关系 　　 　　式中：下标a表示与稀溶液特性有关的参数；下标1表示与浓溶液特性有关的参数；下标g表示与蒸汽特性有关的参数；ξ表示浓度；i表示焓。 　　(2)两相流截面含气率的计算采用米洛保尔斯基提出的利用滑动比S计算垂直上升管中截面含气率α的方法，即 　　 Pc为临界压力，Pa。 　　 　　式中：G为面积质量流速，kg/(m2?s)t0为剪应力，N/m2；g为重力加速度，从式(5)可看出，压降梯度由摩阻、重位和加速压降梯度3部分组成，即式中：xg为提升管出口处的干度，在本文中，可将热虹吸管由起沸点开始往上部分沿管长的热流密度视为常数，热虹吸管起沸点以下部分可按单相流体流动计算其压力降，汽液两相流摩阻压降采用按分相模型整理的两相流摩阻数据经验公式，这里采用罗卡特一马蒂内里(L-M)关系式来计算摩阻压降比较适合实验工况。 　　L－M关系式的基本假定是：气液两相流体在流动时气、液完全分开流动；单位管长上气、液两相的摩擦阻力压力降相等，且等于管子的摩擦阻力压力降；流动时不再考虑两相间的相互作用力，因此，两相流中各相的压降梯度dPa/dZ和dPfg/dZ应等于各相单独流过该相在两相流中所占流通截面时的压 　　 　　 　　2　计算结果及分析 　　 　　2．1　不同结构参数热虹吸泵的特性曲线 　　将计算结果整理成热虹吸特性曲线，如图2～图4所示，其中d1表示弦月型通道内管外径，d2表示弦月型通道外管内径，h表示稀溶液浸没高度，由图可知，在热虹吸提升管结构参数、运行条件一定的情况下，对于确定的热虹吸泵，有一个确保该热虹吸泵正常工作的发生器热源加热负荷范围，这从Sflow曲线与∑△P曲线的交点可以看出，超出这个范围，溶液的流动压头将不足以克服总的流动阻力来使浓溶液提升到热虹吸管的顶部出口，从而使制冷循环中断，而且发生器热源加热负荷过大还可能导致溶液浓度超出极限值，产生结晶现象。 　　可以通过增大Sflow或减小∑△P的各种措施来增大热虹吸泵的正常工作范围，在设计热虹吸泵时，选好弦月型通道的组合管径、严格控制热虹吸管的提升高度及浸没高度等具有重要意义，这些结构和尺寸设计得恰当，既可以保证热虹吸泵正常运行，又可使Sflow曲线与∑△P曲线的交点向两侧移动，扩大泵的正常工作范围。使其在所能提供的加热负荷范围内运行，提高运行的经济性与可靠性。 　　 　　由图2～图4还可看出，d1/d2=16/32的组合套管在浸没高度为57.5cm时可运行到更高的加热负荷，而d1/d2=19／