径向入射激波聚焦实验和数值模拟研究

为得到不同入射激波马赫数对凹腔内径向入射激波聚焦过程的影响,对马赫数为1.25和2.18的径向入射激波在实验段内反射聚焦的过程进行了实验和数值模拟研究。用高速CCD拍摄了凹腔中气流流场的纹影图片,并采用保持强稳定性的龙格—库塔格式、加权基本无振荡算法和块结构动态自适应网格加密对激波反射聚焦的过程进行了数值模拟,实验和数值模拟结果吻合较好。复杂的激波形状及其之间的相互作用(如激波反射,激波聚焦和不同激波的相互作用)在实验和数值模拟中都能很好地观测到。通过比较不同马赫数入射激波的反射聚焦过程,发现较高的马赫数下聚焦后能产生射流,而较低的马赫数下却没有。同时发现,入射激波马赫数能影响二次激波出现的位置和强度,进而影响聚焦过程流场的紊乱程度。     

轴向入射激波反射聚焦的实验和数值模拟

为研究两级脉冲爆震发动机凹腔入口宽度对第2级凹腔中激波聚焦过程的影响，设计轴向入射激波聚焦的实验系统，对不同入口宽度下激波反射聚焦的过程进行了实验和数值模拟。采用纹影系统拍摄凹腔中的流场结构，并测量了凹腔顶点处的动态压力。采用保持强稳定性（SSP）的Runge Kutta格式、weighed essential non oscillation（WENO)格式和块结构网格自适应加密（SAMR）算法对激波反射、衍射和聚焦的过程进行了数值模拟，数值模拟结果和实验吻合较好。通过分析比较相同马赫数激波通过不同宽度凹腔的反射聚焦过程，发现较小宽度凹腔入口下激波反射和聚焦较强，顶点的压力峰值较高，激波聚焦诱导的射流较强。      

凹面腔内激波聚焦起爆爆震波过程的数值模拟

运用CFD方法和基元反应机理对环形向心射流产生的激波在凹面腔内的反射聚焦起爆爆震波过程进行了数值模拟,并根据流场分布及变化情况分析了激波聚焦及其在H2-Air混合物中起爆爆震波的全过程。对于本文建立的模型,起爆点在抛物形壁面的底部顶点处,聚焦起爆爆震波后的瞬间压力达到21.3MPa,温度达到4540K。爆震波在凹面腔内向开口端传播过程中仍会出现两次聚焦,压力达到约18MPa,温度达到4000K左右。研究结果表明:利用环形向心射流产生的激波在凹面腔内聚焦可以成功直接起爆爆震波,是一种有效的爆震直接起爆方法。     

射流马赫数对圆弧形凹面腔内激波聚焦过程的影响

为研究射流马赫数对连续超声速射流对撞流场演化及激波聚焦过程的影响,保持入射导流深度L=0.00mm不变,通过更换不同的Laval喷管,对马赫数分别为Ma=1.2,1.4,1.6和1.8时的实验工况凹面腔内反射聚焦过程进行了实验研究,用高速CCD(Charge coupled device)拍摄了圆弧形凹面腔中气流流场纹影照片,并用动态压力传感器测量了聚焦过程中流场的压力变化,对径向入射激波在凹面腔内的反射聚焦过程进行了描述。通过对比不同射流马赫数下激波反射聚焦过程,发现在低马赫数1.2时,表现出较强的激波完全聚焦特性,即前导激波碰撞形成反射激波,并反射聚焦形成三波点,从而在凹腔底部形成高温高压区触发爆震,前导激波完全聚焦过程在凹面腔内流场演化中占据主导地位。随着马赫数的增加,完全聚焦强度降低,在流场中的主导优势逐渐减弱;其激波聚焦频率受射流马赫数的影响较小,频率差值较小,基本保持一致。     

导流块深度对抛物形凹面腔中激波反射聚焦的影响

为研究导流块深度对抛物形凹面腔内径向入射激波聚焦过程的影响,对马赫数为1.41的径向入射激波在导流块深度分别为0,5,10,15mm的凹面腔内反射聚焦过程进行了实验研究。结合高速CCD拍摄到的凹面腔中气流流场纹影照片和动态压力传感器测得的聚焦过程中流场的压力变化,对径向入射激波在凹面腔内的反射聚焦过程进行了描述。通过比较不同导流块深度下激波反射聚焦过程,发现随着导流块深度的加深,前导激波聚焦和反射激波聚焦的时间差会减小,使激波聚焦的强度增大,当导流块深度从0mm增加到15mm时,激波聚焦所致的峰值压力从0.39MPa增加到0.51MPa。但是,随着导流块深度的加深,压力增益下降,并使排气过程的难度加大,因此导流块深度为10mm左右时能取得较好的聚焦效果。     

二维凹面腔内暂冲式激波聚焦机理的实验

针对两级脉冲爆震发动机中第二级凹面腔内激波聚焦起爆能量和周期的研究,设计了二维凹面腔内暂冲式激波聚焦实验系统。通过分析粒子图像测速系统捕捉的照片和动态压力数据,获得了凹面腔内激波聚焦高压区压力特性和多循环流场演化特征。并且开展了A和B两种型面凹面腔的对比实验,进一步分析了多循环激波聚焦的形成机理。结果表明:入射激波第一次在凹面腔底部聚焦后形成的聚焦反射激波在喷口射流边界上发生反射,产生的反射激波再次在凹面腔底部聚焦,形成了第二次激波聚焦,如此往复便形成了周期性的激波聚焦,聚焦频率高达7~10kHz;并且,速度场和射流强度随入射激波碰撞而减弱,在排气过程中又增强,呈现周期性变化。      

斜激波与内凹半圆柱面无粘相互作用

针对三维内转式进气道流动中激波入射内凹壁面主导的复杂波系干扰问题，采用斜激波入射内凹半圆柱面的简化构型，在来流马赫数为6的条件下，通过无粘数值模拟结合理论分析，研究了激波角βi=14°～29°的斜激波在内凹半圆柱面反射形成的三维流场。结果表明，流场对称面均会出现显著高于二维情况的逆压梯度。当βi≤25°时，从侧壁到对称面，斜激波经历了从马赫反射（MR）到规则反射（RR）的转变，形成了MR-RR型流场，转变点处产生的桥激波向对称面延伸，桥激波在对称面反射后产生的压力峰值高于二维斜激波反射；当βi≥25°时，斜激波在侧壁和对称面均发生马赫反射，形成了MR-MR型流场，两种马赫反射分界点处产生的桥激波向侧壁发展，侧壁气流在对称面相撞后产生的压力峰值高于正激波后的压力；当βi=25°时，流场存在MR-RR型和MR-MR型双解现象。通过降维分析理论，揭示了两类流场中转变点和分界点的形成机制，并厘清了桥激波的产生原因和初期演化特征。当βi≥18°时，无粘激波干扰所主导的侧壁气流加剧向对称面汇聚，并在对称面附近产生流向涡对。无粘分析获得的认识，有助于揭示内转式进气道中流动汇聚和流向涡对等现象的形成机...     

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

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

分隔式射流通道对凹面腔内爆震起爆与反传的影响

针对凹面腔内激波聚焦起爆的爆震波反传现象,采用燃料与氧化剂以分隔式通道进入凹面腔的非预混燃料入射方案,通过数值计算,主要探讨了射流入射通道截面积、各截面积比和通道数量对激波聚焦起爆的影响作用。计算结果表明:燃料/氧化剂射流入射通道截面积比为化学恰当比时,凹面腔底部的掺混程度较好,有利于点火起爆;入射通道截面积比为化学恰当比的多入射通道方案,其凹面腔底部激波聚焦起爆时温度和压力更高,相同时刻的燃烧流场其温度和压力峰值均较高,更利于增强凹面腔内的掺混程度,使燃料浓度分布均匀;小截面积比的射流入射通道,有助于增大波阻提高抑制压力反传效果,减小压力反传强度。     

耦合传热下激波对超声速气膜冷却影响

针对离散孔式超声速平板气膜冷却，在主流区引入楔角形成激波环境，以研究激波与超声速气膜之间的相互作用。通过计算楔角在0°、15°、20°和25°产生的四种激波强度下，超声速气膜与高温壁面的耦合传热。所得结果表明：适当强度的激波能够抑制气膜入射后产生的反向涡旋对，降低主流对气膜的卷吸，增大壁面平均H2摩尔分数并降低壁面温度。对金属层温度场的分析表明，壁面冷却效果随着激波角的增加而先增加后降低，其中楔角为20°时的流场结构最有利于壁面温度保护。小楔角生成的激波在低冷流马赫数下对冷却效果的改善更明显，大楔角则在高冷流马赫数下更明显，热障涂层（TBC）不影响这种变化趋势；激波的存在削弱了TBC的影响范围。可以揭示超声速气膜在耦合传热条件下的传热机理，为超声速气膜冷却的设计提供参考，或为现有超声速气膜冷却结构的优化提供依据。     

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.     

激波对凹腔火焰稳定器流场影响

为揭示激波对超声速流中凹腔流场的影响规律,在超声速流中利用长9 mm,坡度为23°的斜坡产生激波。组合利用高速摄影仪和纹影仪,采用碘钨灯、连续激光和脉冲激光作为光源,摄像曝光时间分别为1 m s,8 us,8 ns,获得了斜坡激波入射在单个凹腔中部、中后部和后部,以及激波入射在不同构型凹腔相同位置的流场纹影图像。结果表明,凹腔自由剪切层在受到激波撞击后,在入射点前会发生弯曲而偏向主流,在入射点后加速发展和破碎;激波入射点靠近前沿则引起凹腔前沿激波的增强,靠近凹腔后沿导致再附激波减弱。     

流向涡与斜激波相互作用的试验研究和数值模拟

为研究流向涡与斜激波相互作用在超声速燃烧中的应用,进行了由翼产生的流向涡与楔块产生的二维斜激波相互作用的燃烧室冷流试验研究.在不同激波强度下,纹影仪捕捉到了强、中、弱不同的涡/波作用现象.仿真与试验结果符合得较好.试验结果表明:对马赫数为2.3流场,在翼攻角10°时,能产生强流向涡,此工况下,锲角越大,涡/波作用越强.仿真结果表明:马赫数对涡/波作用影响较大,总压影响不明显,总温可影响亚声速回流区的尺寸.      

A theoretical method for solving shock relations coupled with chemical equilibrium and its applications

In this study, a theoretical method is proposed to solve shock relations coupled with chemical equilibrium. Not only shock waves in dissociated flows but also detonation waves in combustive mixtures can be solved. The global iterative solving process is specially designed to mimic the physical and chemical process in reactive shock waves to ensure good stability and fast convergence in the proposed method. Within each global step, the single-variable equations of normal and oblique shock relatio...     

面向三维激波问题的装配方法

给出了一种基于非结构动网格技术的三维激波装配方法。在该方法中,三维激波面由被标记为激波属性的网格点连接构成,标记为激波属性的网格点称为激波点。激波点具有两组参数分别代表激波的上下游,利用激波点上下游参数求解R-H关系式获得激波点运动速度。非结构动网格技术的使用允许激波大幅度运动,降低了对初始激波位置的要求。通过引入网格属性定义避免了对计算网格进行分区,增加了装配激波的灵活性。通过球柱体绕流问题验证了该三维装配方法的合理性,针对三维激波装配中比较困难的交点装配问题,通过对三维激波反射以及三维激波相交等算例进行研究找到了可用的三维激波交点运动速度的确定方法,保证了激波运动过程中交点运动与流场求解之间的相容性,获得了相应的装配结果。     

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.     

Ion acceleration at shocks in interplanetary space: A brief review of recent observations

The mechanism by which ions are accelerated near the Earth's bow shock and near shocks propagating outward from the Sun in response to solar activity appears to be essentially the same. For both types of shock the solar wind thermal distribution acts as a seed population. Leaked magnetospheric ions and resident flare ions are additional seed populations for the bow shock and outward propagating shocks respectively. The acceleration of solar wind ions at these shocks begins with either the reflection of ions off the shock or leakage of shocked plasma back through the shock. Interaction with a disruption wave field self-generated by these backstreaming ions is responsible for the remainder of the acceleration at the bow shock. Both the disruption wave field and the ambient interplanetary wave field play important roles in accelerating ions at outward propagating shocks, but on different time scales. The geometry of the shock and the duration of field line connection to the shock play decisive roles in determining what is observed.     

尖劈诱导的斜爆轰波的精细结构及其影响因素

为了揭示斜爆轰的精细结构及其影响因素,基于带化学反应的Euler方程,利用五阶WENO格式,对尖劈诱导的斜爆轰进行了数值模拟。结果表明,在劈面的压缩下,斜爆轰波阵面被分为三种结构,依次是类ZND结构、单三波点阵面结构和双三波点阵面结构。同时,根据三波点的轨迹绘制了胞格结构,其中单三波点的轨迹为平行直线,而双三波点的轨迹为倾斜的蜂窝状结构。讨论了来流马赫数和尖劈角度对劈面压缩性的影响,随着马赫数和劈角的增加,尖劈的压缩更厉害,爆轰波阵面也就越光滑。