Instability theory of shock wave in a shock tube

The instability theory of shock wave in a shock tube including the effects of tube wall and contact surface is studied. The experimental data of unstable shock wave coincide with one of instability criteria derived in the present paper.     

激波的传播与干扰

激波的传播特性既取决于激波的产生条件,也与所处的传播环境密切相关。驱动条件、几何边界、介质的物理化学属性等发生变化时,都会引起激波传播特性的改变,而激波的变化反过来又会对其波及的流场产生影响。尽管激波传播及其干扰现象广泛存在于自然界和人类科技活动之中,其复杂机理的认识、规律的描述乃至应用潜力的挖掘仍有漫长的路要走。本文根据气体中激波传播和干扰现象以及与之相关的理论描述特征,在对激波传播以及反射、折射等基本现象进行简要阐述的基础上,重点围绕目前的热点问题,包括激波/激波干扰、激波/边界层干扰、激波与湍流作用、激波的聚焦与点火以及激波作用下气体界面不稳定性等研究进行了介绍和讨论,旨在对近年来该领域的进展及获得的成果做一个概述和归纳,期望对将来的深入研究有一个鉴借意义。     

Shock wave propagation and interactions

激波的传播特性既取决于激波的产生条件,也与所处的传播环境密切相关.驱动条件、几何边界、介质的物理化学属性等发生变化时,都会引起激波传播特性的改变,而激波的变化反过来又会对其波及的流场产生影响.尽管激波传播及其干扰现象广泛存在于自然界和人类科技活动之中,其复杂机理的认识、规律的描述乃至应用潜力的挖掘仍有漫长的路要走.本文根据气体中激波传播和干扰现象以及与之相关的理论描述特征,在对激波传播以及反射、折射等基本现象进行简要阐述的基础上,重点围绕目前的热点问题,包括激波/激波干扰、激波/边界层干扰、激波与湍流作用、激波的聚焦与点火以及激波作用下气体界面不稳定性等研究进行了介绍和讨论,旨在对近年来该领域的进展及获得的成果做一个概述和归纳,期望对将来的深入研究有一个鉴借意义.     

流场非均匀性对非平面激波诱导的Richtmyer-Meshkov不稳定性影响的数值研究

基于Navier-Stokes方程组,采用可压缩多介质黏性流动和湍流大涡模拟程序MVFT (multi-viscousflow and turbulence),模拟了均匀流场与初始密度呈现高斯函数分布的非均匀流场中马赫数为1.25的非平面激波加载初始扰动air/SF6界面的Richtmyer-Meshkov (RM)不稳定性现象。数值模拟结果表明,初始流场非均匀性将会影响非平面激波诱导的RM不稳定性演化过程。反射激波加载前,非平面激波导致的界面扰动振幅随着流场非均匀性增强而增大;反射激波加载后,非均匀流场与均匀流场条件下的界面扰动振幅差异有所减小。进一步,定量分析流场中环量分布及脉动速度统计量揭示了前述规律的原因。此外,还与平面激波诱导的RM不稳定性进行了简单对比,发现由于非平面激波波阵面区域的涡量与激波冲击界面时产生的涡量的共同作用,使得非平面激波与平面激波诱导的界面失稳过程存在差异。     

Nature of the surface heart transfer fluctuation in a hypersonic separated turbulent flow

This paper presents the results of an experimental study of the unsteady nature of a hypersonic separated turbulent flow. The nomimal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5×10⁷/m. The separated flow was generated using finite span forward facing steps. An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make multi-channel measurements of the fluctuating surface heat trtansfer within the separated flow. Conditional sampling analysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness. The compression waves converge into a single leading shock beyond the boundary layer. The shock structure is unsteady and undergoes large-scale motion in the streamwise direction. The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave. There exists a wide band of frequency of oscillations of the shock system. Most of the frequencies are in the range of 1–3 kHz. The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave. This intermittent phenomenon is considered as the consequence of the large-scale shock system oscillations. Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed. The project supported by National Natural Science Foundation of China     

PSE as applied to problems of secondary instability in supersonic boundary layers

Parabolized stability equations (PSE) approach is used to investigate prob-lems of secondary instability in supersonic boundary layers. The results show that the mechanism of secondary instability does work, whether the 2-D fundamental disturbance is of the first mode or second mode T-S wave. The variation of the growth rates of the 3-D sub-harmonic wave against its span-wise wave number and the amplitude of the 2-D fundamental wave is found to be similar to those found in incompressible boundary layers. But even as the amplitude of the 2-D wave is as large as the order 2%, the maximum growth rate of the 3-D sub-harmonic is still much smaller than the growth rate of the most unstable second mode 2-D T-S wave. Consequently, secondary instability is unlikely the main cause leading to transition in supersonic boundary layers.     

轴向压应力波下圆柱壳弹性动力失稳的判据与机理

基于动力失稳瞬间能量的转换和守恒,推导提出了受轴向力圆柱壳弹性动力失稳的判据和两个临界条件,由第一个临界条件导出的圆柱壳弹性动力失稳的控制方程、边界条件、屈曲变形连续与利用Hamilton原理的结果完全相同,但不足以确定包含在本问题中的两个特定特征参数（临界载荷参数和动力特征参数）,由第二个临界条件导出压缩波前的屈曲变形约束方程,基于控制方程有满足边界条件、变形连续条件和波前约束方程的非平凡解的条件,导出关于两个特征参数的一对特征方程,基于特征方程的解,精确地计算出临界载荷参数和动力特征参数的值以及动力失稳态态,由此建立了轴向力作用下圆桩壳弹性动力失稳的特征值分析方法。     

Impact on the surface of a compressible liquid

The evolution of a weak, nearly plane shock wave produced by the impact on the plane boundary of a compressible liquid is considered. At the initial moment the liquid is at rest and occupies the lower half-plane. Then the points of its boundary get instantly velocities directed into the liquid domain. This leads to the formation of a shock wave the intensity of which is non-uniform due to a non-uniform distribution of the impact velocities. Initially the shock wave is plane but then it bends due to the non-linear effects and can later be focused. To analyze the liquid flow, the method of matched asymptotic expansions is used. For finite times the flow and the evolution of the shock wave are described within the framework of the acoustic approximation. For large times the flow becomes non-linear, and the form of the shock front depends essentially on the characteristics of the liquid flow behind it. If the non-uniformity of the impact velocity distribution is slight then the focusing of the shock wave is shown not to occur. The influence of viscosity of the liquid on the structure of its motion is discussed.     

Formation of blast waves as a result of a low pressure domain decomposition

In this article the formation and propagation of blast wave as a result of a focusing of shock wave in a domain with low pressure and density are examined in the frame of an ideal gas model. We consider the decomposition of a pressure-density discontinuity on the boundary of the spherical (or cylindrical) domain which is filled by a gas whose pressure and density are lower than the pressure and density of gas that filling the external space. At initial moments of this decomposition the rarefaction wave propagates in the external space and the converging shock wave is formed in the domain of low pressure and density. The intensity of the converging shock wave gradually increases, and the wave transitions to the self-similar regime. After implosion, a diverging shock wave is formed which propagates through the disturbed gas. The values of wave intensity, wave impulse and other parameters for some magnitudes of initial parameters have been determined by means of numerical calculations. Received 10 August 1997 / Accepted 13 July 1998     

具有引力场的气体中激波传播的近似理论

本文在三维不定常流体力学特征理论的基础上,运用CCW方法,获得了激波速度的变化和激波阵面曲率、引力的大小和方向之间的解析关系,讨论了引力场对激波的加强和减弱的影响以及超新星爆炸中反弹激波的问题。利用得到的激波面积变化公式,给出了二、三维激波在引力场作用下传播问题的模拟处理方法。作为简单的例子,具体地研究了引力场变密度气体中强激波传播问题。     

气压对压力容器前壁超高速撞击损伤的影响

充气压力容器在超高速撞击下的典型损伤包括穿孔及其边缘的裂纹失稳破坏,会导致气体泄漏或爆炸,内压对容器前壁损伤的影响仍不明确。以不同内压的球形铝合金充气压力容器为研究对象,开展了球形铝合金弹丸超高速撞击实验和数值模拟计算,分析了内充气体压强对前壁穿孔形貌特征、穿孔直径、孔边环向应力等的影响规律和影响机理,讨论了气体冲击波的传播行为及影响前壁穿孔边缘裂纹失稳破坏的机制。结果表明:前壁穿孔边缘内翻边形貌与内压相关,内压越高,弯折程度越轻;穿孔直径与内充气体压强正相关,但气体对孔径的影响远小于容器壁厚及撞击速度的影响;穿孔边缘使裂纹失稳破坏的环向拉应力不仅受到后壁反射冲击波的影响,也与容器壁内应力波的传播有关,与内压成正比。     

燃气射流起始冲击波场形态的实验研究

介绍利用激光摩尔偏折仪和压阻式压力传感器对发动机喷流起始冲击波的基本形态进行的研究． 研究结果验证了在燃气流团前存在着起始冲击波及发展初期具有非球面的特性,并提出了喷口起始冲击波波阵面发展的新结构,对前人有关起始冲击波发展结构的研究结果进行了修正．     

多次激波诱导正弦扰动预混火焰界面失稳的数值研究

激波诱导火焰失稳是实际中常见的现象,为深入研究火焰失稳特性,采用三维单步化学反应的Navier-Stokes方程和9阶weighted essentially non-oscillatory (WENO)的高精度格式,对不同马赫数的入射激波及其反射激波多次诱导正弦型预混火焰界面失稳的现象进行了三维数值模拟,并对计算结果的可靠性进行了验证。研究结果显示,在激波的多次作用下,火焰界面的演变过程主要受Richtmyer-Meshkov (RM)不稳定特性和化学反应特性的双重影响,且随着入射激波强度的增强,上述2种特性均得到进一步强化。为构造体现反应性RM不稳定特性的参数,根据火焰界面混合区平均涡量和化学反应速率,提出了表征界面受不稳定性和化学反应影响的量纲一参数。通过分析发现,在同一入射激波强度下,该参数的对数形式随入射激波和反射激波的多次作用呈基本相同的线性增长趋势;对不同马赫数的入射激波,该参数对数形式的线性增长率也基本一致。这样的变化表明该量纲一参数能够反映反应性RM不稳定过程中火焰界面发展的内在规律。     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

一种三阶有限体积法及其在欠膨胀射流激波结构数值模拟中的应用

以喷管出口欠膨胀射流为研究对象,在Lagrange坐标系下建立欠膨胀射流二维积分形式的流动方程。通过在单元交接面处进行三阶ENO(essentially nonoscillatory)格式插值,构造得到一种适用于求解该方程的三阶ENO有限体积法。采用该格式对一维Sod激波管算例和喷管出口欠膨胀射流进行数值计算。计算结果表明,该方法具有高精度、基本无振荡的特点,能很好地捕捉包含激波、滑移线以及三波交点等复杂流场波系结构。计算得到的波系结构中马赫盘的位置与实验结果吻合很好,相对误差小于1.1%。     

On flows behind shock waves reflected from a solid wall

The flow fields associated with the interaction of a normal shock wave with a plane wall kept at a constant temperature were studied based on kinetic theory which can describe appropriately the shock structure and its reflection process. With the use of a difference scheme, the time developments of the distributions of the fluid dynamic quantities (velocity, temperature, pressure and number density of the gas) were obtained numerically from the BGK model of the Boltzmann equation subject to the condition of diffusive-reflection at the wall for several cases of incident Mach number:M ₁=1.2, 1.5, 2.0, 3.0, 4.0, 5.0 and 6.0. The reflection process of the shocks is shown explicitly together with the resulting formation of the flow fields as time goes on. The nonzero uniform velocity toward the wall occurring between the viscous boundary layer and the reflected shock wave is found to be fairly large, the magnitude of which is of the order of several percent of the velocity induced behind the incident shock, decreasing as the incident Mach number increases. It is also seen that a region of positive velocity (away from the wall) within the viscous boundary layer manifests itself in the immediate vicinity of the wall, which is distinct for larger incident Mach numbers. Some of the calculated density profiles are compared with available experimental data and also with numerical results based on the Navier-Stokes equations. The agreement between the three results is fairly good except in the region close to the wall, where the difference in the conditions of these studies and the inappropriateness of the Navier-Stokes equations manifest themselves greatly in the gas behavior.This article was processed using Springer-Verlag TEX Shock Waves macro package 1990.     

二维会聚波在气体内传播不稳定性的数值模拟

利用CE/SE(conservation element and solution element )格式研究了柱面会聚波在气体中传播时间断面的不稳定问题和波阵面的演变问题,并利用level set函数追踪了驱动气体与低压气体间断面的发展过程。得到了间断面的Rayleigh-Taylor(R-T)和Richtmyer-Meshkov(R-M)不稳定性发展成典型的尖钉和气泡结构的图像,初始正弦扰动下的会聚波产生尖角和尖瓣结构。结果表明,CE/SE格式在涉及会聚波的数值计算中是可行的。     

空化可压缩流动空穴溃灭激波特性研究

为深入研究空化可压缩流动中空泡/空泡团溃灭过程中激波产生、传播及其与空穴相互作用规律,本文采用数值模拟方法对空化可压缩流动空穴溃灭激波特性展开了研究.数值计算基于OpenFOAM开源程序,综合考虑蒸汽相和液相的压缩性,通过在原无相变两相可压缩求解器的控制方程中耦合模拟空化汽液相间质量交换的源项,实现了对空化流动的非定常可压缩计算.利用上述考虑汽/液相可压缩性的空化流动求解器,对周期性云状空化流动进行了数值模拟,并重点研究了空穴溃灭激波特性.结果表明:上述数值计算方法可以准确捕捉到空穴非定常演化过程及大尺度脱落空泡云团溃灭激波现象,大尺度脱落空泡云团溃灭过程分为3个阶段:(1) U型空泡团形成; (2) U型空泡团头部溃灭; (3) U型空泡团腿部溃灭.在U 型空泡团腿部溃灭瞬间,观察到激波产生,并向上游和下游传播,向上游传播的激波与空穴相互作用,导致水翼吸力面新生的附着型片状空穴回缩,直至完全溃灭.并且空穴溃灭激波存在回弹现象, 抑制了下一周期的空化发展.      

一种新的运动激波--非定常边界层相互干扰现象

本文探讨了一种新的激波－非定常边界层相互干扰现象，这种激波－边界层干扰现象既不同于定常激波－边界层干扰现象，又不同于激波在端面反射后与该激波所诱导的边界层之间的干扰现象，而是运动激波与稀疏波和第一激波所诱导的非这常边界层之间的干扰现象，本文对这种现象用微波动力学理论进行分析，并把这种干扰现象看成激波的绕射现象，同时在稀疏波破膜的双驱动激波管中进行实验观察，最后把理论分析与实验观察进行了比较。