欠膨胀椭圆射流的流动结构

为了探讨超声速欠膨胀椭圆射流的流动特征,采用大涡模拟（LES）方法与高精度混合格式对出口压力比N分别为14、24、40的欠膨胀椭圆射流流场结构进行数值模拟。结果清晰描述了欠膨胀椭圆射流的三维(3 D)结构特征与发展规律,并分析了因喷嘴方位曲率不一致而导致长轴与短轴平面上激波结构出现差异的原因。另外,结果还发现:当N为24时,短轴平面上射流域内的激波结构已由正规反射转变为马赫反射,但长轴平面上仍维持正规反射,而当N为40时,长轴与短轴平面上的激波结构均为马赫反射结构,由此可知喷嘴方位曲率变化越平缓,马赫反射形成所需的出口压力比越大。      

马赫数1.95超声速欠膨胀喷流声辐射数值分析

针对马赫数为1.95的欠膨胀超声速喷流声辐射特性,采用高精度计算格式的大涡模拟（LES）方法进行数值研究。通过对喷流平均流及湍流脉动统计结果的对比,确定数值方法的精确性。细致分析了超声速喷流流动特征,特别关注了其中激波胞格和湍流相互作用现象。基于流动与声场辐射间的关联性,分析了超声速喷流马赫波、宽频激波噪声辐射特性及形成机理。数值结果表明超声速喷流的剪切层演化及激波胞格与湍流间的相互作用构成欠膨胀超声速喷流的主要噪声源。     

斜激波极值规律的边界层影响

采用边界层理论与斜激波/膨胀波精确算法，建立一种结合Eckert参考温度法和Illingworth-Stewartson变换法优势的边界层权重算法，用于研究超声速黏性楔面边界层位移厚度对斜激波极值规律的影响。分别应用层流Navier-Stokes方程和湍流Navier-Stokes方程的CFD解算器对边界层新模型进行了算例精度评估。在来流马赫数为1.2~2.4和楔面角为3°~20°的范围内，压强比的相对误差小于0.1%。计入层流与湍流边界层影响的理论模型研究表明，边界层影响使得最优马赫数增加；对于层流边界层，最优马赫数增量约为0.001 5~0.003 3；对于湍流边界层，最优马赫数增量约为0.002 8~0.006 1。     

高度欠膨胀声速射流的自相似特性

构建了高压气体射流模型和高分辨率的计算网格,对喷压比（NPR）为5.60、流动雷诺数Re为105量级的高度欠膨胀射流进行了三维大涡模拟（LES）计算.讨论了时均的射流近场结构,发现大涡模拟成功捕捉到了高度欠膨胀射流近场的典型波系结构,并与文献结果吻合较好.研究着重定量考察流向速度和氮气质量分数的分布规律,以揭示高度欠膨胀射流的自相似特性.结果表明:在射流核心区之后,不同流向位置上流向速度和氮气质量分数沿径向的分布呈现出与亚声速射流类似的自相似特性.但射流流场开始呈现自相似的位置仍然为超声速,这是高度欠膨胀射流自相似特性的独有特征.提出了高斯拟合公式,在靠近和远离射流中心线的区域内均可较好地表征射流的自相似特性.此外,还考察了射流剪切层的发展特征,比较了计算得到的射流剪切层增长率与先前实验测量结果的差异,并分析了其中可能的原因.      

燃气自由射流的正格式数值模拟

采用二阶正格式方法对复杂燃气自由射流进行了数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对不同欠膨胀压力比下的燃气射流进行了数值计算。计算结果表明,该方法能够较好捕捉到波胞、滑移线、射流边界以及三波点等复杂射流流场波系结构,与实验照片反映的流动特征及高精度、高分辨率的三阶间断有限元方法的计算结果吻合较好,马赫盘位置与理论预估值和实验测量结果误差较小。证明该方法对激波具有较强的捕捉能力,为此类流动的数值研究提供了一种新手段。     

计算流体动力学在航天飞行器射流逆流中的应用

超声速射流逆流通常用于导弹、航天飞机、卫星和飞船等飞行器运动状态的控制。欠膨胀超声速射流逆流的流场包含有多激波（如弓形激波和马赫盘）、接触间断和剪切层，其结构非常复杂。本文采用激波高分辨率有限差分（ＴＶＤ）格式，对恒定自由流条件，各种不同射流出出口压比的超声速轴对称逆射流进行了数值模拟，且对各种条件下的物理现象给予了分析。计算的马赫盘和弓形激波位置与实验值相吻合，为此类流动问题提供了一种有效的预测     

不同飞行状态下固体火箭发动机尾喷焰数值研究

为了研究飞行状态对固体火箭发动机尾喷焰的影响,建立了含化学反应项和组分输运项的N-S(Navier-Stokes)控制方程,并采用MUSCL(monotonic upstream-centered scheme for conversation laws) Roe格式的有限体积法进行求解.在此基础上,结合热力计算结果对某固体火箭发动机在不同飞行状态下的化学反应流和单一组分流进行仿真计算.结果表明:初始倾角和第一马赫波节长度随飞行高度的增加而增大,随来流马赫数的升高而减小.同单一组分流相比,化学反应流对第一马赫波节的影响较初始倾角显著.射流轴线温度因飞行状态的不同有较大涨落差异.初始倾角、第一马赫波节由来流总压决定,射流轴线温度由激波结构、复燃效应、来流条件等因素共同决定.      

斜切喷管欠膨胀超声速射流的数值模拟

采用三阶MUSCL TVD格式解三维可压缩平均雷诺纳维尔-斯托克斯方程组，湍流模型为Spalart-Allmaras代数模型，数值模拟了两种斜切角喷管在喷口压比为2.0，2．8，3．4和4．0条件下的欠膨胀超声速射流场，获得了射流场流谱和参数分布，计算结果和实验数据符合良好，并给出了喷管斜切角和喷口压强比与射流偏转角和扩张角的分布曲线。     

静音锥对超声速客机声爆水平的影响

静音锥低声爆是通过在超声速飞机头部加装静音锥将机头产生的强激波转化为一系列互不叠加的弱激波,从而降低声爆。以一种“梭式”布局的超声速客机为基本模型,采用计算流体力学和波形参数法相结合的方法,研究不同参数的单级和多级静音锥对超声速客机声爆水平的影响。结果表明：静音锥的长度可以调节激波的干涉程度;静音锥的直径和圆锥顶角可以改变静音锥的初始超压值;静音锥级数对上升时间影响显著。     

定几何混压式轴对称超声速进气道气动特性数值仿真和实验验证

 针对一种Ma4一级定几何混压式超声速轴对称进气道进行了数值仿真研究,并和风洞实验结果进行对照,验证了本文所采用计算方法的可靠性。利用CFD方法获得了进气道激波系分布、内通道流场分布和沿程静压分布并对Ma4下稳定亚临界状态进行了分析。研究结果表明：（1）超临界状态下,随着进气道出口反压的提高,结尾激波系向喉道方向移动,结尾激波损失减小,总压恢复系数提高;（2）攻角的增加对进气道的迎风侧和背风侧影响增大,结尾激波系由对称分布向一边倾斜的趋势增大,背风侧的承受反压能力下降,总压恢复系数随着下降;（3）随着来流马赫数的增加,激波损失加大,总压恢复系数随之下降,同时由于激波角变小,激波也越靠近外唇罩,溢流减小,流量系数增大,在激波贴口后流量系数基本保持不变;（4）通道内的静压分布曲线清晰的反映了内通道沿程激波系情况;（5）在大于贴口Ma数工作时,结尾激波系被推出唇口的情况下,由于滑流层作用出现一个类似外压缩式的气动通道从而存在稳定的亚临界状态。     

超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

当量比对带凹腔超声速燃烧室流动及燃烧特性的影响

通过调节燃料氢喷射压力参数从而改变油气当量比,分别在冷态喷流与燃烧情况下计算了当量比对流场波系及燃料与空气的掺混、燃烧特性的影响.研究表明:增大当量比会增大燃烧反压,扩大上游边界层分离区,逐渐将激波串推向燃烧室进口;同时发现燃料与空气的混合性能在近场处受当量比影响较大,随当量比增大混合效果明显增强,而在远场处基本不再受当量比变化的影响;增大当量比一定程度增大射流穿透深度,但并不能明显改善燃料分布区不合理问题,导致燃烧效率随当量比单调下降;流场总压损失及凹腔阻力系数均同当量比成正比例关系,而流场摩擦阻力随当量比单调下降,且燃烧使摩擦阻力更小.      

脉冲能量沉积对超声速射流/激波相互作用掺混的控制研究

为探究超声速射流掺混增强的有效控制方法,基于二维非稳态雷诺平均Navier-Stokes方程,研究了脉冲能量沉积对超声速射流与斜激波相互作用后增强掺混的控制机理和规律。宽度为5 mm的低密度平面射流同向完全膨胀射入马赫数为2.5的主流中,并与20°压缩斜面所产生的斜激波相互作用。在流场中引入脉冲能量沉积对射流掺混进行控制,并考虑射流马赫数、射流长度、能量沉积位置及激励频率对掺混效果的影响。对比激励流场与基础流场,结果表明:脉冲能量沉积的加入诱导射流形成大尺度涡结构,显著提高射流的掺混效果;能量沉积位置对于不同马赫数射流的掺混增强具有不同的控制效果;激励前后射流宽度的对比表明,脉冲能量沉积的无量纲激励频率在0.22附近时,可使得射流的掺混效果最佳。     

涡轮叶栅超声速流场流动特征与气膜冷却特性

应用shear strain transport(SST) k-ω 两方程湍流模型,对超声速涡轮叶栅通道内气膜冷却特性进行数值研究,得到不同气膜孔倾角和吹风比下叶栅通道内流场流动特征以及气膜冷却效率的变化规律.在激波入射点附近的气膜射流能够向分离区边界层中补充动量,克服逆压力梯度,有效改善由于激波引起的局部过热.亚声速流动状态下的气膜入射角度对冷却效率的影响能够在较大吹风比下得以体现,而超声速主流状态下,气膜冷却效率与入射角度基本无关,说明亚声速的气膜冷却射流对超声速主流的穿透力要弱于对亚声速主流的穿透力;超声速主流条件下,在激波入射位置的气膜冷却效率要高于激波入射位置下游的气膜冷却效率,这与气膜孔出流在当地的湍流度有关.      

进口斜激波、膨胀波干扰下等直隔离段内的激波串特性

 超燃冲压发动机的隔离段在实际工作中会受到进气道唇罩激波及肩部膨胀扇的显著干扰,本文针对这一特定问题进行了专门研究。提出了唇罩入射激波及肩部膨胀扇的模拟方法,并利用德国Achen的风洞试验对其进行了检验,而后以此研究了入射激波及肩部膨胀扇干扰下隔离段内激波串的基本形态,并分析了出口反压和激波入射位置的影响。仿真结果表明：当激波串在隔离段内不断前移时,受唇罩入射激波及其反射激波的干扰,其高速核心区交替地偏向上下壁面;与无激波入射的情况相比,此时激波串的耐反压能力显著降低,且入射点位置越高,降低幅度越大,管道内的沿程静压分布规律与Waltrup经验公式偏离程度也越来越大。该文结果可为进气道/隔离段的一体化设计提供依据。     

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

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

凹腔前缘角对超声速燃烧室性能的影响

针对带有不同前缘角的凹腔内流动和燃烧过程,分别在冷态和燃烧条件下探讨了前缘角对凹腔内流动损失及阻力特性的影响.研究表明:在壁面垂直喷射的喷口上游和凹腔内部均会形成低速、高温回流区,有利于点火及火焰稳定,燃烧反压通过边界层的亚声速区域上传,形成激波/边界层干扰结构.减小前缘角,可使剪切层分离位置提前,更偏向凹腔内部,导致凹腔后壁面再附激波增强,进而增大了总压损失,降低了总压恢复系数;亦可导致凹腔前、后壁面压差阻力增大,阻力系数上升.进一步认识了凹腔内部流场及稳焰增混机理,进而为优化凹腔结构设计提供依据.      

2级脉冲爆震发动机中超声速射流对撞诱导激波聚焦试验研究

为研究2-St age PDE中连续超声速射流及结构参数对撞诱导激波聚焦的影响规律,在冷态条件下开展了2-St age PDE中连续超声速射流对撞诱导激波聚焦的试验。分析了喷口宽度、导流环深度、凹面腔开口端与喷管间距、尾喷管角度、射流入射压力等参数对凹面腔底部峰值压力的影响。结果表明:喷口宽度、导流环深度、凹面腔开口端与喷口间距、射流入射压力越大,凹面腔底部峰值压力越大,激波聚焦效果越好;尾喷管角度越大,凹面腔底部峰值压力越小,激波聚焦效果越差;喷口宽度、导流环深度、射流入射压力对激波聚焦的影响较大。     

煤油在超声速气流中非定常横向喷射的实验观察

本文目的是研究煤油在超声速气流中非定常横向喷射问题,并与静止和亚声速气流中的结果进行对比。改造普通油嘴并采用单次高压燃油喷射系统,得到不同延时的流场阴影照片。结果表明：在超声速气流中,射流柱破碎是由迎风面的表面波引起的,破碎点位于表面波的波谷。雾化过程过分四个区域。射流激波由近射流柱的迎风面,并在人壁反射,非定常喷射的波系演化复杂。煤油在亚声速气流中喷射,不存在射流激波,破碎点沿射流柱上移,穿透深     

超音速流中球头反向喷流流场的数值模拟

采用CSCM(Conservative Supra-Characteristics Method)格式结合激波捕捉法求解了Navier-Stokes方程。通过算例对喷口总压等于2.08倍到6.54倍波后总压等4种情况进行了计算。得到壁面压力分布、等密度线和典型流场速度矢最图。与实验结果和其它数值计算值进行了比较和分析,得到本文计算出的弓形激波脱体距离在喷口压力低时与实验值符合得很好;压力高时差些。