马赫数2.5超燃室中超/亚声流流向涡掺混冷态流场数值研究

结合实验,对一种带亚声速预燃室和流向涡掺混器的超声速燃烧模型燃烧室,在其进口马赫数为2.5的来流条件下,进行了冷态流场的数值研究.计算与实验得到的燃烧室沿程压力分布相当一致,计算与实验得到的激波结构也基本吻合.计算结果表明:首先,在马赫数2.5的来流条件下,亚燃预燃室易于达到启动状态;其次,流向涡掺混器增强超/亚声速流之间的掺混的效果明显,但其掺混深度尚有限;最后,流向涡掺混的超燃室掺混段有着复杂的激波膨胀波波系,波涡干涉和激波附面层干涉结构.     

吸附式叶栅抽吸流与激波相干性研究

以某高负荷静子叶栅为研究对象,应用数值模拟方法分析了其在不同工况下附面层拓扑结构特点,针对该分离现象实施附面层抽吸,分析了激波与附面层相干过程,探讨了抽吸流与槽道激波的相干作用.结果表明:(1)在有激波、无抽吸条件下,附面层发展一般会经历分离泡产生、破碎、附着叶片表面的过程之后进入大分离状态;(2)吸力面开孔进行附面层抽吸,在一定程度上可以提高静压比,同时损失变化不大;(3)在激波后实施附面层抽吸,会使激波向下游漂移,其后附面层分离更为严重,因此在槽道存在激波时,若实施附面层抽吸,应该预先考虑抽吸气流与激波的相互干涉作用.     

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

非结构网格下可压缩黏性流体的隐式算法

采用牛顿线化和Roe近似Jacobians矩阵隐式算法,实现基于非结构网格的三维定常黏性流动问题的求解.分别对平板、跨声涡轮叶片Mark-Ⅱ和CW-22这三种具有代表性的流动进行数值模拟,并详细分析了流动细节,模拟结果与理论以及实验值吻合良好,准确预测了流场内激波、二次流等重要的流动现象.验证了基于非结构网格的可压缩流体隐式计算方法对于黏性流动问题的计算是正确的.      

多组件结构系统布局拓扑优化中处理组件干涉约束的惩罚函数方法

包含大量组件的多组件结构系统布局拓扑优化设计中存在大量的组件干涉约束,研究了包含大量组件的结构系统整体式拓扑布局优化设计问题,基于有限包络圆方法（FCM）提出了处理组件干涉约束的惩罚函数方法,构造了包含结构刚度和组件之间几何干涉函数的内外混合惩罚函数,应用基于梯度的优化算法对包含数十个组件上百个干涉约束的多组件结构系统进行刚度优化设计,得到了清晰的支撑结构构型和无干涉的组件布局位置,充分体现了提出的混合惩罚函数方法在解决多组件结构系统布局拓扑优化设计中组件干涉问题上的有效性和适用性。     

“咽”式高超进气道流动特性及性能分析

采用数值模拟的方法比较分析了一种矩形型面的内收缩进气道和一种椭圆型面的"咽"式进气道的流动特性和性能.这两种内收缩进气道都是以四道平面斜激波三维流场为基本流场,利用流线追踪技术得到的.研究结果表明,该"咽"式进气道对设计状态下的攻角变化不太敏感;在非设计状态下具有较高的流量捕获和压缩能力;另外,由于其浸湿面积小,进气道内附面层增长缓慢,激波与附面层干扰较弱.因此,这种"咽"式流道可作为吸气式高超声速飞行器进气道的一个有利选择方案.      

粘性高超声速激波—激波碰撞的计算和分析

本文用特片变量型的隐式ＮＮＤ格式计算了五种类型的高超声速激波－激波干扰流场，得到了清晰的流场结构图画。激波和剪切层等分辨得很清楚。计算的物面压力与实验甚为吻合。借助于激波极曲线，对第４类激波干扰结构进行了分析，解释了干扰引起物面压力增高的原因。     

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

光学干涉技术定量显示双马赫反射的密度场

本文介绍了对激波管流动采用两次曝光全息干涉、Ｍ－Ｚ干涉和差分干涉三种技术所获得的双马赫反射二维密度场的定量测量结果，并与数值计算结果作了比较，两者十分吻合。     

刷式密封刷丝摩擦磨损特性实验研究

建立了基于柱面圆周摩擦的直接测力法与间接转矩法两种刷式密封摩擦因数实验识别模型,设计搭建了刷式密封摩擦磨损特性实验装置,设计加工了8种不同结构参数的刷式密封实验件,实验对比分析了两种摩擦因数识别模型的实验结果,研究了结构参数对刷丝与转子表面正压力影响、工况参数对刷丝摩擦因数影响、干涉量对刷丝磨损影响。结果表明：直接测力法与间接转矩法测得的刷丝摩擦因数彼此相差较小,直接测力法稳定性优于间接转矩法。刷式密封刷丝与转子表面静态和动态正压力均随后挡板保护高度的增加而增大,随刷丝径向长度的增加而减小,随刷丝束厚度的增加而增大;相同干涉量下,静态正压力大于动态正压力;刷式密封刷丝与转子表面正压力在干涉加载阶段大于干涉卸载阶段,在转速升高阶段大于转速降低阶段,刷丝表现出滞后效应。刷式密封刷丝摩擦因数随刷丝与转子间干涉量的增加而降低,随转子转速的增加而降低;随摩擦时长增加,摩擦因数在磨损初期先迅速降低,之后基本保持稳定。刷式密封刷丝磨损量随刷丝与转子间干涉量的增加而增大,刷丝与转子间干涉量由0.3 mm增加至0.4 mm时,刷丝磨损量增大了296.66%。     

立楔诱导的高超声速流动分离数值分析

高超声速飞行器进气道等关键部件引起的激波与边界层相互作用将导致流动分离,从而改变当地压力分布与局部受热情况,影响飞行稳定性与飞行安全,因此需要对高超声速流动的分离现象进行细致研究。采用高精度5阶特征型WENO格式与3阶TVD型Runge-Kutta方法,求解三维Navier-Stokes方程,对立楔诱导的高超声速激波与边界层相互作用引起的分离流动流场结构进行了细致的数值模拟与分析。结果表明,5阶特征型WENO格式分辨率远高于类TVD格式;Ma=6时得到清晰的激波结构、分离涡结构及其演化过程和壁面极限流线的拓扑结构,证明了WENO格式应用于高超声速分离流动的可行性与高分辨率;对不同来流Mach数的对比证明Mach数的增大抑制流动分离,导致分离涡减小。     

超声速流中三波结构的参数优化研究

对三波结构的参数进行了探讨，给出了在反射激波和主激波后一些气动力变量的比值存在极值，文中称这种结构为优化结构，当给定气流的马赫数时，理论上就可确定优化三波结构的入射激波强度和存在区域，这可为一些超声速气动装置的设计提供参考依据。     

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

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

爆震发动机超声速进气道的非定常耦合流场数值研究

采用带化学反应的非定常数值仿真方法,对爆震发动机(PDE)内进气道、燃烧室、尾喷管的耦合流场进行了分析,并着重研究了爆震压力波及其反射波对超声速进气道内流动结构的非定常干扰过程,对比了3种进气道/燃烧室耦合方式下进气道流动特性的差异.结果表明:通过在燃烧室内填充以化学恰当比预混的燃料和氧化剂,数值仿真可获得稳定自持的爆震波.在爆震波压强作用下所产生的结尾波系和反射激波会发生耦合干扰作用,它们在流道内的运动传播影响了结尾波系的前传幅度.通过对比不同进气道/燃烧室耦合方式的流场特征发现:采用突扩式进气道/燃烧室连接段可利用燃烧室头部固壁的阻碍以及其对爆震高压波的反射作用将大部分的压缩气体存储于燃烧室内,进而减弱进气道与燃烧室之间的耦合干扰作用,对提高进气道的工作稳定性有利.另外,内凹式燃烧室头部段的引入还为压缩气体提供了额外的存储空间,故可以进一步缓冲爆震高压,提高进气道的稳定工作裕度.      

Effects of sweep angles on turbulent separation behaviors induced by blunt fin

An experimental study was conducted on turbulent separation behaviors induced by blunt fins with different sweep angles at Mach number 6.0. The Nano-particle based Planar Laser Scattering technique(NPLS) was applied to visualize the flowfield, complemented by pressure tests.Sweep angles of the fins were 10°, 20°,...,60°, with the same leading edge diameter of 10 mm. Fine structures of the interference flowfield induced by blunt fins have been obtained, including the shock systems and vortexes. I...     

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。      

等截面矩形隔离段内流场的三维数值模拟

针对三维矩形等截面隔离段中由激波／附面层干扰诱导的复杂流场，采用有限体积方法，OC－TVD格式结合B－L代数湍流模型隐式求解了三维雷诺平均N－S方程，数值模拟了三维流场内复杂的激波串和旋涡结构，计算结果表明所采用的数值方法能够很好地捕捉隔离段内流场的复杂流动现象。     

高超声速多体干扰与分离试验

进行了高超声速飞行器多体系统分离过程中存在的气动力干扰试验研究.试验模型是某构型的可重复使用航天飞行器,由助推器以及再入体两部分组成.研究在FL-31风洞中进行,试验马赫数为Ma=6.97.试验结果表明:分离过程中助推器和再入体之间存在复杂的激波干扰现象,多体系统分离过程中的气动干扰本质上是激波干扰引发的.      

Continuous Measurement of Shock and Contact Discontinuities Velocity

An experimental method applying microwave techniques to obtain continuous measurement of both the shock and contact discontinuities bounding an air plasma generated in a cylindrical hypersonic shock tube is developed. X-band microwave signals excited in the TE11 mode reflect from the moving shock and contact surfaces. The resulting FM interference pattern is demodulated, yielding continuous velocity versus time data. Results depicting subtle detail of velocity behavior, particularly of the shock front, were obtained over a shock Mach number range of 9 to 13.     

Shock interactions in the outer heliosphere

Observations of plasma and magnetic fields by Pioneer 10 and 11 and Voyager 1 and 2 reveal that MHD shocks are an important component of the large-scale solar wind structures in the outer heliosphere. This review discusses recent progress in simulation studies of the nonlinear evolution of the solar wind structures, and in particular concentrates on the theoretical development and applications of the shock interactions model. Various stream propagation models, which do not use the Rankine-Hugoniot relations to calculate the jump conditions at shock crossings, have been used to simulate the essential evolution process of isolated streams and the formation and propagation of corotating and transient shocks. They produce fairly good results in the region up to a few AU. In 1984, the shock interactions model was introduced to study the evolution of large-scale solar wind structures in the region outside 1 AU up to several tens of AU. The model uses the exact Rankine-Hugoniot relations to calculate the shock speed and shock strength at all shock crossings. So that the model can more accurately calculate the shock speeds and the accumulated irreversible shock heating of plasma at several tens of AU. The applications of the shock interactions model are presented in three groups. (a) The first group covers the basic interaction of a shock with the ambient solar wind, the formation and propagation of shock pairs, and the collision and merging of shocks. (b) The second group covers the use of the shock interactions model to simulate the nonlinear evolution of large-scale solar wind structures in the outer heliosphere. These simulation results can provide the detailed evolution process for large-scale solar wind structures in the vast region not directly observed. Two selected studies are reported. (c) Finally, the shock interactions model is applied to studying the heating of the solar wind in the outer heliosphere. The model calculations support shocks being chiefly responsible for the heating of the solar wind plasma in the outer heliosphere at least up to 30 AU.