用激波管研究超音速气固两相流

在激波管里进行了可压缩性气固两相流的实验研究。测量了激波通过颗粒群时的压力的衰减过程。用纹影仪拍摄了激波与颗粒群相互干涉的照片。试验了颗粒群的不同构造对压力衰减的影响。指出了激波反射、聚焦等非线性气动因素是可压缩性气固两相流的关键问题。     

激波与可运动颗粒群相互作用反射与透射机理实验研究

对水平圆柱形激波管内可压缩性气体与颗粒群的相互作用进行了实验研究与理论分析。利用由压力传感器、信号放大器、示波器和计算机组成的压力测量系统对激波与颗粒群作用时的动态压力进行了测量。发现激波管内发生的是一个复杂的过程,包括激波与颗粒群作用时伴随了激波和膨胀波的反射与透射现象、激波和膨胀波与接触面的干涉、以及激波从激波管端壁的反射等现象。当颗粒装载比α=1时,透射激波被直径为6mm的颗粒群衍射并且有膨胀波紧随其后,因此由透射激波引起的第一个压力峰急剧下降。颗粒装载比和颗粒直径影响透射激波的衍射和衰减。     

秦岭公路隧道2号竖井地应力与岩爆分析

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

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

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

激波驱动下固体颗粒抛撒的实验研究

本文使用阴影照相技术、高速摄像技术及压力测试手段,实验记录和研究了激波与固体颗粒群的作用及激波作用后固体颗粒群的抛撒和云团的形成过程.结果表明:在激波与固体颗粒群作用过程中,存在着清晰的激波透射、反射及绕射现象,同时激波强度在作用后有明显的下降趋势;在固体颗粒抛撒及云团形成过程中,实验发现对同一粒径的颗粒抛撒来说,抛撒的颗粒群质量越大,云团形成的均匀性及稳定性越好,而对不同粒径的颗粒群来说,粒径越大,形成的云.团集中性越强.     

收敛柱激波在粉尘气体中的传播

取稀颗粒群气固两相耦合的双流体简化模型，采用高解数值方法研究了收敛柱激波在粉尘气体中的传播和波后流场特性。通过与纯气体情况比较，揭示了固体颗粒出现对收敛柱激波传播特性的影响。     

气体—颗粒激波松弛流动特性研究

提出了一种解气体－颗粒粘性流动方程组的方法。在此方法中，计算颗粒阻力时采用Ｈｅｎｄｅｒｓｏｎ阻力系数公式，得到两相激波松弛流动数值结果，并与实验结果作了比较，详细分析和讨论了计算结果，揭示了两相激波松弛流动特性随颗粒容积比分的变化规律     

超声速双层翼翼型的阻力特性研究

以布兹曼双层翼为基础,采用基于压力梯度自适应的非结构网格求解欧拉方程的计算流体力学(CFD)方法,计算分析了双层翼翼型的厚度和翼面间距对阻力特性的影响。在马赫数为1.7的情况下,由于激波的反射和干涉,超声速双层翼翼型的阻力系数仅为0.00189,为相同厚度菱形翼型的1/15。本文通过进一步的研究发现:减少翼型厚度对于双层翼翼型设计马赫数的阻力系数有一定的影响,且与超声速状态相比,厚度对于亚声速状态的阻力影响更大,厚度减少20%,亚声速状态的阻力系数减少可达60%以上;翼面间距对阻力特性的影响相对复杂,设计马赫数之前的阻力系数与翼面间距成反比,而设计马赫数之后的阻力系数与翼面间距成正比。在此基础上,基于激波的反射及干涉效应,提出了一种双设计状态的双层翼翼型,在最佳设计点之前,双层翼之间的激波/膨胀波会有两次反射,使翼型前后的压力基本相同,阻力系数出现一次下降。随着马赫数的增加马赫角减少,激波经过一次反射就能使翼型前后的压力基本相同,使翼型达到最佳设计状态。计算结果表明,双设计状态双层翼能够使双层翼翼型在两个设计点都具有较低的阻力系数。     

地面效应影响下的高速火箭橇气动特性数值与风洞试验研究

精确预示地面效应下高速火箭橇的气动特性及流场规律对高速火箭橇的设计和评估具有重要意义。本文应用有限体积方法,研究了湍流模型对火箭橇气动特性计算精度的影响,建立了基于realizable k-ε湍流模型的火箭橇气动特性的高精度数值计算方法;结合风洞试验方法,研究了雷诺数和地面效应对高速火箭橇流场流动规律的影响,分析了火箭橇气动特性。结果表明,火箭橇阻力系数随雷诺数增大而减小,升力系数和俯仰力矩系数随雷诺数增大而增大,但雷诺数对高速火箭橇气动特性的影响较小;地面效应会使火箭橇流场发生激波-激波干扰、激波-边界层干扰和激波反复反射等复杂气动现象,大幅提升了火箭橇的升力系数和俯仰力矩系数,但对阻力系数的影响较小。研究为高速火箭橇气动外形的设计及运动稳定性的评估提供依据。     

两相介质中火焰诱导激波现象的研究

基于双流体模型，采用跟踪法处理火焰内边界，ＴＶＤ格式计算气相激波，ＭａｃＣｏ－ｍａｃｋ格式计算颗粒相流场，对两相可燃介质中火焰诱导激波现象进行了理论分析与实验模拟，研究了颗粒相参数对火焰诱导激波现象的影响。     

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

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

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

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

Unsteady reflection of a shock wave from a body with a cylindrical recess

In a number of cases of supersonic flow past bodies with recesses pulsations in the flow arise [1–3]. Experiments [4, 5] indicate that stabilization of the steady supersonic flow past the body with a recess on which a shock wave is incident takes place after a series of oscillations of the bow wave. Numerical calculation of the interaction of a supersonic jet with a cylindrical cavity [6] reveals that damped pressure pulsations arise inside the cavity if the jet is homogeneous, and undamped pulsations it is inhomogeneous. The authors explain the damping of the pulsations by the influence of artificial viscosity. This paper investigates experimentally and theoretically (by numerical methods) the oscillations of the bow shock wave and the parameters of the flow behind it in the case of unsteady reflection of a shock wave from a body with a cylindrical recess turned towards the flow. The problem is posed as follows. A plane shock wave with constant parameters impinges on a cylinder with a cavity. The unsteady flow originating from this interaction is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 199–202, September–October, 1984.     

Study of unsteady flow past bodies behind shock waves

Several theoretical and experimental studies have been devoted to the problem of the nonstationary action of the stream behind a shock wave on bodies of varied shape. In particular, in [1], the pressure and density are calculated for flow about bodies of the more typical shapes in the initial stage of the process. The basic relations which accompany the interaction of shock waves are considered in [2, 3]. The analysis of the phenomena of diffraction of shock waves on the sphere, cylinder, and cone is presented in [4]. Problems of unsteady flow about a wing are examined in [5, 6]. A detailed review of the foreign studies on unsteady flow is given in [7]. Of great practical interest is the question of the time for flow formation and the magnitudes of the unsteady loads during this period. Experimental investigations have been made recently [8, 9] in which some criteria are presented for estimating the bow shock formation time for supersonic flow about the sphere and cylinder with flat blunting. However the question of the formation time of the stationary pressure on the body surface is not referred to in these studies and no relationship is shown between the transient position of the reflected wave and the corresponding unsteady pressure on the surface. Moreover, in [8] the dimensionless time criterion is determined very approximately, independently of the Mach number of the shock wave. The present study was undertaken with the object of determining the basic criteria which characterize unsteady flow about bodies behind a plane shock wave which has time-independent parameters, and clarification of the shock wave reflected from the body and the pressure on the surface of the body during the transient period. The most typical body shapes were studied: 1) a cylinder with flat face aligned with the stream; 2) a spherically-blunted cylinder; and 3) a cylinder transverse to the stream. The experiments were conducted in a conventional shock tube using the single-diaphragm scheme. The measurements of the pressure on the models and the velocity of the incident shock wave were made using the technique analogous to that of [10, 11]. A highspeed movie camera was used to record the pattern of the wave diffraction on the body. The Mach number of the incident shock wave varied in the range from M=1.5 to M≈6.0, which corresponded to a range of Mach numbers M_∞ of the stream behind the shock wave from 0.6 to 2.1. The calculations of the required gas dynamic parameters for high temperatures were made with account for equilibrium dissociation of the air on the basis of the data of [10, 12, 13]. The magnitude of the relative maximal shock wave standoff Δ at the stagnation point obtained in the present experiments was compared with the values of Δ from other studies. In the case of the flat-blunted cylinder it was in good agreement with the results of [8–14], and in the case of the spherically-blunted cylinder and the transverse cylinder it was in agreement with the results of [15].     

Importance of unsteady contributions to force and heating for particles in compressible flows : Part 1: Modeling and analysis for shock–particle interaction

Shock–particle interaction is an important phenomenon. The interaction can be accurately resolved by direct numerical simulations. However, as the length scales of interest are much larger than the particle size in many applications, fully resolving the flow around the particle is impractical. Therefore, rigorous model for momentum and energy exchange in the interaction is very important. Shock–particle interaction is strongly time-dependent, so unsteady mechanisms play important roles in momentum and energy transfer. A model that includes unsteady contributions to force and heating is proposed. The model is used to investigate particle interactions with a planar shock wave and a spherical shock wave. The peak values and the net effects of unsteady contributions are used to measure their importance. The results show the peak values of unsteady contributions are much larger than the quasi-steady ones for a wide range of particle parameters. The net effects of unsteady contributions are important when the particle-to-gas density ratio is small. For the flow behind the spherical shock is unsteady and non-uniform, unsteady contributions have long-time influence on the particle evolution.     

Time-resolved stereo PIV measurements of shock–boundary layer interaction on a supercritical airfoil

Time-resolved stereo particle-image velocimetry (TR-SPIV) and unsteady pressure measurements are used to analyze the unsteady flow over a supercritical DRA-2303 airfoil in transonic flow. The dynamic shock wave–boundary layer interaction is one of the most essential features of this unsteady flow causing a distinct oscillation of the flow field. Results from wind-tunnel experiments with a variation of the freestream Mach number at Reynolds numbers ranging from 2.55 to 2.79 × 10⁶ are analyzed regarding the origin and nature of the unsteady shock–boundary layer interaction. Therefore, the TR-SPIV results are analyzed for three buffet flows. One flow exhibits a sinusoidal streamwise oscillation of the shock wave only due to an acoustic feedback loop formed by the shock wave and the trailing-edge noise. The other two buffet flows have been intentionally influenced by an artificial acoustic source installed downstream of the test section to investigate the behavior of the interaction to upstream-propagating disturbances generated by a defined source of noise. The results show that such upstream-propagating disturbances could be identified to be responsible for the upstream displacement of the shock wave and that the feedback loop is formed by a pulsating separation of the boundary layer dependent on the shock position and the sound pressure level at the shock position. Thereby, the pulsation of the separation could be determined to be a reaction to the shock motion and not vice versa.     

Unsteady oblique shock waves

A spatially and temporally local analysis is provided for unsteady, oblique shock waves, in which the flow is assumed to be two-dimensional or axisymmetric. Three unsteady parameters, in a laboratory frame, are viewed as the known independent variables. These are the upstream Mach number, the shock Mach number, and the angle of the shock relative to the instantaneous upstream velocity. Other steady and unsteady parameters, such as the velocity turn angles and downstream Mach numbers, are evaluated in closed form, in terms of these three quantities. Trends are assessed, and a sensitivity analysis is provided. It is suggested that the theory may find application in converting a shock capturing algorithm, at an early time during the computational process, into a shock fitting algorithm. Received 30 April 1999 / Accepted 29 November 1999     

Reconstruction of velocity fields from wall pressure measurements in a shock wave/turbulent boundary layer interaction

We present here experimental results in a shock wave/turbulent boundary layer interaction at Mach number of 2.3 impinged by an oblique shock wave, with a deflection angle of 9.5°, as installed in the supersonic wind tunnel of the IUSTI laboratory, France. For such a shock intensity, strong unsteadiness are developing inside the separated zone involving very low frequencies associated with reflected shock motions.The present work consists in simultaneous PIV velocity fields and unsteady wall pressure measurements. The wall pressure and PIV measurements were used to characterize the pressure distribution at the wall in an axial direction, and the flow field associated. These results give access for the first time to the spatial-time correlation between wall pressure and velocity in a shock wave turbulent boundary layer interaction and show the feasibility of such coupling techniques in compressible flows. Linear Stochastic Estimation (LSE) coupled with Proper Orthogonal Decomposition (POD) has been applied to these measurements, and first results are presented here, showing the ability of these techniques to reproduce both the unsteady breathing of the recirculating bubble at low frequency and the Kelvin–Helmholtz instabilities developing at moderate frequency.     

The effect of an unsteady drag force on the structure of a non-equilibrium region behind a shock wave in a gas-particle mixture

For numerical analysis of shock wave propagation in gas-particle mixtures, drag coefficients of a sphere in steady flows are generally used. However, it is shown both experimentally and numerically that a shock loaded solid sphere experiences unsteady drag forces. The paper describes a model of unsteady drag force and its effect on the structure of the non-equilibrium region behind a shock front traveling in a dusty gas. The results are compared with those obtained by using a steady drag coefficient and are discussed. It is demonstrated that the large drag force at the early stage of the interaction between shock-wave induced flow and a solid particle affects the flow structure that is obtained with a steady drag force.      

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