双弯曲入射激波的可控中心体内收缩基准流场设计

采用有旋特征线法设计了一种双弯曲入射激波的可控中心体内收缩基准流场,两道入射激波交于中心体起始点,入射激波和反射激波通过给定激波径向总压恢复系数分布进行反设计,壁面通过给定轴向马赫数分布规律进行反设计.该基准流场分为“三波四区”且压缩效率较高.基于该基准流场设计了圆形进口内收缩进气道并进行了黏性修正,数值计算结果表明:内收缩进气道设计点核心区的流场特征和激波形状与基准流场基本一致;在来流马赫数为4.0~7.0时进气道具有较高的压缩效率和良好的流量捕获能力,设计点喉道截面增压比和总压恢复系数分别为17.7和0.729;来流马赫数为5.0~7.0时内部总阻力系数变化平缓,从0.23下降为0.22.      

三维内转式进气道双激波基准流场的设计方法

探索了一种三维内转式进气道基准流场的设计新思路,基准流场由特征线方法设计的曲面压缩系统组成,包含一道入射激波和一道末端激波,消除了激波在内通道的连续反射。通过数值模拟对该设计思路进行了验证,结果表明:该双激波轴对称基准流场,前缘激波和末端激波入射位置与设计吻合,末端激波入射在肩点且完全实现消波;特征线计算获取的外壁面马赫数分布和CFD结果吻合较好;经过设计,在喉部截面上流动参数比较均匀,总压恢复系数达到0.91;无粘条件下流线追踪进气道完全继承了基准流场的流动特征,流量捕获系数0.999,喉道总压恢复0.88,与同设计条件流线追踪Busemann进气道相当。     

考虑进气道喉道非均匀流场影响的隔离段直连试验

为了研究进气道喉道非均匀流场对隔离段流动特性的影响,设计了能模拟隔离段入射激波和非对称附面层的直连试验装置,并进行了马赫数3.0,3.4和3.8的吹风试验,完成了壁面静压、纹影和NPLS测量,获得了有效的试验数据。结果显示,唇口角对隔离段直连台的起动和耐反压能力影响非常大,随着唇口角度的增加,隔离段的极限反压降低,降低幅度达到22%~32%。与均匀来流相比,激波串很难稳在隔离段入口附近,造成激波串极易被推出隔离段。激波串形态受隔离段入射激波和肩部附近附面层状态的影响较大。NPLS测量系统观察到隔离段内部流场精细结构,如下壁面附面层的转捩现象、激波串的马赫盘、分离和滑移线等。     

基于前体激波的内转式进气道一体化设计

在腹部进气的乘波前体/内转式进气道的一体化设计中,为使进气道捕获截面和唇口型线的形状与飞行器前体激波较好匹配,提出一种基于前体激波形状的一体化设计方法。首先,计算乘波前体流场并提取前体激波形状;其次,将进气道捕获型线（ICC）投影在前体激波曲面上,得到可全流量捕获的进气道唇口型线（IFCC）;再次,给定进气道基本流场的中心体轴线位置,确定基本流场的入射激波形状;然后,给定基本流场的沿程压缩规律,应用特征线法确定进气道的基本流场;最后,将ICC顺来流方向投影至进气道入射激波曲面上,经流线追踪和黏性修正得到最终的进气道型面。数值模拟结果表明,对于典型飞行器前体,在设计马赫数为7.0的条件下,应用该方法得到的进气道流量捕获系数达0.976,隔离段出口截面的马赫数、压比和总压恢复系数分别为3.17、38.9和0.487。     

平面激波绕30°出峰后的绕射与反射

本文针对一平面激波S绕过坡角为30°的山峰后与地面作用所产生的波系图案等问题,通过激波管光测模型实验手段进行了分析和研究。实验是在入射激波马赫数Μ_i=1.18,1.32,1.74,2.00和2.29等五种条件下进行的。实验结果表明,在入射激波S和反射激波R之间出现了一接触面C,并与Skews关于平面激波绕凸角壁面时所得到的结果做了比较。与此同时,从我们所拍摄的纹影流场照片上可以看出,平面激波S绕过30°山峰后与地面作用所产生的反射波均属于马赫反射的范围。     

激波管壁面边界层脉动特性的分析与计算

本文分析了激波管壁面湍流边界脉动压力的强度、脉动压力的频率特性及其传播与衰减特性，给出了这些特性的计算方法，通过算例说明了②区和③区壁面脉和压力的强度及其频率随激波马赫数的变化规律，并讨论了壁面脉动对激波管气流品质的影响。     

激波聚焦引燃可燃混合气体的实验研究

在94mm×94mm方截面激波管上,入射激波马赫数为1 60～1 95内开展了汽油/空气可燃混合气体在激波聚焦方式下的点火实验。当汽油/空气混合比处于化学恰当比附近,激波马赫数为1 69～1 95的激波聚焦能够成功地引燃汽油/空气的可燃混合气体,而当入射激波马赫数低于1 69时,同样条件下的汽油/空气混合气体不能被点燃。流场的纹影照片显示,与未燃状态相比,汽油/空气可燃混合气体被引燃后的反射激波波面变宽。     

径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

高超声速二元进气道前体曲线激波逆向设计

为实现高超声速二元进气道前体激波的压缩能力、压缩效率、流量捕获特性和结构长度可控,提出了平面曲线激波逆向设计方法。采用B-Spline曲线控制激波形状。利用有旋特征线法,求解激波的影响域及决定激波的壁面。设计了一凹激波,并对波后流场进行CFD无黏数值模拟,CFD结果和设计结果一致,验证了设计方法的可行性。此外,还设计了直激波、凹激波和凸激波3种激波,并对其在设计点和非设计点处的性能开展了数值研究。在设计点处,分析了激波的压缩比、激波的总压恢复系数、激波压缩区的压升、激波压缩区出口马赫数和流动角随激波控制角的变化规律。在非设计点处,分析了激波压缩区的流量系数和总压恢复系数随攻角和马赫数的变化规律。     

高超声速弯曲激波压缩侧压式进气道数值研究

采用壁面马赫数呈线性分布的曲面压缩系统改进参考侧压式进气道的顶板,得到弯曲激波压缩侧压式进气道,并与参考侧压式进气道进行了比较.数值研究结果表明:设计状态无黏时曲面压缩顶板壁面马赫数分布与给定的马赫数分布基本一致,并且有黏时其壁面压力分布也与二维曲面的基本相同;同参考侧压式进气道相比,顶板采用曲面压缩能够一定程度地改善壁面压力分布,使其末端压力梯度变化平缓;并且非设计状态下的性能也得到有效地改善,特别是来流马赫数为4时,其流量系数提高6.0%、达到0.799,喉道截面总压恢复系数提高1.9%;来流马赫数为5时,其流量系数提高5.2%、达到0.909,喉道截面总压恢复系数提高3.2%.随着攻角增大,该进气道流量捕获能力增强、隔离段出口截面流场畸变减小,但喉道截面总压恢复系数下降剧烈.      

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

非均匀超声来流矩形隔离段内流场实验

针对超燃冲压发动机隔离段的非均匀进口条件设计了隔离段实验风洞，通过测量隔离段壁面压力和拍摄流场纹影照片研究了矩形隔离段内激波/紊流附面层相干流场。研究发现，隔离段进口的非均匀流使隔离段流场压升特征与附面层发展规律与均匀进口的隔离段流动有显著差异。用截面当量直径取代Waltrup公式中的圆管直径可以取得较好的吻合效果。在进口马赫数小于2时，升高同样的压力，非均匀进口隔离段产生的激波串长度比Waltrup公式预测的长度要长。纹影仪观察发现隔离段内激波存在严重的振荡现象。     

具有无源控制空腔时正激波/湍流附面层干扰的数值模拟

采用雷诺平均N-S方程和B/L代数湍流模型计算了具有无源控制空腔时正激波/湍流附面层干扰流场。计算与实验结果的比较表明,本文方法可较准确地预测激波结构、激波与附面层干扰区流动基本特征及波后流动分离状态、激波位置、波前马赫数等参数。      

一种用于RM不稳定性研究的竖直环形激波管的设计与验证

设计并加工了一套竖直环形同轴无膜激波管,可用于环形汇聚激波诱导下的Richtmyer-Meshkov不稳定性实验研究。与前人工作相比,本文在流体界面的形成以及流场的观测方法上做了较大的改进。通过实验和数值方法,对该竖直激波管产生的环形柱状汇聚激波的参数进行测量和分析,验证了同轴激波管形成柱状汇聚激波方法的可行性和可靠性。在界面形成方面,采用细丝约束肥皂膜技术形成正八边形气体界面,并利用数值方法考察了细丝对界面发展的影响。结果表明在界面发展的前期,细丝的影响几乎可以忽略。利用连续激光片光结合高速摄影相机对流场进行观测,获得了正八边形air/SF6气体界面在环形汇聚激波及其反射激波冲击下的演化过程,并与数值结果进行了对比,获得了较好的一致性,进一步验证了汇聚激波的对称性以及细丝约束肥皂膜技术用于形成多边形气体界面的可靠性。     

跨音速喷管流动的数值模拟

数值方法采用的差分格式是具有高分辨率和快速稳定收敛性质的全变差衰减(TVD)格式。喷管壁面形状和反压可任意给定,对于典型轴对称与二维喷管(即壁面由收缩角为45°,扩张角为15°的两条直线用圆弧光滑连接而成。圆弧的无量纲半径等于0.625,其参考长度为喉道半宽),进行了有激波与无激波流动的数值模拟。计算结果与实验结果以及其他计算结果符合良好,此法可推广到非定常跨音速喷管流动的计算,并可用于工程中喷管设计。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

激波/平板湍流边界层干扰的数值模拟

采用理性GAO-YONG可压缩湍流模型,模拟了激波/平板湍流边界层干扰现象,结果表明:计算所得到的壁面压力分布、摩阻系数分布和速度型分布均与实验值吻合很好,并且比较准确地预报出了入射斜激波与平板湍流边界层干扰所引起的边界层分离点和再附点等流动特性.由此表明GAO-YONG可压缩湍流模型能够准确地用来模拟激波/平板湍流边界层的干扰流动.      

超燃冲压发动机隔离段非对称来流下激波串受迫振荡流动研究

用非对称马赫1.98的隔离段直连实验研究了不同波形不同频率的脉动反压对隔离段内流动的影响。实验结果表明:隔离段能有效地隔离不同频率不同波形周期性反压脉动对上游进气道的影响;反压脉动以第二特征波速向上游传播,但在激波振荡区域,压力脉动主要受激波振荡的影响;厚附面层一侧的下壁面激波振荡区域内压力脉动前传时以指数规律衰减,衰减的程度随着频率和波形的改变而改变。而薄附面层的上壁面激波振荡区域内压力脉动前传时不是单调下降,而是波动的。     

基于边界变差最小的高精度有限差分格式构造

从加权紧致非线性（WCNS）高精度有限差分格式出发，在其重构形式的基础上，根据边界变差最小（BVD）原理，遵循单元边界两侧重构物理量值之差最小的准则，在每个单元内通过两步空间重构，构造了一种新的高精度有限差分格式。一般对WCNS等加权非线性格式的改进都是基于改善色散耗散特性、优化非线性权、提高分辨率等单一途径，将它们作为重构候补函数进行结合，既保持了各自优势所在，又能控制格式整体黏性，所得格式具有丰富的应用场景。通过数值实验，将结果与单一格式进行对比，新格式既能在流场光滑区保证设计的精度，对激波等间断附近的振荡也有很好的抑制作用，提高了对高波数区的分辨率，而且在长时间计算后也有较为精确的结果。面向广泛发展的数值格式，还可以构造出其他新方法，对包含强间断和多尺度的流动问题可以获得更好的结果。