粗糙度对大间隙环流偏心转子动特性系数的影响

基于作者建立的大间隙环流中转子运动的理论模型,用摄动法推导了大间隙环流流场非线性控制方程的一阶摄动方程,采用数值方法研究了静子和转子壁面粗糙度对大间隙环流中偏心转子动特性系数的影响。研究结果表明：静子和转子壁面粗糙度对大间隙环流中偏心转子动特性系数有较大影响。所得到的数值结果与已有的解析解和实验结果具有较好的一致性。     

大间隙环流中偏心转子动特性系数的数值分析方法

基于作者用整体流动理论和Moody壁面摩擦系数方程建立的大间隙环流中转子动特性系数数值计算模型,应用摄动方法推导了大间隙环流流场非线性控制方程组的一阶摄动方程,提出了求解大间隙环流中偏心转子动力学特性系数的数值分析方法。用该方法得到的数值结果与已有的解析解和实验结果具有较好的一致性。     

大间隙环流壁面摩擦及偏心转子静特性研究

基于紊流整体流动模型和Ｍｏｏｄｙ壁面摩擦系数方程简化了大间隙环流三维流非线性偏微分控制方程组,推导了零阶摄动方程并采用数值方法对大间环流中偏心转子静特性进行了深入研究。实例计算结果表明：大间隙环流中转子与以灿承和密封为代表的小间隙环流中转子的静特性有很大同;壁面摩擦系数沿周向变化,同时转子和静子壁面粗糙度以及偏心率对大间隙环流中转子静特性有较大的影响。     

环状密封转子动力系数的数值值计算

本文以Ｈｉｒｓ紊流润滑方程为基础，应用摄动分析法及数值计算获得了有限长环状密封的转子动力系数，模型中考虑了入口预旋的影响，通过与实验及用有限差分法获得的结果比较，证明本文提出的方法具有一定的实用性。     

环状密封转子动力系数的数值计算

本文以Ｈｉｒｓ紊流润滑方程为基础，应用摄动分析法及数值计算获得了有限长环状密封的转子动力系数，模型中考虑了人口预旋的影响。通过与实验及用有限差分法获得的结果比较，证明本文提出的方法具有一定的实用性。     

惯性效应对超高速倾斜端面气膜密封稳动态特性影响

研究了惯性效应和端面倾斜对超高速气膜端面密封稳动态特性的影响. 考虑气体惯性效应,建立了气膜端面密封稳动态特性数值分析模型,采用有限差分法求解稳态和微扰雷诺方程,获得端面膜压分布. 数值分析了惯性效应和端面倾斜度对开启力、气膜刚度和泄漏率等稳态性能参数以及刚度系数和阻尼系数等动态特性系数的影响规律,并以获得较大刚度系数为目标,获得了螺旋槽关键几何参数的优选值范围. 结果表明:在超高速条件下,考虑惯性效应后的干气密封泄漏率显著减小,刚漏比明显增大,而开启力、气膜刚度和动特性系数变化不大;倾斜端面气膜密封相较于平行端面气膜密封具有更佳的低频刚度和高频阻尼.      

对流扩散方程的摄动有限体积(PFV)方法及讨论

在有限体积(FV)方法的重构近似中,引入数值摄动处理,即把界面数值通量摄动展开成网格间距的幂级数,并利用积分方程自身的性质求出幂级数的系数,同时获得高精度迎风和中心型摄动有限体积(PFV)格式.对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型PFV格式,这些PFV格式的结构形式及使用基点数与一阶迎风格式完全一致,迎风PFV格式满足对流有界准则;二阶和四阶中心PFV格式对网格Peclet数的任意值均为正型格式,比常用的二阶中心格式优越.用一维标量输运和方腔流动算例说明PFV格式的优良性能,并把PFV方法与性质相近的摄动有限差分(PFD)方法及相关的高精度方法作了对比分析.     

蜂窝密封动力特性计算方法的改进及应用

针对蜂窝密封的转子动力学特性分析,本文由双控制体模型入手,建立相应的控制方程,将其无量纲化,并通过摄动分析分解为零阶及一阶无量纲方程。通过零阶方程求得定常流场,确定流场是否发生阻塞现象,继而求解一阶方程得到密封的动力特性系数。本文的双控制体模型在蜂窝密封静子和转子表面采用Blasius模型的摩擦系数,且舍弃转子轴心传统的圆轨迹小扰动假设,直接设定为椭圆轨迹小扰动,加之矩阵形式的表达,得到了更为清晰的模型及更为简洁的算法。随后与实验结果相比较,验证模型及计算结果的准确性,并比较一般气体密封的计算结果,体现蜂窝密封在稳定性方面的优越性。最后分析了几种密封参数对其动力特性的影响,揭示某些规律并给出一些建议。     

分数阶拟周期Mathieu方程的动力学分析

分数阶微积分有着诸多优异的特点, 目前在动力学领域主要用来提高非线性系统振动特性研究的准确性. 本文在拟周期Mathieu方程的基础上, 引入分数阶微积分理论, 研究了分数阶微分项参数对方程稳定性的影响. 首先, 采用摄动法得到方程稳定区和非稳定区分界线(即过渡曲线)近似表达式, 利用数值方法验证了解析结果的准确性, 图像显示两者吻合较好. 随后, 通过归纳总结不同情况下的过渡曲线近似表达式, 发现在系统中分数阶微分项以等效线性刚度和等效线性阻尼的方式存在. 根据这一特点, 得到了系统等效线性阻尼和等效线性刚度的一般形式, 并且定义了非稳定区域厚度. 最后, 通过数值仿真直观地分析了分数阶微分项参数对方程稳定区域大小和过渡曲线位置的影响. 结果发现, 分数阶微分项不仅具有阻尼特性还具有刚度特性, 并且以等效线性刚度和等效线性阻尼的方式影响着方程稳定区域大小和过渡曲线位置. 合理选择分数阶微分项参数可以使其呈现不同程度的刚度特性或阻尼特性, 方程稳定区域的大小和过渡曲线的位置也因此产生了不同程度的变化.      

局部裂纹损伤简支梁的曲率模态特性

将裂缝损伤简化成矩形凹槽,采用delta函数表示简支梁的裂纹损伤位置,得到了全梁范围内截面转动惯量和单位长度质量的表达式,建立了局部裂缝损伤简支梁的横向自由振动方程.利用摄动方法给出了裂纹摄动项的一般表达式,根据摄动项和完整梁都同时满足边界条件的特点,将一阶和二阶摄动项都表示成完整梁模态的线性组合,结合delta函数的性质,最终获得了受损简支梁的特征值和模态振型的解析表达式.最后,通过数值计算得到结构模态参数,对比了一阶摄动和二阶摄动对计算结果的影响,分析了不同阶固有频率和模态曲率的变动量,为简支梁的损伤监控和检测提供了理论依据.     

轻载径向滑动轴承中Taylor涡动的产生和影响研究

本文用原始变量法直接求解了三维的Ｎ－Ｓ方程，计算分析了高速旋转有限长圆柱轴承中油膜层流失稳产生涡动的临界Ｔａｙｌｏｒ数及流场、压力场和摩擦阻力。轴承端部泄油量等的变化。结果表明，在有限长同心圆柱轴承中，随着轴旋转速度的提高，轴承磨擦阻力线性增大，油膜层流失稳出现的涡动增加轴承摩擦阻力并减少轴承端部泄油量，油膜层流失稳后，轴承长度方向均匀地排列着一些流体涡，涡动的强度从轴承中间截面向轴承端部逐渐减弱     

Performance model of metallic concentric tube recuperator with counter flow arrangement

A performance model for counter flow arrangement in concentric tube recuperator that can be used to utilize the waste heat in the temperature range of 900–1,400°C is presented. The arrangement consists of metallic tubular inner and outer concentric shell with a small annular gap between two concentric shells. Flue gases pass through the inner shell while air passes through the annular gap in the reverse direction (counter flow arrangement). The height of the recuperator is divided into elements and an energy balance is performed on each elemental height. Results give necessary information about surface, gas and air temperature distribution, and the influence of operating conditions on recuperator performance. The recuperative effectiveness is found to be increased with increasing inlet gas temperature and decreased with increasing fuel flow rate. The present model accounts for all heat transfer processes pertinent to a counterflow radiation recuperator and provide a valuable tool for performance considerations.     

Entrance flow in eccentric annular ducts

A control-volume-based solution of the complete set of Navier-Stokes equations for the laminar, three-dimensional developing flow in straight, eccentric, cylindrical annular ducts is described. Numerical results for velocity and pressure development, pressure defect and entrance lengths are presented for a wide range of duct parameters, i.e. relative eccentricity ? and radius ratio γ. The present results match very well with earlier numerical solutions for the limiting cases of developing flow in concentric ducts and fully developed flow in eccentric ducts. Comparison with earlier approximate results for developing flow in eccentric ducts indicates that the approximate model predicts the velocity and pressure development with an error of about 10%. However, the development length predicted by the approximate model is grossly in error. The pressure defect and development length in eccentric ducts are very high compared with their counterparts in concentric ducts. The pressure defect, development length and maximum velocity increase with the radius ratio for eccentric ducts, while the reverse is true for concentric ducts. Also, the apparent friction factor decreases as the eccentricity increases.     

Computation of rotordynamic coefficients associated with leakage steam flow through labyrinth seal

A mathematical model of calculating rotordynamic coefficients associated with leakage steam flow through labyrinth seals was presented. Particular attention was given to incorporating thermal properties of the steam fluid into prediction of leakage flow and subsequent derivation of rotordynamic coefficients, which quantitatively characterize influence of aerodynamic forcing of the leakage steam flow on the rotordynamics. By using perturbation analysis, we determined periodic and analytic solutions of the continuity and circumferential momentum equations for the time-dependent flow induced by non-axisymmetric rotation of the rotor encompassed by a labyrinth seal. Pressure distributions along labyrinth seal cavities and rotordynamic coefficients were compared at the same condition for air and steam flows. Influence of steam flow through the labyrinth seal cavities on rotordynamic coefficients was analyzed in terms of inlet pressure, inlet swirl velocity and rotor speed.     

Pore-Scale Modeling of Viscous Flow and Induced Forces in Dense Sphere Packings

We propose a method for effectively upscaling incompressible viscous flow in large random polydispersed sphere packings: the emphasis of this method is on the determination of the forces applied on the solid particles by the fluid. Pore bodies and their connections are defined locally through a regular Delaunay triangulation of the packings. Viscous flow equations are upscaled at the pore level, and approximated with a finite volume numerical scheme. We compare numerical simulations of the proposed method to detailed finite element simulations of the Stokes equations for assemblies of 8–200 spheres. A good agreement is found both in terms of forces exerted on the solid particles and effective permeability coefficients.     

Velocity distributions in Taylor vortex flow with imposed laminar axial flow and isothermal surface heat transfer

Fully-developed flow in a concentric annulus formed by a stationary outer cylinder, which may be heated isothermally, and a rotatable inner cylinder has been studied experimentally by means of hot-wire anemometry techniques. Velocity profiles for the axial, tangential, and radial directions of flow have been obtained for adiabatic conditions and for wide and narrow annular gaps.It has been shown that the onset of Taylor vortex flow has a pronounced effect on the velocity profiles for all three directions. However, while the profiles for the axial and tangential directions are explicable, those for the radial direction are not so, at present. Also, it was found that heat transfer through the outer annular surface had a greater effect on the radial velocity profile than on the axial or tangential, but in the narrow gap case only.     

STABILITY ISSUES OF CONCENTRIC PIPES CONTAINING STEADY AND PULSATILE FLOWS

In this paper, we study both the static and dynamic instabilities of submerged and inclined concentric pipes conveying fluid. The governing equation for the inner tubular beam is derived under small deformation assumptions. We obtain the discretized dynamical equations using spatial finite-difference schemes. In the case of steady flow, both buckling and flutter instabilities are investigated. In the case of pulsatile flow, we compute the eigenvalues of the monodromy matrix derived from the discretized linear system with periodic coefficients, and deduce the dynamical stability information. In addition, for a special case, in which the concentric pipes have the same length, we compare the dynamic stability results with the corresponding solutions obtained with the Bolotin method.