CN111737836B - 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 - Google Patents
一种动静腔流动微沟槽-超疏水耦合减阻设计方法 Download PDFInfo
- Publication number
- CN111737836B CN111737836B CN202010610401.4A CN202010610401A CN111737836B CN 111737836 B CN111737836 B CN 111737836B CN 202010610401 A CN202010610401 A CN 202010610401A CN 111737836 B CN111737836 B CN 111737836B
- Authority
- CN
- China
- Prior art keywords
- micro
- groove
- dynamic
- drag reduction
- cross
- 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
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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 OR CALCULATING; 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 OR CALCULATING; COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Application Of Or Painting With Fluid Materials (AREA)
Abstract
本发明涉及一种动静腔流动微沟槽‑超疏水耦合减阻设计方法,利用动盘近壁区流动数学模型,构建微沟槽型线控制方程。首先,代入动静腔工作参数,计算微沟槽型线控制方程,得到微沟槽的型线数据。然后,选取微沟槽横截面形状,设计微沟槽横截面高度和弧长,选取微沟槽数量,在动盘上加工出微沟槽。最后,在带微沟槽的动盘上喷涂超疏水材料,保证涂层不淹没微沟槽。本发明通过在动盘上实施微沟槽‑超疏水耦合减阻方法达到节能目的。
Description
技术领域
本发明涉及一种动静腔流动微沟槽-超疏水耦合减阻设计方法,属于旋转机械动静腔流动减阻领域。
背景技术
动静腔是指旋转圆盘与静止圆盘之间的流动区域,如图1所示,既是一个经典的流体力学研究对象,又是一个广泛存在于石油化工、水利水电、海洋工程、航天动力、核电装备等领域旋转机械中的工程问题。在大多旋转机械中,动静腔是不可避免的形成的,动盘与腔内流体的摩擦功率属于能量损失,减小动盘转动阻力(即减小壁面周向切应力)是实现旋转机械节能降耗的一个重要途径。
仿荷叶超疏水壁面是指液体表观接触角大于150°且滚动角较小的壁面。该壁面上的微纳结构及低表面能阻止液体侵入,从而在壁面贮存了一层极薄的气膜,气液界面速度滑移具有减阻功能且减阻效果优异。超疏水壁面可通过喷涂含聚四氟乙烯、氧化铝纳米线等添加剂的超疏水材料制得。但是,在较大雷诺数或压力下,很难维持长期稳定气膜减阻。
仿鲨鱼皮微沟槽在平板流动中的减阻机理可以总结为:(1)粘性底层和缓冲层的厚度增大,对数律区抬升,近壁区平均速度梯度减小;(2)流向沟槽限制了流向涡的展向运动等,湍流猝发强度、雷诺应力等统计量减弱。由于平板上主流是直线型,因而微沟槽型线设计成直线流向型,不需要设计成曲线型。在动静腔中实施微沟槽减阻策略时,显然不能把微沟槽型线设计成平板流动中的直线型。这是因为动盘壁面对流体的剪切作用、腔体对流体的约束等,动盘近壁区流动方向为周向和径向,其流线是曲线型的。因此,为达到良好的减阻效果,动盘壁面上的沟槽应该设计成曲线型的。但是,该曲线设计方法较缺乏。
通过在动静腔动盘壁面上耦合超疏水减阻技术和微沟槽减阻技术,克服各自不足之处,达到两者相互影响、协同作用下的长效增强减阻目的。基本原理是:加工曲线型微沟槽,喷涂超疏水材料,使得微沟槽减小边界层湍流强度、雷诺应力等并在槽谷产生低速流动,削弱对超疏水壁面贮存气体的剪切作用,实现超疏水壁面长效减阻;壁面贮存气体的滑移、刚度和阻尼作用又反过来减小湍流边界层动能损失,实现微沟槽壁面增强减阻。
经检索,CN201810412446.3《一种多功能复合减阻实验测试装置》给出了一种动静腔减阻测试方法,但不涉及微沟槽-超疏水耦合减阻方法。
发明内容
本发明的目的就是针对上述存在的问题,提供一种动静腔流动微沟槽-超疏水耦合减阻设计方法。
本发明的目的是这样实现的,一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,包括以下步骤:
S1、利用动盘近壁区周向和径向速度方程,即二维流线方程,设计出微沟槽型线;
动盘近壁区流体径向速度:
Cr=1.18·(10-5·Rer+2)-0.49(2)
vr是流体径向速度,径向坐标是流体径向速度的自变量;Cr是动盘径向速度系数;r是径向坐标;υ是流体运动粘度;ω是动盘角速度;K是携卷系数,流体周向速度与动盘周向速度之比;z是轴向坐标;δr动盘边界层厚度;Rer是雷诺数;
选取边界层高度系数和携卷系数K,代入径向坐标、动盘转速、携卷系数、边界层高度系数、流体粘度到式(1)-(4),就可计算出[Rs1,Rs2]区间内的径向速度分布;Rs1为微沟槽前缘半径,Rs2为微沟槽尾缘半径,R为动盘半径;
流体角速度:
选取时间步Δt,选取原则是:使得型线光顺;型线计算公式如下:
rn=rn-1+vr n-1·Δt(7)
r0=Rs1(8)
r是极径;是极角;利用式(5),计算出式(6)中的将rn-1代入式(1)-(4),计算出式(7)中的vr n-1,通过式(6)和(7)的不断累加,直至rn≥R,就可得到[Rs1,R]区间内n+1或n个坐标点 描绘出微沟槽的型线;
S2、确定微沟槽其它设计参数的值;
确定微沟槽横截面形状,选取横截面高度系数z根据S1中的δr和计算得到,h是微沟槽横截面高度,选取横截面形状参数s是微沟槽横截面弧长;横截面高度和弧长均随半径变化,在径向上两者比值保持不变,设计微沟槽数量Z;
S3、喷涂超疏水材料;
根据S1和S2方法,设计出微沟槽全参数,在动盘上加工出微沟槽,然后,选用含聚四氟乙烯、氧化铝纳米线添加剂的超疏水材料喷涂;t为微沟槽内最大喷涂厚度,用来表征喷涂尺度。
携卷系数K的取值范围是0.4378≤K≤1。
微沟槽前缘半径Rs1的取值范围是0≤Rs1≤0.9·R。
微沟槽尾缘半径Rs2与动盘半径R相等。
微沟槽横截面形状为V形、U形、矩形。
通过本发明,提供了一种动静腔流动微沟槽-超疏水耦合减阻设计方法,利用动盘近壁区流动数学模型,构建微沟槽型线控制方程。首先,代入动静腔工作参数,计算微沟槽型线控制方程,得到微沟槽的型线数据。然后,选取微沟槽横截面形状,设计微沟槽横截面高度和弧长,选取微沟槽数量,在动盘上加工出微沟槽。最后,在带微沟槽的动盘上喷涂超疏水材料,保证涂层不淹没微沟槽。本发明通过在动盘上实施微沟槽-超疏水耦合减阻方法达到节能目的。
附图说明
图1为动静腔;
图2为微沟槽;
图3a1为微沟槽横截面形状(V形);
图3a2为喷涂超疏水材料后的微沟槽横截面形状(V形);
图3b1为微沟槽横截面形状(U形);
图3b2为喷涂超疏水材料后的微沟槽横截面形状(U形);
图3c1为微沟槽横截面形状(矩形);
图3c2为喷涂超疏水材料后的微沟槽横截面形状(矩形);
图4为带微沟槽-超疏水耦合减阻的动盘;
图中:1动静腔、2动盘、3静盘、h微沟槽横截面高度、R动盘半径、Rs1微沟槽前缘半径、Rs2微沟槽尾缘半径、s微沟槽横截面弧长、t微沟槽内最大喷涂厚度。
具体实施方式
下面结合具体实施方式对本发明作进一步的详细说明。
本发明提出一种用于动静腔的微沟槽-超疏水耦合减阻设计方法,分为以下步骤:
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是雷诺数。
选取边界层高度系数的值,取值范围是选取携卷系数K,取值范围是0.4378≤K≤1;代入径向坐标、动盘转速、携卷系数、边界层高度系数、流体粘度到式(1)-(4),就可计算出[Rs1,Rs2]区间内的径向速度分布。Rs1为微沟槽前缘半径,单位m,取值范围是0≤Rs1≤0.9·R,R是动盘半径,单位m;Rs2为微沟槽尾缘半径,单位m,与动盘半径R相等,即Rs2=R。
流体角速度:
选取时间步Δt,选取原则是:使得型线光顺。型线计算公式如下。
rn=rn-1+vr n-1·Δt(7)
r0=Rs1(8)
r是极径;是极角。利用式(5),计算出式(6)中的将rn-1代入式(1)-(4),计算出式(7)中的vr n-1。通过式(6)和(7)的不断累加,直至rn≥R,就可得到[Rs1,R]区间内n+1或n个坐标点 描绘出微沟槽的型线。
S2确定微沟槽其它设计参数的值。
微沟槽横截面形状为V形、U形、矩形。为横截面高度系数,取值范围是其中,z根据S1中的δr和计算得到。h是微沟槽横截面高度,单位m。横截面形状参数的取值范围为s是微沟槽横截面弧长,单位m;横截面高度和弧长均随半径变化,在径向上两者比值保持不变。设计微沟槽数量Z,使得微沟槽稠密度为
S3喷涂超疏水材料。
根据S1和S2方法,设计出微沟槽全参数,在动盘上加工出微沟槽。然后,选用含聚四氟乙烯、氧化铝纳米线等添加剂的超疏水材料喷涂。t为微沟槽内最大喷涂厚度,单位m,使得微沟槽增厚系数为即保证涂层不淹没微沟槽。
以一动静腔为例,如图1所示,动盘转速1450r/min,动盘半径160mm,流体运动粘度是8.93×10-7m2/s。
S1利用动盘近壁区周向和径向速度方程,即二维流线方程,设计出微沟槽型线。在MATLAB中编写式(1)-(9)的计算程序,选取边界层高度系数选取携卷系数K=0.5,选取微沟槽前缘半径Rs1=50mm(即Rs1=0.3125·R),微沟槽尾缘半径为Rs2=160mm(Rs2=R),选取时间步Δt=0.001s(使得型线光顺),得到234个坐标点,如表1所示。设计出的微沟槽型线如图2所示。
表1微沟槽型线坐标点
S2确定微沟槽其它设计参数的值。
微沟槽横截面形状选取为矩形,如图3c1、图3c2所示。首先,选取横截面高度系数然后,根据S1中的δr和计算得到z。最后,利用和z,计算出微沟槽横截面高度h,如表2所示。选取横截面形状参数利用和h计算出s,如表2所示。选取微沟槽数量为60,计算得稠密度为
表2微沟槽横截面几何参数
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。
Claims (10)
1.一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,包括以下步骤:
S1、利用动盘近壁区周向和径向速度方程,即二维流线方程,设计出微沟槽型线;
动盘近壁区流体径向速度:
Cr=1.18·(10-5·Rer+2)-0.49 (2)
vr是流体径向速度,径向坐标是流体径向速度的自变量;Cr是动盘径向速度系数;r是径向坐标;υ是流体运动粘度;ω是动盘角速度;K是携卷系数,流体周向速度与动盘周向速度之比;z是轴向坐标;δr动盘边界层厚度;Rer是雷诺数;
选取边界层高度系数和携卷系数K,代入径向坐标、动盘角速度、携卷系数、边界层高度系数、流体运动粘度到式(1)-(4),就可计算出[Rs1,Rs2]区间内的径向速度分布;Rs1为微沟槽前缘半径,Rs2为微沟槽尾缘半径,R为动盘半径;
流体角速度:
选取时间步Δt,选取原则是:使得型线光顺;型线计算公式如下:
r0=Rs1 (8)
是极角;利用式(5),计算出式(6)中的将rn-1代入式(1)-(4),计算出式(7)中的vr n-1,通过式(6)和(7)的不断累加,直至rn≥R,就可得到[Rs1,R]区间内n+1或n个坐标点描绘出微沟槽的型线;
S2、确定微沟槽其它设计参数的值;
确定微沟槽横截面形状,选取横截面高度系数z根据S1中的δr和计算得到,h是微沟槽横截面高度,选取横截面形状参数s是微沟槽横截面弧长;横截面高度和弧长均随半径变化,在径向上两者比值保持不变,设计微沟槽数量Z;
S3、喷涂超疏水材料;
根据S1和S2方法,设计出微沟槽全参数,在动盘上加工出微沟槽,然后,选用含聚四氟乙烯、氧化铝纳米线添加剂的超疏水材料喷涂;t为微沟槽内最大喷涂厚度,用来表征喷涂尺度。
3.根据权利要求1所述的一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,携卷系数K的取值范围是0.4378≤K≤1。
4.根据权利要求1所述的一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,微沟槽前缘半径Rs1的取值范围是0≤Rs1≤0.9·R。
5.根据权利要求1所述的一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,微沟槽尾缘半径Rs2与动盘半径R相等。
6.根据权利要求1所述的一种动静腔流动微沟槽-超疏水耦合减阻设计方法,其特征在于,微沟槽横截面形状为V形、U形、矩形。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010610401.4A CN111737836B (zh) | 2020-06-30 | 2020-06-30 | 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010610401.4A CN111737836B (zh) | 2020-06-30 | 2020-06-30 | 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111737836A CN111737836A (zh) | 2020-10-02 |
| CN111737836B true CN111737836B (zh) | 2021-06-22 |
Family
ID=72653610
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010610401.4A Active CN111737836B (zh) | 2020-06-30 | 2020-06-30 | 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111737836B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4313978A (en) * | 1978-12-20 | 1982-02-02 | Minnesota Mining And Manufacturing Company | Antistatic compositions and treatment |
| US7935734B2 (en) * | 2005-07-08 | 2011-05-03 | Anna Lee Tonkovich | Catalytic reaction process using microchannel technology |
| CN103953580A (zh) * | 2014-05-27 | 2014-07-30 | 山东理工大学 | 中大型转子叶片上均匀渐变型高逼真仿鲨鱼沟槽微结构的设计方法 |
| CN108997888A (zh) * | 2018-08-03 | 2018-12-14 | 西北工业大学 | 一种超疏水减阻涂层的制备方法 |
| CN110484151A (zh) * | 2019-08-29 | 2019-11-22 | 浙江大学 | 微沟槽减阻柔性薄膜及其制备方法 |
| CN111723447A (zh) * | 2020-06-30 | 2020-09-29 | 扬州大学 | 一种动静腔流动微沟槽减阻设计方法 |
-
2020
- 2020-06-30 CN CN202010610401.4A patent/CN111737836B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4313978A (en) * | 1978-12-20 | 1982-02-02 | Minnesota Mining And Manufacturing Company | Antistatic compositions and treatment |
| US7935734B2 (en) * | 2005-07-08 | 2011-05-03 | Anna Lee Tonkovich | Catalytic reaction process using microchannel technology |
| CN103953580A (zh) * | 2014-05-27 | 2014-07-30 | 山东理工大学 | 中大型转子叶片上均匀渐变型高逼真仿鲨鱼沟槽微结构的设计方法 |
| CN108997888A (zh) * | 2018-08-03 | 2018-12-14 | 西北工业大学 | 一种超疏水减阻涂层的制备方法 |
| CN110484151A (zh) * | 2019-08-29 | 2019-11-22 | 浙江大学 | 微沟槽减阻柔性薄膜及其制备方法 |
| CN111723447A (zh) * | 2020-06-30 | 2020-09-29 | 扬州大学 | 一种动静腔流动微沟槽减阻设计方法 |
Non-Patent Citations (1)
| Title |
|---|
| 表面活性剂溶液与壁面纵向微沟槽协同减阻研究;魏进家等;《化工学报》;20180131;第69卷(第1期);第472-482页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111737836A (zh) | 2020-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111460699B (zh) | 平壁表面减阻功能微织构的设计方法 | |
| CN102588600B (zh) | 具有回流功能的型槽端面非接触式机械密封 | |
| CN104632893A (zh) | 基于疏水界面的多孔水静压轴承结构 | |
| CN108358154A (zh) | 一种具备减阻特性的梭形仿生微结构 | |
| Monfared et al. | Experimental study on the friction drag reduction of superhydrophobic surfaces in closed channel flow | |
| CN103016743B (zh) | 一种可实现零磨损和零泄漏且防固体颗粒的机械密封 | |
| CN111723447B (zh) | 一种动静腔流动微沟槽减阻设计方法 | |
| CN102266933A (zh) | 一种铝合金熔液变流量电磁低压充型方法 | |
| CN111737836B (zh) | 一种动静腔流动微沟槽-超疏水耦合减阻设计方法 | |
| Sasaki | Surface texturing for friction control: A review on existing technology and prospects | |
| Liu et al. | Design and optimization of bionic Janus blade in hydraulic torque converter for drag reduction | |
| Zhu et al. | Influence of surface slip on hydrodynamics and flow field around a two-dimensional hydrofoil at a moderate Reynolds number | |
| CN117533453A (zh) | 一种聚合物减阻剂注入式水翼 | |
| CN108999846B (zh) | 一种超疏水减阻肋条结构 | |
| Ma et al. | Drag reduction characteristics of bionic mucous membrane acting on the turbulent boundary layer | |
| CN205592513U (zh) | 一种应用于管道上的微沟槽及具有该结构的管道 | |
| CN116926483B (zh) | 提高海水液压泵奥氏体不锈钢摩擦配副表面性能的方法 | |
| Yamagishi et al. | Numerical simulation of flow around a circular cylinder with curved sectional grooves | |
| Tian et al. | Study on the hydrodynamic characteristics of an outboard engine propeller with hydrophobic coating | |
| CN215861731U (zh) | 水滴型微凹坑织构端面机械密封结构 | |
| Sekhar | Heat Transfer Analysis of Second Grade Fluid Over a Stretching Sheet Through Porous Medium Under the Influence of Chemical Reaction Parameter | |
| CN109236727A (zh) | 一种基于仿生学的轴流泵叶片 | |
| CN116484765A (zh) | 一种用于摩擦润滑与磨屑调控的设备及优化方法 | |
| CN115374721A (zh) | 一种凹坑微结构的减阻参数优化方法、系统及存储介质 | |
| Chen et al. | Numerical simulation of the lubrication characteristics of bionic surface textures on mechanical sealing rings |
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 |