基于湍流模型的高速螺旋槽机械密封稳态性能研究

针对高速工况下的液膜润滑螺旋槽端面机械密封，建立了其湍流润滑模型，采用有限单元法结合松弛迭代技术实现了润滑方程和液膜湍流模型的数值求解，对比分析了层流模型和湍流模型下不同螺旋槽几何参数和工况参数对密封性能的影响. 结果表明:液膜湍流效应显著提升了螺旋槽机械密封端面液膜的动力润滑效应，密封的开启力、泄漏率和刚度明显大于层流模型预测值. 在不同条件下，比较而言螺旋槽内产生更加明显的湍流效应，其内液膜流动行为远不同于层流模型. 以开启力为优化目标，湍流模型获得的优化螺旋槽几何参数在螺旋角、槽深明显不同于层流模型. 在高速和低黏度介质下，机械密封的湍流效应不可忽略.      

微纳米结构壁面对流场及动压润滑的影响研究

固体边界具有的微纳米结构将影响流体在近壁面处的流动行为,进而由于尺度效应改变流体在整个微间隙的流动或润滑规律.将壁面可渗透微纳米结构等效为多孔介质薄膜,采用Brinkman方程来描述流体在近壁面边界渗透层内的流动,并将其与自由流动区域的不可压缩流体Navier-Stokes控制方程耦合,在界面处的连续边界条件下求解和分析了速度分布规律和压力变化规律.针对恒定法向承载力的油膜润滑条件,进一步讨论了静止表面或运动表面的微纳米结构对近壁面流动行为的影响;并揭示了考虑壁面微纳米结构的流体动压润滑的油膜厚度和摩擦系数的变化规律.论文结果为具有可渗透微结构表面的微间隙流动与润滑提供了理论参考.     

强吸气旋转圆筒壁面湍流边界层建模及计算

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.      

MODELLING AND CALCULATION OF THE TURBULENT BOUNDARY LAYER ON A ROTATING CYLINDER SURFACE WITH STRONG SUCTION 1)

提出了湍流边界层的一种简单、快速计算方法, 用以求解强吸气作用下旋转圆筒表面边界层流动. 首先, 理论分析了同心圆筒间的旋转流体运动, 外筒静止、内筒旋转且为多孔吸气条件. 强吸气情况下旋转流动主要表现为内筒壁面附近的边界层流动, 基于这一事实得到了周向速度分布的解析表达式. 其次, 通过引入新参数扩展Cebeci-Smith代数湍流模型, 使其能考虑流线曲率、壁面吸气、低Reynolds数效应等因素. 针对这些因素的综合影响, 采用解析修正和经验参数对模型进行调整. 同时, 基于Reynolds应力湍流模型的仿真结果, 校准代数湍流模型中的经验参数. 最后, 给出基于广义Cebeci-Smith湍流模型的旋转壁面边界层流动的迭代算法, 该算法适用于需要特殊迭代过程的轴向及周向流动均匀情况. 计算了不同旋转速度和吸气强度组合工况下的边界层流动, 其周向速度和湍流强度分布与基于Reynolds应力湍流模型的计算结果非常接近. 并且表明, 当Reynolds应力湍流模型数值模拟预测内筒边界层为稳定层流时, 该方法也再现了相同初始条件下的层流边界层.     

Reynolds方程在纹理表面动压润滑计算中的有效性评价

从理论上阐述了纹理表面动压润滑计算中决定Reynolds方程有效性的两个关键因素为油膜厚度与纹理特征长度的比值h/L和缩减的雷诺数re;只有当h/L和re同时趋近于零时Reynolds方程才能够获得准确的结果,并由此在h/L-Re平面上标注了Reynolds方程的适用范围.继而以二维矩形沟槽为实例,采用数值方法计算了h/L和re对Reynolds方程误差的影响规律;分析了Reynolds方程在不同条件下的失效机制;分析了矩形沟槽纹理表面Reynolds方程有效性的评价标准:当缩减的雷诺数re小于0.20,并且h/L小于0.015时能够保证Reynolds方程的误差在10%以下.     

基于润滑特性仿真的燃油泵滑动轴承优化设计

磨损加剧是低介质黏度和自冷却结构燃油泵滑动轴承主要失效机制,为提高轴承使用寿命,论文提出了一种基于滑动轴承润滑特性分布规律的轴承优化设计方法.首先基于油膜动压润滑流动的Reynolds方程和等效黏度润滑流动模型,以绝热流动为假设简化滑动轴承内部流动的能量方程,构建一种联合Reynolds方程和绝热流动能量积分方程的燃油泵滑动轴承热流润滑模型.其次采用CFD数值模拟和有限差分法相结合的混合仿真方法,分别对不同间隙比、偏心率、宽径比条件下的滑动轴承的润滑特性进行了仿真,最后采用基于遗传算法的间接法优化方法进行滑动轴承结构优化设计.研究结果表明:优化后的轴承偏心率为0.822 4,宽径比为1.2,油膜厚度为2.7μm,平均温升为36℃,滑动轴承油膜承载力试验测试值与仿真计算值间的误差最大不超过5%,表明轴承能够很好地承受燃油泵的径向力载荷,从而保证轴承与轴瓦不发生摩擦接触并具有良好的润滑性能.     

材料表面润湿性调控及减阻性能研究

设计合成不同结构的自组装分子,使其可以在不改变表面粗糙度的情况下改变表面的润湿性能;利用低表面能涂层修饰粗糙表面得到超疏水表面.采用流变仪和水洞试验分别在层流和湍流流动状态下测试了具有不同润湿行为的亲、疏水材料的减阻性能.结果表明:在层流流动状态,随着不同表面的接触角从13°增加到45°、113°和161°,减阻率随之从1.8%增大到7.2%、7.9%和14.9%;在湍流流动状态下,自组装涂层接触角为13°、45°和113°的三组模型的平均减阻率为0.8%、1.9%和6.8%,最大减阻率分别可达3.6%、9.2%和18.0%.两种流体流动中均存在材料表面水接触角增加减阻效率增大的行为.     

基于稀薄效应的微气体径向轴承稳态性能

针对微气体轴承给出参考努森数的定义,根据空气不同温度时的黏度,得到参考努森数的分布范围;考虑气体稀薄效应,给出基于Burgdorfer一阶滑移速度边界的微气体径向轴承润滑Reynolds方程的修正形式; 采用有限差分法求解修正的Reynolds方程,得到不同参考努森数$Kn_0$, 轴承数以及轴颈偏心率情况下轴承的压力分布、无量纲承载能力及偏位角. 数值分析表明:随气体稀薄程度的增强,气体径向轴承的压力明显降低,无量纲承载力降低,而轴承偏位角增大. 当偏心率小于0.6时,轴承偏位角变化平缓,受$Kn_0$数的影响不明显. 当轴承数较小时,气体稀薄程度对轴承的无量纲承载力、偏位角影响较小.      

低雷诺数翼型蒙皮主动振动气动特性及流场结构数值研究

针对低雷诺数(Re)翼型气动性能差的特点,文章通过对翼型柔性蒙皮施加主动振动的方法,提高翼型低Re下的气动特性,改善其流场结构.采用带预处理技术的Roe方法求解非定常可压缩Navier-Stokes方程,对NACA4415翼型低Re流动展开数值模拟.通过时均化和非定常方法对比柔性蒙皮固定和振动两种状态下的升阻力气动特性和层流分离流动结构.初步研究工作表明在低Re下柔性蒙皮采用合适的振幅和频率,时均化升阻力特性显著提高,分离泡结构由后缘层流分离泡转变为近似的经典长层流分离泡,分离点后移,分离区缩小.在此基础上,文章更加细致研究了柔性蒙皮两种状态下单周期内的层流分离结构及壁面压力系数分布非定常特性和演化规律.蒙皮固定状态下分离区前部流场结构和压力分布基本保持稳定,表现为近似定常分离,仅在后缘位置出现类似于卡门涡街的非定常流动现象.柔性蒙皮振动时从分离点附近开始便产生分离涡,并不断向下游移动、脱落,表现为非定常分离并出现大范围的压力脉动.蒙皮振动使流体更加靠近壁面运动,大尺度的层流分离现象得到有效抑制.      

限量供油对速度交叉工况下润滑特性的试验研究

速度交叉效应对摩擦副表面润滑剂迁移及润滑状态改善具有重要意义.利用表面速度异向光干涉润滑油膜测量装置，在限量供油条件下，改变两固体表面卷吸速度与滑动速度夹角ε，研究了该角度变化对润滑油膜及油池边界的影响.结果表明，卷吸速度与滑动速度呈一定夹角时，部分润滑剂再分配至接触区参与成膜；由于润滑速度交叉行为导致润滑剂横向运输效果增强，限量供油条件下的润滑状态改善.此外，低供油量(0.1μl)限制了速度交叉效应的润滑增效性，产生的低油池边界减弱了两侧润滑剂的有效交汇.     

Linear stability analysis of flow in a periodically grooved channel

We have conducted the linear stability analysis of flow in a channel with periodically grooved parts by using the spectral element method. The channel is composed of parallel plates with rectangular grooves on one side in a streamwise direction. The flow field is assumed to be two‐dimensional and fully developed. At a relatively small Reynolds number, the flow is in a steady‐state, whereas a self‐sustained oscillatory flow occurs at a critical Reynolds number as a result of Hopf bifurcation due to an oscillatory instability mode. In order to evaluate the critical Reynolds number, the linear stability theory is applied to the complex laminar flow in the periodically grooved channel by constituting the generalized eigenvalue problem of matrix form using a penalty‐function method. The critical Reynolds number can be determined by the sign of a linear growth rate of the eigenvalues. It is found that the bifurcation occurs due to the oscillatory instability mode which has a period two times as long as the channel period. Copyright © 2003 John Wiley & Sons, Ltd.     

Near wall characteristics of flow over grooved circular cylinder

An experimental investigation was made of the flow over a grooved circular cylinder of different aspect ratios. Based on the drag coefficient, Strouhal number and mean shear stress, three flow regimes of subcritical, critical and supercritical were found, all of which are below the subcritical Reynolds number of a smooth cylinder. Boundary layer characteristics within the different flow regimes were measured. The shift of the boundary layer away from the grooves and an estimated change in the virtual origin are used to establish the similarity of the flow characteristics of grooved and smooth cylinders.Now at Development and Planning Division, Hong Kong Electric Co. Ltd., Hong Kong     

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一, 对缓解日益严峻的能源危机极为重要. 在重力式管道实验系统中, 测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态, 对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力, 沟槽结构采用机械方法加工, 并在表面喷涂超疏水涂层. 结果表明, 持续通气能解决超疏水沟槽表面气膜层流失问题, 实现气膜层长时间稳定维持; 恒定雷诺数下, 随通气速率增大, 超疏水沟槽表面气膜铺展更趋均匀, 减阻率上升; 由于通气速率影响气膜横向扩展能力, 致使恒定通气速率下, 减阻率随雷诺数的变化呈现两种模式; 在固定雷诺数及通气速率时, 减阻率随沟槽尺寸的扩大先增后减, $S^{+}\approx 76$时减阻率最大. 分析其原因在于, 沟槽结构增大沾湿面积的同时, 显著提升了通气状态下超疏水表面气膜层的稳定性, 因而展示出与超疏水表面和沟槽表面均不相同的减阻规律, 且效果更佳.      

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.     

中低速下干气密封型槽仿生改进研究

为解决干式气体端面密封在中低速或启动停车阶段密封性能不稳定的问题,在现有螺旋槽端面型槽结构的基础上借鉴仿生学设计方法提出了仿生多流通道螺旋槽干气密封端面结构.基于气体润滑理论模型,采用数值分析方法计算了气体端面气膜压力控制方程-雷诺方程,研究对比了普通螺旋槽和多流通道槽的端面压力分布.在改变操作参数和端面气膜厚度的条件下分析对比了普通螺旋槽和多流通道槽两种干气密封的开启力,泄漏率,刚度和刚漏比的变化规律,并研究了槽的通道数对密封性能的影响规律.结果表明:中低速下多流通道细长螺旋槽具有高的稳定性和好的密封性,汇流槽在增加密封稳定性的同时可降低泄漏率.     

Identification of a laminar-turbulent interface in partially turbulent flow

Viscous incompressible fluid flow along the clearance between two parallel eccentric circular cylinders is studied numerically. The geometric parameters and the Reynolds number are taken so that the turbulent flow regime takes place in a part of flow, namely, in the zone where the clearance is wide, while the low is laminar in the narrow clearance zone. Criteria using which these zones can be distinguished are discussed.     

液膜润滑机械端面密封流体动压沟槽侧壁效应研究

针对非接触机械密封端面开槽后所出现的膜厚不连续处存在的侧壁效应,在沟槽边界处将广义伯努利方程引入传统润滑方程,建立了考虑动压沟槽侧壁效应的液膜润滑螺旋槽端面机械密封数值分析模型. 采用有限单元法结合拉格朗日乘子法求解润滑方程,研究了不同螺旋槽几何参数和工况条件下沟槽侧壁效应对密封性能的影响. 结果表明:数值模型可方便捕捉沟槽边界处的压力跃变,侧壁效应在不同螺旋槽深度下表现出截然不同的影响规律,高转速、大螺旋角和小密封间隙下动压沟槽的侧壁效应较为显著. 理论模型和计算方法可为超高速工况螺旋槽机械密封的设计和局部惯性效应的研究提供指导.      

Effect of the pattern of initial perturbations on the steady pipe flow regime

The influence of the inlet flow formation mode on the steady flow regime in a circular pipe has been investigated experimentally. For a given inlet flow formation mode the Reynolds number Re* at which the transition from laminar to turbulent steady flow occurred was determined. With decrease in the Reynolds number the difference between the resistance coefficients for laminar and turbulent flows decreases. At a Reynolds number approximately equal to 1000 the resistance coefficients calculated from the Hagen-Poiseuille formula for laminar steady flow and from the Prandtl formula for turbulent steady flow are equal. Therefore, we may assume that at Re > 1000 steady pipe flow can only be laminar and in this case it is meaningless to speak of a transition from one steady pipe flow regime to the other. The previously published results [1–9] show that the Reynolds number at which laminar goes over into turbulent steady flow decreases with increase in the intensity of the inlet pulsations. However, at the highest inlet pulsation intensities realized experimentally, turbulent flow was observed only at Reynolds numbers higher than a certain value, which in different experiments varied over the range 1900–2320 [10]. In spite of this scatter, it has been assumed that in the experiments a so-called lower critical Reynolds number was determined, such that at higher Reynolds numbers turbulent flow can be observed and at lower Reynolds numbers for any inlet perturbations only steady laminar flow can be realized. In contrast to the lower critical Reynolds number, the Re* values obtained in the present study, were determined for given (not arbitrary) inlet flow formation modes. In this study, it is experimentally shown that the Re* values depend not only on the pipe inlet pulsation intensity but also on the pulsation flow pattern. This result suggests that in the previous experiments the Re* values were determined and that their scatter is related with the different pulsation flow patterns at the pipe inlet. The experimental data so far obtained are insufficient either to determine the lower critical Reynolds number or even to assert that this number exists for a pipe at all.     

Effects of riblets upon flow stability

The onset of transition in the boundary layer over a grooved surface is studied by introducing the presence of grooves into the standarde ⁹ model with the aid of a previously obtained equivalent boundary condition. Under conditions of self-similarity and smallness of the grooves, Tollmien-Schlichting waves are found to be excited at a slightly smaller Reynolds number and Görtler vortices at a slightly larger one than on a smooth surface.     

On the Lower Critical Reynolds Number for Flow in a Circular Pipe

Flow in a circular pipe is investigated experimentally at Reynolds numbers higher than that at which the resistance coefficients calculated from the Blasius formula for laminar flow and from the Prandtl formula for turbulent flow are equal. The corresponding Reynolds number based on the mean-flow velocity and the pipe diameter is about 1000. The experiments were performed at a high level of inlet pulsations produced by feeding gas into the pipe through a hole with a diameter several times smaller than the pipe diameter. In our experiments the critical Reynolds number was determined as the value, independent of the distance from the inlet, at which the ratio of the axial to the mean-flow velocity as a function of the Reynolds number deviated from 2. At the maximum ratio of the pipe cross-sectional area to the area of the hole through which the gas entered the pipe, equal to 26, the critical Reynolds number was about 2300. After a fivefold increase in the hole area the critical Reynolds number increased by approximately 4%.At Reynolds numbers below 2000, after at a high level of the inlet pulsations an almost laminar flow had developed in the pipe, a perturbation was introduced by inserting a diametrically oriented cylindrical rod with a diameter 10–20 times smaller than the pipe diameter. In these experiments, at Reynolds numbers higher than 1000, at a distance from the rod equal to 50 pipe diameters the axial to mean-flow velocity ratio was less than 2, approaching this value again at large distances from the rod. The insertion of the rod led to a decrease in the critical Reynolds number by approximately 12%.