正交射流煤粉燃烧器冷态流场的数值计算

本文采用ｋ－ε双方程湍流模型对正交射流煤粉燃烧器出口三维湍回流流动进行了数值研究，得到了副射流与主流速度比为２。５时的正交射流尾迹流场持速度分布、回流特性，并比较了不同速度比时的流场计算结果。计算结果与实测数据符合较好，所得结论对分析正交流燃烧器的稳燃机理具有一定作用。     

转轴位置对俯仰动态特性的影响

研究了转轴位置对机翼、机身和三个组合体上翼面的动态特性的影响，对于机翼和机身，当转轴位置后移时其动态迟滞效应更为强烈。当翼面处于转轴之前，上仰时会出现更为的动态迟滞效应。因此鸭翼达到失速状态较晚。而当翼面处于转轴之后时，如尾翼，将会更早地出现动失速，转轴位置对动态迟滞特性的影响可以认为主要是由于俯仰运动时的动态洗流影响所导致的实际有效迎角的变化。由此影响了动态分离涡及其发展。     

大迎角三角翼旋涡运动及其破碎特性的数值研究

从流体力学的基本方程出发，利用Ｈａｌｌ的涡核准柱假设，导出反映涡核运动的Ｎ－Ｓ方程。采用差分方法计算旋涡流场，进而分析三角翼上前缘分离涡的运动特点及其破碎机下。从计算结果可以看出，旋涡的轴向速度向下游逐渐下降，且涡心处于降较快，外缘下降较慢，反映了粘性作用自涡心外缘逐渐下降的特点；涡核外缘的径向速度开始为负，说明开始阶段有流体流入涡核，随着旋涡向下游运动，径向速度有所增加，到一定位置后增加迅速，说     

高速风洞支架干扰数值修正研究

提出了一种计算高速风洞支架系统对飞行器模型纵向气动力干扰量的数值计算方法 ,从跨声速全位势积分方程出发 ,编制了适用于飞行器全机模型及其带支架情况下的跨声速绕流计算程序。通过对双垂尾模型和GBM 0 3模型两个算例的计算 ,讨论了尾支撑位置及其几何外形参数对模型气动力的影响 ,并对GBM 0 3模型带短支杆情况下的纵向实验结果进行了修正。表明该方法对于分析研究风洞模型支架干扰问题并进行支架干扰修正是可行的、有效的 ,可以作为选择尾支撑位置及其几何外形参数和对跨声速风洞纵向实验结果进行支架干扰修正的工具。     

超声速流中激波/湍流附面层干扰数值模拟

采用修正的B／L湍流模型以及多块结构化网格求解了二维N-S方程。分别对超声速流和高超声速流中的激波／湍流附面层干扰进行了数值研究。本文首先研究了进口马赫数为2．96的超声速流。计算结果准确预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征：分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波。计算的壁面压力分布与实验值吻合较好，计算的分离区长度与实验值比较有一定误差。本文还对进口马赫数为9．22的高超声速流中压缩角引起的激波／湍流附面层干扰进行了数值研究。计算结果与实验结果吻合较好。     

Laminar flow control research at TsAGI: Past and present

This paper presents a brief review of activities in laminar flow control being performed at the Central Aerohydrodynamic Institute named after Prof. N.E. Zhukovsky (TsAGI). These efforts are focused on the improvement of the existing laminar flow control methods and on the development of new ones. The investigations have demonstrated the effectiveness of aircraft surface laminarization applications with the aim of friction drag reduction. The opportunity of considerable delaying of laminar-turbulent transition due to special wing profile geometry and using boundary layer suction and surface cooling has been verified at sub- and supersonic speeds through various wind tunnel testing at TsAGI and during flying laboratory experiments at the Flight Research Institute (LII). The investigations on using hybrid laminar flow control systems for friction drag reduction were also carried out. New techniques of laminar flow control were proposed, in particular, the method of local heating of the wing leading edge, boundary layer laminarization by means of receptivity control, and electrohydrodynamic methods of boundary layer stability control.     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

球形粒子在流体中的跟随性

本文从Maxey和Riley的粒子运动基本方程出发,得到了球形固体粒子在相对流动雷诺数很小情形下的分析解。讨论了粒子跟随性对外力、初始条件和流场性质的依赖关系;对均匀湍流场,分析了不同密度比和扰动频率对跟随性的影响