标准孔板流量计内部流场的CFD数值模拟

更新时间:2023-05-26 13:32:01 阅读量: 实用文档 文档下载

说明:文章内容仅供预览,部分内容可能不全。下载后的文档,内容与下面显示的完全一致。下载之前请确认下面内容是否您想要的,是否完整无缺。

 第22卷 第2期

 2008年06月

文章编号:1672-9897(2008)02-0051-05

实验流体力学

JournalofExperimentsinFluidMechanics

Vol.22,No.2 Jun.,2008 

标准孔板流量计内部流场的CFD数值模拟

陈家庆,王 波,吴 波,初庆东

(北京石油化工学院机械工程学院,北京 102617)

摘要:通过在标准孔板流量计中引入CFD数值模拟,为流出系数的获取提供了新途径。对不可压缩流体在不同流量、不同直径比、不同孔板轴向厚度和不同流动介质下的内部流场进行了数值模拟计算,并将计算出的流出系数与根据ISO公式计算出的流出系数进行了分析对比。结果表明,随着结构参数和工作条件的改变,流出系数都会随之发生变化,但ISO公式对孔板厚度的变化不太敏感。CFD数值模拟可以作为标准孔板流量计的辅助设计与标定手段,以进一步提高孔板流量计计量的准确性。  关键词:标准孔板流量计;流场计算;数值模拟;CFD  中图分类号:TE65;V211.3   文献标识码:A

CFDsimulationofflowfieldinstandardorificeplateflowmeter

CHENJia-qing,WANGBo,WUBo,CHUQing-dong

(SchoolofMechanicalEngineering,BeijingInstituteofPetrochemicalTechnology,Beijing 102617,China)  Abstract:TheutilizationofCFDnumericalsimulationinstandardorificeflowmeterprovidesanew

methodontheestimationandoptimizationofdischargecoefficient.Thispapercomplementednumericalsimula-tionontotheflowfieldofincompressibleflowunderdifferentflowrate,diameterratio,orificethicknessandflowmedia,bycomparingandanalyzingthedischargecoefficientachievedfromnumericalsimulationandtheISOempiricalformula.Itisindicatedthat,withthevariationofstructuralparametersandworkingcondition,thedischargecoefficientwillchangecorrespondently,butISOformulaisnotverysensitivetothethicknessoforificeplate.TheCFDnumericalsimulationcanbeadoptedasacomplementarydesignandcertificationtoolinstandardorificeflowmeter,toenhanceitsmeasurementaccuracyfurthermore.  Keywords:standardorificeflowmeter;flowfieldcalculation;numericalsimulation;computationalfluiddynamics(CFD)

0 引 言

孔板流量计因其结构简单、耐用而成为目前国际上标准化程度最高、应用最为广泛的一种流量计,但也存在着流出系数不稳定、线性差、重复性不高等缺点。采用CFD数值模拟来分析研究管内孔板类节流元件的相关流场在国外已有数十年的历史[1-5]。Sheikholeslami等人和Barry等人使用Fluent软件模拟了孔板流量计工作性能随着雷诺数、直径比、管道表面粗糙度、上游旋涡以及上下游流动边界条件的变化情况,在使用二维轴对称模型的情况下,他们认为采用80×60的网格足以得到与经验数据相差在2%以内的流出系数值。美国TexasA&M大学的Morri-son等人采用CFD对孔板流量计进行了参数化研究,

[6-7]

并测量了孔板下游的流场[8-9];后来又采用实验和CFD模拟对计量管内部的粗糙度规格进行评定改进,认为可以通过CFD模拟来获得任意雷诺数和管壁粗

糙度下的管流通用速度分布图[10]。挪威科技大学(NTNU)的Erdal等人采用Phonics软件研究了充分发展流动条件下,单孔孔板下游的流场,并采用二维轴对称模型分析讨论了不同边界条件、差分算法和湍流模型等对孔板前后流场模拟结果的影响[11-12]。悉尼大学的Langrish等人利用CFX软件中的标准湍流模型,模拟了三维轴对称突扩管内雷诺数达到10时的湍流流动情况[13]。2003年,美国福特公司的研究人员采用Fluent6.0软件模拟了汽车传动液压控制系统管路用孔板节流元件的流场,讨论了流出系数与直径比、孔板横截面形状、孔板轴向厚度、孔板出入口几

5

收稿日期:2007-09;修订日期:2007-11-16

作者简介:陈家庆(1970-),男,湖北宜昌人,工学博士,教授.研究方向:计算机辅助工程(CAE)、多相流分离技术等.

52实 验 流 体 力 学               (2008)第22卷 

何形状之间的关系[14]。

近5年来,国内一些单位也开始围绕孔板类节流元件的流场问题进行数值模拟研究[15-17]。但迄今尚无人利用商业CFD软件专门针对孔板流量计的内部

流场进行系统分析,与相关经验公式进行对比讨论,因此开展这方面的工作非常有必要。

孔板流量计的流出系数C;在工程实际应用过程中,只需通过测定实际的■p值,将C、■p值代入式(2),即可得到所关心的实际体积流量值qV。对于不可压缩流体,当采用标准孔板结构时,也可不实测标定,而使用国际标准化组织(ISO)的里德-哈利斯/加拉赫公式确定流出系数,该公式是基于大量实测实验而回归出的一个经验公式。

在已知qV的前提下,可以通过CFD数值模拟得出孔板前后D和D/2截面上的压力差■p,然后将qV、■p代入式(2),求出数值模拟流出系数C′。1.2 建模与求解

运用Gambit直接建立标准孔板流量计D和D/2取压时的三维实体模型,但利用对称的特点沿轴向考虑1/2的实体。上游管段和下游管段直径D取100mm,孔板上游管段长取20D,下游管段长取10D,孔板厚度取3mm。在固-液交界壁面处(图2(a)中的线段AB、HG、IJ、OP处)进行边界层处理,边界层的第一行百分比选用15%,共5层,比例设为1.1。为了准确捕捉孔板前后流场的变化情况,以上、下游直管段内与孔板等孔径的圆柱面为分界面实施Split操作(图2(a)中所示为线段DC、FE、KL、MN),并将边界层作为一种网格加密的技巧在此予以应用:分界

[19]

[18]

1 CFD模型及计算结果

1.1 理论基础

标准孔板流量计有D和D/2取压、角接取压和法兰取压等多种方式,其中D和D/2取压法的结构

如图1所示。对于不可压缩流体的水平管流动,在忽略沿程摩擦阻力损失的情况下,

根据流体流动的伯努

图1 D和D/2取压标准孔板流量计的结构示意图Fig.1 SketchofDandD/2pressuretakinginstandardorifice

plateflowmeter

利方程(能量守恒)和连续性原理,可以得出管道中流体理论体积流量QV的计算公式(1)

QV=

2··41-βρ

(1)

式中,β为孔板流量计的直径比,βD分D,d、别为孔板孔径和上游管道内径;■p′为截面S1-S1、S2-S2上平均流体压力p1、p2之差;ρ为管道内流体介质的密度。

实际上,对于不可压缩流体,下游取压口并非设置在截面S2-S2处,而是在与S2-S2有一定距离的截面S3-S3处。考虑到在截面S1-S1、S3-S3上测取的平均流体压力差■p一定大于■p′,故定义流出系数C来修正上述公式,可得实际体积流量值的计算公式

qV=C·QV=

2

··41-βρ

(2)

图2 标准孔板流量计的CFD数值模拟网格划分示意图

Fig.2 MeshillustrationforCFDnumericalsimulationinstandardori-ficemeter

,

 第2期            陈家庆等:标准孔板流量计内部流场的CFD数值模拟53

圆柱端面向外、向内边界层的第一行百分比选用15%,共10层,比例设为1.1;分别将圆环面和半圆面以Map、Pave的方式进行网格划分。上游直管段的轴向网格密度沿BA、CD、EF、GH方向以1.1的比例由密变疏,下游直管段的轴向网格密度沿IJ、KL、MN、OP方向以1.1的比例由密变疏,孔板轴线方向(CK、EM)上的网格平均分布。最后采用Cooper进行网格划分,最终所得网格划分情况如图2(b)所示。  在初始验证性算例中,取β=0.4、qv=0.5m3/h;以常温下的水作为流体介质,其密度、粘度根据Flu-ent6.2.16数据库中对应的物性参数来选取。进、出口的边界条件分别设置为速度入口、出流(outflow),取流体重力加速度沿着-Y轴方向。此时对缨的入口流速v=0.0176839m/s,则上游直管段内和节流孔处的雷诺数Re分别为1760、4400,可见流体在上游直管段内为层流流动,在节流孔内为湍流流动。为此,Fluent数值模拟时采用3ddp求解器,选择标准k-ε两方程湍流模型和强化壁面处理,离散方程组的压力速度耦合选择SIMPLE算法,动量、湍流动能、湍流耗散率均采用一阶迎风差分格式。图3(a)为内部的速度大小分布等值线图,水流经由孔板节流后,形成了一个对称的速度尖峰,中心轴线上的速度最大,

逐渐向

两边递减;在孔板的下游形成一个尖核状的速度峰,而在下游靠近壁面处形成一个回旋区。图3(b)为湍流动能的分布图,湍流动能在孔板下游区域较强,并在孔板内壁所在面附近形成双峰。  通过在Fluent中读取孔板前后D和D/2轴截面上的平均压力值ph、pl,得出Δp=14.05Pa,进而计算出数值模拟流出系数C′=0.6508;根据ISO经验公式计算出的推荐流出系数C=0.6323,两者的相对误差δ为2.93%,可见C′与C吻合较好。虽然CFD数值模拟与实验实测一样都存在着各种误差影响,但仍足以证明CFD数值模拟模型的正确性。

2 各参数变化对流出系数影响的讨论

为了研究不同流量、直径比、孔板厚度和流体介质对标准孔板流量计流动情况的影响,得出一些具有指导意义的结论,在工况温度均为300K的情况下,每次仅改变其中某个参数,利用Fluent进行数值模拟和相关分析讨论。如不特别指出,所建三维模型边界层的第一行百分比都选用15%。

2.1 流量的影响

以水为流体介质,对β=0.5、E=3mm的标准孔板流量计,根据具体的流动情况在Fluent求解器中选用层流或k-ε湍流模型,但网格划分模式不变(即都采用边界层网格加密处理)。得出不同流量下的流场计算结果如表1所示。

由表中可以看出,在包含层流、过渡流和湍流状态的不同流量下,数值模拟流出系数C′与ISO公式流出系数计算值C均吻合得较好,并且在层流状态下(Re<2100),C′与C的误差维持在2%以内。随着流量的不断增大,C逐渐减小,而C′随着流动状态的不同其变化规律也不同。在层流状态下,随着流量的增大C′逐渐减小;在过渡流状态下(2100<Re<4000),随着流量的增大C′逐渐增大;在湍流状态下(Re>4000),随着流量的增大C′逐渐减小。并且,在湍流状态下,C′始终大于C。2.2 直径比β的影响

标准孔板流量计的直径比β一般在0.2~0.75之间选取。图4为以水为流体介质,在Q=10m/h、E=3mm时,不同直径比β下流动情况的计算结果。

3

图3 主要流动参数在纵向对称轴剖面上的分布图Fig.3 Distributionofsomemainflowparametersonthelongitu-dinalsymmetryplane

由于β较小时边界层第一行百分比过大易使模拟效果失真,故边界层第一行的百分比都选用5%。于

表1 标准孔板流量计在不同流量下的流场计算结果

Table1 Calculationresultsofflowfieldinstandardorificeplateflowmeterunderdifferentflowrate

流量Q/(m3·h-1)上游直管段Re

0.2

0.30.40.50.70.80.91.02.04.06.08.010.0

704

1056140817602464281631683520704014080211202816035198

ph/Pa101324.90101324.80101324.70101324.55101324.25101324.05101323.85101323.60101320.55101310.85101298.80101281.55101260.00

pl/Pa101324.15

101322.95101321.15101318.90101313.00101309.45101303.55101301.30101232.65100958.20100485.6599817.1798957.35

■p/Pa0.75

1.853.555.6511.2514.6018.3022.3087.90352.65813.151464.382302.65

C0.7043

0.67640.66180.65260.64170.63810.63530.63290.62150.61470.61200.61040.6094

C′0.7067

0.67500.64970.64370.63870.64070.64380.64800.65280.65180.64390.63970.6377

相对误差δ0.24%

0.21%1.83%1.36%0.47%0.41%1.34%2.39%5.04%6.04%5.21%4.80%4.64

%

0.02D之间、孔板厚度E应在e与0.05D之间,对应

本文即为0.5mm<e<2mm、e<E<5mm。图5为以水作为流动介质,在Q=10m3/h、β=0.5时,选取不同孔板厚度E的计算结果。为了使不同厚度孔板之间的网格具有可比性,采用按孔板厚度等比例确定网格数量的方式,使不同厚度孔板上的网格尺寸基本相等。由图中可以看出,孔板厚度对推荐流出系数

C

图4 标准孔板流量计在不同直径比β时流出系数的计算结果

Fig.4 Calculationresultsofdischargecoefficientinstandardori-ficeplateflowmeterunderdifferentβvalue

ISO公式流出系数计算值C,且C值始终随着β的增大而增大。在β值较小时,随着β的增大C′有降低的现象;在β=0.55时,C′较为反常,出现了最小值。但总体上来看,C′呈上升趋势。2.3 孔板厚度的影响

根据相关规定,节流孔的厚度e应在0.005D与

图5 标准孔板流量计在不同孔板厚度E时流出系数的计算

结果Fig.5 Calculationresultsofdischargecoefficientinstandardori-ficeplateflowmeterunderdifferentorificethickness

表2 标准孔板流量计在改变流体介质的流场计算结果

Table2 Calculationresultsofflowfieldinstandardorificeplateflowmeterunderdifferentflowmedia

流体介质戊烷己烷庚烷辛烷甲醇甲苯液苯水煤油柴油

密度ρ/(kg·m-3)

626660684720785866875998.2780830

粘度μ×104

/(N·s·m-2)

2.293.24.095.45.4955.865.8910.0324.033.2

ph/Pa101291.65101287.25101283.55101280.10101276.40101271.50101270.90101260.00101250.60101238.05

pl/Pa99829.6599745.0799685.9899607.2199448.4499253.5299232.5798957.3599553.9399431.35

■p/Pa1462.001542.181597.571672.891827.962017.982038.332302.651696.671806.70

C0.60630.60700.60760.60830.60810.60800.60790.60940.61630.6189

C′0.63380.63360.63380.63540.63470.63450.63460.63770.65670.6565

相对误差

δ4.54%4.38%4.31%4.46%4.37%4.36%4.39%4.64%6.56%6.08%

的影响很小,基本维持不变;而随着孔板厚度的增大,数值模拟流出系数C′逐渐增大。当孔板为零厚度时,C′小于C;但随着孔板厚度的增加,C′逐渐超过C,并且其差值逐渐增大。其实不难想象,随着孔板厚度的增大,由于狭长通道对高速水流有一个强制性的收聚汇集作用,下游形成的对称速度尖峰越长,因此下游取压口处的压力值必然变小,致使流出系数相应增大。

2.4 流体介质物性参数的影响

取Q=10m3/h、β=0.5、E=3mm,计算出不同流体介质下标准孔板流量计的流动情况如表2所示。由表中可以看出,流出系数与流体介质的种类密切相关。一般而言,随着流体介质的粘度不断增大,C与C′基本上保持着逐渐增加的趋势,并且C′始终大于C。

[8][7][4]

[5]

[6]

3 结论与展望

(1)采用CFD数值模拟可以有效获取标准孔板流量计内部的流场分布情况,并可根据具体应用场合计算出相应的流出系数,因此能够弥补因受各种局限而不能进行流出系数实测标定的缺憾和不足;

(2)随着流量、直径比、孔板厚度以及流体介质的改变,流出系数也会随之发生变化。一般情况下,ISO公式流出系数计算值小于数值模拟流出系数,且难以准确反映出孔板轴向厚度变化对内部流场的影响;

(3)虽然采用CFD数值模拟与实测标定一样都存在着误差,但应用本文所建立的CFD数值模拟模型可以对孔板上游直角区污物堆积、孔板轴向入口锐角变钝等几何形状变化对流动情况的影响进行研究,还可以直接推广到喷嘴、文丘里管等节流差压式流量计的分析研究中。参考文献:

[1] DAVISRW,MATTINGLYGB.Numericalmodelingofturbu-lentflowthroughthinorificeplate[A].Symp.onFlowinOpen

ChannelsandClosedConditions[C],Gaithersburg,23-25,1977.[2] READER-HARRISMJ,parisonofcompu-tationandLDVmeasurementofflowthroughorificeandperfo-ratedplates,andcomputationoftheeffectofroughpipework

onorificeplates[A].Proc.oftheInt.Symp.onFluidFlowMeasurement[C].Washington,DC,USA,1986.[3] putationofflowthroughorifice

[9]

[10]

[11]

[12]

[13]

[14]

[15]

[16][17]

[18]

[19]

plates[J].NumericalMethodsinLaminarandTurbulent

Flow,6:1907-1917,1989.

 DRUSTF,WANGAB,FOUNTIM.Similarityphenomenaandcomputationsoftheflowthroughanaxisymmetricring-typeobstacleattachedtoapipewall[A].Hydrocomp89Dubrovnik,1989:484. DURSTF,WANGAB.Experimentalandnumericalinvestiga-tionoftheaxisymmetric,turbulentpipeflowoverawall-mountedthinobstacle[A].SeventhSymposiumonTurbulentShearFlows[C].StanfordUniversity,21-23,1989. SHEIKHOLESLAMIMZ,PATELBR,KOTHARIK.Numeri-calmodelingofturbulentflowthroughorificemeters—apara-metricstudy[A].2ndInternationalConferenceonFlowMea-surement[C].London,UK,11-13,1988.

 BARRYJJ,SHEIKOLESLAMIMZ,PATELBR.Numericalsimulationofflowthroughorificemeter[R].GasResearchIn-stitute,GRI-92/0060.1,1992.

 MORRISONGL,PANAKDL,DEOTTERE.Numericalstudyoftheeffectofupstreamflowconditionuponorificeflowmeterperformance[A].OffshoreMechanicsandArcticEngi-neeringProceedings(OMAE)[C],1992.

 MORRISONGL,DEOTTERE,NAILGH,etal.Meanve-locityandturbulencefieldsinsideaβ=0.50orificeflowme-ter[A].AIChESpringNationalMeeting[C].1992. MORROWTB,MORRISONGL.Effectofmetertuberough-th

nessonorificeCd[A].The4InternationalSymposiumonFluidFlowMeasurement[C].Denver,CO,28-30,1999. ERDALA,TORBERGSENLE,RIMESTADS,etal.Evalua-tionofaCFDmodelforsimulationofsimplifiedflowcondi-tioners[A].FluidFlowMeasurement3rdInternationalSympo-sium[C].SanAntonio,1995.

 ERDALA,ANDERSSONTHI.Numericalaspectsofflowcomputationthroughorifices[J].FlowMeas.Instrum.,1997,8(1):27-37.

 BAOYUGUO,TIMAG,DAVIDF.Simulationofprecessioninaxisymmetricsuddenexpansionflows[A].TheSecondIn-ternationalConferenceonCFDintheMineralsandProcessIndustries[C].CSIRO,Melbourne,Australia,1999.

 CHEOLWANKim,DENGFUCHIM-YuanPERNG,ZHANGDengfu,etal.OrificeFlowCharateristics[J].FordTechnicalJournal,2003,6(3):1-15. 程耕,程平,李受人.节流管孔流动参数与雷诺数关系的数值研究[J].计算机工程与设计,2005,26(3):575-607.

 程勇,汪军,蔡小舒.低雷诺数的孔板计量数值模拟及其应用[J].计量学报,2005,26(1):57-59.

 GENGYANFENG,ZHENGJINWU,SHITIANMING.Studyonthemeteringcharacteristicsofaslottedorificeforwetgasflow[J].FlowMeasurementandInstrumentation,2006,17:123-128.

 GB/T2624-93.流量测量节流装置用孔板、喷嘴和文丘里管测量充满圆管的流体流量[S].1993-02-03批准,1993-08-01实施. 孙淮清,王建中.流量测量节流装置设计手册[M].北京:化学工业出版社,2005.

本文来源:https://www.bwwdw.com/article/emo4.html

Top