Propagation of a pressure wave in two-phase flow with very high void fraction

An experimental and theoretical study of a finite amplitude pressure wave propagating through a two-phase media of about 0.9999–0.99999 void fraction is performed. This two-phase media consists of many parallel liquid films in a gas. The films are perpendicular to the wave propagation direction and result in a two-phase fluid of extremely high void fraction. Experiments are done in a vertical shock tube and show that the shock wave is broken down into an initial sharply rising wave and a second gradually rising wave. The velocity of the first wave agrees well with the theoretical prediction assuming an adiabatic thermal equilibrium change, which approaches the gas sonic velocity in the two-phase flow in the low void fraction region. The second wave is caused by the complex reflection and destruction of the waves.     

Analytical Solution for Isothermal Flow in a Shock Tube Containing Rigid Granular Material

Analytical solution of shock wave propagation in pure gas in a shock tube is usually addressed in gas dynamics. However, such a solution for granular media is complex due to the inclusion of parameters relating to particles configuration within the medium, which affect the balance equations. In this article, an analytical solution for isothermal shock wave propagation in an isotropic homogenous rigid granular material is presented, and a closed-form solution is obtained for the case of weak shock waves. Fluid mass and momentum equations are first written in wave and (mathematical) non-conservation forms. Afterwards by redefining the sound speed of the gas flowing inside the pores, an analytical solution is obtained using the classical method of characteristics, followed by Taylor’s series expansion based on the assumption of weak flow which finally led to explicit functions for velocity, density and pressure. The solution enables plotting gas velocity, density and pressure variations in the porous medium, which is of high interest in the design of granular shock isolators.     

Interaction of a shock wave with a loose dusty bulk layer

The interaction of a planar shock wave with a loose dusty bulk layer has been investigated both experimentally and numerically. Experiments were conducted in a shock tube. The incident shock wave velocity and particle diameters were measured with the use of pressure transducers and a Malvern particle sizer, respectively. The flow fields, induced by shock waves, of both gas and granular phase were visualized by means of shadowgraphs and pulsed X-ray radiography with trace particles added. In addition, a two-phase model for granular flow presented by Gidaspow is introduced and is extended to describe such a complex phenomenon. Based on the kinetic theory, such a two-phase model has the advantage of being able to clarify many physical concepts, like particulate viscosity, granular conductivity and solid pressure, and deduce the correlative constitutive equations of the solid phase. The AUSM scheme was employed for the numerical calculation. The flow field behind the shock wave was displayed numerically and agrees well with our corresponding experimental results.      

Breakup of gas bubbles in a liquid by shock waves

The nature of the propagation of shock waves in various media is related to the characteristics of the latter, including their compressibility, thermophysical properties, the presence of multiple phases, etc. The structure of a shock wave varies appreciably as a function of the properties of the medium. The most significant property of a liquid mixture with gas bubbles is the compressibility of the latter under the influence of an externally applied pressure, for example, in a shock wave propagating in the liquid—gas medium. The transfer of momentum and energy between phases and the pressure variation behind the wave depends on the behavior of the gas bubbles behind the shock front.     

Dynamics of shock waves in a liquid containing gas bubbles

Results are presented of a numerical solution of the Korteweg-de Vries-Burgers equation that describes the propagation and establishment process for a stationary structure to a shock wave in a gas-liquid medium. Data are obtained on the time for the establishment of a stationary structure of a shock wave, propagation velocity, and amplitude oscillations in the front of the shock wave. Experiments are discussed on the basis of the results obtained for the study of shock waves in a liquid containing gas bubbles.     

Pressure waves in a separated gas-liquid layer in a horizontal duct with a step

Pressure wave propagation into a separated gas-liquid layer in a horizontal duct with a step is investigated analytically. The linear solution is derived assuming a large density ratio of liquid to gas. The solution can be found first for the gas layer and then for the liquid layer. The linear wave in a liquid layer is valid even for fairly large initial pressure ratios, and clearly exhibits the dispersive characteristics of the pressure wave in a liquid layer. As the initial pressure ratio is increased, the pressure wave in the gas layer becomes a shock wave. Thus, its effect on the wave in a liquid layer can be found analytically by modifying the boundary condition in part. The wave in a liquid layer consists of a main wave, which propagates with the shock speed in gas, and a precursor wave, whose front propagates with the speed of sound in liquid. The precursor wave has an oscillatory structure; its amplitude increases with increasing shock strength and also with liquid layer thickness.     

颗粒材料中致密波结构研究

采用一维两相流模型与相应颗粒构形应力函数,研究了致密波的形成及其结构．用简化两相流模型系统地讨论致密波对有关因素的依赖关系．分析指出：小于基体材料音速的致密波仅能在非理想颗粒材料中存在,从波前到波后,所有状态物理量光滑过渡．大于基体材料音速的致密波,波头可能存在间断．应力函数与致密粘性确定后,致密波速度决定致密波结构、宽度、终态压实度．采用一维两相流模型模拟了活塞驱动颗粒床形成致密波这一动态过程．用线方法（MOL）对该方程组求数值解．计算表明,经过短暂的非稳态过程,颗粒床中形成一稳态致密波．分析了活塞速度与初始孔隙率对致密波结构的影响,并对简化两相流模型与两相流模型的计算结果进行了对比．     

Contributions to numerical developments in shock waves attenuation in porous filters

The paper deals with the numerical method of the compressible gas flow through a porous filter emphasizing the treatment of the interface between a pure gaseous phase and a solid phase. An incident shock wave is initiated in the gaseous phase interacting with a porous filter inducing a transmitted and a reflected wave. To take into account the discontinuity jump in the porosity between the gaseous phase and the porous filter, an approximate Riemann solver is used to compute homogeneous non-conservative Euler equations in porous media using ideal gas state law. The discretization of this problem is based on a finite volume method where the fluxes are evaluated by a “volumes finis Roe” (VFRoe) scheme. A stationary solution is determined with a continuous variable porosity in order to test the numerical scheme. Numerical results are compared with the two-phase shock tube experiments and simulations of a shock wave attenuation and gas filtration in porous filters are presented.      

Formation of shock waves in gas-liquid foams

The purpose of the present investigation is to analyze the phenomenon of shock wave formation in gas-liquid foams and to explain the qualitative differences which are found when comparing results from shock tube experiments performed with foams and bubbly liquids. It is well known that oscillatory pressure waves in bubbly liquids may reach an amplitude twice as large as that of the original pressure impulse. However, experiments showed that pressure disturbances in foams always attenuate without significant change in the wave pressure profile. In the present study this behavior is explained by analyzing shock wave formation using the Burgers equation which is derived from the conservation laws for a bubbly liquid. It is shown that the parameter of non linearity in the Burgers equation describing wave propagation in bubbly liquids is about 40 times higher than in foams. At the same time coefficient of bulk viscosity of a foam is about 10³ times greater than that of a bubbly liquid. This explains why in shock tube experiments with foams shock waves are not detected while they are easily observed when bubbly liquids are used under similar conditions.     

Blast wave propagation in air from a gaseous charge

In the point explosion problem it is assumed that an instantaneous release of finite energy causing shock wave propagation in the ambient gas occurs at a space point. The results of the solution of the problem of such blasts are contained in [1–4]. This point model is applied for the determination of shock wave parameters when the initial pressure in a sphere of finite radius exceeds the ambient air pressure by 2–3 orders of magnitude. The possibility of such a flow simulation at a certain distance from the charge is shown in papers [4, 5] as applied to the blast of a charge of condensed explosive and in [6, 8] as applied to the expansion of a finite volume of strongly compressed hot gas. In certain practical problems the initial pressure in a volume of finite dimensions exceeds atmospheric pressure by a factor 10–15 only. Such cases arise, for example, in the detonation of gaseous fuel-air mixtures. The present paper considers the problem of shock wave propagation in air, caused by explosion of gaseous charge of spherical or cylindrical shape. A numerical solution is obtained in a range of values of the specific energy of the charge characteristic for fuel-air detonation mixtures by means of the method of characteristics without secondary shock wave separation. The influence of the initial conditions of the gas charge explosion (specific energy, nature of initiation, and others) is investigated and compared with the point case with respect to the pressure difference across the shock wave and the positive overpressure pulse.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–118, May–June, 1986.     

当量比对汽油燃料两相旋转爆轰发动机工作特性影响实验研究

为了研究当量比对汽油燃料两相旋转爆轰发动机工作特性的影响,开展了以高总温空气为氧化剂的气液两相旋转爆轰实验研究。旋转爆轰发动机环形燃烧室外径、内径和长度分别为202、166和155 mm。汽油和高温空气采用高压雾化喷嘴与环缝对撞喷注的方式混合,以此提高推进剂的掺混效果与活性,采用预爆轰管作为点火装置。实验通过改变汽油质量流量改变推进剂当量比,并基于燃烧室内测得的高频动态压力和平均静压,对气液两相旋转爆轰波的传播模态和传播特性以及发动机的工作特性进行了详细分析。实验结果表明:在当量比为0.79～1.25时,燃烧室内均实现了旋转爆轰波的连续自持传播,且随着当量比的增加,爆轰波传播模态从双波对撞/单波的混合模态转变为单波模态;降低当量比至0.61～0.66,爆轰波传播稳定性变差,传播模态表现为间断爆轰以及零星爆轰;进一步降低当量比至0.53,爆轰波起爆失败。此外,燃烧室平均绝对压力与爆轰波平均传播频率均随着当量比的增加呈先增大后减小的趋势,极大值出现在当量比1.19附近。在此工况下获得了最佳实验结果,旋转爆轰波的平均传播频率为1 900.9 Hz,平均传播速度为1 110.8 m/s,与高频压力信号经快速傅里叶变换得到的主频基本一致,爆轰波传播速度存在严重亏损。     

瓦斯/空气预混气体爆燃流场测试系统多参数耦合同步控制实验方法

为了更精确地获得爆炸激波管内瓦斯/空气预混气体爆燃过程中，激波形成过程、压力和火焰传播速度以及火焰与惰性阻燃剂相互作用的流场演化图像。通过分析激波管测试系统中多个目标的时间响应特征及控制方式，利用超高速相机、光电倍增管、时间延时器、固态继电器、电荷放大器和数据采集系统等设备，设计实验方案，分别对激波管中瓦斯/空气预混气体爆燃高压点火系统的响应时间和惰性介质阻燃剂喷射系统的响应时间进行测试。实验结果表明电火花点火的响应时间为微秒量级，而阻燃剂喷射系统的响应时间为毫秒量级，以响应时间为依据，通过设置精确的延迟时间实现多目标同步控制，为完成激波管内瓦斯/空气预混气体爆燃过程的微观流场显示奠定基础。     

Wave propagation in an elastic tube: a numerical study

In this paper, a fluid–wall interaction model, called the elastic tube model, is introduced to investigate wave propagation in an elastic tube and the effects of different parameters. The unsteady flow was assumed to be laminar, Newtonian and incompressible, and the vessel wall to be linear-elastic, isotropic and incompressible. A fluid–wall interaction scheme is constructed using a finite element method. The results demonstrate that the elastic tube plays an important role in wave propagation. It is shown that there is a time delay between the velocity waveforms at two different locations and that the peak velocity increases while the low velocity decreases in the elastic tube model, contrary to the rigid tube model where velocity waveforms overlap each other. Compared with the elastic tube model, the increase of the wall thickness makes wave propagation faster and the time delay cannot be observed clearly, however, the velocity amplitude is reduced slightly due to the decrease of the internal radius. The fluid–wall interaction model simulates wave propagation successfully and can be extended to study other mechanical properties considering complicated geometrical and material factors.     

Choking and shock phenomena in a single component two-phase flow including vibrational effects

An analytical and experimental investigation including vibratory effects of flashing flow in a tube with a sharp edged entrance is presented. A free streamline flow model is applied to predict choking in single-component two-phase flow. By identifying three separate regimes (i.e. jet flow, two-phase homogeneous flow, and single-phase liquid flow) in the flashing flow system, an expression is obtained for the prediction of the minimum stagnation pressure loss under choked flow conditions. A normal shock located between the flashing two-phase mixture and the single-phase liquid was experimentally observed. The location of the shock is predicted as a function of the stagnation pressure drop across the tube. The analytical predictions are verified by experimental data.     

Nonlinear vibrations of a vapor film in the presence of intense heat fluxes

Waves propagating along the interface between a thin vapor film and a liquid layer in the presence of a heat flux are investigated. The boundary conditions on the vapor-liquid phase surface take into account the temperature dependence of the pressure and the possibilities of formation of the metastable state of the superheated liquid and mass flow. Variations in the saturation pressure as functions of the temperature and mass flux lead to generation of weakly damped periodic waves of low amplitude whose velocity can be much higher than the velocity of the gravity waves. The waves ensure stability of the vapor film beneath the liquid layer in the gravity field. The finite-amplitude waves on the surface of the vapor film differ from the Stokes surface waves on the free surface of isothermal fluid. Instability regimes related with superheating of the liquid ant its explosive boiling when the amplitude of an initially small wave increases to infinity in a finite time can develop in a certain working-parameter regime.     

Effect of diameter on two-phase pressure drop in narrow tubes

The effect of tube diameter on two-phase frictional pressure drop was investigated in circular tubes with inner diameters of 0.6, 1.2, 1.7, 2.6 and 3.4 mm using air and water. The gas and liquid superficial velocity ranges were 0.01-50 m/s and 0.01-3 m/s, respectively. The gas and liquid flow rates were measured and the two-phase flow pattern images were recorded using high-speed CMOS camera. Unique flow patterns were observed for smaller tube diameters. Pressure drop was measured and compared with various existing models such as homogeneous model and Lockhart-Martinelli model. It appears that the dominant effect of surface tension shrinking the flow stratification in the annular regime is important. It was found that existing models are inadequate in predicting the pressure drop for all the flow regimes visualized. Based on the analysis of present experimental frictional pressure drop data a correlation is proposed for predicting Chisholm parameter “C” in slug annular flow pattern. For all other flow regimes Chisholm’s original correlation appears to be adequate except the bubbly flow regime where homogeneous model works well. The modification results in overall mean deviation of pressure drop within 25% for all tube diameters considered. This approach of flow regime based modification of liquid gas interaction parameter appears to be the key to pressure drop prediction in narrow tubes.     

Enhancement of shock waves in a porous medium saturated with a liquid containing soluble-gas bubbles

Evolution of a moderate-intensity shock wave and its enhancement after reflection from a rigid surface embedded in a porous medium are studied experimentally. The medium is saturated with a liquid that has bubbles of a soluble gas. A physical mechanism of shock wave enhancement in a saturated porous medium is proposed. Experimental data on the amplitude and velocity of reflected waves are compared with results of theoretical modeling. The process of gas bubble dissolution behind a shock wave is studied.     

Fundamentals of rotating detonations

A rotating detonation propagating at nearly Chapman–Jouguet velocity is numerically stabilized on a two-dimensional simple chemistry flow model. Under purely axial injection of a combustible mixture from the head end of a toroidal section of coaxial cylinders, the rotating detonation is proven to give no average angular momentum at any cross section, giving an axial flow. The detonation wavelet connected with an oblique shock wave ensuing to the downstream has a feature of unconfined detonation, causing a deficit in its propagation velocity. Due to Kelvin–Helmholtz instability existing on the interface of an injected combustible, unburnt gas pockets are formed to enter the junction between the detonation and oblique shock waves, generating strong explosions propagating to both directions. Calculated specific impulse is as high as 4,700 s.      

Shock structure in bubbly liquids: comparison of direct numerical simulations and model equations

Shock wave structure in a bubbly mixture composed of a cluster of gas bubbles in a quiescent liquid with initial void fractions around 10% inside a 3D rectangular domain excited by a sudden increase in the pressure at one boundary is investigated using the front tracking/finite volume method. The effects of bubble/bubble interactions and bubble deformations are, therefore, investigated for further modeling. The liquid is taken to be incompressible while the bubbles are assumed to be compressible. The gas pressure inside the bubbles is taken uniform and is assumed to vary isothermally. Results obtained for the pressure distribution at different locations along the direction of propagation show the characteristics of one-dimensional unsteady shock propagation evolving towards steady-state. The steady-state shock structures obtained by the present direct numerical simulations, which show a transition from A-type to C-type steady-state shock structures, are compared with those obtained by the classical Rayleigh–Plesset equation and by a modified Rayleigh–Plesset equation accounting for bubble/bubble interactions in the mean-field theory.      

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。