Hydroshock theory in a two-phase gas-liquid mixture in a slug flow regime

A study is made of the intensity of a hydroshock in a two-phase gas-liquid mixture in a slug flow regime in the case when a pipeline is shut off by a liquid slug. The intensity is studied as a function of the length of the shut-off section of the liquid slug, the content of gas bubbles in the liquid slugs, and the pipeline shut-off law, and with allowance for the shock-wave character of the process [1, 2]. The calculated data using the shock-wave theory agree well with the experimental data of [3] and, unlike the results of the linear theory of [3], make it possible to determine the intensity of the hydroshock not only in the case of weak waves, but also in the case of waves of moderate intensity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 188–190, September–October, 1985.     

Intensification of a detonation wave in the zone of secondary reactions in a two-phase medium

A method is proposed for the numerical calculation of one-dimensional nonsteady-state flows of a mixture of a gas with particles, based on the separation of a system of differential equations for a two-phase medium into two subsystems. The problem is solved concerning the propagation of a plane detonation wave in a mixture of a detonating gas with particles, behind the front of which secondary chemical reactions are taking place between the vapors of the particle material and the detonation products. The velocity profiles of the gas and of the thermodynamic functions behind the detonation wave front are determined, and also the dependence of the detonation velocity on the distance to the point of initiation. The conditions for intensification of the detonation wave are obtained in the zone of secondary reactions.Leningrad. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 92–96, September–October, 1972.     

Thermal mechanism of wave damping in a gas-liquid mixture

The effect of interphase heat transfer on the process of propagation of a wave in a monodisperse gas-liquid mixture is investigated. A three-wave equation of the Boussinesq type is derived and formulas for the dependence of the thermal components of the dissipative coefficients on the thermophysical parameters of the mixture are obtained. The limits of applicability of the short-wave method are determined. The theoretical and experimental values of the phase wave velocity are compared and found to be in good agreement.Yerevan. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 75–83, November–December, 1994.     

Relaxation of overdriven detonation waves with finite reaction rate

A numerical study is made of the interaction of a detonation wave having finite reaction velocity with a rarefaction wave of different intensity which approaches it from the rear, for the Zeldovich-Neumann-Doring (ZND) model with a single irreversible reaction A B. It is found that, for a fixed value of the parameter characterizing the initial supercompression (depending on the activation energy and the heating value of the mixture), the considered interaction leads either to a gradual relaxation of the detonation wave and its transition to the Chapman-Jouguet (CJ) regime, or to the development of undamped oscillations.Interest in the problems of detonation and supersonic combustion has increased in recent years. This is associated with the appearance and development of new experimental and theoretical techniques; it is also associated with the further development of air-breathing reaction engines, and other practical requirements. The present state of detonation theory is reflected in the survey [1].It has been established [2] that the detonation wave in gases nearly always has a complex nonuniform structure. Transverse disturbances are observed under a wide range of conditions and differ both in amplitude and wavelength. At the same time, behind the detonation leading front there is a region of uncompletely burned gas corresponding to the effective ignition induction period [3]. In spinning detonation the induction period is significantly longer than the heat release period and transverse detonation waves traveling in the induction zone of the head wave appear [3, 4]. Such a secondary detonation wave is free of transverse disturbances. The same is true of the detonation waves observed in the wake behind a body moving at high speed in a combustible medium [5] or in a gas which has been preheated by a shock wave [6].Although it is possible, under favorable conditions, to study in detail the system of discontinuities accompanying detonation, information on the extensive zones in which heat release takes place is scarce, the mechanism of detonation wave autonomy (in particular, the role of the rarefaction zone behind the wave) is not entirely clear, and the fact that, in spite of the complex structure, an autonomous detonation propagates with the CJ velocity calculated on the basis of one-dimensional theory has not yet been explained.In studying the nonlinear phenomena associated with the finite reaction rate it is quite acceptable to investigate only the simple one-dimensional detonation model, with which it is convenient to restrict ourselves to a single effective chemical reaction. This model is particularly reasonable since, in certain cases, the real detonation is virtually one-dimensional.The question of the stability of the one-dimensional detonation wave to disturbances of its structure has been examined by several authors [7–13]. The use of computers makes possible the direct computation of flows with heat release and the study of their properties. This method has been used in [11–13] to study the stability problem for a detonation wave with respect to finite disturbances.In the present paper we present a numerical study of the interaction of a detonation wave having finite chemical reaction rate with a rarefaction wave of different intensity approaching it from the rear for the ZND model with a single irreversible reaction A B. It is found that for a fixed value of the parameter characterizing the difference between detonation and the CJ waves, depending on the activation energy E and the mixture heating value Q_m, the interaction in question leads either to a gradual relaxation of the detonation wave and its transition to the CJ regime (this relaxation may be accompanied by decaying oscillations) or to the appearance of undamped oscillations (the unstable regime). The parameters E and Q_m affect the wave stability differently: with increase of Q_m, the wave is stabilized; with increase of E, it is destabilized. The boundary between the stable and unstable detonation wave propagation regimes is found. This boundary has a weak dependence on the rarefaction wave intensity. Estimates and calculated examples show that the amplitude of the unstable wave oscillations is finite and that the average detonation propagation velocity is close to the CJ velocity computed for the given heating value Qm.The author wishes to thank G. G. Chernyi for his guidance and L. A. Chudov for advice on computational questions.     

Wave front in the blunt body immersion problem

The immersion of a three-dimensional blunt convex body in a compressible fluid with nonpositive acceleration is considered in the linear formulation. It is shown that at every instant the perturbation zone will be convex. The fluid particle velocity and pressure are calculated at each point on the wave front. At every instant the wave front is wholly determined by the initial supersonic stage of propagation of the boundary of the body-fluid interaction zone.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.4, pp. 5–11, July–August, 1992.The author wishes to thank A. G. Khovanskii for his constant interest.     

Theory of regular and mach reflection of shock waves in a two-phase gas-liquid medium

An equilibrium model of a two-phase gas-liquid medium, with allowance for the proportion, density, and compressibility of the components, and with a difference from [1] in that the adiabatic velocity of sound is introduced, has been used in order to study the regular and Mach (elementary theory) reflection of a shock wave of moderate intensity from a solid wall throughout the whole range of gas proportions. A complicated nonmonotonic variation has been found for the pressure on the wall behind the reflected wave, the angle of reflection, and the angle of departure of the triple point as functions of the gas proportion, the angle of incidence, and the intensity of the incident wave. In particular, it is shown that oblique reflection for moderate and low gas contents leads to the formation of a stronger reflected shock wave than does normal reflection. The effect of the gas proportion on the position of the boundary between the regions of regular and Mach reflection has already been studied in [2]. The results are described of serial calculations of the parameters of reflection for an air-water mixture, and these results agree fairly well for normal reflection with the known experimental data [3] for low and moderate gas contents. In the limiting case, the results agree with the known results for single-phase media [4, 5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 188–190, March–April, 1985.     

Nonlinear wave motions on liquid surfaces

A study is made of the formation of a shock wave (bore), produced by the movement of an initially weak discontinuity in the spatial derivatives of velocity and liquid depth in an area of stationary current in a channel of constant inclination. The formation of shock waves from compression waves was first studied by Riman [1]. Frictional resistance was considered in the Chezy form. The equations obtained therein for determination of the moment in time and spatial coordinates of the point at which the shock wave is formed, as well as the laws for propagation of shock waves are applicable to the problem of one-dimensional transient motion in a gas, the pressure of which is dependent on density. Instantaneous collapse of waves, as well as formation and movement of bores in rivers for an idealized flow model in a channel with horizontal bottom, neglecting friction, were described by Khristianovich, Mikhlin, and Devison [2], and Stoker [3]. Recently in the work of Sachdev and Bhatnagar [4], using numerical integration of the equation for bore intensity, the problem of shock wave propagation in a channel of constant inclination with consideration of fluid resistance in the Chezy form was studied. Gradual wave collapse and the bore formation mechanism were studied by Stoker [3] on the basis of the shallow-water theory. Neglecting friction on the horizontal channel bottom, he calculated the moment of time and coordinates of the point at which the shock wave is formed in the case where the initial disturbance is sinusoidal. The dependence of these values on wave amplitude for a channel of constant inclination was obtained by Jeffrey [5], who also neglected friction on the channel bottom and considered the initial disturbance to be sinusoidal. Lighthill and Whitham [6] discovered that for Froude numbers greater than two, the linear theory led to unlimited growth in the intensity of the flood wave. We note that the studies of flood-wave motion in the region of the first characteristic, performed in [3, 6], differ only in the forms of the resistance laws and dependences of the unknown functions on the variables. Physical peculiarities of various liquid wave motions were also examined by Lighthill in [7].Saratov. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 62–66, March–April, 1972.     

Propagation of shock waves in free space

The problem of the exit of a shock wave from an axisymmetric channel and its propagation in a free space occupied by an ideal gas is examined. This problem has been studied earlier in [1], in which the shock wave front was considered planar, as well as in [2], in which the wave front was regarded as a surface of an ellipsoid of revolution. The solutions obtained in these studies assumed the presence of two regions in the wave-front surface: the region of the original shock wave and a region stemming from the decomposition of an infinitesimally thin annular discontinuity of the gas parameters, with the wave intensity over the front surface in each region being considered constant, i.e., the wave character of the process over the front was not considered. In this study a solution will be achieved by the method of characteristics [3–5] of the equations of motion of the shock-wave front, as obtained in [6, 7]. Flow fields are determined for the region immediately adjacent to the shock-wave front for a wide range of shock-wave Mach numbers M _a =1.6–20.0 for = 1.4. On the basis of the data obtained, by introduction of variables connected with the length of the undisturbed zone, as calculated from the channel cross-section along the x axis, together with the pressure transition at the wave front, relationships are proposed which approach self-similarity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 163–166, September–October, 1971.In conclusion, the author thanks S. S. Semenov for his valuable advice on this study.     

Flame propagation in the turbulent flow of a burning mixture

The propagation of a flame in the turbulent flow of a burning mixture is analyzed theoretically. An equation is derived for the gas temperature and velocity probability distributions. The solutions of this equation are analyzed and the rate of flame propagation is determined.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 3–15, September–October, 1976.The author wishes to thank V. A. Sabel'nikova for critical comments.     

Theory of flame propagation in a gas and drop mixture

The equations of a reacting multiphase continuous medium [1] are used to investigate the problem of steady-state flame front propagation in a gas mixture with evaporating drops. A simple model for ignition of the liquid drops is proposed which is based on the application of the method of equally accessible surfaces [2] to the heat and mass exchange processes between the microflames surrounding the separate drops, the drops, and the carrying gas medium. The parameter distributions in the macroscopic flame front as well as the dependences of the flame propagation velocity in the gas suspension on a number of parameters governing the process under investigation are represented.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 101–108, July–August, 1973.     

Calculation of a turbulent two-phase isobaric jet

A numerical calculation is carried out by the finite-difference method based on proposed equations for a turbulent submerged jet containing an admixture of solid particles. The relative longitudinal particle velocity and the influence of particles on the turbulence intensity are taken into account. The calculated results adequately agree with available experimental data. A turbulent two-phase jet is examined in [1] on the basis of the theory for a variable density jet, assuming equal mean velocities for the gas and particles and not considering the influence of particles on the turbulence intensity. Particles are analogously taken into account by a noninertial gas mixture in [2, 3], and a particle Schmidt number of 1.1 is assumed in [4]. A model is proposed in [5] which takes into account the influence of particles on the turbulence intensity of the gas phase. Problems concerning the initial and main sections of a submerged jet were solved in [6] by the integral method on the basis of this model and the assumed equality of the mean velocities of the gas and particles. Turbulent mixing of homogeneous two-phase flows with allowance made for dynamic nonequilibrium of the phases is considered in [7]. However, the neglect of turbulent transfer of particle mass and momentum led to a physically unrealistic solution for the particle concentration in the far field of the mixture. A two-phase jet is considered in the present work on the basis of the theory of a two-velocity continuous medium [8, 9] with allowance made for turbulent transfer of particle mass and momentum. The influence of particles on the turbulence intensity of the gas phase is taken into account with the model of [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 57–63, September–October, 1976.The author acknowledges useful comments and discussion.of the work by G. N. Abramovich and participants of his seminar. The author sincerely thanks I. N. Murzinov for scientific supervision of the work.     

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.     

Amplification of shock waves in a bubble-filled liquid with gas concentration gradient

The properties are studied of the propagation of unsteady shock waves in a gas-liquid system of bubble structure in the case when the volume concentration of the gas changes in the direction of motion of the shock wave. It is established that when there is a sufficiently rapid drop in the gas content, an effect of amplification of the shock wave is observed which is due to the deceleration of the medium behind the shock wave. A study is made of the laws of the evolution of long- and short-wave pulsed perturbations in such systems. The authors consider processes of reflection of waves from obstacles and their passage from a gas into a bubble liquid, from a two-phase mixture into a pure liquid. The contribution is determined of nonequilibrium effects to the process of amplification of a wave.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 49–54, January–February, 1988.The authors wish to express gratitude to R. I. Nigmatulin for his interest in the study and for useful discussions.     

Water hammer and propagation of perturbations in elastic fluid-filled pipes

A Korteweg—de Vries (KV) approximation is constructed in this paper for the perturbations being propagated in elastic pipes filled with fluid. On the basis of the approximation constructed and the equation obtained for the perturbations of a finite-amplitude velocity, the water-hammer phenomenon is analyzed in the Zhukovskii formulation, and the water hammer in systems with preliminary longitudinal tension is considered separately. Special attention in the study of the perturbations is paid to the signal structure and evolution in the hydraulic line. Taking account of the inertial properties of the pipe in the approximation mentioned permitted the indication of new effects, in principle, which are essential for applied problems of the propagation of perturbations in elastic hydraulic lines. In particular, it is shown that the initial signal can be doubled in such lines by redistributing its intensity over the frequencies. It is established that the origination of an oscillating forerunner is possible in hydraulic lines with preliminary tension. Starting with [1], the water-hammer phenomenon was investigated in many papers, in [2], for example. The main attention in these papers was paid to the propagation velocity of the water hammer and its intensity. After simplification, the initial system of Zhukovskii equations contains no mechanism hindering the twisting of the wave profile, and, therefore, there is no possibility of stationary shock formation within the framework of this theory. Moreover, the Zhukovskii theory of the water hammer and of propagation of perturbations in elastic pipes results in the conclusion that the wave structure, velocity, and amplitude depend essentially on the characteristics of the initial perturbation and can differ significantly from the water hammer predicted by theories for powerful signals in sufficiently long pipes.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–8, July–August, 1976.     

Asymptotic solution of the Orr-Sommerfeld equation in the outer region of the boundary layer with allowance for nonparallelism

There have been many publications devoted to the investigation of the hydrodynamic stability of nonparallel flows on the basis of the modified Orr-Sommerfeld equation [1–4]. Taking into account the additional terms associated with the presence in the flow of a transverse component of velocity and acceleration can lead not only to a significant quantitative discrepancy as compared with calculations based on the usual Orr-Sommer-feld equation but also to qualitatively new results (nonclosure of the neutral curves for flow on a permeable surface in the presence of strong injection [4]). In this paper an asymptotic solution of the Orr-Sommer-feld equation, valid in the outer region of boundary layer flow, is constructed for self-similar gradient flow over a surface (Falkner-Skan flow). The continuity of the eigenvalue spectrum for an unbounded increase in the perturbation propagation velocity is demonstrated on the basis of the solution obtained. For the ordinary Orr-Sommerfeld equation a continuous transition of the spectrum through the value of the perturbation propagation velocity C_r=1 (which coincides with the velocity of the external flow) is impossible [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 171–173, January–February, 1987.     

Effect of particle size on nonlinear wave propagation in a fluidized bed

The effect of particle size (Archimedes number) on the propagation of a kinematic particle concentration wave in a fluidized bed is investigated. The dependence of the characteristic wave velocity on the porosity of the bed (particle concentration) and the Archimedes number (or the Reynolds number for flow past individual particles of the dispersed phase) is determined. The evolution of a nonlinear perturbation of the bed porosity is investigated and the formation of discontinuities in the concentration of the dispersed phase is studied in relation to the particle size (Archimedes number). It is shown, in particular, that, as distinct from a bed of small particles, in a bed of large particles with quadratic interphase interaction only compression discontinuities can be formed. The results obtained can be used to analyze the formation of inhomogeneities (slugs and bubbles) in a fluidized bed in relation to the particle size.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 96–100, May–June, 1987.     

Population inversion behind a detonation wave propagating in a variable-density medium

The creation of an active medium by means of detonation has been investigated on a number of occasions. It was suggested that one could use the expansion of the detonation products of an acetylene-air mixture in vacuum [1] or the cooling of the detonation products of a mixture of hydrocarbons and air through a nozzle [2, 3]. In [4], the detonation of a solid high explosive was used to produce population inversion in the gas mixture CO₂-N₂-He(H₂O). Stimulated emission from HF molecules was observed in [5] behind the front of an overdriven detonation wave propagating in an F₂-H₂-Ar mixture in a shock tube. Population inversion behind a detonation wave was studied in H₂-F₂-He mixtures in [6–8] and in H₂-Cl₂-He mixtures in [9] with energy release on a plane and on a straight line in a medium with constant density. Similar problems were solved for shock waves propagating in both a homogeneous gaseous medium [7, 10] and in the supersonic part of an expanding nozzle. In the present paper, we study theoretically population inversion behind an overdriven detonation wave propagating in a mixture (fine carbon particles + acetylene + air) which flows through a hypersonic nozzle. The propagation of detonation in media with variable density and initial velocity was considered, for example, in [11, 12]. Analysis of the gas parameters behind a detonation wave propagating in a medium with constant density (for a given fuel) showed that the temperature difference across the detonation front is insufficient to produce population inversion of the vibrational levels of the CO₂ molecule.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 65–71, January–February, 1980.I am grateful to V. P. Korobeinikov for a helpful discussion of the results.     

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.     

A contribution to the theory of streamer breakdown

The development of a systematic theory of streamer breakdown of a gas requires the consideration of the transport of the region of ionization toward the ionized gas in an electric field depending on the form of the streamer, which in turn is determined by the transport mechanisms [1–3]. In this form the problem is very complicated,and the theory takes the path of investigation of different qualitative models of a streamer [4]. It is assumed in [4] that the rates of anode-directed and cathode-directed streamers are determined by the drift velocity of the electrons. The mechanism of propagation of anode-directed streamers is taken to be the development of avalanche from the leading front of the electrons traveling to the anode. On the side of the cathode, electrons before the front of the cathodedirected streamer are produced due to the transport of radiation from the ionized region [1]. It is shown in [5] that direct photo-ionization is ineffective because of the small range of the quantas, and a mechanism of development of cathode-directed streamer related to the associative ionization of excited atoms is proposed. These atoms are formed by long-span resonance photons from the wings of the spectral line. An interesting prediction of the theory [4] was a linear dependence of the velocity of the streamers on their length. This dependence was confirmed in experiments on the study of streamer breakdown initiated at the center of the discharge gap in spark chambers [6, 7]. At the same time, for streamers developing from avalanche initiated at one of the electrodes the velocity of propagation of the breakdown wave remains constant with a good accuracy in gaps having lengths of the order of 1 m. In the present work a qualitative theory is developed which permits one to calculate the velocity of the an ode-directed streamer in the case where it is independent of the length. Since for pressures of the order of atmospheric pressure the diffusion coefficient of excited atoms [8] is comparable with the electron diffusion coefficient, the effect of radiation transport is disregarded. The stability of the front of the streamer to infinitely small perturbations is investigated. It is shown that, when the finite thickness of the front is taken into consideration, the streamer is stable. It is unstable in the approximation of infinitely thin leading fronts.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 56–65, January–February, 1973.The authors thank A. A. Vedenov, E. P. Velikhov, A. P. Napartovich, and O. B. Firsov for valuable discussions.     

Korteweg-de vries-type equations in the theory of surface waves in electrically conducting liquids

Equations are obtained which describe the propagation of long waves of small, but finite amplitude in an ideal weakly conducting liquid and on the basis of these equations the influence of MHD interaction effects on the characteristics of the solitary waves is investigated. The wave equations are derived under less rigorous constraints on the external magnetic field and the MHD interaction parameter than in [1–3]. It is shown that the evolution of the free surface is described by the KdV-Burgers or KdV equations with a dissipative perturbation, and that the propagation velocity of the solitary waves depends on the strength of the external magnetic field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 177–180, November–December, 1989.