CN102865233B - 一种基于全工况内流测量的离心泵滑移系数确定方法 - Google Patents
一种基于全工况内流测量的离心泵滑移系数确定方法 Download PDFInfo
- Publication number
- CN102865233B CN102865233B CN201210382442.8A CN201210382442A CN102865233B CN 102865233 B CN102865233 B CN 102865233B CN 201210382442 A CN201210382442 A CN 201210382442A CN 102865233 B CN102865233 B CN 102865233B
- Authority
- CN
- China
- Prior art keywords
- centrifugal pump
- under
- slip
- operating mode
- different operating
- 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
- Control Of Non-Positive-Displacement Pumps (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种基于全工况内流测量的离心泵滑移系数确定方法。本发明采用PIV实验测量离心泵不同工况下的内部非定常流动规律；基于PIV实验得到的绝对速度分布，采用算术平均法得到离心泵不同工况下叶轮出口处的平均滑移速度，并计算出离心泵不同工况下滑移系数；计算不同工况下的离心泵液流角系数；对不同比转数的离心泵内部流动进行PIV测量，得到不同比转数、不同工况下的液流角系数；采用一元线性回归方法建立不同工况下的液流角系数与比转数之间的关系，从而确定不同工况下的离心泵滑移系数计算公式。本发明不仅能较为准确计算不同工况下离心泵滑移系数，还可以在此基础上建立起离心泵能量性能计算公式并对已有的泵进行节能改造。
Description
技术领域
本发明属于离心泵能量性能计算领域，具体涉及一种基于全工况内流测量的离心泵滑移系数确定方法。
背景技术
目前，普遍采用的滑移系数计算公式主要有斯托道拉公式、威斯奈尔公式、普夫莱德尔公式。斯托道拉公式和威斯奈尔公式仅考虑轴向漩涡对出口相对运动的影响，普夫莱德尔公式则假设叶轮中液体的流动不产生偏离。与叶轮内流体的真实流动相比，这三个公式都没有考虑液流角的影响，计算出来的滑移系数存在一定的误差。随着社会经济的快速发展，很多离心泵在设计时必须要满足多个工况的性能要求，如核电用离心泵和舰艇用离心泵，而滑移系数是基于损失法的离心泵多工况水力设计中的关键技术之一，因此有必要在离心泵真实流动的基础上建立离心泵全工况滑移系数确定方法。
迄今为止，尚未见离心泵全工况滑移系数计算方法的文献报道，仅有一些学者做与离心泵全工况滑移系数计算方法相近的研究工作。山东大学硕士论文《滑移系数对离心泵外特性的影响分析》和山东大学学报（工学版）《考虑边界层堵塞的离心泵滑移系数修正》(2009年第2期)考虑了由粘性效应引起的边界层堵塞的影响，并通过计算得到的边界层位移厚度对叶轮流动滑移量进行了修正。流体机械《离心油泵滑移系数的计算》（2008年第5期）考虑了输送介质的粘度影响，对已有滑移系数计算公式进行了改进，发展了适用于离心油泵的滑移系数计算公式。排灌机械《离心泵滑移系数精度的比较》（2006年第6期）在大量优秀离心泵水力模型的基础上比较了目前常用的斯托道拉、威斯奈和斯基克钦等滑移系数计算公式的精度。水泵技术《离心泵叶轮滑移系数的研究》（2006年第1期）分析了叶片排挤系数对滑移系数的影响并修正了斯托道拉滑移系数计算公式。流体工程《超短叶片离心泵叶轮的滑移系数》（1991年第6期）采用三孔圆柱探针测量了超短叶片离心泵叶轮出口处液体的绝对速度，并在此基础上研究了叶片出口安放角对长短叶片滑移系数的影响，但是探针则影响了叶轮出口处的真实流动。
发明内容
本发明旨在提供一种基于全工况内流测量的离心泵滑移系数确定方法，通过采用PIV内流测试的方法来得到离心泵叶轮出口处液体的真实流动状态，并在此基础上确定离心泵全工况滑移系数计算公式。
为达到以上目的，采用如下技术方案：
通过离心泵内部非定常流动的PIV测量，揭示不同工况下的离心泵叶轮出口处速度分布规律，采用算术平均法得到离心泵不同工况下叶轮出口处的平均滑移速度，并在此基础上，将液流角系数k α 引入斯托道拉滑移系数计算中，从而建立一种离心泵全工况滑移系数确定公式。
其具体步骤如下：
（1）采用PIV实验测量离心泵不同工况下的内部非定常流动规律；
采用相位平均法测量离心泵不同工况、不同相位下的内部非定常流动；采用二维PIV测量的具有时间序列的绝对速度场为：f（X，Y，t i ）、f（X，Y，t i +T）、f（X，Y，t i +2T）、……、f（X，Y，t i +nT），其中X、Y为二维流场的空间坐标，t i 为不同相位的初始时刻，i=0、1、……，T为叶轮旋转周期，n为拍摄流场图像的数量；将n组绝对速度场中对应坐标的瞬时流速进行算术平均即可得到不同工况、不同相位下离心泵内部流动的绝对速度分布；
（2）基于PIV实验得到的绝对速度分布和叶轮出口处的速度三角形，计算出不同工况、不同相位下的离心泵叶轮出口处的圆周分速度v u2，实验；
（3）计算不同工况下的离心泵叶轮出口处的圆周分速度v u2，理论，然后计算出不同工况、不同相位下的离心泵叶轮出口处的滑移速度，滑移速度为Δv u2=v u2，理论-v u2，实验；采用算术平均法得到离心泵不同工况下叶轮出口处的平均滑移速度 ，同时计算出离心泵不同工况下滑移系数，其中u 2为叶轮出口处的圆周速度；
（4）计算不同工况下的离心泵液流角系数k α ，液流角系数为，其中z为叶轮叶片数，β 2为叶片出口安放角；
（5）根据（1）、（2）、（3）和（4），对不同比转数n s的离心泵内部流动进行PIV测量，得到不同比转数、不同工况下的液流角系数；
（6）采用一元线性回归方法对不同比转数、不同工况下的液流角系数进行回归，建立不同工况下的液流角系数k α 与比转数n s之间的关系k α =f(n s/100)，从而确定了全工况下的离心泵滑移系数计算公式。
本发明的优点在于：
（1）基于离心泵内真实流动建立的离心泵全工况滑移系数计算公式，不仅考虑叶轮内轴向漩涡对出口相对运动的影响，还考虑了叶轮出口处液流角变化的影响；
（2）建立的离心泵全工况滑移系数计算公式，适用于所有流量范围内的离心泵滑移系数的计算；
（3）不仅能较为准确地计算离心泵不同工况下的滑移系数，还可以在此基础上建立起离心泵能量性能计算公式并对已有的泵进行多工况优化或节能改造。
附图说明
图1为一种基于全工况内流测量的离心泵滑移系数确定方法的流程图。
图2为本发明实施例的双叶片离心泵PIV测量平面。
图3为本发明实施例的0.8倍设计工况、φ/6相位下双叶片离心泵叶轮出口处的绝对速度分布示意图。
具体实施方式
实施例：
（1）采用PIV实验测量离心泵三个工况下的内部非定常流动规律。
一比转数为115.6的双叶片离心泵，其3个工况下的性能参数为：0.8倍设计工况Q 1=27.56m3/h，H 1=4.96m，η 1=56.32%；设计工况Q 2=34.48m3/h，H 2=4.50m，η 2=58.48%；1.2倍设计工况Q 3=41.40m3/h，H 3=4.03m，η 3=57.89%。
1） 测量区域及示踪粒子的选取。
选择叶轮中间截面为测量平面，如图2所示。叶轮顺时针转动，转速为1000r/min，正对隔舌处叶片吸力面与X轴之间夹角θ=28°。每个工况各测量6个不同相位，分别记为φ/6、φ/3、φ/2、2φ/3、5φ/6、φ，相位差为30°。
选择能随流体同步流动、且能较好地散射激光的空心玻璃球作为示踪粒子。
2） 离心泵内速度场的测量方法。
采用相位平均法测量离心泵3个工况、6个相位下的内部流动。采用二维PIV测量的具有时间序列的绝对速度场为：f（X，Y，t i ）、f（X，Y，t i +T）、f（X，Y，t i +2T）、……、f（X，Y，t i +nT），其中X、Y为二维流场的空间坐标，t i 为不同相位的初始时刻，i=0、1、……、5，T为叶轮旋转周期，n为拍摄流场图像的数量。
将n组绝对速度场中对应坐标的瞬时流速进行算术平均即可得到不同工况、不同相位下离心泵内部流动的绝对速度分布。小流量工况、φ/6相位下叶轮出口处的绝对速度分布如图3所示。
（2）基于PIV实验得到的绝对速度分布和叶轮出口处的速度三角形，计算出3个工况、6个相位下的离心泵叶轮出口处的相对液流角和圆周分速度v u2，实验。
（3）计算3个工况下的离心泵叶轮出口处的圆周分速度v u2，理论，然后计算出3个工况、6个相位下的离心泵叶轮出口处的滑移速度Δv u2:
（1）
式中：u 2为叶轮出口处的圆周速度；v m2为叶轮出口处的轴面速度。
（2）
（3）
（4）
（5）
式中：Q为泵的实际流量；η v为离心泵的容积效率；D 2为叶轮出口直径；b 2为叶轮出口宽度；ψ 2为叶片出口排挤系数；q 1为叶轮密封环的泄漏量；μ为密封环间隙的速度系数，；η为圆角系数，一般取0.5-0.9；λ为水力阻力系数，一般取0.04-0.06；l为间隙长度；f为密封环间隙的过流断面面积，f=2πR m b；R m为密封环半径；b为间隙宽度；ΔH m为间隙两端的压差，n s≤100时，ΔH m=0.6H，n s≥100时，ΔH m=0.7H；z为叶片数；s u2为叶轮出口处叶片的圆周厚度。
采用算术平均法对6个相位下的离心泵叶轮出口处的滑移速度进行平均，得到了离心泵3个工况下叶轮出口处的平均滑移速度，同时计算出离心泵3个工况下滑移系数。
1） 0.8倍设计工况
当相位为φ/6、φ/3、φ/2、2φ/3、5φ/6、φ时，叶轮出口处的滑移速度分别为6.358m/s、6.444m/s、6.391m/s、6.106m/s、5.890m/s、5.795m/s。叶轮出口处的平均滑移速度为6.164m/s。该工况下的滑移系数为0.415。
2） 设计工况
当相位为φ/6、φ/3、φ/2、2φ/3、5φ/6、φ时，叶轮出口处的滑移速度分别为6.182m/s、6.230m/s、6.222m/s、5.978m/s、5.758m/s、5.693m/s。叶轮出口处的平均滑移速度为6.010m/s。该工况下的滑移系数为0.429。
3） 1.2倍设计工况
当相位为φ/6、φ/3、φ/2、2φ/3、5φ/6、φ时，叶轮出口处的滑移速度分别为6.079m/s、6.117m/s、6.081m/s、5.892m/s、5.708m/s、5.611m/s。叶轮出口处的平均滑移速度为5.915m/s。该工况下的滑移系数为0.439。
（4）计算3个工况下的离心泵液流角系数k α ，液流角系数为。
1） 0.8倍设计工况下，k α =0.923。2） 设计工况下，k α =0.900。3） 1.2倍设计工况下，k α =0.886。
（5）根据（1）、（2）、（3）和（4），对其它不同比转数n s的离心泵内部流动进行PIV测量，得到不同比转数、3个相同工况下的液流角系数k α 。
1） 0.8倍设计工况。n s=64.2，k α =0.973；n s=68.8，k α =0.985；n s=85.9，k α =0.974；n s=87，k α =0.968；n s=90.7，k α =0.964；n s=92.8，k α =0.945；n s=125.3，k α =0.941；n s=129.6，k α =0.932。
2） 设计工况。n s=64.2，k α =0.977；n s=68.8，k α =0.981；n s=85.9，k α =0.962；n s=87，k α =0.958；n s=90.7，k α =0.946；n s=92.8，k α =0.961；n s=125.3，k α =0.925；n s=129.6，k α =0.916。
3） 1.2倍设计工况。n s=64.2，k α =0.971；n s=68.8，k α =0.96；n s=85.9，k α =0.957；n s=87，k α =0.942；n s=90.7，k α =0.933；n s=92.8，k α =0.945；n s=125.3，k α =0.909；n s=129.6，k α =0.915。
（6）采用一元线性回归方法对不同比转数、3个相同工况下离心泵的液流角系数进行回归，建立3个工况下的液流角系数k α 与比转数n s之间的关系k α =f(n s/100)。
1） 五台不同比转数的离心泵液流角系数的一元线性回归，比转数n s分别为68.8、85.9、90.7、115.6和125.3。建立的液流角系数k α 与比转数n s之间的关系式依次为：0.8倍设计工况下，k α =1.0548-0.1002（n s/100）；设计工况下，k α =1.0634-0.1240（n s/100）；1.2倍设计工况下，k α =1.048-0.1224（n s/100）。
采用（5）中n s=92.8的离心泵液流角系数值对建立的3个工况下液流角系数k α 与比转数n s之间的关系式进行验证，其计算误差分别为：1.785%、-1.318%、-1.011%。
2） 八台不同比转数的离心泵液流角系数的一元线性回归，比转数n s分别为64.2、68.8、85.9、87、90.7、115.6、125.3和129.6。建立的液流角系数k α 与比转数n s之间的关系式依次为：0.8倍设计工况下，k α =1.036-0.0820（n s/100）；设计工况下，k α =1.0425-0.1004（n s/100）；1.2倍设计工况下，k α =1.0336-0.1037（n s/100）。
采用（5）中n s=92.8的离心泵液流角系数值对建立的3个工况下液流角系数k α 与比转数n s之间的关系式进行验证，其计算误差分别为：1.591%、-1.212%、-0.707%。
3） 离心泵液流角系数的计算误差的比较。从（6）的第1）步和第2）步中的计算误差可以看出：3个工况下，第2）步中的计算误差都比第1）步中的计算误差小。因此，不同比转数的离心泵液流角系数值越多，回归得到的液流角系数k α 与比转数n s之间的关系式的计算误差就越小。
（7）建立3个工况下的离心泵滑移系数计算公式。
Claims (1)
1.一种基于全工况内流测量的离心泵滑移系数确定方法，其特征在于，在离心泵全工况下叶轮出口处真实流动测量的基础上，将液流角系数引入斯托道拉滑移系数计算公式中，从而建立了一种离心泵全工况滑移系数的确定方法；具体步骤以下：
（A）采用PIV实验测量离心泵不同工况下的内部非定常流动规律；
采用相位平均法测量离心泵不同工况、不同相位下的内部非定常流动；采用二维PIV测量的具有时间序列的绝对速度场为：f（X，Y，t i ）、f（X，Y，t i +T）、f（X，Y，t i +2T）、……、f（X，Y，t i +nT），其中X、Y为二维流场的空间坐标，t i 为不同相位的初始时刻，i=0、1、……，T为叶轮旋转周期，n为拍摄流场图像的数量；将n组绝对速度场中对应坐标的瞬时流速进行算术平均即可得到不同工况、不同相位下离心泵内部流动的绝对速度分布；
（B）基于PIV实验得到的绝对速度分布和叶轮出口处的速度三角形，计算出不同工况、不同相位下的离心泵叶轮出口处的圆周分速度v u2，实验；
（C）计算不同工况下的离心泵叶轮出口处的圆周分速度v u2，理论，然后计算出不同工况、不同相位下的离心泵叶轮出口处的滑移速度，滑移速度为Δv u2=v u2，理论-v u2，实验；采用算术平均法得到离心泵不同工况下叶轮出口处的平均滑移速度 ，同时计算出离心泵不同工况下滑移系数，其中u 2为叶轮出口处的圆周速度；
（D）计算不同工况下的离心泵液流角系数k α ，液流角系数为，其中z为叶轮叶片数，β 2为叶片出口安放角；
（E）根据（A）、（B）、（C）和（D），对不同比转数n s的离心泵内部流动进行PIV测量，得到不同比转数、不同工况下的液流角系数；
（F）采用一元线性回归方法对不同比转数、不同工况下的液流角系数进行回归，建立不同工况下的液流角系数k α 与比转数n s之间的关系k α =f(n s/100)，从而确定了全工况下的离心泵滑移系数计算公式。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210382442.8A CN102865233B (zh) | 2012-10-11 | 2012-10-11 | 一种基于全工况内流测量的离心泵滑移系数确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210382442.8A CN102865233B (zh) | 2012-10-11 | 2012-10-11 | 一种基于全工况内流测量的离心泵滑移系数确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102865233A CN102865233A (zh) | 2013-01-09 |
| CN102865233B true CN102865233B (zh) | 2015-04-22 |
Family
ID=47444301
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210382442.8A Active CN102865233B (zh) | 2012-10-11 | 2012-10-11 | 一种基于全工况内流测量的离心泵滑移系数确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102865233B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107423473A (zh) * | 2017-04-26 | 2017-12-01 | 江苏大学 | 一种泵内粗颗粒滑移速度的计算方法 |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103742445B (zh) * | 2013-12-31 | 2017-01-11 | 江苏大学 | 一种核主泵极大流量水力设计方法 |
| CN105022879B (zh) * | 2015-07-21 | 2019-02-19 | 江苏大学 | 一种离心泵非定常特性的联合仿真分析方法 |
| CN107102163B (zh) * | 2017-04-26 | 2019-08-27 | 江苏大学 | 一种基于清水内流计算的泵内颗粒滑移速度计算方法 |
| CN110132452B (zh) * | 2019-05-13 | 2020-10-09 | 常州大学 | 一种利用中性层分叉确定复合梁界面滑移系数的方法 |
| CN111680372B (zh) * | 2020-06-10 | 2023-06-20 | 大连海事大学 | 考虑自然预旋时离心风机叶轮做功能力的一维计算方法 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4965713A (en) * | 1988-08-15 | 1990-10-23 | Viking Pump Inc. | Terminal element |
| EP0719944B1 (en) * | 1994-12-28 | 2002-05-29 | Ebara Corporation | Turbomachinery having a variable angle flow guiding device |
| JP3565938B2 (ja) * | 1995-03-14 | 2004-09-15 | トヨタ自動車株式会社 | 車両の制動力制御装置 |
| CA2184882A1 (en) * | 1995-09-08 | 1997-03-09 | Hideomi Harada | Turbomachinery with variable-angle flow guiding vanes |
| CN101956710B (zh) * | 2010-10-27 | 2012-05-23 | 江苏大学 | 一种基于损失的离心泵多工况水力优化方法 |
-
2012
- 2012-10-11 CN CN201210382442.8A patent/CN102865233B/zh active Active
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107423473A (zh) * | 2017-04-26 | 2017-12-01 | 江苏大学 | 一种泵内粗颗粒滑移速度的计算方法 |
| CN107423473B (zh) * | 2017-04-26 | 2020-09-25 | 江苏大学 | 一种泵内粗颗粒滑移速度的计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102865233A (zh) | 2013-01-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102865233B (zh) | 一种基于全工况内流测量的离心泵滑移系数确定方法 | |
| Wang et al. | The method for determining blade inlet angle of special impeller using in turbine mode of centrifugal pump as turbine | |
| CN103353908B (zh) | 一种基于数值计算的管路阻力系数精确计算方法 | |
| Kang et al. | Experimental study on the three-dimensional flow within a compressor cascade with tip clearance: Part I—Velocity and pressure fields | |
| Štefan et al. | Study of the internal flow field in a pump-as-turbine (PaT): Numerical investigation, overall performance prediction model and velocity vector analysis | |
| Zhang et al. | Experimental study of flow field in interference area between impeller and guide vane of axial flow pump | |
| CN101956710A (zh) | 一种基于损失的离心泵多工况水力优化方法 | |
| Shao et al. | Experimental and numerical study of external performance and internal flow of a molten salt pump that transports fluids with different viscosities | |
| CN103226635A (zh) | 基于三维动网格的旋转叶轮机械非定常流场的计算方法 | |
| Zhang et al. | Investigation of the integrated model of side chamber, wear-rings clearance, and balancing holes for centrifugal pumps | |
| CN111400941A (zh) | 一种叶片泵内部回流及回流旋涡空化的数值预测方法 | |
| CN102878091B (zh) | 一种基于内流测量的离心泵全工况理论扬程确定方法 | |
| Aakti et al. | On the performance of a high head Francis turbine at design and off-design conditions | |
| CN103883556B (zh) | 一种轴流转轮的设计方法 | |
| Chen et al. | Solid-liquid multiphase flow and erosion in the energy storage pump using modified drag model and erosion model | |
| Capurso et al. | Dependency of the slip phenomenon on the inertial forces inside radial runners | |
| CN104234996A (zh) | 多转子螺杆泵间隙泄露量检测方法 | |
| Pospelov et al. | Effect of the parameters of a computational model on the prediction of hydraulic turbine characteristics | |
| CN104179670B (zh) | 应用再制造技术后水泵机组节能率的测量方法 | |
| CN102251955A (zh) | 中低比转速泵/风机变速调节运行工况的测量方法 | |
| 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 | |
| Shi et al. | Optimal design of multi-conditions for axial flow pump | |
| Li et al. | Design and analysis of flow rectifier of gas turbine flowmeter | |
| Talukder et al. | Effects of number of stator blades on the performance of a torque converter | |
| Anagnostopoulos | Numerical calculation of the flow in a centrifugal pump impeller using Cartesian grid |
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 |