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

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

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

为了研究在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。燃烧震荡来源于分离区喷注与燃烧的耦合效应,震荡机制与涡无关,可能是由声振机制引起的凹腔自激震荡。     

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

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

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

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

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

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

无碰撞激波的理想MHD模拟

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

黄道面内行星际激波相互碰撞作用过程研究—间距效应

本文用两维半MHD数值模型,模拟研究了不同间距的两扰动源所产生的激波在行星际空间黄道面内的相互碰撞过程。在内边界(18R_s),两扰动中心的间距取为24°、36°、48°、54°、72°、120°,模拟了两激波速度均为1000km/s和2000km/s两种情形。结果表明:1.速度为1000km/s的两中等激波,相距24°时,在1AU以内强烈相互作用而汇合成一新激波;相距36°时,在1AU附近发生汇合,但尚可分辨;相距48°时,在1AU以内有相互作用,但可近似认为是独立传播的;相距72°时,将独立向前传播。2.速眨为2000km/s的两强激波,相距≤48°时,将强烈汇合成一新激波;相距54°时,在1AU附近产生汇合,但尚可分辨;相距72°时,相互作用不足以汇合成新的激波;相距120°时,其传播相互独立地进行。并给出了上述情形下两激波汇合的时间和位置,详细比较了汇合激波与相同能量的单激波的等离子体特性。     

计算尖锥迎面干扰的一种新方法

本文从运动气流的二维激波动力学方程组出发,提出了一种计算平面运动激波和尖锥头激波迎面干扰问题的新方法。该方法可以计算穿透激波的强度和形状,马赫反射时三波点的轨迹,马赫杆的强度和几何形状,以及锥面峰值压力的大小。计算结果与数值模拟、其它计算方法以及实验数据进行了比较,结果表明本方法不仅精度较高,而且保持了简捷,省机时的优点。     

尾缘厚度对涡轮叶栅性能影响的数值研究

采用数值模拟方法,分别通过改变亚声速和超声速叶型尾缘厚度,研究尾缘厚度变化对涡轮叶栅损失的影响,并在宽广工况范围内探讨尾缘厚度对涡轮叶栅性能影响的敏感性。结果表明:尾缘厚度对亚声速叶型的影响较小,涡轮叶栅损失随尾缘厚度的增大而增大;尾缘厚度对超声速叶型的影响较为明显,随着尾缘厚度的增大,尾缘附近的激波强度增强,叶栅通道中的损失明显增大。对于本文所研究的超声速叶型,尾缘厚度的影响在非设计攻角下不会被放大;但随着马赫数的变化,尾缘厚度的影响规律不同。