[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN111723447B - 一种动静腔流动微沟槽减阻设计方法 - Google Patents

一种动静腔流动微沟槽减阻设计方法 Download PDF

Info

Publication number
CN111723447B
CN111723447B CN202010610372.1A CN202010610372A CN111723447B CN 111723447 B CN111723447 B CN 111723447B CN 202010610372 A CN202010610372 A CN 202010610372A CN 111723447 B CN111723447 B CN 111723447B
Authority
CN
China
Prior art keywords
micro
groove
dynamic
cross
drag reduction
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
CN202010610372.1A
Other languages
English (en)
Other versions
CN111723447A (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.)
Yangzhou University
Original Assignee
Yangzhou University
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 Yangzhou University filed Critical Yangzhou University
Priority to CN202010610372.1A priority Critical patent/CN111723447B/zh
Publication of CN111723447A publication Critical patent/CN111723447A/zh
Application granted granted Critical
Publication of CN111723447B publication Critical patent/CN111723447B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force 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)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Moulds For Moulding Plastics Or The Like (AREA)

Abstract

本发明涉及一种动静腔流动微沟槽减阻设计方法。考虑动盘近壁区周向和径向流动，建立了微沟槽的空间型线方程组。首先，给定动静腔工作参数，采用数值方法求解该方程组，得到微沟槽的空间型线。然后，选取微沟槽横截面形状，设计微沟槽横截面高度和弧长，选取微沟槽数量。最后，结合微沟槽空间型线、横截面几何参数和数量，在动盘上加工出微沟槽。本发明通过在动盘上加工微沟槽达到减阻目的。

Description

一种动静腔流动微沟槽减阻设计方法
技术领域
本发明涉及旋转机械动静腔流动减阻领域，尤其是一种动静腔流动微沟槽减阻设计方法。
背景技术
动静腔是指旋转圆盘与静止圆盘之间的流动区域，如图1所示，既是一个经典的流体力学研究对象，又是一个广泛存在于石油化工、水利水电、海洋工程、航天动力、核电装备等领域旋转机械中的工程问题。在大多旋转机械中，动静腔是不可避免的形成的，动盘与腔内流体的摩擦功率属于能量损失，减小动盘转动阻力(即减小壁面周向切应力)是实现旋转机械节能降耗的一个重要途径。
在平板流动中，仿鲨鱼皮微沟槽减阻机理可以总结为：(1)粘性底层和缓冲层的厚度增大，对数律区抬升，近壁区平均速度梯度减小；(2)流向沟槽限制了流向涡的展向运动等，湍流猝发强度、雷诺应力等统计量减弱。由于平板上主流是直线型，因而沟槽设计成直线流向型，不需要设计成曲线型。根据微沟槽减阻理论，在动静腔中实施微沟槽减阻策略。显然，动静腔流动中的微沟槽不能设计成平板流动中的直线型。这是因为动盘壁面对流体的剪切作用、腔体对流体的约束等，动盘近壁区流动方向为周向和径向，其流线是曲线型的。因此，为达到良好的减阻效果，动盘壁面上的沟槽应该设计成曲线型的，但该曲线设计方法较缺乏。
发明内容
本发明的目的就是针对上述存在的问题，提供一种动静腔流动微沟槽减阻设计方法。
本发明的目的是这样实现的，一种新的微沟槽减阻设计方法，分为以下步骤：
S1、考虑动盘近壁区周向和径向流动，设计微沟槽空间型线。
动盘近壁区流体径向速度：
Cr＝1.18·(10-5·Rer+2)-0.49 (2)
vr是流体径向速度，单位m/s，径向坐标是流体径向速度的自变量；Cr是动盘径向速度系数；r是径向坐标，单位m；υ是流体运动粘度，m2/s；ω是动盘角速度，单位rad/s；K是携卷系数，流体周向速度与动盘周向速度之比；z是轴向坐标，单位m，以动盘表面为基准；δr动盘边界层厚度，单位m；Rer是雷诺数；Cqr是径向流量系数；Q是动静腔通流量，单位m3/s。
给定动盘转速ω和通流量Q；选取边界层高度系数的值，取值范围是/>采用数值方法求解式(1)-(6)，计算出[Rs1,Rs2]区间内的径向速度分布。Rs1为微沟槽前缘半径，单位m，取值范围是0≤Rs1≤0.9·R，R是动盘半径，单位m；Rs2为微沟槽尾缘半径，单位m，与动盘半径R相等，即Rs2＝R。
流体角速度：
是角速度，单位rad/s，径向坐标是流体角速度的自变量。给定动盘转速ω和通流量Q，采用数值方法求解式(3)、(5)、(6)、(7)，计算出[Rs1,R]区间内的角速度分布。
选取时间步Δt，选取原则是：使得型线光顺。型线计算公式如下。
r0＝Rs1 (10)
式中，r是极径；是极角。将rn-1代入式(3)、(5)、(6)、(7)，计算出式(8)中的/>将rn-1代入式(1)-(6)，计算出式(9)中的vr n-1。通过式(8)和(9)的不断累加，直至rn≥R，就可得到[Rs1,R]区间内n+1或n个坐标点(/>r)，描绘出微沟槽的空间型线。
S2、确定微沟槽其它设计参数的值。
微沟槽横截面形状为V形、U形、矩形。为横截面高度系数，取值范围是/>其中，z根据S1中的δr和/>计算得到。h是微沟槽横截面高度，单位m。横截面形状参数/>的取值范围为/>s是微沟槽横截面弧长，单位m；横截面高度和弧长均随半径变化，在径向上两者比值保持不变。选取微沟槽数量Z，使得微沟槽稠密度为/>
通过本发明，提供了一种动静腔流动微沟槽减阻设计方法。考虑动盘近壁区周向和径向流动，建立了微沟槽的空间型线方程组。首先，给定动静腔工作参数，采用数值方法求解该方程组，得到微沟槽的空间型线。然后，选取微沟槽横截面形状，设计微沟槽横截面高度和弧长，选取微沟槽数量。最后，结合微沟槽空间型线、横截面几何参数和数量，在动盘上加工出微沟槽。本发明通过在动盘上加工微沟槽达到减阻目的。
附图说明
图1为动静腔示意图；
图2为微沟槽型线示意图；
图3a为微沟槽横截面形状(V形)；
图3b为微沟槽横截面形状(U形)；
图3c为微沟槽横截面形状(矩形)；
图4为微沟槽示意图；
图5为带微沟槽的动盘示意图；
图中：1动静腔、2动盘、3静盘、h微沟槽横截面高度、R动盘半径、Rs1微沟槽前缘半径、Rs2微沟槽尾缘半径、s微沟槽横截面弧长。
具体实施方式
下面结合具体实施方式对本发明作进一步的详细说明。
一种动静腔流动微沟槽减阻设计方法，包括以下步骤：
S1、微沟槽空间型线方程；
流体径向速度：
Cr＝1.18·(10-5·Rer+2)-0.49 (2)
vr是流体径向速度，径向坐标是其自变量；Cr是动盘径向速度系数；r是径向坐标；ω是动盘角速度；K是携卷系数；z是轴向坐标；δr动盘边界层厚度；Rer是雷诺数；Cqr是径向流量系数；Q是动静腔通流量；
给定动盘转速ω和通流量Q；选取边界层高度系数的值；采用数值方法求解式(1)-(6)，计算出[Rs1,Rs2]区间内的径向速度分布，Rs1为微沟槽前缘半径，Rs2为微沟槽尾缘半径；
流体角速度：
是角速度，径向坐标是其自变量；给定动盘转速ω和通流量Q，采用数值方法求解式(3)、(5)、(6)、(7)，计算出[Rs1,R]区间内的角速度分布，R是动盘半径；
选取时间步Δt，选取原则是：使得型线光顺，型线计算公式如下：
rn＝rn-1+vr n-1·Δt (9)
r0＝Rs1 (10)
式中，r是极径；是极角，将rn-1代入式(3)、(5)、(6)、(7)，计算出式(8)中的/>将rn-1代入式(1)-(6)，计算出式(9)中的vr n-1，通过式(8)和(9)的不断累加，直至rn≥R，就可得到[Rs1,R]区间内n+1或n个坐标点(/>r)，描绘出微沟槽的空间型线；
S2、确定微沟槽其它设计参数的值；
确定微沟槽横截面形状，选取横截面高度系数横截面形状参数/>微沟槽数量Z，z是S1中设置的值，h是微沟槽横截面高度，s是微沟槽横截面弧长，横截面高度和弧长均随半径变化，在径向上两者比值保持不变。
其中，边界层高度系数的取值范围是/>微沟槽前缘半径Rs1的取值范围是0≤Rs1≤0.9·R。微沟槽尾缘半径Rs2与动盘半径R相等。微沟槽横截面形状为V形、U形、矩形。横截面高度系数的取值范围是/>其中，z根据S1中的δr和/>计算得到。横截面形状参数/>的取值范围是/>选取微沟槽数量Z，使得微沟槽稠密度为/>
以一动静腔为例，如图1所示，动盘转速500rpm，动静腔向心通流量0.01m3/s，动盘半径160mm。
S1考虑动盘近壁区周向和径向流动，设计微沟槽空间型线。
在MATLAB中编写式(1)-(11)的求解程序，选取边界层高度系数选取Rs1＝0.5·R，选取时间步Δt＝0.05s，得到136个坐标点(/>r)，如表1所示。设计出的微沟槽的空间型线如图2所示。
表1微沟槽型线坐标点
S2确定微沟槽其它设计参数的取值范围。
微沟槽横截面形状为V形、U形、矩形，如图3a所示，本实施例采用V形。选取横截面高度系数选取/>微沟槽型线各点的横截面高度h和横截面弧长s如表2所示。设计出的单个微沟槽如图4所示。
表2微沟槽横截面几何参数
微沟槽数量设置为10000，则稠密度图5是带微沟槽的动盘示意图。使用ANSYS分别对带光滑动盘的动静腔和带微沟槽动盘的动静腔进行数值模拟。模拟结果显示，带微沟槽的动盘转矩是光滑动盘的98.21％，达到了减阻效果。
所述实施例为本发明的优选的实施方式，但本发明并不限于上述实施方式，在不背离本发明的实质内容的情况下，本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。

Claims (8)

1.一种动静腔流动微沟槽减阻设计方法，其特征在于，包括以下步骤：
S1、微沟槽空间型线方程；
流体径向速度：
Cr＝1.18·(10-5·Rer+2)-0.49 (2)
vr是流体径向速度，径向坐标是其自变量；Cr是动盘径向速度系数；r是径向坐标；ω是动盘角速度；K是携卷系数；z是轴向坐标；δr动盘边界层厚度；Rer是雷诺数；Cqr是径向流量系数；Q是动静腔通流量；
给定动盘转速ω和通流量Q；选取边界层高度系数的值；采用数值方法求解式(1)-(6)，计算出[Rs1,Rs2]区间内的径向速度分布，Rs1为微沟槽前缘半径，Rs2为微沟槽尾缘半径；
流体角速度：
是角速度，径向坐标是其自变量；给定动盘转速ω和通流量Q，采用数值方法求解式(3)、(5)、(6)、(7)，计算出[Rs1,R]区间内的角速度分布，R是动盘半径；
选取时间步Δt，选取原则是：使得型线光顺，型线计算公式如下：
rn＝rn-1+vr n-1·Δt (9)
r0＝Rs1 (10)
式中，r是极径；是极角，将rn-1代入式(3)、(5)、(6)、(7)，计算出式(8)中的/>将rn-1代入式(1)-(6)，计算出式(9)中的vr n-1，通过式(8)和(9)的不断累加，直至rn≥R，就可得到[Rs1,R]区间内n+1或n个坐标点(/>r)，描绘出微沟槽的空间型线；
S2、确定微沟槽其它设计参数的值；
确定微沟槽横截面形状，选取横截面高度系数横截面形状参数/>微沟槽数量Z，z是S1中设置的值，h是微沟槽横截面高度，s是微沟槽横截面弧长，横截面高度和弧长均随半径变化，在径向上两者比值保持不变。
2.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，边界层高度系数的取值范围是/>
3.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，微沟槽前缘半径Rs1的取值范围是0≤Rs1≤0.9·R。
4.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，微沟槽尾缘半径Rs2与动盘半径R相等。
5.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，微沟槽横截面形状为V形、U形、矩形。
6.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，横截面高度系数的取值范围是其中，z根据S1中的δr和/>计算得到。
7.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，横截面形状参数的取值范围是/>
8.根据权利要求1所述的一种动静腔流动微沟槽减阻设计方法，其特征在于，选取微沟槽数量Z，使得微沟槽稠密度为
CN202010610372.1A 2020-06-30 2020-06-30 一种动静腔流动微沟槽减阻设计方法 Active CN111723447B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010610372.1A CN111723447B (zh) 2020-06-30 2020-06-30 一种动静腔流动微沟槽减阻设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010610372.1A CN111723447B (zh) 2020-06-30 2020-06-30 一种动静腔流动微沟槽减阻设计方法

Publications (2)

Publication Number Publication Date
CN111723447A CN111723447A (zh) 2020-09-29
CN111723447B true CN111723447B (zh) 2023-11-14

Family

ID=72571966

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010610372.1A Active CN111723447B (zh) 2020-06-30 2020-06-30 一种动静腔流动微沟槽减阻设计方法

Country Status (1)

Country Link
CN (1) CN111723447B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111737836B (zh) * 2020-06-30 2021-06-22 扬州大学 一种动静腔流动微沟槽-超疏水耦合减阻设计方法
CN112253452B (zh) * 2020-10-16 2022-02-22 扬州大学 一种带螺旋形流道的微型圆盘泵设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5637942A (en) * 1994-10-18 1997-06-10 Arthur D. Little, Inc. Aerodynamic drag reduction arrangement for use with high speed rotating elements
WO2008121418A1 (en) * 2007-03-30 2008-10-09 The University Of Alabama A passive drag modification system
CN107742011A (zh) * 2017-09-26 2018-02-27 南京航空航天大学 叶轮叶片减阻微织构的设计方法
CN110635583A (zh) * 2018-08-31 2019-12-31 北京金风科创风电设备有限公司 电磁装置的铁心及其叠片

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5637942A (en) * 1994-10-18 1997-06-10 Arthur D. Little, Inc. Aerodynamic drag reduction arrangement for use with high speed rotating elements
WO2008121418A1 (en) * 2007-03-30 2008-10-09 The University Of Alabama A passive drag modification system
CN107742011A (zh) * 2017-09-26 2018-02-27 南京航空航天大学 叶轮叶片减阻微织构的设计方法
CN110635583A (zh) * 2018-08-31 2019-12-31 北京金风科创风电设备有限公司 电磁装置的铁心及其叠片

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于RNG紊流模型的立式轴流泵站三维流动数值模拟及性能预测;成立;刘超;汤方平;周济人;;机械工程学报(第03期);全文 *
基于大涡模拟的圆管脉动湍流减阻数值分析;宁涛;顾春伟;;清华大学学报(自然科学版)(第05期);全文 *

Also Published As

Publication number Publication date
CN111723447A (zh) 2020-09-29

Similar Documents

Publication Publication Date Title
Wu et al. Application of CFD to model oil–air flow in a grooved two-disc system
CN111723447B (zh) 一种动静腔流动微沟槽减阻设计方法
Chitrakar et al. Study of the simultaneous effects of secondary flow and sediment erosion in Francis turbines
Guo et al. Effects of the number of inducer blades on the anti-cavitation characteristics and external performance of a centrifugal pump
CN111460699B (zh) 平壁表面减阻功能微织构的设计方法
De Maesschalck et al. Blade tip carving effects on the aerothermal performance of a transonic turbine
Li et al. Improved soft abrasive flow finishing method based on turbulent kinetic energy enhancing
Zhou et al. A novel suction-side winglet design philosophy for high-pressure turbine rotor tips
Zhong et al. Effects of tip gap size on the aerodynamic performance of a cavity-winglet tip in a turbine cascade
Desando et al. Numerical analysis of honeycomb labyrinth seals: Cell geometry and fin tip thickness impact on the discharge coefficient
Chitrakar et al. Particle image velocimetry investigation of the leakage flow through clearance gaps in cambered hydrofoils
Liu et al. Numerical study of scale effects on the open water performance of a rim-driven thruster
Cao et al. Effect of gap flow on the torque for blades in a rim driven thruster without axial pressure gradient
Han et al. Role of wall roughness on interaction of leakage flow and main flow in a mixed flow pump with tip clearance
CN105045987B (zh) 孔板厚度对孔板能量损失系数的影响关系计算方法
An et al. Research on the mechanism of drag reduction and efficiency improvement of hydraulic retarders with bionic non-smooth surface spoilers
Das et al. The effect of midplane guide vanes in a biplane wells turbine
CN111737836B (zh) 一种动静腔流动微沟槽-超疏水耦合减阻设计方法
Ma et al. Experimental study of effects of grooved tip clearances on the flow field in a compressor cascade passage
Yang et al. Effects of Reynolds number and tooth front angle on leakage loss and heat transfer characteristics in a rotating labyrinth seal
CN113044201B (zh) 一种具有主动射流结构的翼型
Myronov et al. Improving the Energy Performance of a High-Head Francis Turbine
Jiao et al. Investigations of an enclosed annular rotor-stator system with LES method
Prasad et al. Poiseuille flow across an eccentrically confined stationary/rotating cylinder
Will et al. Investigation of the fluid flow in a rotor-stator cavity with inward through-flow

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