涡轮叶片外换热系数计算方法和比较

从工程实用观点出发，引入考虑主流紊流度影响的经验关系式，对McNally程序的前驻点、层流区计算方法和转捩模型进行了改进；引入“紊流度粘性”和间歇因子的概念，并综合考虑压力梯度、紊流度和进、出口雷诺数等因素对转捩的影响。对STAN5程序的层流区和过渡区计算方法进行了改进，采用相似解方法计算前驻点区的流动和换热，作为STAN5程序的初始条件。用改进的程序，对C3x、MKⅡ、Turner、Daniels和VKI叶型，在多种工况下的二维外边界层流动和换热进行了计算研究，分析和比较了程序中积分法和差分法的特点，总结出一套适用于工程设计计算的方法和程序。     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

转动园盘上可压缩流动换热的计算研究

本文将边界层方程无因次化,用“BOX”算法求解。研究计算了转动园盘上层流区、过渡区和紊流区的扭矩系数和局部换热系数以及变物性的影响,并研究了适合计算转盘问题的双层紊流模型。     

一种跨音速叶栅流场分析系统

本文建立了一种工程实用的跨音速叶栅流场分析系统,可进行跨音速叶栅无粘绕流和无粘流-边界层迭代计算。叶栅无粘绕流的计算以Denton拟流向“相对差分”方法为基础,y向采用松驰线性插值算法,所形成的解法能保证以较少的网点获得较准确的结果。边界层的计算采用积分方法,并考虑了紊流边界层的分离问题,文中还介绍了无粘流-边界层的耦合方式。应用本文方法对国内外多种压气机和涡轮叶栅进行了计算分析,所得结果与实验的吻合程度较好,因此,所提出的叶栅流场分析系统,可望为工程设计上实用的计算手段。     

拟合带间歇性的人工风速序列

研究生成既满足给定的工程谱特性，又具有间歇性结构的风速序列作为结构气动计算的输入。生成大气边界层紊流的人工风速序列是做结构风振计算的重要步骤，目前的方法都是生成高斯序列；实际的大气边界层紊流具有多尺度多层次的非高斯结构，集中体现在所谓“间歇性”上。利用余振苏与Leveque提出的对数泊松律模型产生各尺度的小波系数，可以在拟合目标风速谱的同时引入湍流间歇性，通过小波逆变换生成单条风速序列，再用拟合周期图的思路扩展为多重相关序列。     

《航空学报》1985年(第六卷)分类索引

流体力学旋涡运动研究的某些进展北京航空学院刘谋信苏文翰(P .1)用一种新的极点配置方程确定双输入多输出系统的控制率南京航空学院卢叔全(P .13)可压紊流边界层掺混方法力学模型的改进中国气动研究与发展中心包涵令(P .188)分离流动的某些进展中国气动研究与发展中心张涵信(P     

中心进气转静盘腔一维流动模型的改进

采用数值模拟方法,对动盘面边界层夹带流量进行分析,以改进转静盘腔一维流动模型的计算精度.结果表明:边界层速度分布的1/7幂律假设低估了自由盘边界层夹带流量,从而影响了一维流动模型对源区边界的计算精度;传统的关联式未能有效地归纳不同工况下转静盘腔内动盘面边界层夹带流量,从而影响了一维流动模型对核心区气体转速的计算.通过归纳数值模拟结果中的自由盘和转静盘腔内动盘面边界层夹带流量,对传统一维流动模型中源区边界及核心区气体转速的计算提出修正,并对修正后的一维流动模型进行实验验证,表明该修正明显改进了一维模型对腔内径向压差的计算精度.      

二维可压湍流边界层改进的积分方法

采用含高阶项的边界层动量积分方程,同时对常规边界层卷吸方程及延迟方程进行相应的高阶影响修正,得到改进的边界层积分方程组。应用此方程组对多种高亚音速及跨音速翼型边界层流动作了计算,并与一阶Ｇｒｅｅｎ方法计算结果及实验结果作了比较。结果表明,在边界层积分方程中保留法向压力梯度项及雷诺法向应力项,明显改进了翼面边界层接近分离区域处参数的计算精度。     

机翼边界层计算及其在飞机设计中的应用

本文给出一种计算实际机翼可压缩层流-紊流三维边界层的新方法。方法使用了Cebeci提出的一种非正交坐标系,并采用一种二次变换以克服在某处横流发生反号时计算遇到的困难。对Lockheed C-5A机翼和F-8超临界机翼给出了一系列计算结果,并同其它方法的结果以及试验数据作了比较。     

A study on turbulence transportation and modification of Spalart–Allmaras model for shock-wave/turbulent boundary layer interaction flow

It is of great significance to improve the accuracy of turbulence models in shock-wave/ boundary layer interaction flow. The relationship between the pressure gradient, as well as the shear layer, and the development of turbulent kinetic energy in impinging shock-wave/turbulent boundary layer interaction flow at Mach 2.25 is analyzed based on the data of direct numerical simulation(DNS). It is found that the turbulent kinetic energy is amplified by strong shear in the separation zone and the adverse pressure gradient near the separation point. The pressure gradient was non-dimensionalised with local density, velocity, and viscosity. Spalart–Allmaras(S–A) model is modified by introducing the non-dimensional pressure gradient into the production term of the eddy viscosity transportation equation. Simulation results show that the production and dissipation of eddy viscosity are strongly enhanced by the modification of S–A model. Compared with DNS and experimental data, the wall pressure and the wall skin friction coefficient as well as the velocity profile of the modified S–A model are obviously improved. Thus it can be concluded that the modification of S–A model with the pressure gradient can improve the predictive accuracy for simulating the shock-wave/turbulent boundary layer interaction.     

使用GAO-YONG方程组对可压湍流边界层的数值模拟

应用SIMPLE/SIMPLER方法及QUICK格式求解GAO-YONG可压湍流方程组,对二维零压力梯度可压平板湍流边界层进行了数值模拟。结果表明,不需要任何经验系数及壁面函数的GAO-YONG可压湍流控制方程组对较广马赫数范围(Ma=0.2～2.4)平板边界层各项参数(表面摩擦系数C_f,对数律,亏损律等)能做出良好的预测。      

ALGEBRAIC TURBULENCE MODEL WITH MEMORY FOR COMPUTATION OF 3-D TURBULENT BOUNDARY LAYERS WITH VALIDATION

Additional equations were found based on experiments for an algebraic turbulence model to improve the prediction of the behavior of three dimensional turbulent boundary layers by taking account of the effects of pressure gradient and the historical variation of eddy viscosity, so the model is with memory. Numerical calculation by solving boundary layer equations was carried.out for the five pressure driven three dimensional turbulent boundary layers developed on flat plates, swept-wing, and prolate spheroid in symmetrical plane. Comparing the computational results with the experimental data, it is obvious that the prediction will be more accurate if the proposed closure equations are used, especially for the turbulent shear stresses.     

使用GAO-YONG方程组对不可压转捩/湍流平板边界层的计算

与以往的湍流模型不同的是GAO-YONG不可压湍流控制方程组不需要任何经验系数及壁面函数。本文运用SIMPLE方法QUICK格式求解GAO-YONG方程组,对二维零压力梯度下平板转捩/湍流边界层进行了数值模拟。结果表明,方程不仅对于边界层流动的各项细节(如表面摩擦系数C_f,对数律,亏损律等)能作出良好的预测,成功地解决了以往一般模型不能同时计算近壁区和远壁区的难题,而且能够预报层流-湍流转捩过程。本文还对机械能方程如何影响边界层近壁区特性进行了数值研究。      

计算实际机翼湍流边界层流动的积分法

本文采用一种拟坐标变换掺混方法,计算任意外形机翼的三维边界层流动,亦适用于边界层外流为位流或非位流情况。对C-5A机翼的计算结果表明,该方法具有同微分方法可以比拟的计算精度。并分析了在亚音速情况下粘性对机翼升力和力矩的干扰影响。     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

跨声速翼型有粘／无粘强干扰计算

本文介绍一种二元跨声速激波-边界层强干扰的计算方法。边界层计算采用湍流边界层积分反方法,它借助Whitfield和Swafford提供的既适合附着流,也适合分离流的速度剖面表达式。跨声速无粘流用全速势方程模拟。通过边界上排溢速度来考虑粘性的影响,用有粘/无粘迭代得到粘性流解。本方法计算的结果与其它方法以及实验的结果进行了比较,证明该方法可以在工程上推广使用。     

三维激波／湍流边界层干扰产生的起始分离的预测

本文对由平板上直立尖劈产生的激波/湍流边界层干扰引起的起始分离进行了研究,提出了一种预测起始分离条件的方法。边界层内的分析采用了Johnston的三角形模型,边界层外主流分析则利用了Prandtl-Meyer函数。预测方法的结果表明与文献中的实验数据符合良好,而且比McCabe理论和Korkegi准则及Lu的半经验关系有更好的物理基础。     

旋涡运动与湍流模型

本文力图将熟知的湍流大涡运动现象与湍流模型理论联系起来以减少模型理论中的经验关系。 在剪切流中以球的流体动力系数,分析确定了涡球的运动,并以此确定剪切流中的湍流剪应力。由此得到了一个二方程模型,方程中仅包含球的流体动力系数。 平板湍流边界层的实例计算表明,计算的速度剖面在整个边界层中都能很好地与实验符合。