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.     

Propagation of pressure waves in two-phase flow

We deal with a pressure wave of finite amplitude propagating in a gas and liquid medium or in the fluid in an elastic tube. We study the effects of pipe elasticity on the propagation velocity of the pressure wave. Pressure waves of finite amplitude progressing in the two-phase flow are treated considering the void fraction change due to pressure rise. The propagation velocity of the two-phase shock wave is also investigated, and the behavior of the reflection of the pressure wave at the rigid wall is analyzed and compared to that in a pure gas or liquid. The results are compared to experimental data of a pressure wave propagating in the two-phase flow in a vertical shock tube.     

Solution to the Riemann problem for drift-flux model with modified Chaplygin two-phase flows

In this paper, we concern about the Riemann problem for compressible no-slip drift-flux model which represents a system of quasi-linear partial differential equations derived by averaging the mass and momentum conservation laws with modified Chaplygin two-phase flows. We obtain the exact solution of Riemann problem by elaborately analyzing characteristic fields and discuss the elementary waves namely, shock wave, rarefaction wave and contact discontinuity wave. By employing the equality of pressure and velocity across the middle characteristic field, two nonlinear algebraic equations with two unknowns as gas density ahead and behind the middle wave are formed. The Newton–Raphson method of two variables is applied to find the unknowns with a series of initial data from the literature. Finally, the exact solution for the physical quantities such as gas density, liquid density, velocity, and pressure are illustrated graphically.     

Numerical calculation of the motion of a gas from a surface explosion

The motion of a gas by the normal impact of a high-speed body at the interface between a dense half-space and a vacuum is investigated numerically. The motion of the shock wave and the shape and distribution of the parameters of the gas dispersing in the vacuum are obtained. The motion is studied during the formation of a region with high pressure at the boundary with the vacuum of a gas occupying the half-space z > 0. The assumption of cylindrical symmetry relative to the z axis enables this three-dimensional nonsteady-state problem in the general case to be solved as a two-dimensional problem. For the corresponding one-dimensional problem, the numerical solution and, for certain gases also, the analytic solutions are well known and are considered in detail in [1]. As a result of solving the two-dimensional problem, profiles of the gasdynamic quantities are obtained which are similar to the solutions in the one-dimensional case and the result of the solution by a self-similar method. The cup-shaped surface of the shock wave front with a pressure gradient on it “focusses” the dispersing gas so that its velocity component normal to the surface z = 0 is greater by an order of magnitude than the component parallel to the surface of separation of the medium, and only at individual points is their ratio close to 0.4. Therefore, the dispersing gas is formed into the shape of a “jet”, the pressure and density profiles on the axis of which have a shape similar to the one-dimensional problem of a brief shock, but in the plane z = 0 the pressure and density distributions are similar to the distributions of these quantities in the case of a powerful point explosion in an unbounded medium. The initial disturbance in the symmetrical problem being considered may be the result of either the normal impact of the body with a high velocity at the surface of the dense medium, or the consequence of the effect of a giant laser pulse, or some other process when a certain volume is formed with a high pressure at the interface between the dense medium and a vacuum, or with another low-density medium.     

Ideal gas flow behind a finite-amplitude shock wave

The equations of one-dimensional (with a plane of symmetry) adiabatic motion of an ideal gas are transformed to a form convenient for studying flows between a moving piston and a shock wave of variable intensity. The solution is found for the equations of a motion containing a shock wave which propagates through a quiescent gas with variable initial density and constant pressure. This solution contains four arbitrary constants and, in a particular case, gives an example of adiabatic shockless compression by a piston of a gas initially at rest.     

Detonation Wave Propagation in Rotational Gas Flows

This paper studies the propagation of detonation and shock waves in vortex gas flows, in which the initial pressure, density, and velocity are generally functions of the coordinate — the distance from the symmetry axis. Rotational axisymmetric flow having a transverse velocity component in addition to a nonuniform longitudinal velocity is considered. The possibility of propagation of Chapman–Jouguet detonation waves in rotating flows is analyzed. A necessary conditions for the existence of a Chapman–Jouguet wave is obtained.     

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

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

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

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

Shock wave in a gas-liquid medium

The propagation of weak shock waves and the conditions for their existence in a gas-liquid medium are studied in [1]. The article [2] is devoted to an examination of powerful shock waves in liquids containing gas bubbles. The possibility of the existence in such a medium of a shock wave having an oscillatory pressure profile at the front is demonstrated in [3] based on the general results of nonlinear wave dynamics. It is shown in [4, 5] that a shock wave in a gas-liquid mixture actually has a profile having an oscillating pressure. The drawback of [3–5] is the necessity of postulating the existence of the shock waves. This is connected with the absence of a direct calculation of the dissipative effects in the fundamental equations. The present article is devoted to the theoretical and experimental study of the structure of a shock wave in a gas-liquid medium. It is shown, within the framework of a homogeneous biphasic model, that the structure of the shock wave can be studied on the basis of the Burgers-Korteweg-de Vries equation. The results of piezoelectric measurements of the pressure profile along the shock wave front agree qualitatively with the theoretical representations of the structure of the shock wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 65–69, May–June, 1973.     

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.     

Shock tube study of the dynamical behavior of granular materials

The purpose of the present study is to clarify both the compression phenomenon and the gas filtration effect that take place inside a granular medium when it is dynamically loaded by a shock wave. In order to measure the pore pressure and the total stress at different locations along the granular medium, pressure transducers were placed along the side-wall and at the end-wall of the shock tube test section, which was filled with the granular material. In order to elucidate the gas filtration effect, the results of two experiments with identical granular media but with and without filtration were compared. The gas filtration was eliminated by means of a thin plastic film, which was placed at the front edge of the granular medium. Based on this comparison quantitative information on the gas filtration and its role in the stress formation inside granular media of different material and length was obtained. Furthermore, curves of the dynamic compression and the Young moduli of the granular medium for the range of the operating conditions were reconstructed.     

Structure and bifurcations of plane shock waves in a vibrationally excited gas with an external pumping source

The structure of the density profiles in stationary plane shock waves in a vibrationally excited gas is investigated. For self-similar solutions a bifurcation diagram is plotted in the parametric “traveling wave velocity—degree of nonequilibrium” plane. The bifurcation boundaries of the domains of existence of the structures of different types are analytically derived. It is shown that weak plane shock waves are unstable, accelerate, and break down into a sequence of pulses or-at a fairly high pumping rate-waves with nonzero asymptotics, whose amplitude and propagation velocity are independent of the initial disturbance and are determined by the parameters of the medium itself.     

The influence of magnetic field upon the collapse of a cylindrical shock wave

In the present paper, the problem of propagation of collapsing cylindrical shock wave in an ideal gas permeated by a transverse magnetic field with infinite electrical conductivity is investigated. Here it is assumed that the medium ahead of the shock front is uniform and at rest. Also, its counter pressure concerning the motion of the wave front is neglected. This problem admits a self similar solution of second kind. The similarity exponent has been computed by solving a nonlinear eigenvalue problem and integrating numerically the self-similar equations for various values of adiabatic heat exponent and Cowling number. Numerical computations have been performed to determine the flow field behind the shock wave. The influence of magnetic field strength and adiabatic heat exponent on the flow parameters for various cases is presented.     

Propagation of discontinuities against a static background

The solution of the ideal gasdynamic equations describing propagation of a shock wave initiated, for example, by the motion of a piston against an inhomogeneous static background is considered. The solution is constructed in the form of Taylor series in a special time variable which is equal to zero on the shock wave. In the case of weak shock waves divergence of the series serves as the constraint for such an approach. Then the solution is constructed by linearizing the equations about the solution with a weak discontinuity. In the case of a given background the last solution can be always found exactly by solving successively a set of transport equations, all these equations are reduced to linear ordinary differential equations. The presentation begins from the one-dimensional solutions with plane waves and ends by discussion of spatial problems.     

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.      

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     

Delta Shock Wave Solution of the Riemann Problem for the Non-homogeneous Modified Chaplygin Gasdynamics

The motivation of the present study is to derive the solution of the Riemann problem for modified Chaplygin gas equations in the presence of constant external force. The analysis leads to the fact that in some special circumstances delta shock appears in the solution of the Riemann problem. Also, the Rankine–Hugoniot relations for delta shock wave which are utilized to determine the strength, position and propagation speed of the delta shocks have been derived. Delta shock wave solution to the Riemann problem for the modified Chaplygin gas equation is obtained. It is found that the external force term, appearing in the governing equations, influences the Riemann solution for the modified Chaplygin gas equation.

     

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

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

驱动气体的密度梯度对弹丸发射速度的影响

为进一步提升轻气炮的发射能力,提出采用梯度气体替代单一氢气或氦气作为驱动气体的方法,通过对等直径发射器进行分析,建立了弹丸在梯度气体驱动下的加速运动模型,对比了氖-氦梯度气体驱动与单一氦气驱动的发射能力差异,分析了梯度气体参数对发射性能的影响。结果表明,与单一氦气驱动相比,氖-氦梯度气体驱动能够提升0.4～1.4 km/s的发射速度或降低0.2～0.9 GPa的发射过载;气体的密度和活塞的运动速度对发射速度和过载的影响最大,气体压力和多方气体指数的影响次之;梯度气体中,高密度气体应选择多方气体指数和密度较高的气体（如氖气、氩气等）;梯度气体界面位置（高密度气体占比）对发射速度的影响不大,但高密度气体占比少有利于降低弹底压力。