CN108412805B - 一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 - Google Patents
一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 Download PDFInfo
- Publication number
- CN108412805B CN108412805B CN201810184905.7A CN201810184905A CN108412805B CN 108412805 B CN108412805 B CN 108412805B CN 201810184905 A CN201810184905 A CN 201810184905A CN 108412805 B CN108412805 B CN 108412805B
- Authority
- CN
- China
- Prior art keywords
- inducer
- coefficient
- cylindrical
- design
- 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
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/18—Rotors
- F04D29/22—Rotors specially for centrifugal pumps
- F04D29/2238—Special flow patterns
- F04D29/225—Channel wheels, e.g. one blade or one flow channel
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mechanical Engineering (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明提供了一种圆柱形诱导轮设计方法,涉及流体机械技术领域。该圆柱形诱导轮设计方法,将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;利用二维截面涡量对对称流动进行修正,从而得到诱导轮的几何参数及诱导轮的无空化特性参数。本发明的圆柱形诱导轮设计方法,基于二元流动模型进行设计,比一元流动理论的分析设计方法具有更高的精度和设计效率。在此基础上,本发明还提供了一种圆柱形诱导轮和圆柱形诱导轮流场。
Description
技术领域
本发明涉及流体机械技术领域,具体而言,涉及一种圆柱形诱导轮设计方法、圆柱形诱导轮及圆柱形诱导轮流场。
背景技术
诱导轮是低温液体火箭发动机中的关键部件,用于输送液氢液氧至燃烧室,其运行转速高,并且时常伴随着复杂的因空化流动导致的不稳定现象。
其中,诱导轮常见的有圆柱形和圆锥形的外形结构。对于圆柱形诱导轮,技术人员一般采用基于无限叶片假设的一元流动理论分析方法进行设计。
基于无限叶片假设的一元流动理论分析方法,只考虑使得主离心泵叶轮不发生空化时的诱导轮扬程,估算诱导轮空化比转速,再通过扬程系数及诱导轮直径确定进口流量系数。此设计方法预估的性能与典型诱导轮的实际性能除在设计点处较为吻合外,其他区域均存在较大的差距,有待进一步改进。另外,上述设计方法除了难以在全结构尺寸上满足性能要求之外,其计算的中间过程和步骤也并不简单。
所以,如何改良圆柱形诱导轮的设计方法,提高设计效率和精度是本领域技术人员亟待解决的问题。
发明内容
本发明的目的在于提供一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场,以解决现有技术中的诱导轮设计方法存在的设计精度差、计算过程和步骤繁复等问题,并获得一种高效的圆柱形诱导轮及其无空化特性的流场。
为了实现上述目的,本发明采用以下技术方案:
本发明提供的一种圆柱形诱导轮设计方法,将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;利用二维截面涡量对对称流动进行修正,得到诱导轮的几何参数及诱导轮的无空化特性参数。
需要说明的是,上述圆柱形诱导轮设计方法与传统的一元流动理论分析方法不同,是基于二元流动模型的给出的一种新型的圆柱形诱导轮的设计方法。
在上述技术方案的基础上,进一步,其已知输入条件包括:叶片数N,转速Ω,叶尖半径rT,轮毂半径rH,进口叶尖安装角γTle,轴向长度ca和设计流量系数φD;设定诱导轮前缘阻塞系数Ble和尾缘阻塞系数Bte均等于1,包括如下步骤:
步骤一:结合已知输入条件、前缘阻塞系数Ble和尾缘阻塞系数Bte,计算出口导程P、叶片弦长c、叶片实度σ、扩散系数D、叶片冲角α以及尾缘阻塞系数估算值B'te;判断设计结果是否满足第一检验条件,如满足,进行步骤二的设计;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至设计结果满足第一检验条件。
步骤二:计算诱导轮总扬程系数ψt和静扬程系数ψ,判断计算结果是否满足第二检验条件,如满足,结束设计过程;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至计算结果满足第二检验条件。
在上述技术方案的基础上,进一步,第一检验条件包括:0<扩散系数D<0.5、叶片实度σ>2、叶片冲角α/叶片安放角βb<0.5。
本发明还提供一种由上述圆柱形诱导轮设计方法设计得到圆柱形诱导轮。
本发明还提供一种由上述圆柱形诱导轮设计方法设计得到圆柱形诱导轮流场。
本发明具有如下有益效果:
1、本发明提供的圆柱形诱导轮设计方法,将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;再利用二维截面涡量对对称流动进行修正,得到圆柱形诱导轮的几何参数及圆柱形诱导轮的无空化特性参数,与现有的基于无限叶片假设的一元流动理论分析方法所得到的数据,更满足实际性能的需求。
2、本发明提供的圆柱形诱导轮设计方法,可经过多次迭代计算,简化了中间流程,提高了设计效率。
本发明的附加技术特征及其优点将在下面的描述内容中阐述地更加明显,或通过本发明的具体实践可以了解到。
附图说明
为了更清楚地说明本发明具体实施方式的技术方案,下面将对具体实施方式描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的圆柱形诱导轮设计方法的流程图;
图2为本发明提供的圆柱形诱导轮设计方法中三种不同叶栅的边界层动量厚度/弦长之比与扩散系数的关系图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
现有技术说明:
在现有技术中,圆柱形诱导轮普遍采用基于无限叶片假设的一元流动理论分析方法进行设计。基于无限叶片假设的一元流动理论分析方法,只考虑使得主离心泵叶轮不发生空化时的诱导轮扬程,估算诱导轮空化比转速,再通过扬程系数及诱导轮直径确定进口流量系数。此设计方法预估的性能与典型诱导轮的实际性能除在设计点处较为吻合外,其他区域均存在较大的差距,有待进一步改进。另外,上述设计方法除了难以在全结构尺寸上满足性能要求之外,其计算的中间过程和步骤也并不简单。
针对上述现有技术方案存在的技术问题,下面结合具体的实施方式对本发明做进一步地解释说明:
本实施例提供了一种圆柱形诱导轮设计方法,其总体设计思路和流程包括:将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;利用二维截面涡量对对称流动进行修正,得到诱导轮的几何参数及诱导轮的无空化特性参数。
需要说明的是,上述圆柱形诱导轮设计方法与传统的一元流动理论分析方法不同,是基于二元流动模型的给出的一种新型的圆柱形诱导轮的设计方法。
本发明提供的圆柱形诱导轮设计方法,能够较好地解决现有技术中的诱导轮设计方法存在的设计精度差、计算过程和步骤繁复等问题,并获得一种高效的圆柱形诱导轮及其无空化特性的流场:首先,本设计方法得到的圆柱形诱导轮的几何参数及圆柱形诱导轮的无空化特性参数,与现有的基于无限叶片假设的一元流动理论分析方法所得到的数据,更满足实际性能的需求。其次,本设计方法可经过多次迭代计算,简化了中间流程,提高了设计效率。
具体地,在本发明提供的圆柱形诱导轮设计方法中,其已知输入条件包括:叶片数N,转速Ω,叶尖半径rT,轮毂半径rH,进口叶尖安装角γTle,轴向长度ca和设计流量系数φD;设定诱导轮前缘阻塞系数Ble和尾缘阻塞系数Bte均等于1。设计方法包括如下步骤:
步骤一:结合已知输入条件、前缘阻塞系数Ble和尾缘阻塞系数Bte,计算出口导程P、叶片弦长c、叶片实度σ、扩散系数D、叶片冲角α以及尾缘阻塞系数估算值B'te;判断设计结果是否满足第一检验条件,如满足,进行步骤二的设计;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至设计结果满足第一检验条件。
步骤二:计算诱导轮总扬程系数ψt和静扬程系数ψ,判断计算结果是否满足第二检验条件,如满足,结束设计过程;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至计算结果满足第二检验条件。
图1为本发明提供的圆柱形诱导轮设计方法的流程图。如图1所示,圆柱形诱导轮设计方法的具体设计流程如下:
步骤1、输入操作参数和几何参数:
包括叶片数N,转速Ω,叶尖半径rT,轮毂半径rH,进口叶尖安装角γTle,轴向长度ca和设计流量系数φD;同时,设定诱导轮前缘阻塞系数Ble和尾缘阻塞系数Bte均等于1,即忽略流道阻塞情况的影响。
步骤2、计算进口叶尖导程和轴向导程变化关系,求得出口导程:
进口叶尖导程:PTle=2πrTtanγTle,
轴向导程关系:P=PTle+P′z,
出口导程:P=PTle+P′×ca=2πrtetanγte;
步骤4、求解进出口平均直径处相对速度:
出口相对速度:
步骤5、求解平均流线上叶片实度:
平均流线为:
其中,初始条件
则叶片弦长的平均值和叶片稠度为:
步骤6、求解扩散系数及其对应阻塞系数:
步骤7、根据第一检验条件判断上述设计结果是否满足以下范围要求:
叶栅扩散系数:0<D<0.5,
叶片稠度:σ<2,
若满足要求,进行步骤8;若不满足要求,回到步骤2重新迭代;
步骤8、尾缘出口流场计算:
尾缘滑移速度的切向分量:
其中对常数d有:
步骤9、流动偏移预估:v2δ°=Ωr-w2tan(β′2+δ°),
其中:
步骤10、流动损失预估:
步骤11、计算诱导轮扬程系数:
判断总扬程系数ψt和静扬程系数ψ是否满足设计要求,若是,结束设计过程;若否,回到步骤2重新迭代,直至得到满足设计要求的数据。
综上所述,在本发明提供的圆柱形诱导轮设计方法中,首先,忽略边界层阻塞的影响,可得轮毂和叶片的初步预估情况。由此产生的流场将用于估计叶栅扩散系数和流动阻塞,以实现圆柱形诱导轮几何参数的修正和流场的计算。应用3D造型可获得大体连续的几何模型。如果求得的圆柱形诱导轮几何结构不能被接受,程序将不断迭代输入数据直至输出最优结果并进行条件检验。
迭代法在求得圆柱形诱导轮的几何参数的同时还可输出泵的非空化特性。若求得的泵特性不满足使用需求,将对操作参数和几何参数进行修改并重复上述流程,直至输出结果能满足要求。
本发明技术方案具有至少如下优点:
1、本发明提供的圆柱形诱导轮设计方法,将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;再利用二维截面涡量对对称流动进行修正,得到诱导轮的几何参数及诱导轮的无空化特性参数,与现有的基于无限叶片假设的一元流动理论分析方法所得到的数据,更满足实际性能的需求。
2、本发明提供的圆柱形诱导轮设计方法,可经过多次迭代计算,简化了中间流程,提高了设计效率。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (2)
1.一种圆柱形诱导轮设计方法,其已知输入条件包括:叶片数N,转速Ω,叶尖半径rT,轮毂半径rH,进口叶尖安装角γTle,轴向长度ca和设计流量系数φD;设定诱导轮前缘阻塞系数Ble和尾缘阻塞系数Bte均等于1,其特征在于,将诱导轮内的三维流动函数分解为对称流动函数和二维截面涡流函数;利用二维截面涡量对对称流动进行修正,得到诱导轮的几何参数及诱导轮的无空化特性参数;
包括如下步骤:
步骤一:结合已知输入条件、前缘阻塞系数Ble和尾缘阻塞系数Bte,计算出口导程P、叶片弦长c、叶片实度σ、扩散系数D、叶片冲角α以及尾缘阻塞系数估算值B'te;判断设计结果是否满足第一检验条件,如满足,进行步骤二的设计;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至设计结果满足第一检验条件;
步骤二:计算诱导轮总扬程系数ψt和静扬程系数ψ,判断计算结果是否满足第二检验条件,如满足,结束设计过程;如不满足,则将尾缘阻塞系数估算值B'te的数值赋值到尾缘阻塞系数Bte中,重新进行步骤一的运算,直至计算结果满足第二检验条件;
第一检验条件包括:0<扩散系数D<0.5、叶片实度σ>2、叶片冲角α/叶片安放角βb<0.5。
2.一种圆柱形诱导轮,其特征在于,由上述权利要求1所述的圆柱形诱导轮设计方法设计得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810184905.7A CN108412805B (zh) | 2018-03-06 | 2018-03-06 | 一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810184905.7A CN108412805B (zh) | 2018-03-06 | 2018-03-06 | 一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108412805A CN108412805A (zh) | 2018-08-17 |
CN108412805B true CN108412805B (zh) | 2020-01-07 |
Family
ID=63130077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810184905.7A Active CN108412805B (zh) | 2018-03-06 | 2018-03-06 | 一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108412805B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102678617A (zh) * | 2012-05-18 | 2012-09-19 | 江苏大学 | 一种基于离心泵的诱导轮设计方法 |
CN102705263A (zh) * | 2012-06-21 | 2012-10-03 | 江苏大学 | 一种离心泵变螺距诱导轮的优化设计方法 |
CN103104544A (zh) * | 2013-03-07 | 2013-05-15 | 江苏大学 | 一种长短叶片诱导轮的变螺距设计方法 |
CN103307008A (zh) * | 2012-03-08 | 2013-09-18 | 江苏大学 | 一种核主泵的带长短叶片的诱导轮设计方法 |
CN104047890A (zh) * | 2014-07-25 | 2014-09-17 | 江苏大学 | 一种轴流式低扬程前置诱导轮的设计方法 |
CN104500438A (zh) * | 2014-11-21 | 2015-04-08 | 江苏国泉泵业制造有限公司 | 一种两相流泵的水力设计方法 |
CN107273570A (zh) * | 2017-05-10 | 2017-10-20 | 江苏大学 | 一种叶片泵空化诱导噪声数值预报方法 |
-
2018
- 2018-03-06 CN CN201810184905.7A patent/CN108412805B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103307008A (zh) * | 2012-03-08 | 2013-09-18 | 江苏大学 | 一种核主泵的带长短叶片的诱导轮设计方法 |
CN102678617A (zh) * | 2012-05-18 | 2012-09-19 | 江苏大学 | 一种基于离心泵的诱导轮设计方法 |
CN102705263A (zh) * | 2012-06-21 | 2012-10-03 | 江苏大学 | 一种离心泵变螺距诱导轮的优化设计方法 |
CN103104544A (zh) * | 2013-03-07 | 2013-05-15 | 江苏大学 | 一种长短叶片诱导轮的变螺距设计方法 |
CN104047890A (zh) * | 2014-07-25 | 2014-09-17 | 江苏大学 | 一种轴流式低扬程前置诱导轮的设计方法 |
CN104500438A (zh) * | 2014-11-21 | 2015-04-08 | 江苏国泉泵业制造有限公司 | 一种两相流泵的水力设计方法 |
CN107273570A (zh) * | 2017-05-10 | 2017-10-20 | 江苏大学 | 一种叶片泵空化诱导噪声数值预报方法 |
Non-Patent Citations (5)
Title |
---|
"一体式诱导轮与叶轮航空离心泵汽蚀特性研究";李嘉等;《推进技术》;20150731;第36卷(第7期);第1005-1012页 * |
"多工况高抗汽蚀性能的诱导轮设计";李永鹏等;《排灌机械工程学报》;20141130;第32卷(第11期);第948-954页 * |
"液氢诱导轮内部流场数值研究";衣同训等;《推进技术》;20000830;第21卷(第4期);第9-11页 * |
"离心泵诱导轮的数值模拟";刘德民等;《水泵技术》;20080630(第6期);第28-32页 * |
"诱导轮三维设计方法及其性能验证";吴苏青;《通用机械》;20180131;第61-64页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108412805A (zh) | 2018-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lee et al. | Impeller design of a centrifugal fan with blade optimization | |
CN107693868B (zh) | 心脏泵叶轮设计方法和叶轮 | |
CN109625226B (zh) | 一种轴流式高功率密度喷水推进泵设计方法 | |
CN104895832B (zh) | 一种高粘度离心泵叶轮的水力设计方法 | |
CN109595179B (zh) | 一种具有小轮毂比叶轮的排水泵 | |
CN101245711A (zh) | 轴流式叶轮机械的设计方法及叶片 | |
CN110321660B (zh) | 一种径向出水大型混流泵叶轮设计方法 | |
CN110617238A (zh) | 一种离心泵叶轮的优化设计方法 | |
CN108412805B (zh) | 一种圆柱形诱导轮设计方法、圆柱形诱导轮及其流场 | |
CN105673558A (zh) | 一种基于载荷法设计的离心通风机叶片 | |
CN104533828B (zh) | 一种双向轴流泵水力设计方法 | |
CN104047890B (zh) | 一种轴流式低扬程前置诱导轮的设计方法 | |
CN108205607A (zh) | 一种高比转速离心泵叶轮的水力设计方法 | |
CN110608191B (zh) | 一种基于奥森涡的叶片设计方法及其设计的叶片泵 | |
KR101808076B1 (ko) | 효율 및 흡입성능을 동시에 최적화하는 원심 및 사류 펌프 임펠러 최적화 설계 방법, 이에 의하여 설계된 원심 및 사류 펌프 임펠러 및 원심 및 사류 펌프 | |
WO2020134126A1 (zh) | 一种小轮毂比叶轮的设计方法及使用该方法获得的轮缘泵 | |
KR101984022B1 (ko) | 콤팩트하고 단순화한 사류 펌프의 디퓨저 설계 방법, 이에 의하여 설계된 디퓨저 및 이를 구비한 사류 펌프 | |
KR102613148B1 (ko) | 축류펌프 임펠러 설계 방법, 이에 의하여 설계된 임펠러 및 펌프 | |
Shin et al. | CFD-Based optimization of rudder bulb systems | |
CN111079367B (zh) | 适用于不同进水锥管的轴流泵辐条优化设计方法 | |
CN112761969A (zh) | 一种离心泵结构参数的优化方法 | |
CN113250996A (zh) | 一种诱导轮的制作方法 | |
CN113202787A (zh) | 一种蜗壳式离心泵必需汽蚀余量数值模拟预测方法 | |
CN212508963U (zh) | 一种混流泵或轴流泵 | |
Zhang et al. | Blade thickness effect on impeller slip factor |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |