高速三维边界层稳定性问题的数值方法

本文把四阶精度紧致差分格式用于求解声速、高超声速三维边界层稳定性问题。提出了适用地求解高速边界层特征值问题的Ｃｈｅｂｙｓｈｅｖ多项式配置点方法,并将各种数值方法进行了比较。数值结果表明上述两种方法是求解高速三维边界层稳定性问题的有效数值方法。     

边界层转换阶段直接数值模拟的影响矩阵法

提出一种求解时间相关、三维不可压完全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的影响矩阵法。考虑流动区域为上半空间的平极边界层流动,在水平方向,利用周期性条件,采用Ｆｏｕｒｉｅｒ离散,在半无穷方向,定义了一种特殊映像并利用Ｃｈｅｂｙｓｈｅｖ多项式展开。由此将Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程转化为以水平波数（ｋ１,ｋ２）为参数的一组法向Ｈｅｌｍｈｏｌｔｚ方程组,同时可更有效地求解压力。本文讨论了连续问题的影响矩阵法与原问题的等价性。影响矩阵法的离散形式,以及影响矩阵法的算法实现等。最后给出了一维Ｈｅｌｍｈｏｌｔｙ问题的数值实验结果。     

战术导弹零升阻力雷诺数效应及修正方法

叙述了雷诺数对战术导弹零升阻力的影响,并给出了适合战术导弹的零升阻力系数雷诺数效应修正方法,即变雷诺数试验外推修正方法及工程计算方法。修正结果表明,修正方法是可行的,变雷诺数试验外推法得到的修正量比工程计算得到的修正量更为合理,对于外形简单的战术导弹工程计算仍有较好的精确度。     

细长锥小迎角气动特性的Re效应

本文介绍了 10°尖锥的气动特性在M =5、α在 0°附近的Re效应的研究成果。气动特性包括法向力系导数、俯仰力矩系导数、压力中心系数、前体轴向力系数、底部压力系数、俯仰阻尼导数、俯仰气动刚度导数、边界层转捩位置系数、转捩雷诺数、转捩诱导法向力中心系数、动 /静稳定性相关因子、气动滞后时间等 12个参数。研究成果表明 ,上述气动特性除前体轴向力系数外 ,Re效应是明显的 ,主要是边界层的不对称转捩影响 ,并有确定的规律性     

用粘性／非粘性相互作用法求解机翼跨声速分离流

提供了一种计算机翼的跨声速绕流的粘性／非粘性相互作用的计算方法,包括无粘流场计算,混合边界层计算及两者之间的相互作用,其中三维混合边界的计算包括了层流边界层,转捩区,湍流边界层和分离流的积分方法计算了,特别是在靠近分离的区域采用边界层反方法计算,无粘流场由全速势方程计算得到,通过粘流无粘流耦合迭代求得了Ｍ６机翼跨声速绕流的收敛解,与实验结果比较,吻合得较好,本方法能够计算出激波诱导分离泡和后缘分离     

高超声速边界层转捩的若干问题及工程应用研究进展综述

高超声速飞行器绕流存在着激波、边界层、流动分离、稀薄气体效应和高温气体效应等多种复杂流动现象的空气动力学问题,其中高超声速边界层转捩既是空气动力学的基础问题,也是高超声速流动研究的热点和难点。若能对边界层转捩进行准确预示及有效控制,则可以实现对飞行器气动力热特性的精细设计,改进飞行器性能,提高任务执行能力。文章针对工程中具有复杂外形飞行器存在的典型失稳特征进行了研究进展回顾,提出了工程实际中亟需解决的复杂边界层转捩问题,明确了高超声速边界层转捩研究的工程应用方向。文章最后还对高超声速边界层的流动控制进行了回顾,以期在今后高超声速飞行器设计中实现对边界层的流动控制,提高飞行器的飞行性能。     

V形前缘对激波入射边界层流动影响的数值模拟与分析

为指导V形溢流唇口下游的进气道内部流动分析,采用数值模拟开展V形尖前缘对二维斜激波入射平板边界层流动的影响研究。以气流偏转角6°的二元楔面为基准激波发生器,设计了展向气流收缩角α(0°~60°, 0°对应二元构型)的V形前缘构型,开展对比研究。结果表明,V形前缘构型使得激波入射位置沿展向不均匀、流动具有明显三维特征,并且干扰区壁面压强上升、分离区尺度明显增大。在α=0°~60°范围内,干扰区流动的不均匀程度、分离区尺度随α增大单调增加。进一步分析表明,V形前缘构型干扰具有中间平直、侧边斜掠的耦合入射特性,体现为对称面壁面压强符合自由干扰理论,侧边斜掠入射区参数符合斜掠干扰的锥形流特征。对比二元与α=45°构型的无粘模拟结果,V形前缘会诱导展向两侧对称的斜掠激波、并在对称面相互干扰产生平直的“桥”激波,这使得激波入射位置沿展向不均匀并偏向下游。其中对称面处平直入射激波压升比(2.49)高于二元构型结果(2.24),侧边斜掠激波强度与二元构型基本一致。这些因素综合导致V形前缘构型的分离尺度增大。     

马赫数6柱-裙激波/边界层干扰直接模拟

高超声速激波/边界层干扰比超声速工况下具有更强的可压缩效应,再附之后会形成极高的局部力/热载荷,严重影响飞行器飞行安全。而激波/湍流边界层干扰区附近流动的三维特性使得流动更加复杂而难以预测。采用直接数值模拟对高超声速条件下的柱-裙激波/湍流边界层干扰进行了详细研究,特别是对Görtler涡结构对分离泡、物面压力和热流造成的展向差异开展了定性和定量分析。研究发现,干扰区附近的分离泡结构呈现出明显的三维效应,且Görtler涡展向分离位置截面的分离泡要明显小于再附位置,而这两个截面上分离泡的运动基本同步,没有明显的延迟或超前现象。物面压力和热流在展向出现显著的不均匀性,展向再附位置的平均压力和热流要比展向分离位置分别高13%和16.2%,脉动压力和热流比展向分离位置分别高28%和20%。截面流向速度特征正交分解结果显示两个位置上的能量都集中在剪切层附近,并且展向再附位置上低频模态占有更高的能量。在采用模态重构流场分析分离区面积发现,展向分离位置重构误差更小,而展向再附位置上的重构误差收敛更快。     

基于大涡模拟的压缩拐角激波/边界层干扰研究

为了提高对激波/边界层相互干扰的基本理解,采用大涡模拟(LES)对来流马赫数Ma=2.9,转角为24°的压缩拐角激波与湍流边界层相互干扰进行了研究。采用回收/调节方法作为入口湍流生成技术,并在超声速平板湍流边界层上进行了验证。采用涡识别方法和数值纹影图等流动显示方法,研究了干扰区内激波与边界层相互干扰的结构变化特征。通过对比分析湍动能和雷诺正应力在不同流向位置的分布规律,研究表明:经过激波干扰后湍动能主要集中在边界层的外层,并在拐角附近形成一个低湍动能区;雷诺正应力流向分量和法向分量在边界层内的分布呈现为单峰结构,而展向分量呈现为双峰结构。运用间歇因子对分离激波的大尺度流向运动进行研究,发现激波围绕着平均分离点作前后运动,运动的尺度等于进口湍流边界层厚度的72%。证实了拐角下游G?rtler流向涡对的存在,并对其展向分布和空间演化特性进行了详细研究。     

红外热像边界层转捩探测的飞行试验应用研究

针对飞行试验层流测试,开展基于红外热像的边界层转捩探测技术应用研究。分析了高空飞行条件对红外测试的影响因素,包括光路布局、蒙皮表面处理工艺、蒙皮加热等,有针对性地提出红外测试方案。在此基础上,以某民用飞机为试验平台,开展飞行演示验证试验,飞行海拔高度范围5km至7km,马赫数范围0.5至0.65。进行了有/无蒙皮加热、有/无太阳等状态的机翼表面热图采集,对比分析蒙皮加热、太阳辐射等因素对红外转捩探测信噪比的影响特性,通过多个架次试验观测分析机翼前缘污染对转捩的影响特性。为消除由前缘污染物诱发湍流楔等因素导致转捩位置误判,提出了一种基于统计分析的机翼转捩位置判定方法,提高了转捩探测结果的可靠性。本文提供了三组真实飞行条件下的边界层转捩位置数据,雷诺数达到1.5×107。演示验证飞行试验结果表明,通过蒙皮内部加热方式可有效提高红外热图信噪比,利用固定转捩方法验证探测结果,使用本文方法测得的转捩位置偏差量不超过弦长0.5%。