尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

新型畸变流场模拟装置

阐明了发动机进口畸变流场模拟方面技术革新的重要性与迫切性，提出一种新型畸变流场模拟装置，给出了它的设计思想、原理结构、工作特点，试验结果。     

超声速喷流混合流场大涡模拟

以光学窗口外冷喷流为研究背景，采用大涡模拟方法对后台阶外形切向喷流混合流场进行了研究。数值方法基于隐式亚格子模型，采用高精度WENO格式进行空间离散，并通过超声速平面混合层流动对数值方法进行了考核验证。喷流混合流场计算模型与试验一致，来流和喷流马赫数分别为3.4和2.5。数值模拟清晰地捕捉到了流场波系以及混合剪切层、壁面边界层等典型流场结构，并精细预测了混合层发生失稳、转捩及发展为充分发展湍流的时空发展过程。数值模拟得到的湍流大尺度结构的位置和形态与实验图像一致。通过对瞬时流场、统计平均流场和脉动参数的分析，揭示了流场结构特征及其时空演化规律，并获得了流场密度脉动特性。      

细长体大迎角湍流流场的数值模拟

通过数值方法对大迎角细长体超声速流场的模拟，探求正确模拟导弹大迎角绕流湍流流场的简单而有效的方法。数值模拟求解一般曲线坐标系下的三维可压缩Navier-Stokes方程，时间离散采用Euler向后差分，无粘项的空间离散采用二阶TVD格式，分别研究了Bald-win-Lomax代数模型及其修正形式(BLDS)对大迎角分离流动的模拟能力。数值试验表明修正的Baldwin-Lomax模型更精确地预测了流场的旋涡与分离情况，给出合理的表面压力分布，因而更准确地模拟了存在横向分离的导弹流场。     

模型超音速燃烧室流场和性能的数值模拟

研究了简化模型超音速燃烧室流场和性能的数值模拟方法，该燃烧室在支板后缘设置一缝隙式喷嘴，平行于超音速空气流喷射氢气。用椭圆型偏微分方程数学模型，ＭａｃＣｏｒｍａｃｋ差分格式成功地计算了两种不同后绿尺寸的流场及其性能，为了对照比较，还计算了另外两种流场。计算结果表明，支板后的流场存在回流区，因喷氢的压力高于超音速空气流的压力和壁面附面层的影响，燃烧室内将出现斜冲波和膨胀波，压力沿横向变化明显，由此证明流场与“边界层流动”有性质上区别。提出了全流场按性质不同分段进行数值模拟的方法，利用计算得到的流场节点状态参数（如温度、压力等），积分获得燃烧室的性能参数，和文献的实验数据对比，计算结果合理、可靠。     

涡扇发动机燃烧室三维冷态流场的计算

论述了某型号涡扇发动机带甩油盘燃烧室从建模到流场计算的整个过程，对半折流式环形燃烧室进行了冷态流场三维流动数值模拟，运用先进的数值计算软件，模拟甩油盘旋转甩油的过程，数值计算得到了燃烧室的速度场、总压恢复系数的分布以及燃油分布情况。为燃烧室的热态流场计算和燃烧室的设计、反设计及改型奠定了基础。     

NASA TrapWing高升力标模数值模拟研究

针对NASA TrapWing高升力全展襟翼构型，采用计算流体力学（ CFD）方法进行三维复杂流场仿真模拟，考察网格尺度和湍流模型对高升力模型气动特性的影响。采用“超立方体”概念，生成绕TrapWing模型的不同网格密度高质量多块结构网格，通过求解雷诺平均Navier-Stokes方程，研究网格尺度对高升力模型气动特性的影响。在此基础上选取中等规模计算网格，考察Spalart Allmaras和Menter k-ωSST湍流模型对高升力全展构型流场模拟能力，分析湍流模型对气动特性的影响。研究结果表明：网格策略具有较好的网格收敛性；湍流模型对翼稍附近位置上翼面压力系数的预测稍有影响，SA湍流模型预测的压力系数较SST更接近实验值。确认研究工作为大型飞机增升装置数值模拟提供了一定的参考。     

高过载条件下绝热层烧蚀实验方法研究 (Ⅰ)方案论证及数值模拟

分析了高过载对固体火箭发动机流场和绝热层烧蚀的影响规律，提出一种新的研究思路：采用数值模拟方法来预示发动机三维两相流场，建立高过载流场模拟实验装置，开展绝热层烧蚀实验，建立高过载条件下的绝热层烧蚀模型，在此基础上发展高过载发动机绝热层设计和烧蚀预示方法。其中关键技术是高过载流场的模拟，对粒子加入法和弯管分离法两种方案进行了论证，认为弯管分离法原理上是可行的。为了验证这种方案的可行性，开展了弯管通道两相流的数值模拟研究，计算结果表明弯管装置具有使凝相粒子聚集形成高浓度粒子流的功能。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

一种鼻锥钝化高超声速轴对称进气道流动特性实验

前缘钝化尺度是高超声速进气道设计中的关键参数。针对一种前体锥加弯曲压缩面的高超声速轴对称进气道,选取最大尺度为3.2mm(5%唇缘半径)的几种典型鼻锥钝化半径,在马赫数Ma=6来流,及模型安装攻角为0°、4°、7°的条件下开展鼻锥钝化尺度对进气道流动性能影响的实验研究。采用纹影拍摄及压力测量记录各来流条件下进气道前体流场结构及壁面压强分布,并在无攻角来流条件下利用微型扰流器进行边界层强制转捩研究。结果表明,对无攻角来流而言,即使是尺度高达3.2mm的钝化半径对进气道前体流场结构及壁面静压分布也基本没有影响。此来流条件下,几种不同鼻锥钝化半径的前体压缩面均出现小范围流动分离,而添加扰流器后该分离区均消失。钝化尺度的影响随着攻角的增加而显现,尽管不同鼻锥钝化尺度下迎风面流场及壁面压强分布几乎没有差别,但背风面随钝化尺度增大表现为边界层明显增厚、流动趋于不稳定。其中最大钝化尺度R=3.2mm的构型在4°攻角来流时背风面即出现明显的分离区,而7°攻角来流时背风面更是出现大范围流动分离、进气道背风侧不起动,并导致进气道内部壁面压强显著下降。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

端壁等离子体气动激励抑制高负荷压气机叶栅角区流动分离实验

为揭示端壁等离子体气动激励抑制高负荷压气机叶栅角区流动分离的影响规律和流场特征,在不同流场参数和激励条件下分别开展了微秒脉冲和纳秒脉冲等离子体气动激励抑制叶栅流动分离的实验研究.结果表明:端壁等离子体气动激励可以有效抑制叶栅角区的流动分离,其作用效果在攻角为3°时最佳,随攻角的增大逐渐下降;微秒脉冲激励的流动控制效果随来流速度的增大而降低,随激励电压和占空比的增大而提高,最佳非定常脉冲频率为500Hz;在较高来流速度下,微秒脉冲激励的作用效果十分微弱,但纳秒脉冲激励能够有效抑制角区流动分离;纳秒脉冲激励的流动控制效果随激励电压增大而提高,激励频率对控制效果至关重要,作用效果随激励频率的增大而不断增强,但当激励频率为5kHz时,作用效果有所下降.      

翼型大攻角绕流的数值模拟

以求解二维N-S方程数值模拟NACA0012翼型大攻角状态的可压绕流特性;N-S方程是在贴体坐标系中给出的,以代数方法生成C型网格系统。采用LU-ADI格式计算,为提高格式的稳定性在隐式和显式部分分别添加了2阶和4阶人工粘性项。应用BaldwinLomax湍流二层代数模型模拟了大攻角时前缘分离涡的形成,旋涡对流及其非定常现象。在某些Mach数和攻角下NACA0012翼型的湍流解具有周期性。通过比较,本文数值计算结果同实验及国外相应的数值计算结果基本吻合。     

多段翼型局部主动变形流动控制的非定常数值模拟

对30P30N三段翼型失速攻角附近的分离流动进行了数值模拟研究。为了抑制大攻角时背风区的流动分离,在主翼段上表面引入了行波壁变形模型和抛物型局部主动振动模型,利用作者以往发展的动态混合网格技术和相应的非定常计算方法,对变形过程中的非定常分离流动进行了数值模拟,分析了各种变形参数对流动分离的影响。计算结果表明,在适当的条件下,局部主动变形能够抑制翼型背风区的分离,由此可以起到增升减阻的作用,改善翼型的气动性能。     

带支板的冲压燃烧室的燃烧性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了燃料注入角度对燃烧及发动机性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的AUSMPW格式,对扩散项采用了二阶中心差分离散.研究结果表明,燃料有角度注入改善了燃料的混合性能,强化了燃料的燃烧效率,提高了发动机燃烧室的动量增量.此外,燃烧室的性能研究需考虑燃烧效率、总压损失及动量增量等因素.     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

双管头部进气旋流-突扩燃烧室冷态流场研究

本文对双管头部进气旋流-突扩燃烧室模型进行了冷态流场试验研究,以探索进气方式(全旋型或部分旋型)、旋流强度(旋流数S或旋流角(?))、旋流室出口扩张角2α、旋流室长径比l/d、内通道相对面积F等参数对燃烧室流场结构的影响规律.结果表明,适当选择燃烧室结构参数可以在燃烧室中形成稳定的中心回流区和头部旋涡回流区.当(?)=45°、α=15°、l/d=1.3、F=0.41时,除了形成旋流室回流区外,在其尾部还形成了一个较大的、切向分速较低的中心回流区,两者“联串”在一起.部分旋的中心回流区长度与相同旋流角全旋进气时的回流区长度几乎相等,但总压损失却降低63.4%.     

涡扇发动机带动力数值模拟

基于等熵流动原理,通过在排气口处设定总温比／总压比,进气口处设定质量流量,提出了一种新型的发动机进排气边界条件。通过与N-S方程耦合求解,数值模拟了轴对称超高涵道比涡扇模拟器CURF和单独涡轮动力模拟器进排气流场,数值结果表明：计算结果和实验值吻合较好,验证了提出的发动机进排气边界条件的正确性,同时也验证了涡扇发动机带动力流场计算方法的可靠性。