典型激波针减阻降热特性及流动机理

为研究激波针对超/高超声速钝头飞行器进行减阻降热的相关特性,采用数值模拟方法对6种典型激波针构型的流动特征开展系统研究。给出了激波针长度、来流马赫数对流动特征的影响规律,并对其形成机理进行了探讨。结果显示:马赫数较低时,头部有扰流物的5种激波针在回流区即将分裂时减阻率最大;马赫数较高时,减阻率在回流区分裂前后出现局部峰值,但最大减阻率将出现在回流区分裂后更长的激波针长度下。马赫数为3时,6种构型减阻率达最大时的相对长度在0.8~1.2之间,相比而言,球型、半球型和双锥型的减阻效果最好,最大减阻率为45%~50%;圆锥型最差,为20%~25%,明显低于头部有扰流物的构型。相同的激波针长度下,头部有扰流物构型的减阻率随马赫数增大而增大,圆锥型则相反。流场回流区、分离激波、弓形激波、局部膨胀流动等导致的压力分布变化是构型整体阻力变化的主要成因。      

高超声速双模块内转式进气道的流动特性研究-Part II 攻角影响

为了研究飞行攻角对高超声速双模块内转式进气道流动的影响,本文通过试验和仿真方法,获得了0°,4°和6°攻角条件下进气道模块内的流动结构。结果表明：在本文研究的攻角范围内,进气道均可起动,进气道压缩面侧的压力变化体现了基准流场的流动特性。在耦合作用下进气道模块间压缩面诱导的激波形态沿流向由弓形逐渐发展为钟形,并且在外压缩激波的扫掠影响下进气道的三个角区出现了强度不同的旋涡结构。进气道压缩面侧的角区旋涡随着攻角的增加而逐渐增强,而进气道出口截面上低能流区域随着攻角的增加而逐渐减小。低能流区域内的二次旋涡呈现不同的变化趋势,位于上半截面的旋涡随着攻角的增加其逐渐向上移动,而位于下半截面的旋涡位置基本保持不变。     

弹体超声速绕流与底喷流干扰流场的数值模拟

利用ＴＶＤ有限体积格式数值模拟了弹体超声速绕流与底喷注相互干扰的复杂流场。采用四步Ｒｕｎｇ－Ｋｕｔｔａ时间推进多步法及“冻结粘性”技术,加速了迭代收敛,节省了计算机时。     

气相爆轰波马赫反射非自相似性特征的实验

为了研究在CJ(Chapman-Jouguet)状态下气相爆轰波马赫反射的非自相似性特征,搭建了由驱动段、传播段以及观测段组成的矩形管道;利用纹影和烟膜实验方法分别对3种预混气发生马赫反射现象进行了实验研究.实验结果表明:对稳态或非稳态气体,测量了马赫反射三波点位置高度,验证了马赫反射三波点轨迹线为波动的曲线,即自相似性对爆轰波马赫反射失效;非自相似性的重要特征是三波点轨迹线在楔角初始位置时基本遵循无反应冲击波理论计算结果,随后逐渐偏离并向有反应冲击波转变,之后接近平行于反应冲击波理论计算并稳定地处于两个理论值之间.最后给出在初始压力为10kPa和楔角为30°的条件下,对于所研究的3种气体,其马赫反射三波点发生转变的位置高度分别约为0.8,1.05cm以及0.5cm,可见相同初始条件下非稳态气体的马赫反射发生转变比稳态气体提前.      

单边扩张燃烧室燃烧非对称及非稳态现象研究

为了探究单边扩张超声速燃烧室在不同当量比下出现的多种模态及其内在机理,通过实验和数值模拟,对该构型燃烧室在并联凹腔对称喷注的工况下进行了系统研究。实验中,隔离段来流入口马赫数为3.46,来流总温为1486K;采用乙烯作为燃料,当量比范围从0.26到0.51。结果表明:燃烧室工况随着当量比的上升而改变,且可分为三个模态。当量比较低时(0.26),形成只有剪切层内有火焰的对称燃烧。当量比较高时(0.51),形成上凹腔为凹腔内燃烧、下凹腔为剪切层燃烧的非对称燃烧。当量比适中时(0.4),火焰在上下凹腔间来回震荡。在非稳态过程中,预燃激波串沿燃烧室轴向做周期性往复运动,定量分析表明火焰震荡的频率约为170Hz。燃烧震荡来源于分离区喷注与燃烧的耦合效应,震荡机制与涡无关,可能是由声振机制引起的凹腔自激震荡。     

含摩擦和加热的跨声速流动数值模拟研究

喷管中跨声速流动的反应流体动力学是实验和理论研究的一个重要领域.本文主要由两部分构成,在第一部分中对文献[2]中的Murman格式[1]的一般形式进行了简化,得到了守恒形式的欧拉方程,并对其求解方法进行了讨论,同时针对Chattot等的工作展开了分析,并给出了带激波的流动计算结果;而文章的第二部分则研究了包含摩擦、加热以及耗散的跨声速流动,并重点关注了一步Arrhenius化学反应模型下的爆震现象.本文还探究了摩擦力和带燃烧的跨声速流动之间的相互作用,同时详细给出了新的欧拉方程数值求解算法并同已有文献算例进行了对比验证.     

物体穿越斜激波干扰特性的数值研究

当物体以超声速穿越同向运动的斜激波时，它的飞行Mach数和攻角都会发生较大的变化。穿越过程中斜激波与物体脱体激波发生相互作用，形成复杂的激波干扰波系结构，基本的干扰形式为复合Mach反射结构。不同的穿越位置时，由于不同的来流条件，使物体上尾部凹坑内的流场和波系发生很大的变化，随着物体的穿越过程，当头部和凹坑下游壁面与翼面穿越激波时，阻力系数有两次剧烈的变化，法向力基本呈线性下降，压力中心先大幅度后移，物体逐渐穿过激波后又快速前移。因此穿越过程使物体加速，给物体的姿态带来很大的扰动。翼面穿过激波也造成阻力、法向力下降，压心后移。如果对有控制的飞行器，穿越激波时气动力特性和稳定性的急剧变化，将给控制系统带来很大的困难。     

激波和边界层相互作用区域局部烧蚀的理论计算方法

激波和边界层相互作用区域的局部烧蚀计算是极为复杂的非定常湍流分离流动条件下的烧蚀分析问题，本文将这一复杂的物理化学和气体动力学过程简化为二维前向台阶分离区的准定常烧蚀问题来处理，将偏微分的控制方程组简化为超越代数方程组，提供了一个有一定精度的理论分析方法，理论计算结果与实验数据相当吻合。     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用,重点研究了波涡相互作用引起的激波结构的变化过程.应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Euler方程,对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究.漩涡与激波间相对强度的大小对激波结构的影响明显,弱相互作用诱导激波变形,强相互作用诱导激波变形并产生分叉.计算结果表明,应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息.     

无碰撞激波的理想MHD模拟

采用了WENO格式数值求解一维理想MHD方程组，模拟了行星际无碰撞激波，研究了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程，并与粒子模拟的结果进行比对，两者的结果非常类似．模拟结果表明，对涉及无碰撞激波的大部分现象中理想MHD模拟是准确且可行的，同时相对于粒子模拟又有很好的计算效率，便于扩展至二维或三维的情形．