用CFD研究高压比离心叶轮内的二次流动

高压比离心压缩机正在喷气发动机和燃机中得到越来越多的应用。然而,人们对其流道内三维粘性流动的认识却非常有限。本文基于已经过实验数据确认的CFD结果,详细地分析和讨论了一种高压比离心叶轮内在近设计点的三维流场,以揭示其中激波与边界层相互作用和二次流动结构等的特征。      

亚、跨、超声速下空腔流场特性实验研究

在FL-21跨、超声速风洞中对空腔流场特性进行了较为深入的实验研究.空腔的长深比L/D为7、8、9、12、15,试验马赫数M范围为0.6、0.9、1.2、1.5、1.8,试验雷诺数Re为1.23~2.26×107.研究表明:空腔流场特性主要由空腔长深比L/D和来流马赫数所决定,这两个参数的变化均能引起空腔的流场类型发生变化.另外,边界层厚度涠钥涨坏撞康难沽Ψ植家灿凶畔灾挠跋?随着边界层厚度与空腔深度之比?D的减小,空腔的流场类型有向闭式穴流动转变的趋势.     

PSE方法在后部加载叶栅绕流稳定性分析及转捩预测中的应用

采用数值求解与实验测量相结合的方法, 分析后加载叶栅绕流的稳定性并预测转捩.推导出正交曲线坐标系下三维扰动波的抛物化稳定性方程, 在风洞实验中测量叶栅叶片表面的静压分布, 并结合测量结果对设计叶栅的流动稳定性进行特征值分析.分析结果表明:(1)后加载叶型能有效推迟转捩, 提高叶栅的变冲角特性;(2)当冲角变化时叶顶流动稳定性的变化最显著, 负冲角流动稳定, 正冲角流动最不稳定.      

高精度数值方法在DNS中的应用

利用不同格式模拟扰动波在平板边界层中的传播,考察了格式对T-S波传播过程计算结果的影响。研究了不同计算格式对自由剪切层扰动波增长的影响,结果表明:对于二维扰动波模型,采用二阶精度格式,扰动波在平板边界层中只能出现衰减的结果;在自由剪切层中,虽然能模拟出扰动波在传播过程中幅值的增长,但却在相当长的计算范围内看不见三维波的出现。但对于同样的扰动模型,如果采用高精度格式,在边界层中能出现增长的T-S波,并与线性理论非常吻合;在自由剪切层中,扰动波不仅迅速增长,并能在不很长的距离内,明显的看到三维波的出现,并激发出了三维结构。      

有二次涡的激波边界层干扰流动的DPNS和NS方程计算

本文针对具有二次涡复杂分离再附现象的激波边界层干扰流动,数值地考察了扩散抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＤＰＮＳ）方程组的适用情况。壁面摩阻和压力、主涡和二次涡的涡高和涡长、分离再附位置以及流线图等特性的计算表明：ＤＰＮＳ方程组的数值结果均与ＮＳ方程组的数值结果很好相符。     

超声速边界层内小激波的数值模拟

本文通过有限差分计算求解了NS方程,对超声速二维平板边界层的非线性稳定性问题进行了模拟.流场的来流马赫数为4.5.平板从尖前缘开始,流场在前缘处产生一道前缘激波.数值计算模拟了边界层对成倍频关系的双频扰动波的吸收问题.通过计算发现,在一定频率下边界层内可出现小激波.计算给出了小激波的发展过程.     

具有复杂边界的火箭冲压发动机流场数值计算

针对实际复杂的补燃结构和计算边界,以块隐式法计算流场,用代数法生成计算网格,通过改变离散方程来处理三维定常Ｎ－Ｓ方程边界条件。算例表明,该方法正确可靠,迭代过程能够快速收敛。     

带轴对称台阶的圆锥高超声速边界层转捩试验

在Ma 6风洞内,通过高频脉动压力测试技术和基于纳米粒子示踪的平面激光散射(NPLS)技术,分别对带前向、后向轴对称台阶的圆锥高超声速边界层转捩进行了试验研究。采用功率谱密度分析和互相关计算等方法对脉动压力数据进行分析,得到了边界层中扰动波的发展规律,定量分析了第二模态波的相关参数。结果显示：两种模型中第二模态波在沿流向向下游发展的过程中,其幅值均先增大再衰减、特征频率均逐渐减小;特征频率和传播速度整体上均随雷诺数的增大而增加(后台阶模型中特征频率由100 kHz增至196 kHz,前台阶中则由 97 kHz 增至174 kHz)、波长变化规律则与之相反(后台阶中由6.35 mm降至4.54 mm,前台阶由7.35 mm降至 4.66 mm );后台阶模型中第二模态波初次出现位置比前台阶中更靠近上游,边界层转捩位置较前台阶前移。将NPLS结果与高频脉动压力测试结果进行对比,两者吻合较好。     

壁面温度分布的高超声速平板边界层稳定性分析及转捩预测

针对来流马赫数为4.5、6.0和7.0的高超声速平板边界层,取30km高空处的气体参数,壁面为等温、绝热和温度分布等3种不同条件,采用eN方法进行转捩预测。其中,壁面温度分布条件下,在等温壁（冷壁）和绝热壁之间,给出4种流向分布进行分析。取扰动的初始幅值为0.3‰,以幅值达到1.5%作为转捩判断依据。结果表明:当温度为来流温度时,等温壁面条件的转捩位置最靠前,并随马赫数增大更加靠前;绝热壁面条件的转捩位置随马赫数增大而推后;壁面温度分布条件下,在相同时刻,马赫数越大,转捩位置越靠前。相同马赫数下,壁面温度较高者转捩位置较靠后（马赫数为7.0时,不完全满足此规律）。在马赫数为4.5和6.0时,绝热壁面条件转捩由第一模态主导,其余情况主导转捩的都是第二模态。     

非均匀地貌的平屋面建筑风荷载特性研究

通过风洞测压试验,研究单个平屋面建筑物在两类均匀地貌,以及非均匀地貌边界层内的屋面风压分布变化规律。研究结果表明:均匀粗糙地貌下的屋面平均风压系数和脉动风压系数大于均匀平坦地貌,且脉动风压系数差别更显著;从粗糙地貌变化为平坦地貌的非均匀地貌中,用建筑物位置处动压无量纲化时,随着距地貌变化点的距离增大,平均风压系数变化较小,脉动风压系数逐渐减小且变化显著;用上游粗糙地貌动压无量纲化时,随着距地貌变化点的距离增大,屋面迎风分离区平均风压变化较小,迎风下游风压幅值增大,脉动风压幅值略有减小,但在过渡边界层范围内均变化较小;屋面整体平均风荷载的主要影响因素是来流动压变化。