用数值方法分析绝热毛细管的流量特性.docx

用数值方法分析绝热毛细管的流量特性高原田怀璋曾艳袁秀玲( 西安交通大学 , 陕西西安 710049)摘 要 : 根据均相流假设 , 用两相流基本方程建立了毛细管的数学模型 , 采用一种快速 、有效的计算方法对绝热毛细管的流量特性进行数值模拟 , 并把计算结果与实验数据进行了对比分析 。关键词 : 毛细管 ; 两相流 ; 均相流 ; 数值模拟中图分类号 : TB657文献标识码 : ANumerical Simulation of Two2pha se Flo w in Ca pillary2tubesGao Yuan Tian Huaizhang Zeng Yan Yuan XiulingAbstract : Experimental studies and numerical simulation have been performed in order to obtain optimum performance . A mathemat2i cal model based on fundamental governing equations of two2phase flow with homogeneous assumption is presented and is calculated with a more quickly and effective method. The results of numerical simulation to adiabatic capillary tube agree reasonably well with exper2i mental data .Key words : capillary2tube ;two2phase flow ; homogeneous flow ; numerial simulation符号下标A ———截面积cp ———定压比热d ———毛细管内径 e ———比能 λ———摩擦阻力系数 h ———比焓Δhf g ———汽化潜热L ———毛细管长度1m ———质量流量G ———质量流率n ———控制体数目p ———压力P ———周长1q ———热流密度T ———温度 v ———速度 xg ———干度z ———坐标轴ε———收敛率εg ———空泡系数 μ———粘度 Φ ———通用变量g ———气体l ———液体 i ———进口 o ———出口sat —饱和态w ———壁面1引言在家用冰箱 、房间空调器等小型制冷装置中 ,广泛使用毛细管作为节流元件 , 它具有结构简单 、 成本低的优点 , 但是不能象膨胀阀那样在较宽的 工况范围内工作 , 所以毛细管与装置的匹配是系 统优化的重要方面 。研究发现 , 以下因素影响毛 细管的制冷剂流量 :(1) 几何尺寸 :内径 , 长度 , 接管直径 ;(2) 进口压力 , 出口压力 , 进口液体过冷度或 两相混合物干度 ;(3) 内壁面相对粗糙度 ;(4) 壁面温度或热流密度 ; ( 5) 工质的热物理性质 。收稿日期 : 2001 —11 —14流体机械2002 年第 30 卷第 7 期4811m i = m o为确定它们之间的定量关系 , 从 20 世纪 50年代开始 , Bolstand 、Jordan 、Whitesel 等人在实验和 理论方面开展了广泛的研究[1 、2 ] , 提出了多个经 验关系式 。1988 年 , ASHRAE 根据这些研究成果 制定了一系列图表 , 成为工程上选用毛细管的依据 。1990 年 , Kuehl 和 Goldschmidt 在文献 [ 3 ] 中指 出 ,ASHRAE 发表的曲线与实验数据有很大的偏 差 。此后又出现了制冷剂替代问题 , 所以各国学 者对 毛 细 管 的 流 量 特 性 进 行 了 更 深 入 的 研 究[4 、5 ] 。随着计算机技术的发展 , 将流动方程做适当的简化和假设 , 建立数学模型 , 再进行数值模 拟 , 已成为一个重要研究手段 。本文把制冷剂在毛细管内的流动看作一维流 动 , 应用两相流基本方程 , 建立了毛细管的一维数 学模型 , 结合两相流摩擦压降和空泡系数的经验公式 , 对毛细管的流量特性进行了数值模拟 。在 数值模拟过程中 , 对每个控制体应用隐式格式迭 代求解 , 沿长度方向逐个推进 , 避免了全场求解时 差分格式稳定性的问题 , 保证了每一计算步长控 制体内流体均能满足能量方程 、动量方程和连续性方程 , 并将计算结果同实验数据进行比较 。( 1)d p11[ ( pi -po) -Δz ] A = mv- mvg , o g , o g , i g , id z11+ m l , ovl