在正交流线坐标系中数值计算二维升力翼型跨音速绕流

本文将文献[1]提出的流线迭代法进一步应用于计算正交流线坐标系中有攻角时的二维跨音速绕流。用该算法计算了α=1°情况下的6％双圆弧翼型和NACA0012翼型,计算结果与实验数据符合较好。     

D φ—ψ流线坐标系中求解跨音速绕流问题

本文提出了在一种新的流线坐标系，位函数一流函数坐系系中求解跨音速流问题的数值方法。由于跨音速流动对边界条件的变化比较敏感，现在数值计算跨音速流动已大量采用边界拟合坐标系，流线会标系也是一种边界拟合坐标系，且由于流线随着计算的进而而自动调整。因此，本方法与近来发展较快的自适应网格技术有许多相似之处。数值计算表明，本方法在稳定性和收敛性方面都比较好，与实验结果及其它数值计算结果相比，相互吻合程度令人满     

高分辨率TVD格式的应用

在Jameson格式的基础上，利用修正通量法，构筑TVD（Total Variation Diminishing)格式。通过三个典型的一维和二维算例的计算结果表明，TVD格式对激波具有较高的分辨率，它为具有复杂激波系的喷管、叶栅等跨音速流场计算提供了一种有力的工具。     

欧拉方程组的隐式近似因子分解差分方法的研究

本文对求解二维欧拉方程组的隐式近似因子分解差分方法进行了简要的分析,并利用它计算了简单直扩张通道内的亚音速、超音速和跨音速管流、绕无限长拐角和对称棱形翼型的超音速绕流以及绕NACA 0012翼型的亚音速、跨音速和超音速绕流等典型定常流动。计算结果表明,它具有良好的适应性和精度,可以正确地捕获强激波和自动满足Kutta条件,Courant数可以取大于1的数值,从而加速收敛。分析和计算都是在贴体坐标下进行的。方法本身可以推广到三维流的计算。     

弹性机翼的跨音速压强分布计算

本文采用迭代求解的方法计算弹性机翼的跨音速压强分布。用跨音速小扰动势流的二级近似方法求解气动力;结构变形则采用一维简单梁理论或二维平面结构的柔度影响系数矩阵进行计算。对M6弹性机翼所作的示范性计算表明,计算弹性翼的迭代次数与刚性翼计算相当,弹性变形对冀面气动载荷分布和激波强度均有明显影响。     

超声速流场的流线-特征线坐标变换与应用

为了求解二维平面弯曲激波波后流场,讨论并发展了一系列基于流线-特征线坐标系变换的流场代数计算方法。该系列方法依据二维平面流场中气流角、静压沿特征线近似不变的特点,可快速求解平面弯曲激波波后流场。其中,等气流角近似法适合模拟流线轨迹,等静压近似法适合求解波后流场参数。在此基础上,又提出了一种改进的等气流角-等静压混合方法用于计算弯曲激波波后流场。和特征线法对比,等气流角-等静压混合法计算得到的流场基本特征与特征线法得到的结果相同,在来流马赫数Ma=6和Ma=4情况下误差分别仅为0.5%和0.15%,证实了该方法在求解二维平面弯曲激波流场中的适用性。     

旋转冲压压缩转子二维进气流道数值研究

借鉴二维超声速进气道的设计方法,设计了一种旋转冲压压缩转子的二维进气流道,并采用二维雷诺平均N-S方程和Spalart-Allmaras湍流模型对其流场进行了数值仿真,研究了转速、背压对二维进气流道中波系结构、内部流动特性和性能的影响.计算结果表明:所提出的旋转冲压压缩转子二维进气流道设计方法是可行的;二维进气流道中产生的激波系与二维超声速进气道中产生的激波系相类似,所不同的是二维进气流道中产生的激波为弯曲激波,而二维超声速进气道中产生的激波为平直激波;转速和背压对二维进气流道的性能有较大影响.      

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

包括粘性影响的跨音速纵向大扰动非定常翼型绕流计算

本文提出一种计算效率高、并改进小扰动理论的二维跨音速定常和非定常流的计算方法——非定常纵向大扰动流速势方程和边界条件的数值解。本方法还考虑了包括边界层位移厚度以及激波-边界层干扰的粘性影响。文中给出了NACA 0012翼型和NLR 7301超临界翼型绕流的算例,计算结果与实验作了比较。     

正交流线坐标系中计算二维跨音速流的流线迭代法

本文提出了流线迭代法,用以实现文献[7]所提出的正交流线坐标系中的二维定常无粘气体动力学方程组的数值求解。 此方法的基本做法是:从物面起到流场远边界逐条流线地求出流动速度值,同时求出下一条流线的形状;按此方式对全流场进行迭代。迭代的收敛性很好。 用本文方法在零攻角和跨音速绕流的条件下计算了双圆弧翼型、NACA 0012翼型和R.A.E.101翼型。双圆弧翼型结果与实验数据符合较好。钝头翼型结果与其它方法结果相比是令人满意的。     

二元全速势跨音速方程的数值解

本报告提供了二元全速势跨音速方程的一种数值解法。本方法的特点是放弃了繁杂的曲线座标系,而在简洁直观的直角座标系中,利用Jameson的旋转差分格式,并在物面边界上采用速势内嵌场,建立全流场的差分方程。二元翼型的样例计算表明本方法和其他全速势方程的数值解法具有同等的精度。对物面外形作一次附面层位移厚度修正后的算例表明计算结果和试验结果符合较好。 本方法比起采用曲线座标系的方法,更易推广到三元流。     

非结构网格下Euler方程的高分辨率高精度解

提出了一种非结构网格下求解Ｅｕｌｅｒ方程的高分辨率高精度迎风格式。以Ｒｏｅ的矢通量差分分裂为基础,吸收了ＮＮＤ格式的优点,使其具有捕捉激波和滑移线的良好性能;在时间方向上采用Ｊａｍｅｓｏｎ的Ｒｕｎｇｅ－Ｋｕｔａ积分,并结合局部最大时间步长和残差光滑技术加速收敛。最后成功地完成了二维平板激波反射、跨音速Ｌａｖａｌ喷管内的流动和ＧＡＭＭ超音速前台阶绕流等算例,显示了该方法的有效性     

高超声速二维前体进气道一体化优化设计研究

在飞行器前体进气道的一体化优化设计中,最大总压恢复系数是一个必须考虑的参数.本文从二维高超声速进气道的最大总压恢复系数入手,通过理论分析给出了高超声速飞行器从2波系到6波系的二维高超声速飞行器前体进气道的一体化优化设计计算模型.采用拉格朗日乘子法和序列二次规划法(SQP)分别计算了进气道内一道内激波和两道内激波时的情况,给出了进气道的最大总压恢复系数、进气道内马赫数、激波偏转角和激波强度随来流马赫数的变化关系.比较两种方法的计算结果可知采用的计算方法是合理的.     

二维高超声速压缩面斜激波系的优化配置

选择进气道外压缩段的气流转折角为设计变量,以进气道的压升比为约束条件,分别以总压恢复系数和阻力系数为优化目标,基于一维气体动力学,在飞行马赫数6,飞行高度25 km,进气道出口气流方向与自由来流夹角等于0,°2,°4°和6°的条件下,对二维混压式高超声速进气道进行了优化设计,得到了优化的斜激波配置。优化结果表明,在所讨论的各种情况下,以阻力系数最小为目标得到的外压缩段斜激波系与以总压恢复系数最大为目标得到的外压缩段斜激波系基本相同,都近似为等强度斜激波系。运用数值模拟手段验证了等强度斜激波系配置原则,本文提供的方法可以用于二维混压式高超声速进气道的初步设计。     

振动翼型跨音速绕流的数值模拟

从积分形式的非定常Euler方程出发,在固联于振动冀型的贴体坐标系下,采用有限体积法进行离散,并根据对称型TVD格式构造相应的数值通量,数值模拟了振动翼型的非定常跨音速绕流的变化过程。计算结果表明,该方法能较好地捕捉流场中的运动激波,准确地定出激波位置与强度;并且易于推广到求解任意振动方式的跨音速非定常绕流问题。     

基于结构网格的操纵面嗡鸣数值研究

利用无限插值法生成结构化动网格,耦合求解欧拉方程与结构方程,分别模拟了二维带操纵面翼型以及三维NASP机翼周围跨音速流动。对二维带操纵面翼型嗡鸣发生时的流场进行了分析,发现操纵面上的激波后会产生巨大的分离区,且分离区的变化过程滞后于激波变化过程,认为缝隙的存在对嗡鸣发生有促进作用;对三维NASP机翼的研究发现操纵面质量比的增加可以延缓嗡鸣发展。     

底排装置空腔干扰下的火箭燃气流场数值模拟

从二维轴对称Euler方程出发,应用隐式有限体积TVD格式和分区技术数值模拟来流马赫数Ma=2 0,攻角α=0°时超声速绕流与底排装置空腔内的火箭燃气流相互干扰的复杂流场,计算得到了清晰的流场波系结构,包括底排装置空腔中的激波,火箭燃气流在底排装置空腔外底部与外来绕流作用后形成的相交激波、反射激波等,其结果可为底排火箭复合增程弹丸的结构设计和气动特性研究提供理论依据。     

基于二维结构化网格的间断Galerkin方法研究

在二维结构化网格上,对求解Euler方程的间断Galerkin方法进行了研究。应用二阶、三阶精度的问断Galerkin方法对绕NACA0012翼型的跨音速无粘流动进行了数值模拟。在整体逼近精度基本保持不变的前提下,采用moment限制器以消除数值解在激波附近的伪振荡。数值结果表明,间断Galerkin方法能够很好地模拟流场,准确捕捉激波;moment限制器成功抑制了激波附近的伪振荡;高阶格式具有比低阶格式更强的激波捕捉能力。     

前体激波系重构技术及其在高超声速进气道上的应用

提出了一种激波系整体重构的激波控制技术,并对其流动机理和控制规律进行了仿真分析,继而探索了其在定几何可调进气道上的应用.研究结果表明:该激波控制方法可对前体第1级激波进行有效推动,并可减弱第2级激波强度甚至使其完全蜕化,故实现了对前体激波系的重新构造.随着二次流注入缝的角度增大,将前体激波推动至贴口状态所需消耗的二次流流量逐步减少.并且,最多消耗占主流2.04%的二次流,便可使该进气道在马赫数为5.0~6.0范围内保持贴口状态.与已有各级激波独立控制的可调进气道相比,该进气道在调节状态下的总压恢复较高、消耗的二次流较少,具有明显的优势.      

跨音速翼型和机翼的反设计计算方法

给出一种跨音速翼型和机翼的反设计计算方法。对所应用的积分方程反方法引入人工粘性项;采用Riegels因子法消除前缘奇性;对强激波问题采用光滑-松弛过程;并将方程中的系数积分成解析形式;对二维翼型反设计计算还提出了一种封闭形式的正则化条件。算例结果表明,该方法对跨音速翼型和机翼设计是一种有效的工具。