绝热毛细管分布参数模型的仿真.doc

绝热毛细管分布参数模型的仿真 摘? 要? 建立了绝热毛细管长度仿真模型，对比分析目前流行的绝热毛细管内摩阻系数和动力粘度模型，得到了计算准确、稳定的分步参数模型，为整个制冷装置的仿真提供了可靠保证，并在制冷装置设计和实际应用中，对毛细管的选择提供重要依据。关键词? 绝热毛细管 仿真 分布参数模型 1? 引言毛细管作为一种简单、工作稳定的节流元件得到了广泛的使用。对毛细管长度和内径的选择计算前人作了大量的工作，随着计算机技术的发展，模型的建立成为精确计算的重要环节。本文建立了绝热毛细管长度仿真模型，进行计算机模拟，对绝热毛细管内摩阻系数分别以Stoecker方程和Churchill方程进行计算，对制冷剂动力粘度分别以McAdams模型、Cicchitti模型、Dukler模型进行计算，对比分析以上模型，得到准确实用的绝热毛细管分布参数模型。2?毛细管长度模型2.1? 毛细管长度物理模型为了简化毛细管的数学模型，需要在正确了解毛细管内所发生的物理现象的基础上，作以下的物理假设：1)?假设工质流经毛细管是一维绝热均相流动，流动处于紊流区，且制冷剂在流动过程中状态变化是连续的。2)?不考虑毛细管的内径不一致与表面粗糙度的影响。3)?忽略亚稳态流动。建立用于分析的物理模型如图1所示。2.2? 毛细管长度数学模型? 对于稳态流动，连续性方程：??????????????????????????????????? （1）能量方程： ????????????????????????????? （2）动量方程：????????????????????? （3）制冷剂流体进入毛细管，开始处于过冷区，比容不变，焓值也不变，由式（3）积分得：??????????????????????????? （4）随着制冷剂在毛细管中继续流动时，压力逐渐下降，焓保持不变，流动变为两相流，根据假设对上述方程离散化，则：连续性方程：????????????????????????????? （5a）或?????????????????????????????????????????? （5b）能量方程：?????????????????? （6）动量方程：??????????????? （7）?为常数，则：? 由（7）式得：????????????? (8)其中焓、比容由下式确定：? ? ?速度取1、2点的平均值：? 内摩阻系数 采用Stoecker[1]模型：??????????????????????????????????? （9）?取1、2点的平均值： 动力粘度采用Cicchitti[2]模型????????????????????????? （10）若1点参数已知，由以上方程可计算出2点参数,由式（8）计算绝热毛细管长度 ，进而计算出毛细管管长。以R22为例，表1给出当冷凝温度40℃，过冷度0℃，毛细管内径1.63mm，制冷剂流量10g/s时的计算结果。表1：Stoecker模型的毛细管计算值序号?温度(℃)?压力（bar）?干度?焓值(kJ/kg)?速度(m/s)?⊿L(m)?累计长度(m)1?40?15.27?0.0000?249.01?4.23?0.00000?0.000002?39?14.90?0.0075?249.01?4.74?0.22825?0.228253?38?14.54?0.0150?249.01?5.28?0.19769?0.425944?37?14.18?0.0223?249.01?5.85?0.17750?0.603445-34?—?—?—?—?—?—?—35?6?6.02?0.2463?256.40?49.16?0.00262?2.0522136?5?5.84?0.2563?257.43?52.49?0.00119?2.0534037?4?5.66?0.2671?258.60?56.15?0.00025?2.0536538?3?5.48?0.2784?259.94?60.10?-0.00050?2.05315由表1看出，当计算温度到３℃时，计算值 为-0.00050m，与实际不符，这说明制冷剂在毛细管中的流动发生了壅塞现象，则可取上一个计算点的值2.054m作为该条件下的毛细管管长。当冷凝温度40℃，蒸发温度5℃，绝热毛细管管长计算值2.053m,与ASHRAE[3]在相同条件下的推荐值仅相差1.6%，因此该模型精度是较高的，而且计算时间非常短，瞬间完成。2.3? 内摩阻系数和动力粘度模型比较目前常用的绝热毛细管内内摩阻系数的计算方程Chruchill[4]方程：????????????????????????? （1