机身前体横截面形状对三角翼涡破裂位置的影响

通过求解三维定常雷诺时均方程,采用剪切应力输运(SST)湍流模型,在亚声速范围内,分别对融合体前体-三角翼组合体和旋成体前体-三角翼组合体流场中翼涡破裂现象进行了数值模拟.模拟结果显示:与旋成体前体相比,在中、大迎角时,融合体前体的分离涡,涡量集中、强度高,进入机翼上方流场后,能够与翼涡密切耦合,彻底地改变翼涡强度的分...     

细长翼身组合体前体非对称涡特性研究

通过低速和高速风洞试验对翼身组合体的前体非对称分离涡气动特性的研究，以及对旋成体非对称涡进行了大量的资料研究，结果表明：本专题所研制的细长翼身组合体的前体在较大迎角下有多个非对称涡；迎角、旋成体的外形，尤其是头部的几何形状是细长前体出现非对称涡的关键因素。     

细长翼分离涡的布涡计算法

 细长翼在迎角稍大时,前缘卷起螺旋状分离涡,使上表面压力降低,升力增加。涡襟翼技术也是利用前缘涡的这一特性提高升阻比的。为计算有分离涡的机翼特性,须研究分离涡层的卷起和涡层之间相互干扰的计算方法。早期Brown和Michael,Smith等在锥形流假设下,应用细长体理论计算过三角翼的气动特性。Sack和尹协远等放弃锥形流假设,用离散涡代替脱体涡层,但仍用保角转绘法处理横流面内绕翼面流动。这类方法对横截面形状较复杂的细长翼（如带涡襟翼的机翼）,因转绘函数复杂,计算困难。本文为避免转绘带来的困难,采用直接布涡法计算有分离涡的机翼气动特性。     

组合体大迎角侧向气动特性研究

本文通过实验和理论分析，集中研究小展弦细长翼的翼身组合体的大迎角横向气动特性。研究表明，在大迎角定常非对称涡的范围内，由于翼身组合段对后柱体的边界层分离起遮蔽作用，大大削弱了非对称头涡在后柱体上诱导的侧力。实验证实，平置式翼身组合体的侧力要比单独体的侧力大；带两对弹翼的一般翼身组合体，它的侧力主要由前体以及弹翼组成，如果前体涡在弹翼上诱导的侧力与前体的侧力同向，则该侧力要比平置式布局“-О-”的侧力大得多。     

利用涡的有利干扰推迟翼涡破裂

本文通过测力和水槽流态观察试验研究了战斗机和导弹式的翼体组合体翼涡破裂的推迟措施。利用安置于机翼(弹翼)前方和机体两侧的大后掠、小面积的机体边条所产生的边条涡的有利干扰,可以有效地推迟翼涡的破裂,从而达到提高最大升力系数和临界迎角的目的,试验表明,安置在不同位置的机体边条均可不同程度地提高最大升力系数C_(Lmax),在适当位置时,可提高临界迎角α_(kp)达2°～3°。     

有翼飞行器大攻角分离涡模拟

基于势流方程,采用涡格法和三维离散涡相结合,模拟具有升力面分离涡的大攻角弹冀、弹体、尾翼组合体绕流。该方法适用于高亚音速、小攻角或大攻角有翼飞行器的压力分布、法向力和压力中心等气动特性的数值计算,并具有显示三维涡流场的能力。     

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。     

非线性涡格法应用于大攻角战术弹新进展

应用非线性涡格法（ＮＶＬＭ）计算了大攻角战术弹的气动性能。通过迭代方法、分离涡理论模型的改进，以及采用一系列算法技巧，克服了收敛性上的困难。成功地模拟带弹身脱体涡的翼体组合和翼体尾组合各种分离涡和非线性气动性能；其中包括更为复杂的“××”布局。大大增加了ＮＶＬＭ在工程实用中的适应性和灵活性。     

小展弦比弹翼和翼-身组合体大攻角气动特性的数值计算

大攻角时,绕尖前缘细长弹翼上的附面层从前缘分离而形成脱体涡,对于翼-身组合体,除了弹翼上产生脱体涡以外,在弹身背风面上的气流亦发生分离而形成脱体涡(如图1所示),在这些涡系的作用下,使弹翼或翼-身组合体的气动特性随攻角呈非线性变化。     

激光片光技术在非对称分离实验研究中的应用

主要讨论激光片光技术在细长翼身组合体前体非对称分离实验研究中的应用。通过激光的片光显示非对称分离涡，使研究人员能够清晰地看到细长翼身组合体前体的非对称涡的形成和形态。实验采用了四种前体的翼身组合体模型。通过片光显示，得出了每一种模型出现非对称涡的起始迎角。研究表明：l：3尖锥、l：3尖拱、l：3顶点倒圆的尖拱、l：l顶点倒圆的尖拱四种头部出现非对称涡的起始迎角各不相同，它与各自的分离特性关系密切。     

用全位势方程计算侧滑机翼的跨声速绕流特性

本文对大后掠小展弦比细长机翼，由垂直于机翼纵轴的每一横流截面生成Ｏ型网格，从而形成对机翼流场的Ｈ－Ｏ型网格，用守恒型全位势方程，差分方法和隐式近似因式分解ＡＦ２迭代算法计算侧滑机翼的流场。     

低速大迎角气动力快速计算新方法及工程应用研究

以势流理论为基础,采用工程方法并参考试验结果确定分离涡的起始分离位置,用不可压流涡段所产生扰动流场的线性迭加,计算涡分离流动流场。单独翼面、组合翼面标模和FC-1飞机计算结果与实验结果的对比表明,在失速迎角范围内,本文方法是一种有效的低速大迎角气动特性快速计算方法。     

具有后挡板的喷流流场的数值研究

本文对发射台底面与后挡板之间的喷流流场进行了数值研究。采用迎风的Ｐｅｔｒｏｖ－Ｇａｌｅｒｋｉｎ有限元方法求解欧拉方程。描述了喷流出口处及挡板附近的流场细节。计算结果已成功地应用于导弹发射台的设计。     

一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

不同马赫数的无粘和粘性流动高阶精度隐式计算方法

 应用隐式时间推进法对不同马赫数的无粘和粘性流动进行数值分析,给出了基于预处理方法的高阶精度隐式求解方法。方程离散采用改进的高阶精度对流迎风分裂格式。该方法通过求解曲线坐标系可压缩Euler和 Navier-Stokes方程,对低速到超音速范围内的无粘和粘性流动的典型问题进行了数值分析。计算结果与文献的数值结果或实验数据对比分析表明,该方法对低速到超音速范围内的无粘和粘性流动问题进行数值分析是可行而有效的。     

梯度加权局部分析解差分格式及其应用

提出了一种新的针对Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程符合含有梯度项的对流扩散方程的梯度加权局部分析解的有限差分格式，且给出了通用的系数表达式。ＧＷ－ＬＡＤ格式是基于积分型的有限体积法并且通过附加源项的方法来实现。能够适用于粘性和无粘性及不可压缩和可压缩的压力驱动流动。通过对叶栅流道中二维湍流流的计算表明ＧＷＬＡＤ格式与ＥＤ格式相比，具有较高的精度和很快的收敛速度并且得到更好的压力场和速度场。     

基于CFD/CSD耦合的叶轮机叶片失速颤振计算

为了研究叶轮机叶片的失速颤振特性,发展了一种计算流体力学与计算结构力学(CFD/CSD)时域耦合方法。该方法通过每一物理时刻CFD和CSD的循环迭代实现了耦合计算。在CFD分析中,采用鲁棒性较好的空间离散格式AUSM+-UP,并基于延迟脱体涡模型(DDES)模拟了带分离流动。在结构分析中,通过模态法构建了旋转叶片动力学方程并运用杂交多步方法进行求解。以孤立转子Rotor37为例,计算了不同工况下流场总体与细节参数,与实验结果的对比验证了CFD算法的精度。对某转子叶片进行了颤振特性研究,计算所得的广义位移时间响应曲线表明该叶片在近失速工况下会发生失速颤振,其表现形式为一阶弯曲模态发散且各阶模态之间不耦合。分析表明,流场不稳定和非定常效应是引起失速颤振的关键因素,同时折合频率的降低也会导致原本气动弹性稳定的叶片发生失速颤振。     

应用TVD格式分离求解不可压N-S流动

提出一种新的将连续方程与动量方程分离求解的TVD引用方法。对动量方程中的对流项采用三阶精度、迎风TVD格式离散, 粘性项和压力项采用中心差分格式离散。压力采用压力替代法求解。计算所用的网格是具有交错网格主要特点的修正非交错网格, 并采用一种压力修正法和G-S中心对称迭代法加速收敛。给出了二维管流和二维空穴流动在不同状态下的计算结果, 并与其它解法的计算结果或精确解进行比较。结果表明, TVD格式不仅能用于不可压流数值计算, 而且还能得到比非TVD差分法精度更高的结果。      

旋成体大迎角分离流的Euler方程数值模拟

用三维Euler方程计算旋成体的大迎角绕流,我们发现:亚声速自由流时,计算得到的是无分离流动,与实际的分离流不相符合,因而计算的法向力大大小于有分离的实验值;而当自由流为超声速,且横向自由流M数较高时,由于背风表面出现激波,无粘Euler方程的数值计算能自动捕捉到背风面上的分离涡,计算法向力与实验值接近。为了对亚声速自由流用Euler方程数值模拟大迎角旋成体背风面上分离涡,根据问题的物理图像,在假设的分离线上,提出强加Kutta条件的方法,成功地计算出有分离流动,使计算的法向力与实验值接近。     

Improvements on computations of high speed propeller unsteady aerodynamics

This paper presents the application of the CANARI flow solver to the computation of unsteady effects in the aerodynamic interaction of a high speed propeller with the aircraft. The method is first validated on the APIAN isolated propeller test case by comparison with experiment at M=0.7. The method is then applied to the time accurate 3D Euler computation of a generic transport aircraft at M=0.68. Analysis of the results shows significant unsteady effects both on the propeller forces and on the wing aerodynamic flows, by comparison with steady computations.