反射激波作用下三维凹气柱界面演化的数值研究

激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例. 单次激波与二维气柱相互作用已得到广泛关注, 但是反射激波再次冲击气柱 (尤其是三维气柱) 的研究较少, 相关演化规律和机理尚不清楚. 反射激波再次冲击演化中的气柱界面会产生新的斜压涡量, 影响涡量的输运和分布, 从而影响界面的演化. 本文采用自主开发的HOWD (high-order WENO and double-flux methods) 程序, 研究了马赫数为1.29的平面激波冲击N激波与气柱相互作用是Richtmyer-Meshkov不稳定性研究的经典案例.单次激波与二维气柱相互作用已得到广泛关注,但是反射激波再次冲击气柱(尤其是三维气柱)的研究较少,相关演化规律和机理尚不清楚.反射激波再次冲击演化中的气柱界面会产生新的斜压涡量,影响涡量的输运和分布,从而影响界面的演化.本文采用自主开发的HOWD (high-order WENO and double-flux methods)程序,研究了马赫数为1.29的平面激波冲击N_2气柱(气柱外为SF_6)的演化过程,并考察了反射激波对二维和三维凹气柱界面演化的影响规律.在数值模拟中,选取了不同的反射距离(定义为气柱和反射边界的距离),得到了二维和三维凹气柱在反射激波冲击前后的完整演化图像,提取了气柱上特征点位置随时间变化的定量数据,重点分析了不同演化阶段气柱几何特征及斜压涡量分布的变化趋势.研究表明,反射距离决定着反射激波作用气柱时的激波形状和气柱形态,从而影响斜压涡量的生成和分布,进而改变气柱的不稳定性演化过程.对于三维气柱,不同高度截面上的斜压涡量分布不同,从而诱导出复杂的三维演化结构.     

非定常激波三维双楔面反射的数值研究

采用基于MUSCL-Hancock插值的有限体积方法, 在非结构自适应网格上求解三维Euler方程, 研究了非定常激波在三维垂直双楔面上的反射现象. 研究结果表明, 由于三维效应的影响, 通过二维非定常激波反射理论预测的三维马赫反射区域范围存在着一定的局限, 在该区域范围内出现了一种非典型的三维突起结构, 即第二类三维马赫干. 另外, 对于不同的激波马赫数和楔面倾角组合, 非定常激波在三维双楔面上会形成具有四波结构的三维马赫反射或具有三波结构的三维规则反射.      

激波在不同三角楔面内传播特性的数值研究

激波在收缩管内的反射与聚焦会形成高温高压区,点燃可燃混合气并诱导爆轰,因此对爆轰发动机的点火具有重要意义。本文基于二维N-S方程,结合五阶WENO格式,对马赫数为6的正激波在三角形楔面内的反射与聚焦现象进行了数值研究。结果表明,楔面顶角的变化对激波的反射类型以及聚焦均有明显的影响:随着顶角的增加,激波的反射类型从马赫反射向过渡马赫反射和双马赫反射转变,且壁面上的前向射流更加明显;三波点第一次碰撞产生的高温高压区足够满足可燃混合气体的点火条件,且其温度与压力值随顶角的增加而增大;当激波在楔面上发生临界双马赫反射时,温度与压力达到最大;当顶角增加到一定值时,激波在楔面反射转变为常规反射,不会产生激波对碰,因而没有高温高压区。     

不同状态方程对固体介质中斜激波反射的影响

相比气体,固体介质在高压下的状态方程更为复杂,形式也多种多样.现有关于固体介质中激波反射的理论研究,一般直接采用某种状态方程,缺乏对采用不同状态方程得到的结果的对比.本项工作采用激波极曲线的理论分析方法,选择4种不同组合形式的状态方程(一次冲击激波采用线性的冲击波速度与粒子速度关系式,二次冲击激波采用Gr(u|¨)neisen状态方程;一次冲击和二次冲击激波均采用冲击波速度与粒子速度关系式:一次冲击激波采用线性冲击波速度与粒子速度关系式,二次冲击激波采用刚性气体状态方程;以及一次冲击激波和二次冲击激波均采用刚性气体状态方程),研究固体介质中的斜激波反射,比较了采用不同组合形式的状态方程对反射激波波后压力的影响.利用量纲分析方法讨论了简化状态方程达到较高精度的条件.此外,用ANSYS/LS-DYNA软件,对激波极曲线理论给出的结果进行了验证.本项工作可为固体介质中激波反射问题状态方程的选取提供一定的指导.     

反射激波冲击重气柱的RM不稳定性数值研究

数值研究了二维气柱在入射激波以及反射激波作用下的Richtmyer-Meshkov(RM)不稳定性发展规律, 采用有限体积法结合网格自 适应技术的VAS2D程序, 精确刻画激波和界面的演化. 入射平面激波的马赫数为1.2, 气柱界面内气体为六氟化硫(SF₆), 环境气体为空气, 激波管的尾端为固壁. 通过改变气柱与尾端之间的距离调节反射激波再次作用已经变形的气柱的时间, 获得不同时刻下已经变形的气柱形态、界面尺寸以及环量演化受到反射激波的影响. 结果表明, 反射激波再次作用气柱时, 气柱所处发展阶段不同, 界面演化规律以及环量随时间的变化也不相同, 反射激波与气柱相互作用过程中的涡量产生和分布与无反射情况差异较大, 揭示了不同情况下界面演化的物理机理.     

准定常强激波反射马赫杆突出变形准则的探讨

分析了滑移线延长线与楔面交点附近的流动特点,研究了准定常强激波反射中马赫杆的变形,建立了马赫杆突出变形消失的过渡准则,并用激波极曲线方法进行了求解。对马赫杆的变形过程进行了描述,研究了比热比、马赫数的变化对马赫杆突出变形消失条件的影响,并对突出变形消失区域与终点双马赫反射区域进行了比较。结果表明:马赫杆突出变形的消失是低比热比介质中出现的一种激波反射现象。马赫杆突出变形的消失,导致出现一种新的无射流、直马赫杆的双马赫反射结构。     

纳米示踪平面激光散射技术在激波复杂流场测量中的应用

在激波以及激波边界层相互作用这类含激波的复杂流场中,流场结构具有明显的三维特征.研究这类流场,采用纹影、阴影和干涉等传统流动显示技术空间分辨率较低,难以分辨流场的三维特性.基于纳米示踪的平面激光散射技术(nano-tracer planar laser scattering,NPLS),是作者近年来开发的一种新的研究超声速流场的测试与显示技术,可对超声速复杂三维流场进行高时空分辨率流动显示与测量.NPLS技术的特点使其成为测量激波复杂流场的有力手段.近年来,作者以NPLS技术为主要手段,对航空航天领域典型的激波复杂流场进行了试验研究,包括超声速弹头绕流、超声速混合层、超声速边界层,以及激波边界层相互作用流场,显示出NPLS技术在激波复杂流场精细测试与流动显示中优势.本文简要介绍NPLS技术在激波复杂流场测量中应用的研究进展.      

平面激波冲击并排液滴的三维数值模拟研究

采用三维守恒清晰界面数值方法, 研究平面激波冲击并排液滴的动力学过程. 研究的焦点在于激波接触液滴后的复杂波系结构生成, 以及并排液滴相互耦合作用诱导的单个液滴非对称界面演化. 首先, 分析并排液滴之间界面通道内的波系结构发展, 发现在冲击初期由于反射激波相交而形成新的反射激波以及马赫杆; 这些流动现象与液滴另外一侧 (非通道侧) 由激波反射所形成的弯曲波阵面截然不同, 而且所导致的液滴横向两侧流场差异是中后期冲击过程液滴两侧界面非对称演化的主要原因. 其次, 研究冲击中期时, 特别是入射激波已运动至液滴下游并远离并排液滴, 界面形态的演化过程和规律, 揭示通道下游出口处由于气流膨胀导致的界面闭合、以及随后气流阻塞导致的界面破碎等新的流动现象. 最后, 研究液滴间距对并排液滴相互作用的影响规律, 发现液滴间距大小与通道内压力峰值具有明显的关联关系. 研究表明, 更小的液滴间距不仅带来更大的压力峰值, 而且使得峰值出现的时间更早.      

浅层水中爆炸冲击波切断现象浅探

探讨了浅层水中冲击波的切断现象。主要分析了切断现象的特性、产生的机理以及对冲击波参数的影响。实验并分析了水面切断现象起始、终止时间及切断线的斜率等表征切断现象特性的参量。统计分析了水面切断现象对正压作用时间、比冲量等冲击波参数的影响。结合理论和实验数据 ,分析了有两类切断现象的波形 ,对其性质进行了判断。探讨了冲击波在水底反射产生稀疏波的机理 ,简析了水底稀疏波对波形的影响。研究表明 :装药浅层水中爆炸 ,界面对激波的削弱影响中 ,水面反射切断现象作用较强 ,水底反射切断现象的作用较弱。     

激波驱动的气固两相流的动态压力测量

利用水平圆柱形激波管对激波驱动的可压缩性气固两相流进行了试验研究.利用压电式压力传感器、电荷放大器、示波器及计算机组成的压力信号测试系统, 对激波与颗粒作用前后的气相参数进行测量及分析. 试验中测得了激波在管中的传播速度, 波后气流的压力, 反射激波、透射激波的压力和速度等. 分别考察颗粒、装载比、驱动气源以及入射激波马赫数等因素的差异对气相参数的影响.试验结果表明: 激波与颗粒群相互作用时, 会产生反射激波和透射激波,其强度与驱动气源、颗粒大小、颗粒装载比等参数有关;激波衰减率随着装载比、马赫数的增大而减小. 研究指出,在颗粒群被激波加速的初始阶段, 颗粒间的弹性碰撞起着重要的作用.      

Existence criteria for reflection configurations in shock–vortex interactions

The length scale criteria is widely accepted as an explanation for transition and hence existence of different shock wave reflection configurations in pseudo-steady flows. However, there has not been any attempt to validate this criteria using information obtained from a time-dependent numerical simulation. A high resolution time-dependent numerical simulation in pseudo-steady flow is carried out in the present work. Time-dependent numerical data is used to calculate flow features in a laboratory frame of reference to verify validity of the length scale criteria for existence of different shock wave reflection configurations in pseudo-steady flow. This analysis is then extended to the study of unsteady shock wave reflection configurations in shock–vortex interactions. It is shown that the existence of regular reflection (RR) and Mach reflection (MR) configurations in an unsteady flowfield resulting from shock–vortex interactions can also be explained locally based on limiting conditions similar to that prescribed by the length scale criteria for pseudo-steady flow.

---

     

Shock wave propagation into a surface depression

An aircraft travelling at supersonic speeds close to the ground generates a bow wave which is reflected off the ground surface. If a valley is traversed a complex reflection pattern will be generated. Similar patterns will evolve with a plane wave traversing a depression on a surface or structure. To simulate the process a planar shock wave, generated in a shock tube, is moved over several notched wedge configurations. Schlieren photographs were produced to assist in identifying the resulting complex three-dimensional wave structures and then verified and extended by three- dimensional computations. The valley geometries investigated are rectangular, triangular, parabolic and conical with a number of valley floor inclinations. The main features are extracted in surface models to demonstrate the complexity of the flow, and in particular in the case where the reflection is regular on the ground plane and Mach reflection in the valley.      

Numerical confirmation of the hysteresis phenomenon in the regular to the Mach reflection transition in steady flows

Numerical calculations based on the Navier-Stokes equations are carried out to investigate the reflection of shock waves over straight reflecting surfaces in steady flows. The results for a flow Mach number of M₀=4.96 confirm the recent experimental findings of Chpoun et al. (1995) concerning the transition from regular to Mach reflection. Numerical calculations as well as experimental results show a hysteresis phenomenon during this transition and the regular reflection is found to be stable in the dual-solution domain in which theoretically both regular and Mach reflection wave configurations are possible.     

Numerical investigations of transition between regular and Mach reflections caused by free-stream disturbances

Numerical simulations have been performed to study the influence of the free-stream disturbances on the alternation of the steady shock wave reflection configurations in the dual solution domain. Different types of disturbances have been considered. The analysis of interaction between disturbances and the incident shock wave can be substantially simplified for the localized density disturbances. It is shown that such disturbances can indeed cause the transition from regular reflection to Mach reflection and back, so that within a certain range of angles of incidence the shock wave reflection configuration can be considered as a bi-stable system. The threshold amplitude of the localized density disturbance, able to induce the transition, has been estimated theoretically. The results of numerical computations convince of higher stability of the Mach reflection in the dual solution domain compared to the regular reflection, which is in accordance with available experimental data. Received 10 May 2001 / Accepted 15 November 2001 Published online 8 July 2002     

Nonstationary shock reflection over nonstraight surfaces: an approach with a method of multiple steps

Shock reflection phenomena over nonstraight surfaces have been investigated. The models used in this experiment are ordinary circular cylindrical concave and convex wedges and step-like wedges which simulate the former. The step-like wedges were used to investigate the process of reflected-wave formation over circular cylindrical wedges (method of multiple steps). The reflected-wave structure has been photographed with a schlieren apparatus. The formation of the reflected wave over circular cylindrical wedges is physically well understood by comparing it with shock reflection over step-like wedges. In particular, the reason why the reflected wave over a concave circular cylindrical wedge is very weak away from the reflection point is elucidated. Moreover, the structure and the formation mechanism of the so-called transitioned regular reflection (TRR) are illustrated in detail. As a by-product, based on acoustic theory, analytical formulae for the transition wedge angle are found. They are in good agreement with experiments. Received Received 28 February 1996 / Accepted 7 August 1996     

The phenomena of shock wave reflection — a review of unsolved problems and future research needs

Although the phenomenon of shock wave reflection was discovered more than a hundred years ago, active research related to this phenomenon still goes on in many countries in the world (e.g., Australia, Canada, China, Germany, Israel, Japan, Poland, Russia and United States of America). As a matter of fact the research activity increased so drastically in the past decade and a half that a special scientific meeting dedicated to better understanding the reflection phenomena of shock waves, namely The International Mach Reflection Symposium was initiated in 1981 and was held since then in the major research centers actively involved in the research of shock wave reflections. In the present paper the status of the research of the phenomenon of shock wave reflection will be discussed in general, and unresolved problems and future research needs will be pointed out.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Effects of a Single-pulse Energy Deposition on Steady Shock Wave Reflection

The effects of energy deposition in the free stream on steady regular and Mach shock wave reflections are studied numerically. A short-duration laser pulse is focused upstream of the incident shock waves. It causes formation of the expanding blast wave and the residual hot-spot interacting in a complex way with the steady shock wave reflection. It was found that the laser energy addition in the free stream may force the transition from regular to Mach reflection in the dual solution domain. In contrast to previously reported numerical results, the transition from Mach to regular reflection has not been reproduced in our refined computations since the Mach reflection is restored after the flow perturbation.     

Analytical and numerical investigations of the reflection of asymmetric nonstationary shock waves

The reflection of asymmetric nonstationary shock waves is analytically and numerically studied in this paper. An analytical approach, which is a combination of the shock dynamic and shock polar methods, is advanced to predict the reflection wave configurations. The numerical simulations are performed by the finite volume method based on the second-order MUSCL-Hancock scheme and the HLLC approximate Riemann solver, with the self-adaptive unstructured mesh. It is found that the transition between the overall regular reflection and overall Mach reflection in the asymmetric nonstationary reflection agrees with the detachment criterion, which is analogous to the reflection in pseudo-steady flows (i.e. shock reflection over a wedge). Some special reflection wave configurations, which have never been observed in steady or nonstationary shock reflections so far, are found to exist in this asymmetric reflection. Furthermore, the domains and boundaries of various overall reflection wave configurations are analytically predicted, and the effect of mis-synchronization is also discussed.     

Steady shock wave reflections in thermochemical nonequilibrium flows

The shock wave reflection phenomena in hypersonic steady air flows, including thermochemical nonequilibrium effects, are investigated. The main objectives are to study the influence of these effects on the two shock wave reflections (regular and Mach reflections), on the Mach stem height and on the hysteresis behavior. The air computations are performed using a multi-block MUSCL-TVD finite-volume scheme. The computational results with and without thermochemical effects in the air mixture flow at an upstream Mach number equal to 7 are compared. The comparison reveals a strong dependence of the transition angles, of the height and location of the Mach stem on the physical modeling of the gas flow. Received 17 February 2000 / Accepted 30 August 2000