CN103353908B - 一种基于数值计算的管路阻力系数精确计算方法 - Google Patents
一种基于数值计算的管路阻力系数精确计算方法 Download PDFInfo
- Publication number
- CN103353908B CN103353908B CN201310245810.9A CN201310245810A CN103353908B CN 103353908 B CN103353908 B CN 103353908B CN 201310245810 A CN201310245810 A CN 201310245810A CN 103353908 B CN103353908 B CN 103353908B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- pipe
- numerical computations
- pipeline
- line system
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
采用ANSYS‑CFX,对管路进行数值模拟,并将数值模拟中设定的温度和计算得到的单位质量流体能量损失与实验测量温度和单位质量流体能量损失进行数据对比;通过对比,修正了数值模拟中设定的温度和动力粘度,并最终选用最优湍流模型,对管路在不同粗糙度、不同流量情况下进行三维粘性定常数值模拟。将模拟得到的各过流截面总压带入伯努利方程,利用能量损失反求管路不同情况下的沿程阻力系数和局部阻力系数;将数据导入TableCurve3D中进行三维数据曲面拟合,最终得到了沿程和局部阻力系数关于相对粗糙度和雷诺数的最优表达式;本发明不仅能较为准确地计算管路局部和沿程阻力系数,还可以在此基础上建立起管路系统能量性能计算模型并对已有的管路系统进行节能改造。
Description
技术领域
本发明属于管路系统能量性能计算领域,具体涉及一种基于数值计算的管路阻力系数精确计算方法。
背景技术
目前管路系统阻力系数取值仍然是按照上世纪四十年代人们总结的经验公式计算求得。其中沿程阻力系数的计算公式主要是尼古拉兹实验总结的计算公式,尼古拉兹在实验中人为的在管壁上粘结颗粒砂粒,做成人工粗糙管,从而对不同管径、不同流量的管路进行了实验得出沿程阻力系数的经验计算公式。但这种人工粗糙管与工程实际中的粗糙管路系统有一定差别,而且管路不同,粗糙度也会不同,选择尼古拉兹经验公式计算沿程阻力系数肯定会存在一定的误差。而对于计算弯头局部阻力系数的经验公式直接忽略了管路工况和粗糙度的影响。利用经验公式计算的阻力系数引入管路系统能量性能计算,必然也会造成一定误差,这对泵选型和系统节能改造都会带来不利的影响。因此建立一种基于数值计算的管路阻力系数精确计算方法是很有必要的。
迄今为止,尚未见基于数值计算的管路系统阻力系数精确计算方法的文献报道,仅有一些学者做与管路沿程和局部阻力系数计算相关的研究工作。黑龙江水利科技《关于管道沿程阻力系数问题的分析》(2012年第10期)提出利用EXCEL软件来求解柯列勃洛克公式,从而为管道沿程阻力系数计算提供了一种简便方法,即利用Excel中求解勃洛克公式的简便方法,这种方法只是简单地利用EXCEL对已有勃洛克公式进行求解,并没有提出自己有关阻力系数的精确计算公式。系统仿真学报《以分布参数法确定低雷诺数局部阻力系数》(2009年第9期)以分布参数法建立起局部损失的数学模型,实例计算了低雷诺数下直角转向的局部阻力系数,然而此研究只针对于低雷诺数局部流态,并不适用于其他工况。华北水利水电学院学报《灰色线性回归组合模型在计算给水管道阻力系数中的应用》(2009年第1期)运用灰色线性回归组合模型进行求解,得出了计算管道阻力系数关于管径和服役时间的数学模型,但没有涉及到管道粗糙度和不同工况的影响。
发明内容
本发明目的提供一种基于数值计算的管路阻力系数精确计算方法,通过数值计算的方法来模拟管路系统的流动状态,并在此基础上建立管路阻力系数精确计算方法。
为达到以上目的,采用如下技术方案:
将数值计算设定温度和单位质量流体能量损失与实验测得温度和单位质量流体能量损失进行对比,进一步调整计算模型,然后基于调整的计算模型,对不同粗糙度、不同工况管路进行数值计算,通过数值计算得到管路系统各过流截面总压p,并引入伯努利能量方程,利用能量损失反求管路不同情况下的沿程阻力系数λ和局部阻力系数ζ,得到数据样本。基于数据处理软件对数据进行三维数据曲面拟合,数据处理软件可为TableCurve3D,从而建立一种管路局部和沿程阻力系数的精确表示方法。
其具体步骤如下:
(1)对管路系统进行数值计算,计算出管路系统单位质量流体能量损失。在相同工况下,对管路系统进行外特性实验,测出实验中流体温度和管路系统单位质量流体能量损失,通过对比数值计算中设定的温度和实验测量温度,修正数值计算中设定的温度和动力粘度,再通过对比数值计算和实验测量的单位质量流体能量损失,调整数值计算湍流模型,最终选择计算能量损失与实验测量能量损失之间误差最小时的湍流模型为最优湍流模型。
(2)利用(1)中调整后的计算模型对管路系统进行变工况、变粗糙度数值计算。
(3)由数值计算得到两组过流截面(即四个面)各自总压p,将两组过流截面总压引入伯努利能量方程,联合建立伯努利方程组,即:
其中,z各截面中心线物理高程;α为动能修正系数,对于一般工业管道,取α=1;v为管路流速;ghw表示单位质量流体流经两截面途中的能量损失,包括沿程损失和局部损失。则由能量损失即可计算出不同工况、不同粗糙度管路系统的局部阻力系数ζ和沿程阻力系数λ。
(4)由量纲分析法可知,阻力系数只依变于无量纲数雷诺数Re和相对粗糙度Δ/d,利用TableCurve3D,将计算得到的局部和沿程阻力系数分别与雷诺数Re和相对粗糙度Δ/d进行三维曲面拟合,从而得到管路局部和沿程阻力系数的精确表示方法:
本发明的优点在于:
(1)过程中选用的管路系统既考虑了沿程损失,又考虑了局部损失,并同时计算出沿程和局部阻力系数,更具有工程实际意义。
(2)建立的管路阻力系数精确计算方法,适用于不同管径、不同粗糙度、不同工况管路系统阻力系数的精确计算。
(3)不仅能较为准确地计算管路局部和沿程阻力系数,还可以在此基础上建立起管路系统能量性能计算模型并对已有的管路系统进行节能改造。
附图说明
图1为一种基于数值计算的管路阻力系数精确计算方法的流程图
图2为本发明实施例的管径为450mm的管路系统
图3为本发明实施例的管路局部阻力系数关于相对粗糙度和雷诺数的三维曲面拟合图
图4为本发明实施例的管路沿程阻力系数关于相对粗糙度和雷诺数的三维曲面拟合图
具体实施方式
实施例:
结合图1的流程图,管径d=450mm的管路系统,由直管段和局部直角弯头组成,其进、出口处于同一物理高程,如图2所示。
(1)对管路系统进行数值计算和外特性实验,将计算设定温度与实验测定温度进行对比,修正数值计算温度和动力粘度,将计算得到的单位质量水能量损失与实验测得的单位质量水能量损失进行对比,选取最优湍流模型。
在1200m3/h工况下,对管路系统进行数值计算,初步选取计算流体为25°C清水,动力粘度μ=8.899×10-4pa·s,湍流模型为k-ω,计算出管路系统单位质量水能量损失为4.545J。在相同工况下,对管路系统进行外特性实验,测得实验中水的温度为20℃,故将数值计算中水温度修正为20℃,动力粘度修正为μ=1.005×10-3pa·s。实验测得管路系统单位质量水能量损失为4.518J,通过对比数值计算和实验测量的能量损失,调整数值计算湍流模型并重新计算,最终选择计算能量损失与实验测量能量损失之间误差最小时的k-ε湍流模型为最优湍流模型,此湍流模型下计算得到的单位质量水能量损失为4.525J,与实验结果误差为1.5%。
(2)利用调整后的计算模型对管路系统进行变工况、变粗糙度数值计算。
进口边界条件设置为质量流量进口边界,流动方向为边界法向,分别对管路系统400m3/h、800m3/h、1200m3/h、1600m3/h、2000m3/h五个流量进行数值计算。出口边界条件设置为静压出口,相对压力设定为105pa。
固壁上满足无滑移条件,即速度u,v,w=0;压力取为第二类边界条件,即对应上面五个不同流量工况,壁面粗糙度都分别选取30μm、50μm、70μm、90μm、120μm五个等级。即总共计算了25个算例。
(3)分别选取两组过流截面,一组为管路进、出口截面,另一组为第4段直管段进口截面和第7个直角弯头出口截面,由数值计算得到两组过流截面各自总压p,将两组过流截面总压引入伯努利能量方程,联合建立伯努利方程组,即:
其中,两组过流截面各自处于同一物理高程,即z1=z2,z3=z4;同一管路系统,各截面流速相等, hw为水头损失, hf与hj分别为沿程水头损失和局部水头损失。则由能量损失计算出不同工况、不同粗糙度管路系统的局部阻力系数ζ和沿程阻力系数λ。
(4)由量纲分析法可知,阻力系数只依变于无量纲数雷诺数Re和相对粗糙度Δ/d,利用TableCurve3D,将计算得到的局部和沿程阻力系数分别与雷诺数Re和相对粗糙度Δ/d进行三维曲面拟合(如图3、4所示),从而得到管路局部和沿程阻力系数的精确表示方法,局部和沿程阻力系数的表达式均为:
其中z为局部或沿程阻力系数,x为相对粗糙度Δ/d,y为雷诺数Re,a~j为局部和沿程阻力系数对应定常系数,两阻力系数公式中取值不同,这里:
①局部阻力系数表达式中系数为:
a=0.149139371,b=-2.6094e-05,c=35157.48679,d=2.71544e-09,e=7.30891e+08
f=-6.58571761,g=-7.6315e-14,h=-1.2302e+15,i=1.18235e+06,j=5.79135e-05
②沿程阻力系数表达式中系数为:
a=0.014049359,b=3.10644e-06,c=-8386.32301,d=-4.4003e-10,e=-5.5905e+08
f=2.221554936,g=1.22038e-14,h=4.13643e+14,i=-369966.137,j=-2.2668e-05
(5)对拟合得到的局部和沿程阻力系数表达式进行实验验证。
对管径d=450mm、粗糙度为50μm的管路系统进行外特性实验,实验中分别选取数值计算中的400m3/h、1200m3/h、2000m3/h三个工况点。通过对比发现,局部阻力系数表达式计算结果与实验结果最大偏差为1.6%,沿程阻力系数表达式计算结果与实验结果最大偏差为1.2%,计算精度较高,具有一定工程实用价值。
实验数据与阻力系数表达式计算结果对比如下:
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内,本发明要求保护范围由所附的权利要求书其等效物界定。
Claims (4)
1.一种基于数值计算的管路阻力系数精确计算方法,其特征在于,将数值计算设定温度和单位质量流体能量损失与实验测得温度和单位质量流体能量损失进行对比,进而修正数值计算中设定的温度和动力粘度值,调整数值计算湍流模型,最终选择计算能量损失与实验测量能量损失之间误差最小时的湍流模型为最优湍流模型,然后基于调整的CFD计算模型,对不同粗糙度、不同工况管路进行数值计算,通过数值计算得到管路系统各过流截面总压p,并引入伯努利能量方程,联合建立伯努利方程组,利用能量损失反求管路不同情况下的沿程阻力系数λ和局部阻力系数ζ,得到数据样本,基于数据处理软件对数据进行三维数据曲面拟合,从而分别建立了沿程阻力系数λ和局部阻力系数ζ与无量纲数雷诺数Re和相对粗糙度Δ/d的精确关系式:其中f1和f2表示关系函数,即一种管路局部和沿程阻力系数的精确表示方法。
2.根据权利要求1所述的基于数值计算的管路阻力系数精确计算方法,其特征在于,伯努利方程组为:
其中,z各截面中心线物理高程;α为动能修正系数,对于一般工业管道,取α=1;v为管路流速;ghw表示单位质量流体流经两截面途中的能量损失,包括沿程损失和局部损失。
3.根据权利要求1所述的基于数值计算的管路阻力系数精确计算方法,其特征在于,数据处理软件为TableCurve 3D。
4.根据权利要求1所述的基于数值计算的管路阻力系数精确计算方法,其特征在于,计算方法具体步骤如下:
(1)对管路系统进行数值计算,计算出管路系统单位质量流体能量损失,在相同工况下,对管路系统进行外特性实验,测出实验中流体温度和管路系统单位质量流体能量损失,通过对比数值计算中设定的温度和实验测量温度,修正数值计算中设定的温度和动力粘度,再通过对比数值计算和实验测量的单位质量流体能量损失,调整数值计算湍流模型,最终选择计算能量损失与实验测量能量损失之间误差最小时的湍流模型为最优湍流模型;
(2)利用(1)中调整后的计算模型对管路系统进行变工况、变粗糙度数值计算;
(3)由数值计算得到两组过流截面,各自总压p,将两组过流截面总压引入伯努利能量方程,联合建立伯努利方程组,则由能量损失即可计算出不同工况、不同粗糙度管路系统的局部阻力系数ζ和沿程阻力系数λ;
(4)利用TableCurve 3D,将计算得到的局部和沿程阻力系数分别与雷诺数Re和相对粗糙度Δ/d进行三维曲面拟合,从而得到管路局部和沿程阻力系数的精确表示方法:其中f1和f2表示关系函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310245810.9A CN103353908B (zh) | 2013-06-20 | 2013-06-20 | 一种基于数值计算的管路阻力系数精确计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310245810.9A CN103353908B (zh) | 2013-06-20 | 2013-06-20 | 一种基于数值计算的管路阻力系数精确计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103353908A CN103353908A (zh) | 2013-10-16 |
CN103353908B true CN103353908B (zh) | 2016-08-31 |
Family
ID=49310280
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310245810.9A Active CN103353908B (zh) | 2013-06-20 | 2013-06-20 | 一种基于数值计算的管路阻力系数精确计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103353908B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104331603B (zh) * | 2014-10-13 | 2018-01-02 | 南京航空航天大学 | 一种任意组合弯管阻力系数的计算方法 |
CN104462777B (zh) * | 2014-11-13 | 2017-12-29 | 浙江工业大学 | 一种流体管路最小可行能量需求的测算方法 |
CN107808021B (zh) * | 2016-09-09 | 2021-06-15 | 国核示范电站有限责任公司 | 基于cfd的流体设备的阻力计算方法 |
CN106528994B (zh) * | 2016-10-27 | 2017-09-29 | 武汉大学 | 一种基于气液交界面耦合的调压室通气洞风速模拟方法 |
CN106910006A (zh) * | 2017-01-17 | 2017-06-30 | 西南石油大学 | 一种浅层稠油油藏两相流螺杆泵举升的预判方法 |
CN107122516A (zh) * | 2017-03-17 | 2017-09-01 | 浙江工业大学 | 一种灭火系统自段管道的沿程压力损失确定方法 |
CN107291984B (zh) * | 2017-05-22 | 2020-11-17 | 广州市市政工程设计研究总院有限公司 | 一种综合管廊通风系统当量摩擦阻力系数确定方法 |
CN108254122B (zh) * | 2018-01-19 | 2019-08-06 | 山东省环科院环境工程有限公司荣成分公司 | 一种地下自来水管网中管道阻力系数的测量方法 |
CN109215810A (zh) * | 2018-09-17 | 2019-01-15 | 中国核动力研究设计院 | 一种燃料组件定位格架局部阻力系数的精确实验获取方法 |
CN109933861B (zh) * | 2019-02-25 | 2023-05-26 | 中交上海航道局有限公司 | 一种拟合挖泥船排泥管线泥浆阻力特性公式的方法 |
CN110069891B (zh) * | 2019-05-07 | 2023-05-12 | 中建安装集团有限公司 | 一种基于bim的机电设计与复核方法 |
CN111982213B (zh) * | 2020-08-20 | 2022-02-11 | 中国核动力研究设计院 | 一种核反应堆模拟燃料组件流量测量方法及系统 |
CN113033112A (zh) * | 2020-12-30 | 2021-06-25 | 苏州水木科能科技有限公司 | 洁净室风系统建模方法及设备 |
CN113790867B (zh) * | 2021-07-30 | 2022-10-18 | 湖南大学 | 一种梗丝分选通道内流场特性检测方法 |
CN113792443B (zh) * | 2021-10-09 | 2024-05-07 | 中国核电工程有限公司 | 通风系统管网水力平衡的分析方法、管网设计和调节方法 |
-
2013
- 2013-06-20 CN CN201310245810.9A patent/CN103353908B/zh active Active
Non-Patent Citations (7)
Title |
---|
三种空化模型在离心泵空化流计算中的应用评价;刘厚林等;《农业工程学报》;20120815;全文 * |
基于CFD的管道局部阻力的数值模拟;赵月;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120115;全文 * |
弯管阻力损失的理论及实验研究;郭曹虹;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20120715;全文 * |
离心泵内流数值计算的湍流模型研究;舒敏骅;《万方数据》;20110824;全文 * |
离心泵内部空化特性的CFD模拟;王勇等;《排灌机械工程学报》;20110804;第29卷(第2期);全文 * |
离心泵多工况水力性能优化设计方法;王凯;《排灌机械工程学报》;20120515;第1卷(第30期);全文 * |
粗糙度对离心泵性能数值预测的影响;谈明高等;《中国农村水利水电》;20110504(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103353908A (zh) | 2013-10-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103353908B (zh) | 一种基于数值计算的管路阻力系数精确计算方法 | |
Reader-Harris | Orifice plates and venturi tubes | |
Wan et al. | Water hammer simulation of a series pipe system using the MacCormack time marching scheme | |
Sun et al. | Numerical investigation of the effect of surface roughness on the flow coefficient of an eccentric butterfly valve | |
CN102494742B (zh) | 一种大口径流量计在线校准的方法 | |
Streeter | Unsteady flow calculations by numerical methods | |
CN109977615B (zh) | 一种缩进式电磁流量传感器的测量管道结构设计方法 | |
CN104458107B (zh) | 一种简易的差压装置的检测方法 | |
CN101493861A (zh) | 一种自然循环热水锅炉水动力数值计算方法 | |
Merzkirch et al. | Fully developed turbulent pipe flow | |
CN106320257A (zh) | 基于水文观测的湖库槽蓄曲线确定方法 | |
Adamkowski et al. | Uncertainty analysis of liquid flow rate measurement with the pressure–time method | |
Li et al. | 1D-3D coupling investigation of hydraulic transient for power-supply failure of centrifugal pump-pipe system | |
Wang et al. | Turbulent flow in the entry region of a rough pipe | |
Adamkowski et al. | The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks | |
CN108280300A (zh) | 基于计算流体力学的大数据流量计开发方法 | |
Jing et al. | Numerical modeling of flow in continuous bends from Daliushu to Shapotou in Yellow River | |
Krzemianowski et al. | Use of non-uniform rational b-splines for discharge calculation in the velocity area method | |
Yang et al. | Optimize the Rectifier Structure to Improve the Accuracy of Gas Ultrasonic Flowmeter Under Low Flow Conditions | |
CN106500805A (zh) | 一种计量表稳流器整流效果的评估方法 | |
CN103196505B (zh) | 一种三相流体比例和流量测试装置及计算方法 | |
Lv et al. | Numerical simulation of 3D turbulent bend flow based on unstructured grids | |
Huber et al. | Hydraulic optimization of a T-junction hydraulic model tests and CFD-simulation | |
Meng et al. | The development of a multiphase flow meter without separation based on sloped open channel dynamics | |
JP2007107948A (ja) | 流量測定方法、流量測定装置およびこの流量測定装置を備えた水力機械 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |