Calculation of the flow of a two-phase mixture in a Laval nozzle,allowing for the turbulent diffusion of the particles

The flow of a mixture of gas and condensed particles in an axisymmetric Laval nozzle is considered. The motion of the particles is calculated in a specified field of gas flow, with due allowance for their turbulent diffusion. The results of calculations indicating the necessity of allowing for this phenomenon when considering the motion of particles toward the wall of a profiled nozzle are presented.Translated from Izvestiya Akademii Nauk SSSR, No. 2, pp. 161–165, March–April, 1973.     

Flow in the peripheral zone of a freely expanding axisymmetric ideal-gas jet

When an ideal gas flows from a nozzle into a vacuum a substantial part of the jet is occupied by the peripheral zone in which the angle of inclination of the velocity vector W to the axis of symmetry is close to or exceeds /2. The known solutions [1–4] for the far field are unsuitable for describing the jet, since they are valid only at relatively small values of . In this study the author obtains an analytic solution describing in explicit form the shape of the streamlines and the distribution of the parameters in the peripheral zone of a jet flowing into a vacuum from a nozzle with an arbitrary parameter distribution in the exit section. At the nozzle edge the solution describes Prandtl-Mayer flow. As the radial coordinate tends to infinity, the streamlines tend to asymptotes whose angle of inclination depends on the distribution of the parameters in the nozzle exit section, and the local Mach numbers increase without bound.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 145–153, May–June, 1987.     

Investigation of isothermal flow of a mixture of gas and particles in a channel of variable section

Isothermal flow of a gas with particles is investigated analytically, which makes it possible to analyze all possible flow regimes in channels of different shapes. It is shown that in a channel of constant section there are two possibilities: either an equilibrium regime is established with constant flow parameters, or the gas reaches the velocity of sound, and then further flow in the channel is impossible (blocking of the channel). In a contracting nozzle, blocking also occurs if the channel is sufficiently long. In an expanding nozzle when there are particles in the gas with a velocity lower than the gas velocity, it is possible to have flow regimes with transition through the velocity of sound: a subsonic flow goes over into a supersonic flow and, conversely, it is also possible to have a flow in which there is blocking of the channel, which is quite different from the flow of a pure gas in an expanding nozzle and is due to the influence of interphase friction on the flow. The variation of the pressure along the flow can be nonmonotonic with points of local maximum or minimum which do not coincide with the singular point at which the gas velocity reaches the velocity of sound. In the case of nonequilibrium gas flows with particles in a Laval nozzle, the velocity of the gas may become equal to the isothermal velocity of sound not only in the exit section of the nozzle or in its expanding part, as noted in [4–6], but also at the minimal section, since it is possible to have flows for which the velocities of the phases are equalized at this section.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 61–68, October–December, 1981.     

Numerical investigation of nonequilibrium flows of mixtures of fusible metal particles and gases in a laval nozzle

There are numerous papers [1–11] on the determination of the parameters of condensed oxide particles which are formed during combustion of metallized fuels. The ambiguity, and sometimes the contradictoriness, of the test results obtained [3–5, 9–11] indicate the difficulties in conducting correct experimental investigations. In this connection, numerical studies using mixtures of calibrated liquid-metal particles and different gases are of practical interest. Different probes can be calibrated by using calibrated two-phase flows, the two-phase flow around models and probes can be studied, as can the interaction between liquid-metal particles and the front of an aerodynamic compression shock, their intrusion in different entraining media, the interaction between fine particles (particle-projectiles) and large size particles (particle-targets), etc. In many cases, the prehistory of the flow and the parameters of the gas mixture with the particles in the area of the nozzle exit section must be known to investigate the above-mentioned phenomena. The parameters of different nonequilibrium flows of mixtures of gallium particles and gases in a Laval nozzle are investigated numerically in this paper; the maximum diameter (upper boundary of the spectrum) of the particles (d_s = 30 ) which are not destroyed in the nozzle under the effect of the aerodynamic forces and are suitable for use in a calibrated two-phase stream is determined. The computations were carried out in a one-dimensional approximation according to [12–14].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 86–91, March–April, 1976.The authors are grateful to V. K. Starkov and U. G. Pirumov for discussing the results of the research and to N. M. Alekseev for aid in constructing the graphs.     

Investigation of a Gas Flow with Solid Particles in a Supersonic Nozzle

A two-phase flow with high Reynolds numbers in the subsonic, transonic, and supersonic parts of the nozzle is considered within the framework of the Prandtl model, i.e., the flow is divided into an inviscid core and a thin boundary layer. Mutual influence of the gas and solid particles is taken into account. The Euler equations are solved for the gas in the flow core, and the boundary-layer equations are used in the near-wall region. The particle motion in the inviscid region is described by the Lagrangian approach, and trajectories and temperatures of particle packets are tracked. The behavior of particles in the boundary layer is described by the Euler equations for volume-averaged parameters of particles. The computed particle-velocity distributions are compared with experiments in a plane nozzle. It is noted that particles inserted in the subsonic part of the nozzle are focused at the nozzle centerline, which leads to substantial flow deceleration in the supersonic part of the nozzle. The effect of various boundary conditions for the flow of particles in the inviscid region is considered. For an axisymmetric nozzle, the influence of the contour of the subsonic part of the nozzle, the loading ratio, and the particle diameter on the particle-flow parameters in the inviscid region and in the boundary layer is studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 65–77, November–December, 2005.     

Extremal nozzle contours for gas flows with particle lag

The determination of the extremal nozzle contour for gas flow without foreign particles has been carried out in several studies [1–6], based on the calculation of the flow field using the method of characteristics.In [7, 8] the equations are derived for the characteristics and the relations along the streamlines which are required for calculating two-dimensional gas flow with foreign particles. The variational problem for two-phase flow in the two-dimensional formulation may be solved by the method of Guderley and Armitage [9] with the use of equations given in [7] or [8]; however this method is very tedious, even with the use of high-speed computers.In [10, 11] studies are made of two-phase one-dimensional flows by expanding the unknown functions in series in a small parameter, defined by the particle dimensions. In [12] a solution is given for the variational problem (in the one-dimensional formulation) of designing the contour of a nozzle with maximal impulse. However that study does not take account of the static term appearing in the impulse and the solution is obtained in relative cumbersome form. Moreover, the question of account for the losses due to nonparallelism and nonuniformity of the discharge was not considered.The present paper considers in the one-dimensional formulation the flow of a two-phase medium in a Laval nozzle with small particle lags (in velocity and temperature). The variational problem of determining the maximal nozzle impulse is formulated along the nozzle contour for fixed geometric expansion ratio. The impulse losses due to nonparallelism of the discharge are simulated by a function which depends on the ordinates which are variable along the contour and on the slope of the tangent to the contour.The author wishes to thank Yu. D. Shmyglevskii and A. N. Kraiko for helpful discussions and V. K. Starkov for carrying out the calculations on the computer.     

Some peculiarities of two-phase flow in nozzles

Calculations are conducted for unidimensional two-phase flow in nozzles for a wide range of particle concentrations and dimensions. It is established that there exists a maximum in loss of specific momentum due to a lag in particle velocity and temperature relative to the gas. The results obtained are compared with calculations using linearized theory as well as with experimental data. The agreement between calculation and experiment is noted. Equilibrium flow of a two-phase mixture with solidification of liquid particles is considered. The presence of an anomalous flow region is established, where in the model of an ideal unidimensional equilibrium flow in a nozzle with discharge into a vacuum the presence of two successively located minimum sections is necessary.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–57, May–June, 1973.     

Theory of a laval nozzle for a two-phase mixture containing particles of small lag

A two-velocity and two-temperature model is considered for a continuous medium in relation to the flow of a mixture of gas and particles in the subsonic, transsonic, and supersonic parts of a Laval nozzle. It is assumed that the particles are small, and hence that the coefficients f and ^q, which define the interaction with the gas, are large (these coefficients are inversely proportional to the square of the particle radius for a Stokes mode of flow). This means that the velocity or thermal lag of the particles relative to the gas is small. The solution is sought as expansions with respect to the small parameters ₁=1/f and ₂=1/^q.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 89–100, July–August, 1973.     

Experimental investigation of a turbulent jet carrying heavy particles of a disperse phase

The results are given of an experimental investigation of an immersed jet carrying heavy spherical particles in the case of uniform velocity fields of the particles and the gas and uniform concentration of the particles at the nozzle exit. The laser-optical method was used to measure the velocity fields of the gas and the particles and the concentration of the latter in the flow-rate concentration range of the additive from 0 to 1.5 and particle diameters from 35 to 67 m. It was found that the flow in the jet is essentially nonequilibrium, the disequilibrium and intensity of the mixing in the jet depending on the concentration of the additive and the diameter of the particles. The investigation of the initial section revealed a somewhat anomalous behavior of the jet at a comparatively small (₀=0.5) flow-rate concentration of the additive, and this indicates a need to take into account the prehistory of the flow in the construction of a method for calculating the initial section of a two-phase jet.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 26–31, September–October, 1981.     

Dependence of the stagnation temperature of a jet of gas containing solid particles on the concentration of particles and their velocity

This article gives the results of experiments on the measurement of the stagnation temperature of a two-phase jet, issuing from a nozzle. The experiments were made using a mixture of air and aluminum oxide (particle diameter 50) with a ratio of the mass flow rate of the solid phase to the mass flow rate of the gas equal to 0.3–2.5, and at initial temperatures of the mixture of 150–450°C. It follows from the results of the experiments that the stagnation temperature of a two-phase flow considerably exceeds the temperature of the mixture at the inlet of the nozzle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 175–176, July–August, 1970.     

Influence of the mixing of two flows with different total parameters in a laval nozzle on its integral characteristics

A calculation is made of the turbulent zone of mixing of two flows of viscous and heat conducting gas in a Laval nozzle. For such a nozzle of given geometry, a comparison is made of calculations of the integrated characteristics of flows that are nonuniform with respect to the total parameters in the framework of various models: laminar hydraulics, viscous laminar hydraulics, and total mixing without hydraulic losses. The calculations are made for a stationary, nonswlrling flow of a viscous heat conducting gas with nearly discontinuous step distribution of the total parameters at the entrance to an axisymmetric Laval nozzle of given geometry. In this situation, the gas flows with different total parameters at the entrance to the nozzle are separated by a surface near which the profiles of the flow parameters are specified on the basis of boundary-layer theory. In the blocked regime investigated here, the flow in the part where the nozzle becomes narrower and at least at the beginning of the expanding part does not depend on the pressure of the surrounding medium. The integrated characteristics of the nozzle (gas flow rate G, impulse I, specific impulse i = I/G, etc.) depend on the parameter distributions at the entrance to the nozzle, and also on the turbulent mixing of the flows in the mixing zone. To analyze the dependence of the integrated characteristics on the turbulent mixing, the values of these characteristics calculated in the framework of the three models are compared. The model of mixing without hydraulic losses presupposes complete equalization of the parameters of the original inhomogeneous flow in the constant-area chamber in front of the nozzle with conservation of the mass, energy, and momentum fluxes. The model of laminar hydraulics is described in detail in [1, 2]. The model of viscous laminar hydraulics will be described in Sec. 1.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 114–119, July–August, 1979.I thank A. N. Kraiko for supervising the work, A. N. Sekundov for helpful discussions, and I. P. Smirnova and A. B. Lebedev for making available the computer program.     

The structure of particle-laden,underexpanded free jets

Underexpanded, supersonic gas-particle jets were experimentally studied using the shadowgraph technique in order to examine the influence of the dispersed particles on the shape of the free jet and the structure of the imbedded shock waves. The particle mass loading at the nozzle exit was varied between zero and one, and two sizes of particles (i.e. spherical glass beads) with mean number diameters of 26 and 45 m were used. It was found that the Mach-disc moves upstream towards the orifice with increasing particle loading. The laser light sheet technique was also used to visualize the particle concentration distribution within the particle jet and the spreading rate of the particle jet. Furthermore, the particle velocity along the jet centerline was measured with a modified laser-Doppler anemometer. These measurements revealed that the particles move considerably slower than the gas flow at the nozzle exit. This is mainly the result of the particle inertia, whereby the particles are not accelerated to sonic speed in the converging part of the nozzle.In order to further explore the particle behavior in the free jet, numerical studies were performed by a combined Eulerian/Lagrangian approach for the gas and particle phases, including full coupling between the two phases. The numerical results showed that the application of different particle velocities at the nozzle exit as the inlet conditions, which were below the sonic speed of the gas phase has a significant influence on the free jet shape and the configuration of the shock waves. These results demonstrate that the assumption of equilibrium flow (i.e. zero slip between the phases) at the nozzle exit which has been applied in most of the previous numerical studies is not justified in most cases. Furthermore, the numerical calculations of the free jet shape and the particle velocity along the jet axis were compared with the measurements. Although correlations for rarefaction and compressibility effects in the drag coefficient were taken into account, the particle velocity along the center line was considerably overpredicted.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Two-dimensional transonic vortex flows of an ideal gas

Equations are obtained for two-dimensional transonic adiabatic (nonisoenergetic and nonisoentropic) vortex flows of an ideal gas, using the natural coordinates (=const is the family of streamlines, and =const is the family of lines orthogonal to them). It is not required that the transonic gas flow be close to a uniform sonic flow (the derivation is given without estimates). Solutions are found for equations describing vortex flows inside a Laval nozzle and near the sonic boundary of a free stream.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 105–109, September–October, 1973.     

Hypersonic flow past a supersonic two-phase source

The interaction of a uniform hypersonic gas flow with a supersonic two-phase gas-particle source is considered. In the symmetry-axis neighborhood between the bow and termination shock waves, an approximate analytical solution for the carrier-phase parameters is found. On the basis of parametric numerical calculations, the behavior of the particle trajectories and the concentration distribution in the shock layers are studied for both continuum and free-molecule flow regimes around the particles. The appearance of regions with multiple intersections of the particle trajectories and the formation of "layer structures" in the particle concentration distributions (particle accumulation regions near the envelopes of the particle trajectories) are indicated. The dependence of the number of the high concentration layers on the governing parameters is studied. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 134–147, May–June, 1998. The work received financial support from the Russian Foundation for Basic Research (project No. 96-01-00313) and the National Foundation for Natural Sciences of China (joint RFBR-NFNS grant No.96-01-00017c).     

Influence of the viscous boundary layer on the deposition of particles from a gas suspension flowing past a sphere

The flow of incompressible gas containing particles past bodies of simple shapes at moderate and high velocities is investigated in [1–5], in which the flow of the carrier medium is assumed to be irrotational. The estimates made in [3] for the neighborhood of the stagnation point show that it is necessary to take into account the viscous boundary layer in the case of fine particles. In the present paper, the viscous flow of a gas suspension over the front surface of a sphere at Reynolds numbers R = 10³–10⁷ is considered. It is assumed that the carrier gas is incompressible and the particle concentra ion negligibly small. The influence of the boundary layer on the particle trajectories and the deposition of the disperse phase on the surface of the sphere is investigated. It is shown that there is a wide range of flow parameters for the gas suspension in which the influence of the boundary layer is important. The limits of this range are established.Translated from Izvestiya Akademli Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 59–66, January–February, 1982.I thank Yu. P. Savel'ev for a helpful discussion of the work.     

Interaction of the vortex structures formed in parallel pulsed jets

The experimental investigation of the lateral interaction of the heads of pulsed jets and primary shock waves at various nozzle spacings and pressure ratio numbers is described. The various stages of formation of a composite pulsed jet issuing from a multinozzle block are classified and the flow development mechanisms are explored. It is shown for both a block and a single nozzle the shock wave travels with almost the same velocity, whereas the jet front formed at the exit from a single nozzle moves much more slowly than the jet front formed beyond a nozzle block. Long-lived lateral bursts of gas, whose dimensions are an order greater than those of the jet bursts, are detected. Their long period of existence considerably increases the stabilization time of the steady-state structure and parameters as compared with a single pulsed jet with the same flow rate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 153–159, November–December, 1987.     

Investigation of flow beyond a supersonic nozzle by an electron-beam technique

The results of an experimental study of a flow of rarefied gas of density 10^–5 g/cm³ beyond the cutoff of a hypersonic nozzle (M11) by means of an electron beam with energy up to 43 keV are presented. The density and velocity fields at different distances from the nozzle and various receiver pressures were measured using this method and the static and total pressure fields were also measured with the help of a Pitot tube. The flow parameters beyond the nozzle were calculated for two limiting cases: with equilibrium condensation and without condensation. This calculation is compared with the experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–117, January–February, 1976.The authors thank S. N. Romanenko for help with the electron-beam experiments.     

“Lateral” interaction of a supersonic underexpanded ideal-gas jet with surfaces of different shape

A numerical investigation is made of the interaction of an underexpanded jet of an inviscid and nonheat-conducting gas issuing from an axisymmetric conical nozzle with plane, cylindrical, and spherical surfaces. It is assumed that the flow turning angle for flow about a barrier is smaller than the critical angle, and subsonic regions are absent in the flow field studied. The effect of the characteristic parameters (Mach number at the nozzle exit, jet underexpansion) on the flow pattern and jet forces is analyzed. The results of numerical calculations are compared to the results of approximate theories and experimental data. A theoretical solution of the problem of the effect of a supersonic jet on a surface of given shape, even in the approximation of an inviscid, nonheat-conducting gas, is quite difficult. A reason for this is that the flow region contains shock waves interacting with each other, contact discontinuities, and zones of mixed sub-and supersonic flow. As far as is known to the authors, the results obtained for three-dimensional problems for the interaction of supersonic jets with each other or with barriers are primarily experimental (for example, [1–6]). A numerical analysis of the interaction of axisymmetric ideal-gas jets was carried out in [7–10]. In [7] a three-dimensional form of the method of characteristics was used to calculate the initial interaction region for two supersonic cylindrical jets (with Mach number M=10) intersecting at an angle of 60. The interaction of several jets has been considered in [8, 9], where the solution was obtained according to the Lax—Wendroff method without elimination of the discontinuity lines of flow parameters. In [10] the lateral interaction of axisymmetric supersonic jets with each other and with a plate is investigated by means of a straight-through calculationTranslated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 3–8, November–December, 1974.The authors thank A. N. Kraiko for useful discussions of the results, and A. L. Isakov and É. N. Gasparyan for kindly providing the experimental data.     

Groundwater flow in a medium with periodic inclusions

A study is made of the problem of two-dimensional confined steady flow through a porous reservoir whose percolation coefficient is a step function, while the homogeneity zones are rectangular inclusions (blocks) and a matrix (fractures). The velocity distributions, streamlines, isochrones of the motion of marked particles, and their time of passage through the elementary cell are found for this doubly periodic structure under conditions of continuity of the fluid head and the normal component of the flow on the phase contact boundary. Conclusions concerning the convective component of the longitudinal and transverse dispersion are drawn on the basis of these essentially two-dimensional hydrodynamic characteristics. The exact solutions obtained are compared with the results of numerical calculations carried out by the finite-difference method.Kazan'. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 139–148, September–October, 1995.     

Stationary efflux into a vacuum by a dual-temperature fully ionized plasma

It is shown that the macroscopic process of plasma discharge from an expanding nozzle is determined, when the thermal conductivity of electrons and heat transfer between the components are taken into account, by a unique dimensionless parameter: the adiabaticity parameter characterizing the transition from adiabatic flow of a dense plasma to the flow of comparatively rarefied plasma when the free path length of the particles is commensurate with the characteristic dimension of the nozzle. A numerical method is used to find the distribution of gas-dynamic and electrical parameters of the plasma stream, and the relationship between the generalized output parameters. It is shown that the energy associated with the ions at infinity, in the latter case, can be tens of times greater than the energy in adiabatic efflux, because of the high thermal conductivity with respect to electrons, but unrealistically large expansion of the nozzle is needed in order to attain it. Singular flow patterns occurring when stationary discharge of plasma at infinity is calculated are also discussed.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 10–20, November–December, 1971.In conclusion, the author thanks I. K. Fetisov for helpful discussion.