仪器仪表直供热线 18952302362
          三畅创新研发仪表、助力工业自动化升级

          差压式流量计测量脉动流量方法研究

          来源: 北京长城计量测试技术研究所 作者:张永胜, 张毅治, 刘彦军 发布日期:2020-06-09 10:40:48

          摘要:利用脉动流动工况下差压式流量计理论模型,给出了通过双路差压计算瞬时脉动流量的算法,并利用FLUENT软件对文丘里管内流场进行仿真计算,对脉动流量计算方法进行验证。仿真结果表明,单路压差计算脉动流量,由于遗漏流量导数项可产生很大测量误差,非常大误差甚至可达100%;而通过双路差压计算脉动流量,可有效消除流量导数项影响,算例中非常大误差降低至5%左右。

           
          1引言
                  差压式流量计[1]一般以孔板、文丘里管、喷嘴或其它形式的节流元件形成差压,进而计算流体流量。差压式流量计由于其结构简单、无可动部件等特点而广泛应用于测量封闭管道中单相稳定流体的流量,常用于测量空气、天然气、水等流体的流量。差压式流量计是非常为传统且应用广泛一种流量测量装置,该类流量计的准确测量是以平稳流动为基础,当测量流体压力、速度随时间呈周期变化情形时,以稳态流的计算方法直接求取脉动流量将因遗漏流量导数项而引起较大测量误差[2~4]。而脉动流广泛存在于工业管路中,由于旋转式、往复式或可运动传送设备的使用,管道中流量控制设备的周期振荡,以及流程中工艺设备的特定要求,众多因素都会导致流体产生脉动流[5~8]。ISO5167中明确规定差压式流量计不适用于脉动流的测量。
           
                  本文以经典文丘里管[9]为研究对象,通过FLUENT软件[10]对文丘里管在正弦脉动流下的流场进行仿真计算,分析了单路差压计算脉动流量时由导数项引入的测量误差,提出了基于双路差压的脉动流量测量方法,并利用仿真结果证明该方法的有效性。
           
          2脉动流量计算方法
          2.1脉动流理论模型
                  在脉动流状态下,假设流体流过节流件为一维流动[11,12],流量计内流动运动方程为:
          流量计内流动运动方程
                  由于质量流量qm可由式(2)表示,则式(1)可进一步表示为式(3)。
          质量流量 qm公式
                  式中Ax为x处管道横截面积。
           
                  假设节流件上游取压孔处于截面1位置,坐标为x1;下游取压孔处于截面2位置,坐标为x2;在截面1和截面2之间沿流线对式(3)进行积分,可得:
          公式4
           
                  式中:A1和A2分别是截面1和截面2的截面积;Δp为截面1和截面2之间压差;C为节流件流出系数。
          公式5
          2.2双路差压流量计算方法
                  在实际情况中流量导数项无法确定,传统的通过两截面之间单路差压计算流量的方法无法计算瞬时脉动流量,可通过测量两对截面之间双路差压,消除流量导数项的影响,进而计算瞬时脉动流量。在收缩段前后各取2个截面,分别定义为截面Ⅰ、截面Ⅱ、截面Ⅲ和截面Ⅳ,如图1所示。4截面沿轴向的坐标分别是xⅠ、xⅡ、xⅢ和xⅣ。截面Ⅰ与截面Ⅲ内径相同,截面Ⅱ和截面Ⅳ内径相同,截面Ⅰ和截面Ⅱ之间流出系数为C1,截面Ⅲ和截面Ⅳ之间流出系数为C2。
          文丘里管几何模型
                  在截面Ⅰ和截面Ⅱ之间,截面Ⅲ和截面Ⅳ之间具有如式(6)和式(7)所示关系。
          截面Ⅲ和截面Ⅳ之间公式20200609111509.jpg
                  联立式(6)和式(7),消除流量导数项,可得质量流量与双路差压的函数:
           

          质 量流量与双路差压的函数:

           
          3、流体仿真验证
          3.1几何模型
                  本文以文丘里管为对象,利用FLUENT软件进行流体仿真,获得两对截面间双路压差,对脉动流量计算方法进行验证。图1所示的文丘里管进出口内径D=30mm,喉径d=19.5mm,直径比β=0.65,以进口截面中心为坐标原点,截面Ⅰ、截面Ⅱ、截面Ⅲ和截面Ⅳ这4截面坐标分别是xⅠ=25mm、xⅡ=90mm、xⅢ=46mm和xⅣ=78mm。
           
          3.2流体仿真模型
                  利用ICEM软件进行几何建模,忽略重力因素影响,基于文丘里管内流场轴对称特性,建立二维轴对称几何模型并进行网格划分。文丘里管上游设置10倍长直管段,下游设置5倍直管段。仿真模型壁面设置5层边界层网格,中心区域为三角非结构化网格。
           
                  仿真中选用标准k-ε湍流模型,选用液体水作为流体介质。设置下游直管段出口为自由出流条件,直管段及文丘里管壁为静止壁面;上游直管段入口为速度入口,通过UDF函数定义,可计算得到入口流速。
           入口流速公式
                  式中:qm为瞬时质量流量,kg/s;q0=2kg/s为质量流量平均值;A为脉动幅度;f为脉动频率;t为运行时间,s;Vt为瞬时入口速度平均值,m/s;ρ为流体密度,kg/m3;Ain为上游直管段入口截面积,m2。
           
          3.3单路差压流量计算误差
                  由稳态流动下仿真结果,可获得截面Ⅰ和截面Ⅱ之间流出系数C1=0.9815,截面Ⅲ和截面Ⅳ之间流出系数C2=1.007。进一步可获得k1,Ⅰ-Ⅱ=141.75,k2,Ⅰ-Ⅱ=4788.9,k1,Ⅱ-Ⅲ=71.856,k2,Ⅲ-Ⅳ=4549.4。
           
                  稳态流动下,质量流量qmc可由式(11)计算:
          公式11
                  正弦脉动流情况下,流量导数可由式(12)表示,若直接通过式(11)计算瞬时质量流量引入的导数项误差通过式(13)计算:
          算瞬时质量流量引入的导数项误差
                  表1所示为7个工况单路差压计算流量误差数据。其中Ep-max表示误差上限,En-max表示误差下限。
          单路差压计算流量误差表
                  图2是工况1和工况7单路差压流量误差曲线。表1中“×”位置和图2(b)中曲线中断处是由于ΔpⅠ-Ⅱ为负值造成。通过ΔpⅠ-Ⅱ和ΔpⅢ-Ⅳ计算瞬时流量,对应时刻截面Ⅰ-Ⅱ和Ⅲ-Ⅳ误差绝对值只有2个时刻相等,其余前者均大于后者。通过ΔpⅠ-Ⅱ计算脉动流量误差上限可达60%,误差下限非常大可达-100%,甚至差压为负值导致无法计算;通过ΔpⅢ-Ⅳ计算脉动流量误差上限可达34%,误差下限非常大可达-65%。因此减小k1,可有效降低由于遗漏导数项引入的脉动流量计算误差。
          单路差压流量计算误差曲线
          3.4双路差压流量计算
                  基于式(8),利用ΔpⅠ-Ⅱ和ΔpⅢ-Ⅳ计算瞬时脉动流量,表2是双路差压脉动流量误差数据。图3是工况1和工况7双路差压流量计算误差曲线。图4是脉动流量计算误差带分布。
          双路差压计量流量误差数据表
                  通过双路差压计算脉动流量,可明显降低误差上限,提升误差下限,使得误差带变窄,在高频率大幅值下尤其显著。相对于双路差压计算流量,通过ΔpⅠ-Ⅱ计算流量的误差带宽度非常小,为单路的1/24;通过ΔpⅢ-Ⅳ计算流量的误差带宽度非常小,为单路的1/8。
          脉动流量计算误差分布
          4结论
                  研究中针对差压式流量计测量脉动流量时由导数项引入的计算误差,提出了采用双路差压消除导数项的计算方法,并以文丘里管为研究对象利用流体仿真进行验证分析,得到以下几方面结论:
                  1)通过单路差压采用稳态流量公式计算脉动流量,由于遗漏流量导数项分量,造成很大影响,误差极值可达100%,测量误差带覆盖-100%~+60%,甚至由于差压为负值而无法进行计算。
                  2)通过双路差压建立的微分方程组可有效消除导数项影响,脉动流量计算误差极值为4.7%,误差带分布覆盖-4.7%~0.6%,误差带明显收拢。
                  3)建议对双路差压计算脉动流量方法在理论和试验方面进行进一步深入研究,验证该方法的有效性和可行性。
           
          友情提示:感谢本文的作者分享,本文一切观点都来自于本文作者。

          咨询热线

          0517-86998326
          彩票平台排名