CN113688526A - 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法 - Google Patents

一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法 Download PDF

Info

Publication number
CN113688526A
CN113688526A CN202110995829.XA CN202110995829A CN113688526A CN 113688526 A CN113688526 A CN 113688526A CN 202110995829 A CN202110995829 A CN 202110995829A CN 113688526 A CN113688526 A CN 113688526A
Authority
CN
China
Prior art keywords
water pump
prototype
model
pump
speed
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.)
Pending
Application number
CN202110995829.XA
Other languages
English (en)
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.)
Chongqing West Water Resources Development Co ltd
Hadong National Hydroelectric Power Equipment Engineering Technology Research Central Co ltd
Changjiang Institute of Survey Planning Design and Research Co Ltd
Harbin Electric Machinery Co Ltd
Original Assignee
Chongqing West Water Resources Development Co ltd
Hadong National Hydroelectric Power Equipment Engineering Technology Research Central Co ltd
Changjiang Institute of Survey Planning Design and Research Co Ltd
Harbin Electric Machinery Co Ltd
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 Chongqing West Water Resources Development Co ltd, Hadong National Hydroelectric Power Equipment Engineering Technology Research Central Co ltd, Changjiang Institute of Survey Planning Design and Research Co Ltd, Harbin Electric Machinery Co Ltd filed Critical Chongqing West Water Resources Development Co ltd
Priority to CN202110995829.XA priority Critical patent/CN113688526A/zh
Publication of CN113688526A publication Critical patent/CN113688526A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明公开一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法,在模型水泵研发阶段,就能准确的预估原型水泵的扬程。该方法通过水泵模型试验来预测原型水泵扬程时,除了考虑模型和原型间的相似定律外,还需考虑模型水泵和原型水泵由于雷诺数和粗糙度的差异引起的比尺效应。根据统计规律,将过流部件的速度因子和可换算水力比能损失指数表示成最优点比转速的函数,再将模型水泵试验工况的扬程换算到模型水泵参考条件下的扬程,最后将模型水泵参考条件先的扬程换算到原型水泵运行条件下。该方法首次在世界范围建立了雷诺数和粗糙度对原型水泵扬程的影响。不论是对原型水泵的运行还是新建泵站的研发,都有重要意义。

Description

一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算 方法
技术领域
本发明涉及水力机械领域,尤其涉及一种水泵扬程的换算方法。
背景技术
我国是一个严重干旱缺水严重的国家,再加上全国降水量的地 区分布很不均匀,造成了全国水土资源不平衡现象,长江流域和长 江以南耕地只占全国的36%,而水资源量却占全国的80%;黄、淮、 海三大流域,水资源量只占全国的8%,而耕地却占全国的40%,水 土资源相差十分悬殊。在水资源调配过程中,泵站起着及其重要左 右,尤其水泵的流量、扬程、空化性能和稳定运行范围,直接决定 调水工程的经济效益和社会效益。受自然条件限制、再加上不同调 水工程需求差异较大,大多数调水工程的泵站需要量身定制。对于水泵的研发,业内通行的做法的方法是通过在与原型水泵模拟的模 型水泵上进行试验,并根据试验结果来预测原型水泵性能。
对于水泵扬程的换算,目前行业内仅考虑了模型水泵和原型水 泵的相似定律,相似定律的假设条件是模型水泵和原型水泵均为理 想流动。模型水泵及其几何相似的原型水泵在实际运行过程中都会 产生各种水力损失。对于水泵扬程,可以分为局部损失和摩擦损失。 根据水力学原理,在相似工况,模型水泵及相似原型水泵的局部损 失比例相同;水泵的摩擦损失与过流部件的摩擦损失系数相关。根 据Nichtawitz公式,过流部件的摩擦损失系数与过流部件的雷诺数 和表面粗糙度的相关。故目前行业内水泵扬程的换算没有考虑由于 模型和原型间由于雷诺数和粗糙度的差异引起的比尺效应。对于新 建泵站,原型水泵的换算相对较低,而对于表面粗糙度很差的改造 泵站,导致泵的换算扬程偏高,严重的时候会限制水泵的运行范围。
发明内容
有鉴于此,本发明的目的是能够在泵站的模型水泵研发阶段,就 能根据模型和原型水泵的雷诺数和粗糙度准确的预估原型水泵的扬 程,降低原型水泵的运行风险。
本发明在进行水泵扬程换算时,充分考虑水泵各过流部件雷诺 数和粗糙度的影响,首先将模型试验数据换算到参考统一条件下, 得到水泵扬程换算的第一步比尺效应。之后将参考条件下的水泵扬 程换算到原型水泵运行条件下,得到第二步的比尺效应。最后结合 水泵扬程换算的相似定律和比尺效应,得到原型水泵的扬程。
具体的技术方案如下:
一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法 包括如下步骤:
步骤一:确定模型水泵的流量扬程曲线HM~QM和流量效率曲 线ηhM~QM
在水力机械模型试验台进行水泵特性试验,试验转速恒定为nM, 通过阀门依次降低水泵的过流流量,得到一系列工况下水泵性能, 需要测量的物理量有:扬程HM、流量QM、转速nM、力矩TM和水 温TwM,并计算每个工况下模型水泵的效率ηhM和雷洛数ReM,ηhM和ReM分别定义为:
Figure BDA0003234090290000031
Figure BDA0003234090290000032
其中水密度ρM和粘度υM通过水温TwM计算,gM通过试验所在地查取 获得,uM为转轮进口旋转速度,通过以上的试验数据,得到试验水 泵的HM~QM和ηhM~QM关系曲线;
步骤二:测量模型水泵过流部件的粗糙度:模型水泵的特征过 流部件有水泵进水管、泵轮、导叶和蜗型扩散管,分别用粗糙度仪 对这些表面的平均粗糙度进行测量,测量结果分别表示为RaDT,M、 RaRU,M、RaSV,M、和RaSP,M
步骤三:确定模型水泵最优效率和比转速:模型水泵最优效率 定义为在模型试验转速nM时各流量下效率ηhM的最高值ηhMopt,模 型水泵最优点的比转速定义为:
Figure BDA0003234090290000041
式中,HM,opt和QM,opt分别为最优点所对应的扬程和流量;
步骤四:确定模型试验各过流部件的速度因子:模型水泵按最 优工况进行设计,故过流部件尺寸与水泵最优点比转速相关,表征 过流部件速度的无量纲特征速度因子为κu,为过流部件绝对速度vm与转轮旋转速度um的比值,根据水泵最优点的比转速,每个过流部 件速度因子κu与比转速的关系如下所示:
进水管:κuSP=0.27;
泵轮:κuRU=-1.30×NQE+0.79;
导叶:κuSV=-1.40×NQE+0.53;
蜗壳:κuSP=-0.50×NQE+0.31;
其中κuSP、κuRU、κuSV、κuSP分别表示进水管、转轮、导叶和蜗型扩 散管的无量纲特征速度因子;
步骤五:确定每个过流部件的可换算水力比能损失:水泵每个 部件可换算损失指数dEref与水泵的给个过流部件的流速和尺寸相关, 根据水泵最优点的比转速,每个过流部件的可换算损失指数dEref与比 转速的关系如下所示:
进水管:dE,SPref=0.5×NQE+0.05;
泵轮:dE,RUref=3.4×NQE+1.55;
导叶:dE,SVref=-NQE+0.50;
蜗壳:dE,SPref=0.45×NQE
其中dE,SPref、dE,RUref、dE,SVref、dE,SPref分别表示进水管、转轮、导叶 和蜗型扩散管的可换算损失指数;
步骤六:计算模型水泵参考条件下的扬程的比尺效应:模型水 泵参考雷洛数Reref=7×106,过流部件参考雷洛数分别为: RaDT,Mref=0.8μm、RaRU,Mref=0.4μm、RaSV,Mref=0.8μm和 RaSP,Mref=0.8μm,根据经典流体力学的Nichtawitz损失系数λ公式
Figure BDA0003234090290000051
其中λ0为常数0.0085,ks为砂砾粗糙度,dh为流道的水力直径,Reref为常数7×106,Red为流道的雷诺数,流部件由于雷诺数和粗糙度的 差异引起的水泵扬程的比尺效应:
Figure BDA0003234090290000052
其中ΔECO为部件比能的比尺效应,dE,COref为部件可换算损失指数,κuCO为部件无量纲特征速度因子,RaCO,M为模型水泵部件粗糙度,RaCO,ref为水泵部件参考粗糙度,ReM为模型水泵雷诺数,DM为模型水泵直 径;分别将CO等于DT、RU、SV和SP代入上式,分别得到进水 管的参考条件下的比尺效应
Figure BDA0003234090290000053
泵轮的
Figure BDA0003234090290000054
导叶
Figure BDA0003234090290000055
的和蜗壳的
Figure BDA0003234090290000056
将以上四项相加,得到模型水泵在 参考条件下的比尺效应
Figure BDA0003234090290000061
Figure BDA0003234090290000062
Figure BDA0003234090290000063
为从试验条件下向参考条件下的换算;
步骤七:模型水泵参考条件到原型水泵运行条件下扬程的比尺 效应:涉及原型水泵比尺效应的运行条件包括原型机的雷洛数和原 型机过流部件的粗糙度,原型水泵的雷诺数由泵站运行时的水温、 水泵转速和水泵直径决定,其计算公式为:
Figure BDA0003234090290000064
其中DP为原型泵转轮直径,为给定值;nP为原型泵转轮直径,为给 定值;uP为原型泵旋转速度,由转速nP决定;υP为泵站水的粘度, 由水温TwP决定;
原型水泵的过流部件的粗糙度包括进水管、泵轮、导叶和蜗壳, 分别表示为RaDT,P、RaRU,P、RaSV,P、和RaSP,P,在模型研发阶段, 原型机粗糙度由制造厂给出,在原型泵加工过程中进行检查,同理 根据以下公式计算由于原型水泵雷洛数和粗糙度的差异引起过流部 件压力损失的比尺效应:
Figure BDA0003234090290000065
其中RaCO,P为原型水泵部件粗糙度,ReP为原型水泵雷诺数,DP为原 型水泵直径;CO为别为DT、RU、SV和SP,可以计算原型水泵每 个过流部件的比尺效应
Figure BDA0003234090290000071
Figure BDA0003234090290000072
将以上四项相加,得到原型水泵的比尺效应
Figure BDA0003234090290000073
Figure BDA0003234090290000079
Figure BDA0003234090290000074
为从参考条件下向原型水泵运行条件下的换算;
步骤八:计算原型水泵扬程:原型机水泵的扬程由水泵运行的 相似定律和模型到原型机的比尺效应决定,原型水泵的扬程HP为:
Figure BDA0003234090290000075
其中gM为模型试验所在地重力加速度,gP为原型泵站所在地重力加 速度,nP为原型水泵转速,nM为模型水泵转速,DP为原型水泵直径, DM为模型水泵直径,
Figure BDA0003234090290000076
为试验条件下向参考条件下的比能换 算,
Figure BDA0003234090290000077
为参考条件下向原型水泵运行条件下的比能换算。
在上述一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换 算方法中,所述步骤四中通过理论分析和统计分析,得到水泵每个 过流部件速度因子κu与比转速NQE的关系:
速度因子κu为过流部件绝对速度vm与转轮旋转速度um的比值, 即:
Figure BDA0003234090290000078
水泵比转速的定义为:
Figure BDA0003234090290000081
其中,nM时为模型试验转速,gM为试验所在地重力加速度, HM,opt和QM,opt分别为最优点所对应的扬程和流量,通过试验获得;
通过设计数据统计,水泵不同过流部件速度因子与比转速的关 系为:
进水管:κuSP=0.27;
泵轮:κuRU=-1.30×NQE+0.79;
导叶:κuSv=-1.40×NQE+0.53;
蜗壳:κuSP=-0.50×NQE+0.31。
在上述一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的 换算方法中,所述步骤五中,通过统计规律得到水泵每个过流部件 可换算损失指数dEref和速度因子κu与比转速的关系:
进水管:dE,SPref=0.5×NQE+0.05;
泵轮:dE,RUref=3.4×NQE+1.55;
导叶:dE,SVref=-NQE+0.50;
蜗壳:dE,SPref=0.45×NQE
与现有技术相比,本发明的有益效果是:
1.简单易行,根据统计规律将水泵的速度因子和可换算损失指 数表示为比转速的函数;
2.适用于各种比转速的立式离心泵;
3.能在模型水泵的研发阶段,就能准确预估各种运行条件原型 水泵的扬程,保证了泵站的运行范围。
附图说明
图1为模型水泵扬程流量曲线。
图2为模型水泵效率流量曲线。
图3为模型试验条件到模型参考条件的比尺效应。
图4为原型水泵流量扬程曲线。
具体实施方式
下面结合附图对本申请作进一步描述。以下实施例仅用于更加 清楚地说明本发明的技术方案,而不能以此来限制本申请的保护范 围。
一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法 包括如下步骤:
1)通过模型试验确定水泵扬程流量曲线HM~QM和效率流量曲 线ηhM~QM
水泵模型试验依据GB/T15613和IEC60193进行,模型试验转 速恒定为nM,通过阀门依次降低水泵的过流流量,得到一系列工况 下水泵性能,需要测量或者记录的物理量有:扬程HM、流量QM、 转速nM、力矩TM和水温TwM。并计算每个工况下模型水泵的效率ηhM和雷洛数ReM,ηhM和ReM分别定义为:
Figure BDA0003234090290000101
Figure BDA0003234090290000102
其中水密度ρM和粘度vM通过水温TwM计算,gM通过试验所在地查取 获得,uM为转轮进口旋转速度。通过以上的试验数据,将每个工况 下的模型效率ηhM与模型试验流量QM通过样条曲线拟合可以得到该 水泵的效率流量曲线ηhM~QM,如图1所示;将每个工况下的模型试 验扬程HM与模型试验流量QM通过样条曲线拟合可以得到该水泵的 HM流量曲线HM~QM,如图2所示。
2)测量模型水泵过流部件的粗糙度
使用数字接触式粗糙度仪分别测量模型水泵进水管、泵轮、导 叶和蜗型扩散管的粗糙度,测量结果分别表示为RaDT,M、RaRU,M、 RaSV,M、和RaSP,M
3)计算模型水泵最优点的比转速
模型水泵最优效率定义为在模型试验转速nM时各流量下效率 ηhM的最高值ηhMopt,如图1中的点A所示,模型水泵最优点的比转 速定义为:
Figure BDA0003234090290000111
式中,HM,opt和QM,opt分别为最优点所对应的扬程和流量。
4)根据统计规律公式分别计算模型试验各过流部件的速度因子κu和每个过流部件的可换算水力比能损失指数dEref
根据水泵最优点的比转速,每个过流部件的可换算损失指数dEref和速度因子κu通过以下公式计算:
进水管:κuSP=0.27,dE,SPref=0.5×NQE+0.05;
泵轮:κuRU=-1.30×NQE+0.79,dE,RUref=3.4×NQE+1.55;
导叶:κuSV=-1.40×NQE+0.53,dE,SVref=-NQE+0.50;
蜗壳:κusP=-0.50×NQE+0.31,dE,SPref=0.45×NQE
5)计算模型水泵参考条件下的扬程的比尺效应
Figure BDA0003234090290000113
模型水泵参考条件定义为:Reref=7×106、RaDT,Mref=0.8μm、 RaRU,Mref=0.4μm、RaSV,Mref=0.8μm和RaSP,Mref=0.8μm。
根据经典流体力学的Nichtawitz损失系数公式
Figure BDA0003234090290000112
可以得出每个过流部件由于雷诺数和粗糙度的差异引起的水泵扬程 的比尺效应公式:
Figure BDA0003234090290000121
分别将CO等于DT、RU、SV和SP代入上式,可以分别得到 参考条件下进水管的比尺效应
Figure BDA0003234090290000122
泵轮的
Figure BDA0003234090290000123
导叶
Figure BDA0003234090290000124
的和蜗壳的
Figure BDA0003234090290000125
将以上四项相加,可以得到模型水 泵在从试验条件到参考条件下的比尺效应
Figure BDA0003234090290000126
Figure BDA0003234090290000127
Figure BDA0003234090290000128
为从试验条件下向参考条件下的换算,
Figure BDA0003234090290000129
的计 算结果如图3所示。
6)计算模型水泵参考条件到原型水泵运行条件下扬程的比尺效应
Figure BDA00032340902900001210
原型水泵运行条件包括原型机的雷洛数和原型水泵过流部件的 粗糙度,原型水泵的雷诺数由泵站运行时的水温、水泵转速和水泵 直径决定,其计算公式为:
Figure BDA00032340902900001211
υP为泵站水的粘度,由水温TwP决定。
原型水泵的过流部件的粗糙度包括进水管、泵轮、导叶和蜗壳, 分别表示为RaDT,P、RaRU,P、RaSV,P、和RaSP,P。在模型研发阶段, 原型机粗糙度一般由制造厂给出,在原型泵加工过程中进行检查。 同理可以根据以下公式计算由于原型水泵雷洛数和粗糙度的差异引 起过流部件压力损失的比尺效应:
Figure BDA0003234090290000131
CO为别为DT、RU、SV和SP,可以计算原型水泵每个过流部 件的比尺效应
Figure BDA0003234090290000132
Figure BDA0003234090290000133
将以 上四项相加,可以得到原型水泵的比尺效应
Figure BDA0003234090290000134
Figure BDA0003234090290000135
Figure BDA0003234090290000136
为从参考条件下向原型水泵运行条件下的换算。由于 每个工况,原型水泵运行条件和模型参考条件为已知常数,故每 个工况点
Figure BDA0003234090290000137
的值为常数。
7)计算原型水泵扬程HP的计算
如果按照目前行业所采取的方法,原型水泵的扬程和流量仅考 虑模型和扬程的相似定律,计算公式如下:
Figure BDA0003234090290000138
Figure BDA0003234090290000139
其计算结果见图4中的虚线B。
在进行原型机水泵的扬程计算时需同时考虑模型和原型间的相 似定律和比尺效应,其计算公式如下所示:
Figure BDA0003234090290000141
根据模型试验结果和上述公式,可以计算得到原型水泵扬程与 原型水泵流量之间的关系,如图4中实线C所示。
本发明只是对本发明的示例性说明,并不限定它的保护范围, 本领域技术人员还可以对其局部进行改变,只要没有超出本发明的 精神实质,都在本发明的保护范围内。

Claims (3)

1.一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法,其特征包括如下步骤:
步骤一:确定模型水泵的流量扬程曲线HM~QM和流量效率曲线ηhM~QM
在水力机械模型试验台进行水泵特性试验,试验转速恒定为nM,通过阀门依次降低水泵的过流流量,得到一系列工况下水泵性能,需要测量的物理量有:扬程HM、流量QM、转速nM、力矩TM和水温TwM,并计算每个工况下模型水泵的效率ηhM和雷洛数ReM,ηhM和ReM分别定义为:
Figure FDA0003234090280000011
Figure FDA0003234090280000012
其中水密度ρM和粘度υM通过水温TwM计算,gM通过试验所在地查取获得,uM为转轮进口旋转速度,通过以上的试验数据,得到试验水泵的HM~QM和ηhM~QM关系曲线;
步骤二:测量模型水泵过流部件的粗糙度:模型水泵的特征过流部件有水泵进水管、泵轮、导叶和蜗型扩散管,分别用粗糙度仪对这些表面的平均粗糙度进行测量,测量结果分别表示为RaDT,M、RaRU,M、RaSV,M、和RaSP,M
步骤三:确定模型水泵最优效率和比转速:模型水泵最优效率定义为在模型试验转速nM时各流量下效率ηhM的最高值ηhMopt,模型水泵最优点的比转速定义为:
Figure FDA0003234090280000021
式中,HM,opt和QM,opt分别为最优点所对应的扬程和流量;
步骤四:确定模型试验各过流部件的速度因子:模型水泵按最优工况进行设计,故过流部件尺寸与水泵最优点比转速相关,表征过流部件速度的无量纲特征速度因子为κu,为过流部件绝对速度vm与转轮旋转速度um的比值,根据水泵最优点的比转速,每个过流部件速度因子κu与比转速的关系如下所示:
进水管:κuSP=0.27;
泵轮:κuRU=-1.30×NQE+0.79;
导叶:κuSV=-1.40×NQE+0.53;
蜗壳:κuSP=-0.50×NQE+0.31;
其中κuSP、κuRU、κuSV、κuSP分别表示进水管、转轮、导叶和蜗型扩散管的无量纲特征速度因子;
步骤五:确定每个过流部件的可换算水力比能损失:水泵每个部件可换算损失指数dEref与水泵的给个过流部件的流速和尺寸相关,根据水泵最优点的比转速,每个过流部件的可换算损失指数dEref与比转速的关系如下所示:
进水管:dE,SPref=0.5×NQE+0.05;
泵轮:dE,RUref=3.4×NQE+1.55;
导叶:dE,SVref=-NQE+0.50;
蜗壳:dE,SPref=0.45×NQE
其中dE,SPref、dE,RUref、dE,SVref、dE,SPref分别表示进水管、转轮、导叶和蜗型扩散管的可换算损失指数;
步骤六:计算模型水泵参考条件下的扬程的比尺效应:模型水泵参考雷洛数Reref=7×106,过流部件参考雷洛数分别为:RaDT,Mref=0.8μm、RaRU,Mref=0.4μm、RaSV,Mref=0.8μm和RaSP,Mref=0.8μm,根据经典流体力学的Nichtawitz损失系数λ公式
Figure FDA0003234090280000031
其中λ0为常数0.0085,ks为砂砾粗糙度,dh为流道的水力直径,Reref为常数7×106,Red为流道的雷诺数,流部件由于雷诺数和粗糙度的差异引起的水泵扬程的比尺效应:
Figure FDA0003234090280000032
其中ΔECO为部件比能的比尺效应,dE,COref为部件可换算损失指数,κuCO为部件无量纲特征速度因子,RaCO,M为模型水泵部件粗糙度,RaCO,ref为水泵部件参考粗糙度,ReM为模型水泵雷诺数,DM为模型水泵直径;分别将CO等于DT、RU、SV和SP代入上式,分别得到进水管的参考条件下的比尺效应
Figure FDA0003234090280000041
泵轮的
Figure FDA0003234090280000042
导叶
Figure FDA0003234090280000043
的和蜗壳的
Figure FDA0003234090280000044
将以上四项相加,得到模型水泵在参考条件下的比尺效应
Figure FDA0003234090280000045
Figure FDA0003234090280000046
Figure FDA0003234090280000047
为从试验条件下向参考条件下的换算;
步骤七:模型水泵参考条件到原型水泵运行条件下扬程的比尺效应:涉及原型水泵比尺效应的运行条件包括原型机的雷洛数和原型机过流部件的粗糙度,原型水泵的雷诺数由泵站运行时的水温、水泵转速和水泵直径决定,其计算公式为:
Figure FDA0003234090280000048
其中DP为原型泵转轮直径,为给定值;nP为原型泵转轮直径,为给定值;uP为原型泵旋转速度,由转速nP决定;υP为泵站水的粘度,由水温TwP决定;
原型水泵的过流部件的粗糙度包括进水管、泵轮、导叶和蜗壳,分别表示为RaDT,P、RaRU,P、RaSV,P、和RaSP,P,在模型研发阶段,原型机粗糙度由制造厂给出,在原型泵加工过程中进行检查,同理根据以下公式计算由于原型水泵雷洛数和粗糙度的差异引起过流部件压力损失的比尺效应:
Figure FDA0003234090280000051
其中RaCO,P为原型水泵部件粗糙度,ReP为原型水泵雷诺数,DP为原型水泵直径;CO为别为DT、RU、SV和SP,可以计算原型水泵每个过流部件的比尺效应
Figure FDA0003234090280000052
Figure FDA0003234090280000053
将以上四项相加,得到原型水泵的比尺效应
Figure FDA0003234090280000054
Figure FDA0003234090280000055
Figure FDA0003234090280000056
为从参考条件下向原型水泵运行条件下的换算;
步骤八:计算原型水泵扬程:原型机水泵的扬程由水泵运行的相似定律和模型到原型机的比尺效应决定,原型水泵的扬程HP为:
Figure FDA0003234090280000057
其中gM为模型试验所在地重力加速度,gP为原型泵站所在地重力加速度,nP为原型水泵转速,nM为模型水泵转速,DP为原型水泵直径,DM为模型水泵直径,
Figure FDA0003234090280000058
为试验条件下向参考条件下的比能换算,
Figure FDA0003234090280000059
为参考条件下向原型水泵运行条件下的比能换算。
2.根据权利要求1所述的一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法,其特征在于:所述步骤四中通过理论分析和统计分析,得到水泵每个过流部件速度因子κu与比转速NQE的关系:
速度因子κu为过流部件绝对速度vm与转轮旋转速度um的比值,即:
Figure FDA0003234090280000061
水泵比转速的定义为:
Figure FDA0003234090280000062
其中,nM时为模型试验转速,gM为试验所在地重力加速度,HM,opt和QM,opt分别为最优点所对应的扬程和流量,通过试验获得;
通过设计数据统计,水泵不同过流部件速度因子与比转速的关系为:
进水管:κuSP=0.27;
泵轮:κuRU=-1.30×NQE+0.79;
导叶:κuSV=-1.40×NQE+0.53;
蜗壳:κuSP=-0.50×NQE+0.31。
3.根据权利要求1所述的一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法,其特征在于:所述步骤五中,通过统计规律得到水泵每个过流部件可换算损失指数dEref和速度因子κu与比转速的关系:
进水管:dE,SPref=0.5×NQE+0.05;
泵轮:dE,RUref=3.4×NQE+1.55;
导叶:dE,SVref=-NQE+0.50;
蜗壳:dE,SPref=0.45×NQE
CN202110995829.XA 2021-08-27 2021-08-27 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法 Pending CN113688526A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110995829.XA CN113688526A (zh) 2021-08-27 2021-08-27 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110995829.XA CN113688526A (zh) 2021-08-27 2021-08-27 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法

Publications (1)

Publication Number Publication Date
CN113688526A true CN113688526A (zh) 2021-11-23

Family

ID=78583386

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110995829.XA Pending CN113688526A (zh) 2021-08-27 2021-08-27 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法

Country Status (1)

Country Link
CN (1) CN113688526A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115525996A (zh) * 2022-09-21 2022-12-27 中国航发沈阳发动机研究所 一种涡轮工作叶片流量特性旋转修正方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102024069A (zh) * 2009-09-21 2011-04-20 中国轻工业长沙工程有限公司 用于确定泵扬程的计算机辅助设计系统和方法
CN105909537A (zh) * 2016-05-13 2016-08-31 南京工业大学 一种高温熔盐泵的常温模化测试方法
JP6436408B1 (ja) * 2018-02-15 2018-12-12 有限会社北沢技術事務所 ポンプ流量計測装置
CN109271699A (zh) * 2018-09-10 2019-01-25 江苏大学镇江流体工程装备技术研究院 一种大型泵装置性能曲线的计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102024069A (zh) * 2009-09-21 2011-04-20 中国轻工业长沙工程有限公司 用于确定泵扬程的计算机辅助设计系统和方法
CN105909537A (zh) * 2016-05-13 2016-08-31 南京工业大学 一种高温熔盐泵的常温模化测试方法
JP6436408B1 (ja) * 2018-02-15 2018-12-12 有限会社北沢技術事務所 ポンプ流量計測装置
CN109271699A (zh) * 2018-09-10 2019-01-25 江苏大学镇江流体工程装备技术研究院 一种大型泵装置性能曲线的计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
周君亮;: "原型及模型泵水力装置参数换算", 排灌机械, no. 05, pages 273 - 280 *
李彦军;严登丰;袁寿其;: "大型低扬程泵与泵装置特性预测", 农业机械学报, no. 10, pages 36 - 40 *
李龙;杨雪林;李丹;: "考虑粗糙度影响的水泵原模型效率换算", 河海大学学报(自然科学版), no. 03, pages 327 - 331 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115525996A (zh) * 2022-09-21 2022-12-27 中国航发沈阳发动机研究所 一种涡轮工作叶片流量特性旋转修正方法及系统
CN115525996B (zh) * 2022-09-21 2024-03-19 中国航发沈阳发动机研究所 一种涡轮工作叶片流量特性旋转修正方法及系统

Similar Documents

Publication Publication Date Title
Wang et al. The method for determining blade inlet angle of special impeller using in turbine mode of centrifugal pump as turbine
CN103353908B (zh) 一种基于数值计算的管路阻力系数精确计算方法
WO2019136921A1 (zh) 轴流式pat发电模式最优效率点参数和性能曲线预估方法
Quaranta et al. Optimization of breastshot water wheels performance using different inflow configurations
Chakraborty et al. Numerical Studies on Effects of Blade Number Variationson Performance of Centrifugal Pumps at 4000 RPM
CN104454564A (zh) 基于试验的轴流泵装置导叶体水力优化方法
Adhikari et al. A study on developing pico propeller turbine for low head micro hydropower plants in Nepal
CN100348863C (zh) 利用管路特性曲线测量泵特性曲线的方法
CN113688526A (zh) 一种考虑水泵过流部件雷诺数和粗糙度的水泵扬程的换算方法
Gautam et al. Numerical investigation of the effects of leakage flow from guide vanes of Francis turbines using alternative clearance gap method
Li et al. Stability optimization and analysis of a bidirectional shaft extension pump
CN204025064U (zh) 带前导叶测流的新型高效立式轴流泵装置
Thapa et al. Effects of turbine design parameters on sediment erosion of Francis runner
Xue et al. Research and optimization of performances of a pump turbine in pump mode
Miao et al. Study on energy conversion characteristics in volute of pump as turbine
Ge et al. Sediment erosion on pelton turbines: a review
CN116107236A (zh) 水电站过渡过程的半物理仿真模型试验方法
Ida Analysis of scale effects on performance characteristics of hydraulic turbines: Part 1: Scale formulae of hydraulic performance and loss distribution coefficients in model Francis turbines and pump-turbines
CN104791172A (zh) 双向分流的混流式水泵水轮机叶片空间位置的确定方法
Chen et al. Effect of rotational speed on performance of mixed flow pump as turbine
Li et al. Load rejection transient with joint closing law of ball-valve and guide vane for two units in pumped storage power station
CN108915938B (zh) 一种减轻抽水蓄能电站机组和厂房振动的电站运行方法
ENOMOTO et al. Numerical simulation of turbulence flow in a Bulb turbine
Benigni et al. Simulation of damages due to cavitation in non-rotating components in a Kaplan turbine
Mohammadi et al. Analyzing Mathematical and Software Methods for Selecting and Designing Francis Turbine in Hydropower Plants

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination