斜切反喷管性能计算

本文提出了斜切反喷管流场计算方法,内部点用MacCormack二步显格式数值解控制方程,入口边界和固壁边界的参数用物理边界条件和可应用的沿双特征线的相容性方程计算.固壁边界计算中采用了激波装配技术和开关格式,有效地捕获了流场内的激波系列,计算中得到的激波情况和物理分析一致.本文在流场计算的基础上,进行了反喷管的性能计算.本文提出的方法可应用于反喷管的性能预估.     

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

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

运动激波绕射的数值模拟

用高分辨率的ＴＶＤ格式对复杂的激波反射过程做数值模拟。激波斜壁反射计算结果与实验干涉照片、压力分布、密度分布均符合得很好。通过对运动激波绕射凸起物的数值模拟,得到非定常反射流场发展过程及复杂流场间断的清晰图谱。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

高超声速马赫反射

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

反射激波冲击单模界面的不稳定性实验研究

激波诱导流体界面失稳问题广泛存在于惯性约束核聚变、超燃冲压发动机、武器内爆等工程应用中，相关研究具有重要意义。本文在改进的水平激波管中开展入射激波及其反射激波诱导单模气体界面失稳的实验研究，采用线约束肥皂膜技术生成较为理想的空气/六氟化硫（air/SF₆）单模气体界面，借助高速纹影技术捕捉激波冲击界面后的详细不稳定性演化过程，重点关注反射距离对不稳定性发展的影响（反射距离定义为初始界面和激波管尾端固壁的距离）。研究发现，在一定反射距离范围内，反射激波作用后的扰动增长率几乎为一个恒定值，与反射距离无关（即与反射激波作用前的扰动振幅及增长率无关），随着反射距离的进一步增大，扰动增长率降低。通过实验测量值与理论预测值的对比，发现Mikaelian模型和Charakhch’an模型（采用合适的经验系数）均能对反射激波作用界面后的扰动增长率给予有效预测，且这两种模型的经验系数都依赖于反射激波作用前的界面演化状态，如果反射距离的改变引起/（不引起）界面演化状态的变化，则需要/（不需要）改变经验系数的值。     

激波与壁面粉尘相互作用的数值研究

本文采用厚边界层理论建立了考虑两相耦合效应的气固两相湍流边界层方程用以研究激波与壁面粉尘的相互作用。方程中考虑了压力梯度，Ｓａｆｆｍａｎ升力和湍流扩散等因素对颗粒飞行轨迹的影响，并用Ｍｉｒｅｌｓ方程作为边界条件使方程组封闭可解。数值计算时，经适当变换的气相方程组用ＢＯＸ格式，拉格朗日坐标中的颗粒相方程用Ｒｕｎｇｅ－Ｋｕｔｔａ法，而源项则用颗粒源法和求质量加权的平均值获得。计算结果与某些实验结果较为     

深水强爆炸的数值模拟

本文彩用ＮＮＤ格式求解Ｅｕｌｅｒ方程组，模拟深层水域强爆炸以及冲击波固壁反射后的流场。     

考虑粘性时激波反射及壁面热流计算

本文采用显式ＴＶＤ格式，求解薄层Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，计算了激波与楔面相互干扰的流场情况，得出了马赫反射的各种等强度级及固壁上详细的物理分布，根据计算得出的流场，计算了壁面摩擦阻力及热流分布。     

真实气体效应下高马赫数内转进气道特性研究

为初步研究高马赫数内转进气道在真实气体效应下的工作特性,首先设计额定工作状态Ma=12的高超声速内转进气道,再结合不同气体模型对其进行数值模拟。研究结果表明:化学非平衡气体在流场结构、工作性能和气动加热方面与热完全气体较为相近,与热化学非平衡气体存在一定差别。离解反应发生在边界层内和低速涡流区内,热化学非平衡气体的离解反应程度比化学非平衡气体大。在隔离段内激波反射处,相比完全气体,化学反应气体的静温降低了2000～2500K。高热流区在上壁面喉道位置与下壁面激波反射点位置附近,温度较高的等温壁面、热化学非平衡气体均可降低壁面热流密度,不同壁面条件对隔离段出口性能参数影响较为明显。真实气体效应、壁面温度对隔离段涡流区的影响较为复杂,有待进一步研究。     

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

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

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

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.     

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.     

光船对气液两相介质激光支持爆轰波影响的数值模拟

建立了气液两相介质单脉冲激光支持爆轰波的数值模型,结合光船边界条件,数值模拟爆轰波及其周围流场的演化过程,分析温度、压力、速度场在爆轰波发生时的变化特性和演化形态。研究表明:有无光船边界在30μs前爆轰波发展趋势基本一致。激光作用结束后,焦点范围处温度压力均达到最大,速度继续上升达到最大后逐渐衰减。温度最大值一直在焦点附近,压力最大值则是在激光作用时出现在中心处,激光作用结束后,转移到波阵面上,而速度最大值发生在冲击波波阵面上。光船模型中,冲击波会与壁面发生耦合作用,冲击波的反射与叠加作用使最大温度中心会向外偏移,压力和速度会在壁面与冲击波边沿接触面的两侧处达到最大值,向内侧逐渐减小。      

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

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

激波/湍流边界层干扰物面剪切应力统计特性

为了揭示激波/湍流边界层干扰区内物面剪切应力统计特性的演化规律,采用直接数值模拟方法对来流马赫数2.9、12°激波角的入射激波与平板湍流边界层相互作用问题进行了研究。通过与风洞试验数据的比较分析,验证了计算结果的可靠性。系统地探究了干扰区内物面剪切应力的典型特征,如预乘谱、概率密度分布和相干结构等。研究结果表明,分离激波的低频振荡运动对流向及展向分量的预乘谱均没有实质影响,其脉动仍以高频特征为主,低频能量变化较小。干扰区内流向剪切应力概率密度函数分布变化剧烈,分离泡内对数正态分布规律不再适用,而展向剪切应力在干扰区内与正态分布较为接近。相较于上游湍流边界层,分离泡内物面剪切应力矢量夹角与幅值的联合概率密度变化显著,峰值概率降低,峰值范围增大。此外,流向剪切应力脉动场的本征正交分解分析指出,主能量模态与分离激波的低频振荡以及下游再附边界层内的Görtler-like流向涡结构密切相关。     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.