欠膨胀火箭多股喷流干扰流场的数值模拟

本文利用高精度的TVD有限体积差分格式,对欠膨胀火箭多股自由喷流的干扰流场进行了数值研究,给出了两股自由射流的干扰流场的复杂激波系结构,包括内、外层反射激波、相交激波等波系结构。数值结果与实验结果较为吻合。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

准定常强激波马赫反射波形的数值研究

应用DCD频散控制激波捕捉格式,求解二维、多组分、带有化学反应的Euler方程组,数值模拟了准定常强激波的马赫反射问题.研究结果表明:与经典马赫反射理论相比,在强激波条件下,激波诱导的气体分子振动激发和化学反应使马赫反射的三波点轨迹角变小、马赫杆高度变低、楔顶附体激波倾角变小;马赫杆的相对突出量随入射激波马赫数和楔角的增大而增大,而气体分子的振动、离解等真实气体效应能进一步加剧马赫杆的向前突出.     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

高超声速马赫反射

通过实验和数值计算模拟了激波绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项，粘 英和化学源项分别采用迎风ＴＶＤ格式，中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

高分辨率TVD格式的应用

在Jameson格式的基础上，利用修正通量法，构筑TVD（Total Variation Diminishing)格式。通过三个典型的一维和二维算例的计算结果表明，TVD格式对激波具有较高的分辨率，它为具有复杂激波系的喷管、叶栅等跨音速流场计算提供了一种有力的工具。     

冲压发动机进气道掺混段与弹体内外流场一体化数值模拟

采用二阶迎风隐式TVD格式分区求解可压缩N-S方程,数值模拟了导弹、进气道和掺混段内外流一体化通气模型的复杂流场.本文给出了掺混段出口压力与来流压力比为P/P∞=3.5,4.2和4.9情况的数值模拟结果.并分析了流场特性及激波波系结构,总压恢复系数σ和流量系数ψ.同风洞试验结果进行比较,表明计算结果是合理的.文中还对该CFD分区计算内边界处理所采用的数据结构进行了阐述.     

展向宽度对激波入射平板反射类型的影响

激波干扰是高超声速飞行器气动布局和超燃冲压发动机设计中需重点考虑的局部干扰现象,当该现象发生时会产生复杂的波系结构,对流场行为特征产生影响,进而对飞行器物面载荷及发动机性能产生显著影响。采用数值计算方法针对斜激波入射平板问题,在固定来流马赫数5、单位雷诺数7.12×10⁶ m^-1不变的条件下,通过改变上下平板展向宽度研究了6个不同状态下的激波反射类型。结果表明在受限空间内必须考虑侧向溢流对激波反射类型的影响。随平板展向宽度增加,激波反射类型从正则反射逐渐过渡至马赫反射,且马赫杆长度变长并逐渐前移,直至导致内流道堵塞形成脱体激波。采用激波极曲线方法在激波入射角度固定的条件下对两种激波干扰类型的产生机制进行了分析,发现随激波强度增加两束透射激波极曲线上移缩小,进而造成波后流场参数匹配需求,激波反射类型逐渐从正则反射向马赫反射转变,得到了与数值计算结果一致的结论。     

二维折角管道内超声速流动的多重网格计算

本文应用Euler模型和隐式TVD格式计算了二维折角管道内的超声速流动问题,为了有效地模拟激波和膨胀波系的相互干扰及壁面的反射,本文试将多重网格法和TVD格式结合起来进行数值求解。考虑到目前绝大部分文章采用V-循环多重网格法,我们对多重网格法中的几种多层次循环形式(V-循环、W-循环)进行了比较性的数值计算,旨在从中发现合适的多层次循环求解的形式,以期充分发挥多重网格法在Euler方程求解过程中的效益。同时,本文也是对多重网格法和TVD格式结合求解的一次初步尝试。     

Reflection of a moving shock wave over an oblique shock wave

The reflection of a moving shock wave over a wedge immersed in a still gas and the reflection of a wedge induced steady shock wave over symmetrical and asymmetrical reflecting surfaces have received intensive considerations since more than 70 years ago. Here we consider a different shock reflection problem—reflection of a moving shock wave over an initially steady oblique shock wave induced by a wedge immersed in supersonic flow. For the flow condition we considered, five moving triple points, with each connecting an incident shock wave, a reflected shock wave and a Mach stem, are identified. By using the reference frame co-moving with each triple point, the type of each shock wave of this triple point is clarified. The present study is significant in that it treats a new shock reflection problem leading to a new shock reflection configuration and showing potential applications in supersonic flow with unsteady shock interaction.     

平面激波绕矩形障碍物后的绕射与反射

本文通过激波管二维光测模型实验和采用矩形网格的FLIC方法进行数值计算两种手段针对平面激波绕矩形障碍物后的绕射与反射的问题进行了分析和研究。文中给出了绕射波的波形以及平面激波绕过矩形障碍物后与地面相互作用而产生的从规则反射到马赫反射的波系图案。本文考虑了入射激波马赫数M_i及不同矩形障碍物高度H的影响,并将二维光测模型实验同数值计算的结果进行了比较。     

动态间断装配法模拟斜激波壁面反射

基于非结构动网格技术和边界装配思想提出了动态间断装配法,该方法能够应用于求解含有间断的流动问题。无论入射激波还是反射激波都是作为边界进行处理,激波运动速度由兰金-许贡纽(Rankine-Hugoniot)关系确定。激波作为动网格的一部分,其运动由动网格技术实现。采用该方法模拟了超声速二维流场中激波与壁面相交问题,并且与捕捉法进行比较,二者的流场结构符合良好,但是在细节上还是存在明显差异。通过对流动结构的分析,得出采用装配方法得到的流场要优于捕捉方法的结论。激波壁面反射的问题模拟,也说明了边界激波装配方法对于复杂的激波相交问题是具有处理能力的。     

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.     

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

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

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

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

Shock wave configurations and reflection hysteresis outside a planar Laval nozzle

When the pressure ratio increases from the perfectly expanded condition to the third limited condition in which a normal shock is located on the exit plane, shock wave configurations outside the nozzle can be further assorted as no shock wave on the perfectly expanded condition, weak oblique shock reflection in the regular reflection(RR) pressure ratio condition, shock reflection hysteresis in the dual-solution domain of pressure ratio condition, Mach disk configurations in the Mach reflection(MR) pressure ratio condition, the strong oblique shock wave configurations in the corresponding condition, and a normal shock forms on the exit plane in the third limited condition. Every critical pressure ratio, especially under regular reflection and Mach reflection pressure ratio conditions, is deduced in the paper according to shock wave reflection theory. A hysteresis phenomenon is also theoretically possible in the dual-solution domain. For a planar Laval nozzle with the cross-section area ratio being 5, different critical pressure ratios are counted in these conditions, and numerical simulations are made to demonstrate these various shock wave configurations outside the nozzle. Theoretical analysis and numerical simulations are made to get a more detailed understanding about the shock wave structures outside a Laval nozzle and the RRMMR transition in the dual-solution domain.     

振动激发对激波反射的影响

定常激波反射分为规则反射和马赫反射,在不同条件下2种反射结构之间会相互转变。高超声速流动中的激波反射问题常面临高温气体效应,随着温度逐渐升高,最先出现的是空气分子振动激发。通过理论分析和定量计算,研究了振动激发对激波反射及其转变规律的影响。首先给出考虑振动激发的空气热力学模型,并分析其与量热完全气体的差异以及对激波关系的影响;接着分析在规则反射和马赫反射中,振动激发对激波反射流场的影响规律;最后讨论振动激发对激波反射转变2个准则点的影响。研究结果表明,振动激发使激波极线的整体轮廓变大,且这种差异在经过一次激波反射之后被明显放大,会对激波反射的流场产生重大影响;对于激波反射的转变,振动激发使转变的2个准则点都变大,且对规则反射向马赫反射转变的脱体准则影响更大。     

拉瓦尔喷管外发生激波反射工况详细分析

平面拉瓦尔喷管的环境压力与入口总压的比值从设计工况到第三极限工况(管口处为正激波)逐渐增大时,喷管外的激波现象不同.根据已知的激波反射理论,将设计工况到第三极限工况间的工作状况进一步细分为强斜激波工况、马赫反射工况、双解区和正规反射工况、并将由气流偏转角所表示的发生正规反射和马赫反射的条件用环境压力与总压的比值来表示,找出了强、弱斜激波的环境压力与总压比值的临界值,得到这些临界值随着喷管面积比的变化趋势.最后针对面积比为5的平面拉瓦尔喷管,由理论分析得到其管外发生正规反射、马赫反射、强斜激波反射时对应的压比分别为:0.0841,0.0959和0.2005,并用数值模拟来验证理论分析的正确性.      

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。