吸气式高超声速飞行器气动热试验研究

为获得吸气式高超声速飞行器气动热环境的数据,开展了气动热试验研究。在激波风洞中,来流马赫数Ma=6.12,来流单位雷诺数Re/L=1.37×107(1/m)试验条件下,对吸气式高超声速飞行器1/4缩比模型进行了表面气动热的测量。试验获得了小攻角变化范围内的飞行器头部前缘、头部上下交线、机身上下表面中心线、机身横截面周向、平尾垂尾前缘、发动机唇口等位置的热流率分布。研究结果表明,吸气式高超声速飞行器头部前缘、前体进气道壁面、发动机唇口、平尾垂尾前缘气动加热最为严重,另外乘波体外形的设计与布局影响热流的分布。     

一种高阶精度大涡模拟技术用于激波/火焰相互作用

基于Smagorinsky亚格子模型和增厚火焰技术,开发了一种高阶精度的反应流大涡模拟求解器,运用于数值研究边界层对激波与火焰相互作用的影响。该求解器的核心在于加入了基于超黏性的激波捕捉技术和时空三阶精度的两步Taylor-Galerkin紧致（TTGC）有限元格式,并通过对一维Shu-Osher问题和二维激波/气泡相互作用问题的计算,验证了求解器对激波、接触间断和湍流脉动等流动细节的捕捉精度,计算结果与实验数据吻合良好。通过对激波管内激波、火焰与边界层相互作用问题的数值模拟,发现由于激波与边界层的相互作用会产生不稳定的激波分叉现象,激波三分叉点传播速度的发展经历了水平匀速运动、小斜率线性增长和大斜率迅速增长3个阶段,由此揭示了激波分叉促进火焰加速的机理。当火焰面传入激波分叉区后,流场不均匀的回流区起到了稳定火焰的作用,一方面分叉结构内火焰面能够为激波的运动持续地供应热量,另一方面局部超声速区域为火焰的快速传播提供动力,使其能紧跟激波。通过对比相同条件下甲烷和乙烯燃烧的数值结果,发现两者爆震点触发的位置都出现在马赫杆后面,热量释放率的变化趋势也大致相同,但乙烯出现爆震的时间比甲烷早。     

向阳面磁顶区可压缩涡旋诱发重联模拟研究:通量传输事件的形成和结构

本文用二维可压缩MHD模型,模拟研究了向阳面磁顶区的涡旋诱发重联过程。重联过程的结果形成了同心的磁流体涡旋。同时还对可压缩模型与不可压缩模型的异同点作了对比研究。结果表明,在同样条件下,可压缩情况的涡旋发展的增长率比不可压缩情况慢;结构特性无显著变化。本文指出,涡旋诱发重联是产生通量传输事件的一种重要机制。通量传输管可看作为磁流体涡旋管。对通量传输事件诸物理量(磁场、速度、密度、压力和温度)的分布特性作了模拟研究,可解释FTEs的主要观测特性。      

不同波系配置的鼓包压缩面流动特性实验研究

为探究波系配置对鼓包压缩面流动特性的影响,通过纹影和基于纳米粒子示踪的平面激光散射系统(NPLS)观测以及压力测量等实验手段,研究了三种不同波系配置的鼓包压缩面诱导的复杂流动,获得了鼓包压缩面上的压力分布、边界层分布特性以及鼓包压缩面上的精细流场结构。结果表明:三种不同波系配置的鼓包压缩面均具有边界层排移能力,在鼓包对称面上,沿流向边界层不断降低,边界层不断向两侧排移,在鼓包侧面形成堆积。尽管如此,压缩过程的不同,边界层的排移过程存在一定差异,对于单锥鼓包,其最先将低能边界层排移掉,而双锥和等熵锥鼓包略为迟缓,但三种鼓包末端的边界层厚度大体一致。此外,三种鼓包压缩面对称面上前中部压力梯度存在较大差异,但鼓包末端具有相同的压力分布。因此,在保证总偏转角一致的情况下,具体选用哪种鼓包配波方式对于边界层排移、增压能力而言没有显著的差异。     

局部凹槽对高超声速边界层中第二模态扰动演化的影响

采用直接数值模拟方法研究了局部矩形凹槽对来流马赫数为6.0的高超声速平板边界层中第二模态扰动演化的影响。引入了透射系数(定义为凹槽下游与上游扰动渐近幅值的比率)来量化凹槽的局部散射效应。数值结果表明:对于较浅的凹槽,频率较低的第二模态扰动被促进,而高频扰动的规律相反。对于大多数情况,凹槽深度的增加使得透射系数减小,这表明凹槽对扰动的促进作用减弱或抑制作用增强。当凹槽的深度超过某一临界值时,透射系数与凹槽深度的关系变为正相关,这表明了另一机制的出现。临界深度受扰动的频率影响:频率越低,临界深度越大。      

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

模拟弹塑性碰撞的修正SPH方法

从基本的无网格光滑粒子法SPH(Smoothed Particle Hydrodynamics)近似出发,修正了模拟固体力学中大变形弹塑性碰撞的SPH方法.在边界处采用修正的边界条件,弹塑性分析过程中采用增量理论计算应力,迭代过程中用守恒光滑法进行滤波修正消除拉力不稳定.对SPH方法进行了程序实现,给出了杆弹塑性碰撞的算例.计算分析表明,SPH方法节点的影响域较大、精度较相同节点间距有限元法的结果有一定差距,但是通过增加粒子数量可以提高SPH的精度,保持了其简单性和计算大变形的特性.      

粉末高温合金FGH96原始颗粒边界及高温原位高周疲劳研究

利用扫描电镜原位观察的方法研究了粉末高温合金FGH96中不同级别的原始颗粒边界(PPB)在550℃下对合金高周疲劳力学行为的影响。结果表明:采用等离子旋转电极(PREP)制粉+热等静压(HIP)工艺制备的FGH96合金中PPB主要由大尺寸γ'相和碳化物组成;不同级别的PPB对高周疲劳裂纹萌生和扩展均无显著影响,裂纹萌生于晶粒内部,裂纹扩展受晶界与应力轴角度影响,穿晶或沿晶扩展;在裂纹快速扩展区和瞬断区,PPB级别严重的FGH96合金断口呈现穿晶和沿PPB断裂的形貌。     

非定常同向激波边界层干扰的实验和数值模拟

运动激波在波前同向气流中传播时与壁面附近边界层发生的相互作用,既不同于定常的激波边界层干扰,也与普通激波管中端壁反射激波与边界层的作用有很大差异。通过实验和数值模拟对这一问题进行研究表明,激波与边界层作用后,波面在壁面附近发生了明显弯曲,且触及壁面;弯曲的激波与壁面进一步作用后发生反射,由于波前气流速度不同,这种反射存在规则反射和马赫反射两种类型;由于激波反射,波后壁面附近形成了一个高压区;激波与边界层作用后,波后边界层内出现了一系列高低压区相间的流场结构。     

波瓣喷管混合流动的数值计算

借助Ｎ－Ｓ方程对波瓣喷管混合装置的流场进行了数值分析。文中采用三维贴体曲线坐标,非交错网格法和代数湍流模型,对波瓣喷管内外两侧的气流通道及波瓣喷管下游的混合段波场采用分块计算方法进行了数值计算。