摘要:为深入研究金属管浮子流量计的工作机理,采用计算流体力学和流体力学中湍流模式理论,对金属管流量计流场进行了数值仿真研究,获得了低黏度流体介质中浮子受力及浮子在受力平衡下的流量.数值模拟与物理.实验标定数据的对比表明,模拟计算的流量最大满度误差为5.4695%,平均满度误差为2.4731%,说明数值仿真模型能满足金属管浮子流量计设计的需要.
　　金属管浮子流量计是一种传统的变截面流量计，具有结构简单、工作可靠、压力损失小且稳定、可测低速流体介质等诸多优点,广泛应用于测量高温、高压及腐蚀性流体介质。目前，金属管浮子流量计设计采用经典浮子流量公式计算,其中流量系数a受到浮子形状和来流雷诺数等多种因素的影响,需要实验标定来确定.该设计方法需进行大量的实验,费用昂贵,周期很长，且获取的数据有限.为深人研究金属管浮子流量计的工作机理,笔者利用计算流体力学方法对金属管浮子流量计进行数值仿真,通过构造仿真模型为研究金属管浮子流量计的机理建立数值实验平台,优化传感器的结构;同时为金属管浮子流量计提供了低成本、短周期的设计方法.
1基本原理
1.1 金属管浮子流量计的工作原理

　　图1为金属管浮子流量计的工作原理示意.在垂直的锥形管中，放置-.阻力件,即浮子.当流体自下而上穿过锥管时,受到浮子迎流体积阻挡而产生一个压差,从而对浮子形成向上的作用力,同时由于流体本身的黏性,对浮子产生黏性力,当这两个力的合力大于浮子本身的重力时,浮子就会向上升,同时浮子与锥形管间的环形流通面积增大,流速减低,此时浮子对流体阻力作用减小.当浮子受到的力达到平衡时,浮子就会停留在某一高度.传统的金属管浮子流量计设计采用经典浮子流量公式2计算,即.

　　式中:q1为浮子流量计的体积流量;α为流量系数;D0为浮子最大迎流面的直径;h为浮子在锥管中的垂直位置;?为锥形管夹角;Vf为浮子体积;ρf为浮子材料密度;ρ为流体密度;A1为浮子垂直于流向的最大截面面积.
1.2计算流体力学及其控制方程
　　计算流体力学(computationalfluiddynamics,CFD)是门利用计算机求解描述流体流动规律的控制方程组技术，涉及到流体力学.计算方法和计算机图形处理等技术.为了简便,以不可压缩湍流流动为例写出笔者使用的k-?模型通用形式的流体控制方程.在直角坐标系中.流动可由连续性方程和雷诺时均N-S方程描述,即
连续方程

k和ε由各自的输运方程得到，对高Re数问题有

　　笔者的仿真模型建立在已有仿真实验的基础上”,并针对以往模型中存在的人口速度剖面为等值面这一不合理设置,给出合理的圆管人口速度剖面;同时按照实际流量计的构造,在流量计的人口与出口处加入导流架，从而提高了仿真结果的准确性.
2模型建立和数值计算
　　针对25mm1径(DN25)浮子流量计内部机械结构进行流场计算和分析.金属管浮子流量计结构剖面如图2(a)所示.金属管浮子流量计采用1Cr18Ni9Ti材质.流量计设计流量范围为0.4~4.0m3,量程比为10:1,为了检测锥管中浮子的位置.在浮子中安装有水久磁铁,同时利用网格生成器建立了浮子流量计内部结构仿真模型,如图2(b)所示实验所用的低黏度流体介质是在常温下密度为998.0kg/m3、黏度为0.001002kg/(m.s)的水.金属管浮子流量计内部流场是高雷诺数完全发展湍流流动,故采用湍流模式理论中标准k-?模型来计算.

2.1入口速度剖面设定
人口速度由经验性光滑圆管湍流速度分布指数公式表示,即

为使构建模型能满足不同口径金属管浮子流量计
的仿真结果,引人流量修正因子Cq,通过修正人口流
量来适应多种口径金属管浮子流量的仿真模型.
流量的修正公式如下:
qc=Cqq(9)
式中:q为实际流量;qc为修正流量;Cq为流量修正因子,对于25mm口径的金属管浮子流量计,Cq=1.13.
　　根据式(8)~式(11)编写速度分布函数.图3是q=4.8m3/h时金属管浮子流量计人口速度剖面的仿真.图中色标颜色由冷色调到暖色调表示速度由小到大,从图中可以清楚地看到从边壁到中心的速度是由小到大的非线性分布.


3仿真结果及分析
　　图4给出了不同人口流量下,仿真计算得到的金属管浮子流量计内部截面上的速度伪色图.可以清楚地看到金属管浮子流量计中流体在浮子周围以及出入口的速度分布.随着流量的增加,浮子在管中的位置上升,浮子与管道之间环隙变大，流体在管中的速度分布也随之发生明显变化.通过观察伪色图中的速度分布,并根据流体力学基本原理,可以初步判断出计算所得结果是合理的.

　　根据物理实验结果,仿真计算了浮子在锥管中某一高度处流体对它表面的压力、黏性力和浮子重力Gf3项的合力Ff仿真结果中浮子所受合力结果见表2,合力向上为正方向.
　　由表2可知计算所得浮子在各位置的合力接近于零,趋于平衡,为验证仿真模型所能达到的仿真精度，根据浮子只有在所受合力F,为零时才能平衡,又计算了浮子在各个高度保持平衡,即合力为零时,所对应的人口流量.表2还列出了数值计算与物理实验所得流量的比较.并计算出仿真实验的满度误差δs,即

式中;qs为数值模拟流量,m3/h;qp为物理实验流量.

　　中表2可知，数值计算与物理实验所得流量的最大满度误差为大满度误差为5.4695%,平均满度误差为2.4731%.利用本模型对50mm、80mm口径金属管浮子流量计仿真数据与物理实验数据的比较结果也相当接近.表明笔者构建的模型得到令人满意的结果.
4结语
　　通过计算流体力学方法对金属管浮子流量计进行仿真,并对流场进行了定量分析,计算出浮子在金属管中不同垂直位置的受力大小及受力平衡时所对应流量.仿真实验结果与物理实验数据相比,最大满度误差为5.4695%,平均误差为2.4731%,表明笔者构建的金属管浮子流量计数值模型能够满足低黏度介质金属管浮子流量计设计要求,为金属管浮子流量计传感器结构进一步优化提供了可靠的数值仿真平台。

以上内容源于网络，如有侵权联系即删除!