CN102087130B - 基于cfd技术的多声路超声流量计弯管安装的声路优化方法 - Google Patents

基于cfd技术的多声路超声流量计弯管安装的声路优化方法 Download PDF

Info

Publication number
CN102087130B
CN102087130B CN201010549707XA CN201010549707A CN102087130B CN 102087130 B CN102087130 B CN 102087130B CN 201010549707X A CN201010549707X A CN 201010549707XA CN 201010549707 A CN201010549707 A CN 201010549707A CN 102087130 B CN102087130 B CN 102087130B
Authority
CN
China
Prior art keywords
sound travel
pipeline
section
flowmeter
flow
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201010549707XA
Other languages
English (en)
Other versions
CN102087130A (zh
Inventor
何存富
刘飞
吴斌
焦敬品
张力新
樊尚春
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huizhong Instrumentation Co ltd
Beijing University of Technology
Original Assignee
TANGSHAN HUIZHONG INSTRUMENT CO Ltd
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by TANGSHAN HUIZHONG INSTRUMENT CO Ltd, Beijing University of Technology filed Critical TANGSHAN HUIZHONG INSTRUMENT CO Ltd
Priority to CN201010549707XA priority Critical patent/CN102087130B/zh
Publication of CN102087130A publication Critical patent/CN102087130A/zh
Application granted granted Critical
Publication of CN102087130B publication Critical patent/CN102087130B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measuring Volume Flow (AREA)

Abstract

一种基于CFD技术的多声路超声流量计弯管安装声路的优化方法,主要包括,由管路情况计算雷诺数、迪安数及边界层厚度;利用三维造型软件对管路及流量计建模,并根据不同的安装角对管路模型与流量计模型进行装配;选取CFD软件,将装配好的模型导入并进行参数设置;模拟不同流速的流场分布,记录入口平均流量,计算声路断面流量;计算测得流量误差;对计算的各组流量误差进行线性拟合,并根据拟合斜率修正各组误差,筛选在该流速范围内误差值波动范围最小的一组,以该组的声路断面、声路数和安装角度作为该管路特征条件下的最优声路布置方案。本发明可以指导在不具备直管段安装条件下,合理选取多声路超声流量计的声路数量、测量断面及安装角度。

Description

基于CFD技术的多声路超声流量计弯管安装的声路优化方法
技术领域
本发明涉及一种基于CFD技术多声路超声流量计在弯管段安装的声路优化方法,尤指利用CFD软件模拟弯管中的流场分布,调整多声路超声流量计的各个参数,使得多声路超声流量计的测试精度能够适应流场的分布。 
背景技术
随着能源和水资源的全球性匮乏,一批关系国计民生的大型水利工程和引水调水工程在我国迅速发展,如三峡水利枢纽、南水北调工程等。这些工程项目中经常包含一些口径和流量都很巨大的管路,如水电站机组进水管路等,常规流量计无法适应。近年开发应用的多声路超声流量计,较好地解决了大口径水流量测量的技术难题,流量计制造不受管路口径的限制,多声路配置可以适应较为复杂的管路结构和流态分布,故超声流量计已成为大口径水流量测量的最佳技术选择。对于多声路超声流量计测量标准直管段的流量,已有较成型的模式,而对于不具备直管段安装要求,需要在弯管处安装时,如何保证测量精度的研究几乎处于空白。由于弯管段流态的扭曲,声路数量、布置方式及安装角度对多声路超声流量计测量精度影响显著,如何根据流场的分布情况配置最优声路布置方案,对于测量结果显得尤为重要。 
发明内容
本发明的目的在于,通过提供一种基于CFD技术多声路超声流量计在弯管段安装的声路优化方法,对不同流场分布下多声路超声流量计的安装效应给出最优声路布置方案,进而满足测试精度的需求。 
该优化方法基于对湍流模型的理论分析,根据实际管路边界条件进行建模,并利用CFD软件分析管路内流量计安装部位的流场分布。在流场分布确定的情况下,对不同的声路数量、测量断面数和安装角度进行数值模拟,得出各个声路的速度平均值,并利用Gauss-Jacobi数值积分法去计算管路截面的平均流量。将不同进口流速计算的管路截面流量进行线性拟合,利用拟合得到的斜率进行误差修正,确定在该条件下误差的波动范围,选取误差波动范围最小的声路数量、测量断面数和安装角度作为该管路特征条件下流量计安装部位的最优声路布置方案。 
本发明基于CFD技术多声路超声流量计在弯管段安装的声路优化方法包括以下步骤: 
1)根据流体性质、管路环境、流速范围等,利用湍流理论计算雷诺数 
Figure BSA00000351150100021
弯管处的迪安数 边界层厚度 
Figure BSA00000351150100023
2)对整体管路结构进行简化,保留弯头前直管段10D以上,弯头后直管段5D以上(D-管路直径),并利用三维造型软件对简化后的管路建模; 
3)根据欲安装的多声路超声流量计干流标定的各项数据,利用三维造型软件对不同声路数和断面数的流量计进行建模,使模型与流量计实物各项参数均相同; 
4)选取合适的CFD软件,并对管路模型和流量计模型根据不同的安装角进行装配,导入CFD软件进行参数设置; 
5)利用CFD软件对不同的进口流速进行模拟仿真,记录入口平均流量,并提取流量计声路断面的流场分布; 
6)在该流场分布条件下,利用数值积分法计算各声路的轴向速度和横流速度, 并根据式V=Vaxial±Vctgφ计算流量计实测声路速度,其中:Vaxial为声路的轴向速度,Vc为声路的横流速度,φ为声路夹角;再由Gauss-Jacobi数值积分法计算该组声路断面的流量 
Figure BSA00000351150100031
其中:R为管路半径,wi为i声路的权系数,Vi(ri)为i声路的安装位置; 
7)按照步骤6)的方法逐一计算每组声路的截面平均流速,进而可以计算每组声路的流量误差 其中:qmeas为声路断面平均流量,qact为入口平均流量,Vmeas为声路断面平均流速,Vact为入口平均流速; 
8)对在相同安装条件下的各组误差值,以不同的流速值作为横坐标进行线性拟合,计算拟合斜率; 
9)根据拟合斜率修正各组误差,筛选在该流速范围内误差值波动范围最小的一组,以该组的声路断面、声路数和安装角度作为该管路情况的最优声路布置方案。 
与现有多声路超声流量计声路布置优化方法相比,本发明具有以下优点:可以对任意管路流场进行快速、有效地模拟计算,特别是弥补了当前弯管段流量计安装的声路优化方法的空白;可以节省大量的实验,有效地指导流量计出厂时的标定工作。 
附图说明
图1为基于CFD技术声路优化方法步骤框图示意图; 
图2为多声路超声流量计弯管安装实验模型示意图; 
图3为弯管简化模型示意图; 
图4为18声路超声流量计实物示意图; 
图5为18声路超声流量计仿真计算模型; 
图6为0度安装角各组声路仿真线性拟合结果示意图; 
图7为45度安装角各组声路仿真线性拟合结果示意图; 
图8为90度安装角各组声路仿真线性拟合结果示意图。 
具体实施方式
以下结合说明书附图对本发明结合现场实验的仿真优化实例进行说明,具体步骤如图1所示: 
1)根据实验现场提供的数据,管路内径为d=399.489mm,弯管曲率半径为440mm入口流速为0.3m/s~4m/s,流体为水,环境温度为19.7°;经查表计算该温度下水的运动黏度v=1.007×10-6m2/s,由雷诺数计算公式 
Figure BSA00000351150100041
可知雷诺数在1.19×105-1.59×106之间;由迪安数计算公式 
Figure BSA00000351150100042
可知迪安数在1.13×105-1.52×106之间;由于试验管路均为光滑管,按照修正后的普朗特光滑管公式 
Figure BSA00000351150100043
计算沿程阻力系数λ在0.013-0.018之间,沿程阻力系数随着雷诺数的增大而减小,由边界层厚度公式 可得0.076mm≤δ≤0.821mm。 
2)图2为试验所用弯管实物模型图,根据弯管中安装流量计的位置对其进行简化,保留弯头前直管段10D,弯头后直管段5D(D-管路直径),图3为利用Pro/E所建的弯管模型。 
3)试验所用超声流量计为18声路双断面,每个断面9声路,两个断面对称分布,流量计各声路的实测声路长度及声路夹角如表1。图4为流量计实物图,图5为利用Pro/E所建流量计模型。 
表1声路长度及声路夹角 
4)选用Comsol Multiphysics 3.5a化学工程模块湍流分析中的经典k-ε湍流模型。在Pro/E中对0°、45°和90°三个安装角进行装配,并分别导入Comsol中按照步骤(1)中计算的参数进行设置。 
5)首先对在0°安装角下18声路双断面进行了仿真计算,仿真结果能充分地展现管路内各处的流速分布情况;在经过弯管后,流态会发生很大的变化,且在弯管部分会产生局部流态的滞留现象;流量计管段内的流场受探头插入的影响较大,从仿真结果能清楚地看到在探头附近流速很小,甚至个别地方会出现倒流的现象;同时从仿真的后处理从还可以看到横流的分布情况,证实横流对于实测流速的影响是不可忽略的。利用同样的方法,对45°安装角和90°安装角18声路双断面进行仿真计算,从不同安装角流场分布的对比来看,安装角对管路整体的流场分布影响不大,但对流量计管段的流场分布影响显著。 
6)在后处理中调用线积分计算各声路的轴向速度和横流速度,并根据表1中 对应的声路角,利用V=Vaxial±Vctgφ计算流量计声路速度;根据表2所示流量计各个声路对应的权系数,由Gauss-Jacobi数值积分式 
Figure BSA00000351150100061
计算A断面与B断面的截面平均流速,表3为0°安装角仿真数据结果。 
表2声路分布位置与所对应的权系数 
Figure BSA00000351150100062
表3仿真数据 
Figure BSA00000351150100063
7)由流量误差计算式 
Figure BSA00000351150100071
计算各组截面实测流速与标准流速的相对误差。 
8)对在同一安装角,声路断面数相同计算的各组误差值,以不同的流速值作为横坐标进行线性拟合,表4为各组声路线性拟合结果。 
表4各声路线性拟合结果 
Figure BSA00000351150100072
9) 
Figure BSA00000351150100073
修正各组流速,其中V′meas为修正后的流速,a为拟合斜率,b为截距。图6、图7和图8分别为各个安装角修正前后的流速计算的误差值。 
结果表明:(1)仿真测得的流速经修正后,平均误差明显的减少,误差波动幅度得到了显著的改善。(2)在固定安装角的情况下,声路数量对误差产生的影响并不明显,使用A、B双断面测量与单断面测量的结果相差无几。(3)随着安装角度 的增大,误差有明显的上升趋势。在安装角为0°时,所有声路的误差都可控制在±0.5%之内;在安装角为90°时,多数声路在低流速下的误差接近±3%,此时声路数及断面数的增加并不能提高整体的精度,可见安装角度是影响低流速测量精度的主要因素。当流速高于1m/s时,安装角度对其的影响不再明显,除90°安装角在流速为2m/s时误差接近-1%外,其余声路的误差均可控制在±0.5%附近,此时,8声路双断面的效果要优于9声路单断面。 
因此可以通过仿真来优化声路的布置,进而根据不同的精度要求选用最佳的声路布置方案。 
最后应说明的是:以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。 

Claims (2)

1.基于CFD技术的多声路超声流量计弯管安装的声路优化方法,对不同的声路数量、测量断面数和安装角度进行数值模拟,计算管路截面的平均流量;将不同进口流速计算的管路截面流量进行线性拟合,进行误差修正,确定在声路数量、测量断面数和安装角度不同组合条件下的误差波动范围,选取误差波动范围最小的声路数量、测量断面数和安装角度作为具有与模拟管路特征条件相同的实际管路中流量计安装部位的最优声路布置方案;其特征在于包括以下步骤:
1)由管路情况计算雷诺数、迪安数及边界层厚度;
2)利用三维造型软件对管路及流量计建模,并根据不同的安装角对管路模型与流量计模型进行装配;
3)选取CFD软件,将装配好的模型导入并进行参数设置;
4)模拟不同流速的流场分布,记录入口平均流量,并由Gauss-Jacobi数值积分法计算声路断面流量;
5)将计算的各组声路断面流量与入口平均流量代入误差计算式
e = ( q meas - q act ) q act × 100 % = ( V meas - V act ) V act × 100 % , 计算测得流量误差;
其中:qmeas为声路断面平均流量,qact为入口平均流量,Vmeas为声路断面平均流速,Vact为入口平均流速;
6)对计算的各组流量误差进行线性拟合,并根据拟合斜率修正各组误差,筛选在计算流速范围内误差值波动范围最小的一组,以该组的声路断面、声路数和安装角度作为该管路情况的优化声路布置方案。
2.根据权利要求1所述基于CFD技术的多声路超声流量计弯管安装的声路优化方法,其特征在于:所述流量计模型数据全部来源于实测管路标定数据和流量计干式标定数据。
CN201010549707XA 2010-11-19 2010-11-19 基于cfd技术的多声路超声流量计弯管安装的声路优化方法 Active CN102087130B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201010549707XA CN102087130B (zh) 2010-11-19 2010-11-19 基于cfd技术的多声路超声流量计弯管安装的声路优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201010549707XA CN102087130B (zh) 2010-11-19 2010-11-19 基于cfd技术的多声路超声流量计弯管安装的声路优化方法

Publications (2)

Publication Number Publication Date
CN102087130A CN102087130A (zh) 2011-06-08
CN102087130B true CN102087130B (zh) 2012-02-15

Family

ID=44099031

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010549707XA Active CN102087130B (zh) 2010-11-19 2010-11-19 基于cfd技术的多声路超声流量计弯管安装的声路优化方法

Country Status (1)

Country Link
CN (1) CN102087130B (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102538912B (zh) * 2011-11-17 2014-04-16 中国计量科学研究院 超声流量计流场附加误差分析方法
CN102750426B (zh) * 2012-07-19 2014-10-22 华南理工大学 基于cfd技术的流动调整器整流效果评判方法
CN103940495B (zh) * 2014-04-14 2017-03-08 重庆大学 基于流线的小流量超声流量计误差估算方法
CN105698973B (zh) * 2016-02-21 2018-12-07 上海大学 一种超声波热量表信号转换器的校验装置
CN105890684A (zh) * 2016-02-29 2016-08-24 上海安钧电子科技有限公司 一种采用高斯-雅可比多项式确定声道位置的设置新方法
DE102016013328A1 (de) * 2016-11-10 2018-05-17 Nivus Gmbh Verfahren zur Verbesserung der Messgenauigkeit einer Durchfluss-Messanordnung
CN107014449B (zh) * 2017-04-21 2019-05-31 中国农业大学 修正泵站流量测量结果的方法
US10795054B2 (en) * 2018-03-20 2020-10-06 Mitsubishi Electric Research Laboratories, Inc. System and method for sensing wind flow passing over complex terrain
CN112504365B (zh) * 2020-11-25 2022-05-20 合肥工业大学 一种电磁流量传感器的磁路结构优化设计方法
CN112857491A (zh) * 2021-02-10 2021-05-28 北京市南水北调南干渠管理处 一种异型输水管涵的流量测量方法及系统
CN113642274B (zh) * 2021-10-14 2022-01-14 四川大学 一种基于流场模型的河流流量计算方法
CN116698141B (zh) * 2023-07-28 2023-10-27 山东大学 一种不同工况下超声波流量计的测速误差修正方法及系统
CN117664261B (zh) * 2024-02-01 2024-04-05 安徽天康(集团)股份有限公司 一种弯管侧置式气体流量仪

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7299140B2 (en) * 2005-12-14 2007-11-20 Thermo Fisher Scientific Inc. Method and system for multi-path ultrasonic flow measurement of partially developed flow profiles
KR100689771B1 (ko) * 2006-05-10 2007-03-12 주식회사 수인테크 다회선초음파유량계의 제조방법
JP4898724B2 (ja) * 2008-03-06 2012-03-21 パナソニック株式会社 超音波式流体計測装置の多層流路部材の製造方法

Also Published As

Publication number Publication date
CN102087130A (zh) 2011-06-08

Similar Documents

Publication Publication Date Title
CN102087130B (zh) 基于cfd技术的多声路超声流量计弯管安装的声路优化方法
CN109344436B (zh) 一种大型复杂天然气管网系统在线仿真方法
CN103940495B (zh) 基于流线的小流量超声流量计误差估算方法
CN102750426B (zh) 基于cfd技术的流动调整器整流效果评判方法
CN103455725B (zh) 管网系统非恒定流模拟方法
Zheng et al. Study of acoustic transducer protrusion and recess effects on ultrasonic flowmeter measurement by numerical simulation
CN103074873A (zh) 一种非理想流场下的多声道超声流量计声道布置方法
CN106874595B (zh) 基于节点参数化技术的输水管网计算方法
CN107014449B (zh) 修正泵站流量测量结果的方法
CN111784536A (zh) 一种明渠中依据实测水位变化情况估计水位超限时间的方法
CN107014451A (zh) 基于广义回归神经网络推测超声波流量传感器系数的方法
CN109163159A (zh) 用于变截面弯头的导流部件及其制作方法
CN203672427U (zh) 一种航空发动机台架燃油流量测量装置
CN102630426A (zh) 文丘里施肥器的结构优化方法
CN117592257A (zh) 一种闸控水网数字孪生模型的构建方法及系统
CN202255483U (zh) 一种均衡式流量计
CN108506622B (zh) 一种基于弧形导流片的低阻力三通构件
CN108280300A (zh) 基于计算流体力学的大数据流量计开发方法
CN110487346B (zh) 一种大流量低温推进剂供应管路用整流孔板及其设计方法
CN108416161A (zh) 一种基于三维模型的管道水力计算方法及装置
CN105444828A (zh) 超声波流量计的整流器的设计方法
JP2007107948A (ja) 流量測定方法、流量測定装置およびこの流量測定装置を備えた水力機械
Fei et al. Numerical simulation of multi-path ultrasonic flowmeter: Ultrasonic path error analysis
Krzemianowski et al. Use of non-uniform rational b-splines for discharge calculation in the velocity area method
CN103033221A (zh) 方管式平衡流量计

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: TANGSHAN HUIZHONG INSTRUMENTATION CO., LTD.

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20111103

Address after: 100124 Chaoyang District, Beijing Ping Park, No. 100

Applicant after: Beijing University of Technology

Co-applicant after: Tangshan Huizhong Instrumentation Co.,Ltd.

Address before: 100124 Chaoyang District, Beijing Ping Park, No. 100

Applicant before: Beijing University of Technology

C14 Grant of patent or utility model
GR01 Patent grant
CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee after: Beijing University of Technology

Country or region after: China

Patentee after: HUIZHONG INSTRUMENTATION CO.,LTD.

Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee before: Beijing University of Technology

Country or region before: China

Patentee before: Tangshan Huizhong Instrumentation Co.,Ltd.

TR01 Transfer of patent right

Effective date of registration: 20240124

Address after: 063000 No.126, Gaoxin West Road, Tangshan hi tech Industrial Development Zone, Hebei Province

Patentee after: HUIZHONG INSTRUMENTATION CO.,LTD.

Country or region after: China

Address before: 100124 No. 100 Chaoyang District Ping Tian Park, Beijing

Patentee before: Beijing University of Technology

Country or region before: China

Patentee before: HUIZHONG INSTRUMENTATION CO.,LTD.

TR01 Transfer of patent right