资源描述：
《Z型并联管道流量分配的数值研究.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、中国工程热物理学会传热传质学学术会议论文编号：133427Z型并联管道流量分配的数值研究何嘉1，王良璧2(兰州交通大学机电工程学院铁道车辆热工教育部重点实验室甘肃兰州730070)（Tel：0931-4956556，Email：lj811719@126.com）摘要：通过建立简化的流动模型和采用数值分析方法对Z型并联管道层流条件下的流动特性进行了研究。获得了流体在并联管道内的流动及其流量分布特性。研究结果表明：靠近主管道来流入口附近的支管流量较小，离主管道出口近的支管流量较大；所研究的Re范围内，流量在各支管的不均匀性为0.88~1.20，且随着Re的增大，流动
2、不均匀性增加。该工作为准确确定Z型并联管道系统流量分配偏差提供了一种有效的方法。关键词：并联管道；层流；流量分配；数值计算0前言并联管道在工程中，尤其在锅炉、热交换器、太阳能、机车冷却系统等设备中有着广泛的应用。流体在并联管道内的流量分配关系到设备运行的经济性和工作可靠性，研究此类模型的流动分布特性具有相当高的工程实际价值。文献[1~3]采用二维数值方法对层流条件下影响电子元件中冷却剂流量分配的因素进行了分析，重点论述了Re数、模型几何参数的作用，分别提出了获得较好流量分配均匀性的方法。文献[4~6]用不同的方法对流体在太阳能系统中的流动特性进行了研究，文献[4
3、]通过求解代数方程组对充分发展层流条件下的流量分配和压降进行了研究，并考虑了流体等温与不等温的情况。文献[5]采用离散的流体动力学计算模型对影响流量分配的多种因素进行了对比。文献[6]通过实验的方法研究了层流条件下太阳能系统中的流量分配问题。文献[7]通过采用离散模型求解集管内的压力分布对平行流换热器系统中的流量分配进行预测，用实验结果进行验证。文献[8]采用实验和数值模拟方法考察了影响矩形集管平行流换热器流量分配的因素。文献[9]在质量和动量守恒的理论基础上了发展了解决Z型结构的燃料电池流动分配和压降的方法，对描述流动性能，系统几何结构以及流动条件的5个特征参
4、数对流量分配和压降产生的影响进行了研究。文献[10]采用非线性最小二乘法(LMM)与CFD-ACE+软件程序对三维Z型并联管道如何获得流量分配均匀性的问题进行了研究。并联管道的流量分配问题在数学描述和求解方面均比较繁琐复杂，其微分方程是二阶非线性的，一般很难求得其精确解，实验研究通常需要耗费大量的人力和财力。因此为正确认识Z型并联管道中流体的流动特性及分配机理，找到准确性高，省时省力的研究方法对工程实践具有一定的指导意义。文采用三维流体通道数值模型对层流条件下Z型并联管道中的流量分配进行数值计算，得到流体的速度和压力分布，阐明流量分配特性。数值方法可更加全面的了
5、解管道内每个位置的速度和压力分布情况。该数值计算方法可直接用于求解三维资助项目：国家自然科学基金（No.51206073）Z型并联管道内的流动和压力分布。主要讨论了不同Re数（指来流集管流体Re数，以下不再单独说明）对Z型并联管道流动特性的影响。1研究模型及数学公式1.1问题几何结构本文的研究物理模型源于机车冷却系统，通过对实际系统的简化建立如图1(a)所示的物理模型。该模型中有10个支管，支管均为等长、等径的直圆管。其几何尺寸如图1(a)所示。流体从来流主管进入，并行流经各支管后进入出流主管。在计算中，把整个图1(a)所示的区域分成12个计算区域，所生成的网格
6、如图1(b)所示。图1物理模型与网格划分：(a)物理模型几何尺寸；(b)网格划分1.2控制方程及边界条件为便于数值分析，作以下假设：(1)流体物性为常数；(2)流动是定常的，处于层流状态；(3)流体为黏性不可压缩流体；(4)不考虑作用在流体上的体积力。根据以上假设，计算区域的流体流动控制方程为：连续性方程(1)运动方程,(j=1,2,3)(2)式中，ρ为流体的密度，μ为动力黏度。通道壁面上的边界条件,,(3)来流集管进口边界条件,,(4)出流集管出口边界条件为局部单向化条件,,(5)各支管与集管连接处流动边界条件1）通过插值法计算来流集管每个出口的边界速度，将此
7、速度作为对应每个支管的进口边界条件；，，，(i=1)，(j=1,2,…,10)(6)式中，i表示来流集管；j表示各并联支管。2）通过插值法计算出流集管每个进口的边界速度，将此速度作为对应每个支管的出口边界条件；，，，(i=2)，(j=1,2,…,10)(7)式中，i表示出流集管；j表示各并联支管。3）各支管进出口流量修正；因适体坐标下的网格单元各边长均为单元长度，其网格单元在各速度方向上的面积为1，即q=vA=v。所以在适体坐标下，对流量的修正可认为是对速度的修正。，，，(i=1)，(j=1,2,…,10)(8)，，，(i=2)，(j=1,2,…,10)(9)式
8、中，α=uin/uout
