曲壁气膜冷却数值计算的研究

 本文在一般二维边界层流动和换热计算方法的基础上,引入湍流混合长度的修正模型,用以考虑壁面曲率和注入冷气对换热的影响,研究计算了壁面的弯曲部分和曲率恢复部分换热的变化,并对曲面上的气膜冷却问题作了计算,取得较满意的结果。     

基于S1,S2流面理论的液体火箭发动机涡轮内部流场计算

应用S1,S2流面理论对某大型液氢/液氧火箭发动机涡轮泵的两级轴流式涡轮内部跨声速流场进行计算和分析。流面的求解采用流线曲率法。将混合平面法的思想用于流线曲率法的求解过程中,对动、静叶排之间的相互干扰进行处理,将复杂的非定常问题简化为定常问题,从而简化计算。粘性的影响通过几种损失系数进行修正。计算结果较好地描述了涡轮内部的流动情况,为进一步完善涡轮设计、提高涡轮性能提供了理论依据。     

基于曲率分布控制的叶型前缘设计方法

为了实现对压气机叶片的优化,提出了一种基于曲率分布控制的前缘造型方法,实现了对叶型前缘曲率的直接、精确控制。将该造型方法应用于某工业级压气机的可控扩散叶型（CDA）上,通过数值仿真方法计算了叶型在设计来流马赫数下的全攻角工况性能。结果显示增加前缘曲率能有效拓宽许用攻角范围,减小尖峰扩散因子,在相同攻角下能削弱前缘吸力峰,抑制甚至消除前缘分离泡,避免提前转捩的发生。同时,调整曲率分布使其在靠近前缘点处尽可能"饱满"、减缓曲率下降速度,也有同样的效果。理论分析发现前缘曲率通过调整静压分布影响边界层发展起始流态,从而影响叶型性能。设计前缘几何形状时需要确保曲率连续性,调整曲率分布以减小前缘吸力峰的强度,避免分离诱导转捩的出现。     

用粘流－无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进。改进的半反－局部联立耦合方法，考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好。算例表明，翼型高升力状态计算必须包括尾流的作用，也应当对势流计算压力与实际压力之间的差别进行修正。     

二维旋转湍流边界层内气固两相流动的数值模拟

本文用修正的ＤＯＲＯＤ－ＦＥＭ方法对不可压流体二维旋转湍流边界层进行了数值模拟，得到了与实验数据较符合的边界层形状因子和壁面摩擦系数的计算结果。采用单向耦合（Ｏｎｅ－ｗａｙ Ｃｏｕｐｌｉｎｇ）方法，计算了固粒在湍流边界层和无粘流中的运动轨迹，给出了壁面被固粒撞击后的冲蚀磨损分布，分析了近壁流场性态与壁面冲蚀磨损的关系。     

计算三维湍流边界层的积分反方法

本文通过分析、计算比较了现有计算三维湍流边界层反方法几种主要版本。内容包括小横向流和大横向流两部分。在此基础上提出一较好的反方法计算版本。所得结论对于如何选取一合适方法有参考意义。     

用粘流—无粘流相互作用方法计算翼型的分离和失速

本文用势流－边界层相互作用方法计算低速翼型的分离和失速。势流用对称面元法。边界层用改进的滞后掺混法，考虑了高阶项影响，适用于计算分离。文中对粘流－无粘流耦合方法作了改进，改进的半反－局部联立耦合方法。考虑了相邻点之间的作用，收敛性较好。计算了ＮＡＣＡ４４１２翼型在不同迎角下的压力分布和气动力。计算结果与实验符合良好，算例表明，翼型高升力状态计算必须包括尾流的作用，也应用对势流计算压力与实际压力之间     

活塞理论气动力静压修正方法及其在曲壁板颤振分析中的应用

提出了一种采用计算流体力学(CFD)计算的压力分布对活塞理论气动力进行静压修正的方法,将该方法应用到曲壁板的静气动弹性变形及颤振稳定性分析中,并与采用曲率修正活塞理论气动力的计算结果进行了对比.分析结果表明,采用本文提出的活塞理论气动力静压修正方法进行曲壁板的气动弹性分析,在圆柱曲壁板曲率较小的情况下,与采用曲率修正活塞理论气动力方法得到的静气动弹性变形、稳定性边界差别不大;而在曲率较大时,采用本文方法计算得到的曲壁板静气动弹性变形,其曲壁板靠近前缘部分被压的更低,而曲壁板的颤振稳定性边界更小,且这种差别随着圆柱曲壁板曲率的增加而不断增大.该方法突破了曲率修正活塞理论的小曲率限制,扩大了活塞理论气动力在曲壁板颤振分析中的适用范围.     

反映流动曲率影响的非线性湍流模式

本文从一般曲线坐标系下的雷诺应力模式出发，在雷诺应力模式向代数应力模式简化过程中，保留工率项对雷诺应力输运的影响，得到了曲率修正的代数应力模式和非线性κ-ε模式。并且该模式对一典型曲率流动－－U型管道内的湍流流动进行了数值模拟与分析。计算结果显示,率修正的代数应力模式较好地模拟出了凸曲率地雷诺剪切应力的阻滞作用和凹曲率对剪切应力的增强作用，捕捉到了曲率流动的主要特征，是对雷诺应力模式的一个合理而有效的简化。     

弦向缝隙叶栅对边界层分离的控制

本文提出了一种在Ｓ１流面上进行势流与粘性流迭代求解弦向缝隙叶栅流场的计算方法，将Ｓ１流面上势流方程组的解与经Ｉ１１ｉｎｇｗｏｒｔｈ－Ｓｔｅｗａｒｔｓｏｎ变换后的边界层方程结合起来联立求解弦向缝隙叶栅流场。计算值同实验结果相当吻合，在正大攻角时，向缝隙叶栅能效地控制边界层分离。     

跨音速弯曲壁面正激波与湍流边界层干涉的计算结果

本文介绍弯曲壁面上湍流边界层与正激波干涉的某些计算结果。壁面曲率对激波前后的压力比有重要影响。在干涉区中,边界层位移、动量和能量厚度都有显著的增加。这种增加主要决定于激波前后的压力比。(TBTF)     

非惯性非正交曲线坐标系中的边界层方程

应用张量分析方法将 N-S方程写于作非惯性运动的非正交曲线坐标系中 ,然后用匹配渐近展开法导出此坐标系中的边界层方程。消去压强项后 ,与惯性坐标系中的边界层方程相比 ,非惯性系中的边界层方程在 x1和 x2 方向的动量方程中分别只多了一哥氏力项 ,它们是由物体转动角速度的物面法向分量产生的 ,其它惯性力项均不进入边界层方程。方程中所含曲线坐标系的拉梅系数及其导数只需在物面上取值 ,物面曲率影响不进入边界层方程     

航空发动机加力燃油总管的数值计算分析

采用商业软件对航空发动机加力燃油总管进行了数值计算,研究了喷油孔的附面层与质量流量的关系,同时分别对喷油孔允许的最大和最小值进行数值模拟,就其流量值与设计值进行了对比。结果表明:喷油孔的加工质量对加力燃油总管流量的影响很大,喷油孔必须按照公差等级Ⅰ的规定进行加工,以保证其质量流量满足设计要求;对于小孔径结构的数值计算,只有充分考虑附面层的因素,才能提高数值计算的准确性。     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

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

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

考虑流向曲率和压力梯度的可压缩边界层稳定性分析

以圆弧,平板和尖锥为研究对象,基于线性稳定性理论,重点讨论了曲率和压力梯度对可压缩边界层稳定性的影响。     

三维跨声速粘性-无粘相互作用的半逆方法

本文为计算三维跨声速粘性-无粘相互作用提供一半逆方法版本。湍流边界层采用积分反方法。证明了对于跨声速流,选取外部流线角α和等价的不可压缩形参作为输入量,边界层积分方程总是双曲型的。无粘外流采用FL027程序。通过粘性-无粘迭代求得小分离区情况下的与实验符合的收敛解。     

超音速细长梯形翼背风面流型

给出了前缘后掠６５°、双弧形剖面的细长梯形翼背风面流动显示结果。实验Ｍａｃｈ数为１．１０,１．５３,２．５３,３．０１和４．０１,攻角范围为５°～２５°。应用蒸汽屏、纹影和油流技术拍摄了空间和表面流型照片。蒸汽屏显示表明：在机翼背风面三角形区域的空间流型随法向攻角αＮ（在垂直于前缘的平面内流速与弦线间的夹角）和法向Ｍａｃｈ数ＭａＮ（来流Ｍａｃｈ数在垂直于前缘平面内的分量）变化,并可在αＮ和ＭａＮ为坐标的平面上划分出７种流型存在的区域。侧缘区有侧缘分离涡形成;后缘有尾涡拖出。从纹影照片与横截面上的蒸汽屏照片对照可获得机翼锥面激波位置随Ｍａｃｈ数的变化;以及激波－诱导分离线位置随Ｍａｃｈ数和攻角变化曲线。机翼表面油流谱显示出了主再附线、二次分离线、二次再附线和侧缘涡区。显示出的流型与其他有关实验和数值计算结果比较符合得很好