Quasi-three-dimensional calculation of flow in blade passages of turbines

The flow in turbomachines is currently calculated either on the basis of a single successive solution of an axisymmetric problem (see, for example, [1-A]) and the problem of flow past cascades of blades in a layer of variable thickness [1, 5], or by solution of a quasi-three-dimensional problem [6–8], or on the basis of three-dimensional models of the motion [9–11]. In this paper, we derive equations of a three-dimensional model of the flow of an ideal incompressible fluid for an arbitrary curvilinear system of coordinates based on averaging the equations of motion in the Gromek–Lamb form in the azimuthal direction; the pulsation terms are taken into account in the equations of the quasi-three-dimensional motion. An algorithm for numerical solution of the problem is described. The results of calculations are given and compared with experimental data for flows in the blade passages of an axial pump and a rotating-blade turbine. The obtained results are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 69–76, March–April, 1991.I thank A. I. Kuzin and A. V. Gol'din for supplying the results of the experimental investigations.     

Flow in the neighborhood of the stagnation line on a blunt body traveling at supersonic speed through a zone with elevated temperature and vibrational energy of the molecules

The unsteady flow in the neighborhood of the stagnation line on a sphere traveling at supersonic speed through a plane layer of diatomic gas with elevated temperature and nonequilibrium excitation of the molecular vibrations is investigated. (The source of the inhomogeneity could be a gas discharge [1].) The problem is solved using the viscous shock layer model which makes it possible to take molecular transport processes into account and analyze the unsteady heat transfer. Such flows were previously calculated in [2] within the framework of the inviscid gas model.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i GazNo. 3, pp. 183–185, May–June, 1990.     

Determination of relative permeabilities in a three-phase flow in a porous medium

A study is made of the problem of determining the parameters of flow described by the Buckley-Leverett system of equations by using functions that admit direct measurement. The well-known solution to the analogous problem for two-phase flow [1–3] is generalized. In contrast to [4], the general case is considered, when the fractions of the phases in the flow and the phase permeabilities depend on two variables.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 187–189, September–October, 1984.The author wishes to thank B. V. Shalimov for his helpful advice.     

Boundary conditions on a rigid surface in a two-phase flow

A study is made of the boundary conditions on a rigid surface in a two-component disperse flow. Appropriate boundary conditions are obtained for the kinetic equation and macroscopic equations of a pseudogas of solid particles proposed in [1–3]. The reasons for the occurrence of bubbles in two-phase systems are discussed. On the basis of the similitude parameters of the kinetic equation of the pseudogas, disperse systems are classified generally on the basis of the concentration of solid particles and their diameters.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 46–51, May–June, 1980.I thank V. P. Myasnikov for suggesting the problem and for a helpful discussion.     

The boundary layers of a fully ionized two-temperature plasma with given component temperatures at an electrode

The flow of a plasma with different component temperatures in the boundary layers at the electrodes of an MHD channel is investigated without any assumptions as to self-similarity. For the calculation of the electron temperature, the full energy equation for an electron gas [1] is solved with allowance for the estimates given in [2]. In contrast to [3, 4], the calculation includes the change in temperature of electrons and ions along the channel caused by the collective transport of energy, the work done by the partial pressure forces, and the Joule heating and the energy exchange between the components. The problem of the boundary layers in the flow of a two-temperature, partially ionized plasma past an electrode is solved in simplified form by the local similarity method in [5–7]. In these papers, either the Kerrebrock equation is used [5, 6] or the collective terms are omitted from the electron energy equation [7].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 3–10, September–October, 1972.The author thanks V. V. Gogosov and A. E. Yakubenko for interest in this work.     

Flow and Heat Transfer in Underexpanded Air Jets Issuing from the Sonic Nozzle of a Plasma Generator

We present the results of an experimental investigation and numerical simulation of the gasdynamic structure of underexpanded dissociated-air jets and the heat transfer in these strongly nonequilibrium flows under the test conditions realized in the 100-kW electrodeless VGU-4 plasma generator of the Institute for Problems in Mechanics of the Russian Academy of Sciences (IPM RAS). The flow and heat transfer analysis is carried out on the basis of measurements of the static pressure in the plenum chamber, at the sonic nozzle exit, and on the low-pressure chamber wall, the stagnation pressure on the jet axis using a Pitot tube, and the heat transfer at the stagnation points of water-cooled models placed along the jet axis. The numerical simulation, based on complete Navier-Stokes equations, includes the calculation of (1) equilibrium air plasma flows in the discharge channel of the VGU-4 plasma generator; (2) underexpanded nonequilibrium dissociated-air jet outflow into the ambient space; and (3) axisymmetric jet flow past cylindrical models.     

Investigation of hypersonic three-dimensional flow past blunt bodies within the framework of the parabolized Navier-Stokes equations

The flow of a homogeneous gas in a three-dimensional hypersonic viscous shock layer, which includes the shock wave structure, is examined within the framework of the parabolic approximation of the Navier-Stokes equations. The Navier-Stokes equations are simplified on the basis of the asymptotic analysis carried out in [1], are written in variables of the Dorodnitsyn type [2] and are solved by the method proposed in [3, 4] extended to the case of three-dimensional flows. The flow at zero angle of attack past elliptic paraboloids, two-sheeted hyperboloids and triaxial ellipsoids is calculated. Some results of investigating the flow past such bodies are presented. Flow past a sphere in the analogous approximation was considered in [5], where a comparison was also made with the solution of the complete Navier-Stokes equations [6, 7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 134–142, July–August, 1987.In conclusion, the authors wish to express their warm thanks to V. V. Lunev and G. A. Tirskii for useful discussion and valuable comments.     

Calculation of two-dimensional dispersion in a gas in unsteady flow in a porous medium

A two-dimensional problem of the flow of a gas containing an impurity through a porous medium is considered. At the initial time, the gas containing a uniformly distributed impurity is at a high pressure in a spherical cavity in a porous medium at a certain distance from a flat surface. It is assumed that for t > the motion of the carrier gas is described by the system of equations for flow in a porous medium and the dispersion of the impurity is described by the equations of convective diffusion and nonequilibrium adsorption. A numerical method for solving the problem is discussed. Some results of calculations are given. The influence of the flat surface on the flow of the gas and the dispersion of the impurity is analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 61–67, September–October, 1982.We thank V. N. Nikolaevskii for comments which permitted a significant improvement in the paper.     

The steady-state flow of a rarefied gas from a spherical source or sink

The method of characteristics is used to solve problems in the steady-state flows of a rarefied gas on the basis of approximating the kinetic equations. Numerical results are given for the solution of the problem of the flow from a spherical source or sink using the generalized Kruk equation [1] and approximating the Boltzmann equation by the method proposed by the author [2, 3], Various flow conditions are discussed: flow into a vacuum, flow into a flooded volume, flow from a sink.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 58–66, March–April, 1971.     

Chemical relaxation in a viscous shock

We consider the flow of a nonequilibrium dissociating diatomic gas in a normal compression shock with account for viscosity and heat conductivity. The distribution of gas parameters in the flow is found by numerically solving the Navier-Stokes and chemical kinetics equations. The greatest difficulty in numerical integration comes from the singular points of this system at which the initial conditions are given. These points lead to instability of the numerical results when the problem is solved by standard numerical methods. An integration method is proposed that yields stable numerical results-continuous profiles of the distribution of the basic gas parameters in the shock are obtained.We consider steady one-dimensional flow in which the gas passes from equilibrium state 1 to another equilibrium state 2, which has higher values for temperature, density, and pressure. Such a flow is termed a normal compression shock.The parameter distribution in normal shock for nonequilibrium chemical processes has usually been calculated [1–3] without account for the transport phenomena (viscosity, heat conduction, and diffusion). The presence of an infinitely thin shock front perpendicular to the flow velocity direction was postulated. It was assumed that the flow is undisturbed ahead of the shock front. The gas parameters (velocity, density, and temperature) change discontinuously across the shock front, but the gas composition does not change. The composition change due to reactions takes place behind the shock front. The gas parameter distribution behind the front was calculated by solving the system of gasdynamic and chemical kinetics equations using the initial values determined from the Hugoniot conditions at the front to state 2 far downstream.Several studies (for example, [4, 5]) do account for transport phenomena in calculating parameter distribution in a compression shock, but not for nonequilibrium chemical reactions. These problems are solved by integrating the Navier-Stokes equations continuously from state 1 in the oncoming flow to state 2 downstream.We present a solution to the problem of normal compression shock in nonequilibrium dissociating oxygen with account for viscosity and heat conduction using the Navier-Stokes equations.     

Eruption dynamics of high-viscosity gas-saturated magmas

A model of eruption, which is a variant of that described in [4] and takes into account the disequilibrium of the pressure in the bubble and in the liquid in the absence of total solidification is proposed. The fragmentation zone is simulated by a disintegration wave with allowance for the velocity and temperature nonequilibrium of the particles of the gas suspension formed and its polydispersity. On the basis of the model constructed steady-state magma flow calculations are made for a given pressure difference and channel length. The results of the calculations show that taking pressure nonequilibrium into account leads to a qualitatively new dependence of the flow rate on the governing parameters and makes it possible to propose a catastrophic eruption intensification mechanism different from that proposed in [3].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 49–60, March–April, 1993.     

Nonuniqueness of the solution of the problem of viscous interaction on an axisymmetric body

The article discusses solutions of the equations of the hypersonic boundary layer on an axisymmetric offset slender body (with a power exponent equal to 3/4), taking account of interactions with a nonviscous flow. It is shown that, in this case, the equations of the boundary layer have solutions differing from the self-similar solution corresponding to flow around a semi-infinite body. The solutions obtained are analogous to solutions for a strong interaction on a plate with slipping and triangular vanes [1–4], but are obtained over a wide range of values of the parameter of viscous interaction. An asymptotic solution is given to the problem with the approach to zero of the interaction parameter.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 41–47, September–October, 1973.The authors thank V. V. Mikhailova for discussion of the work and useful advice.     

Calculation of unsteady supersonic flow past a two-dimensional cascade of plates ween vorticity inhomogeneities of the flow act on it

In the framework of the linear theory of small perturbations the problem of unsteady subsonic flow past a two-dimensional cascade of plates has been considered in a number of papers. Thus, the unsteady aerodynamic characteristics of a cascade of vibrating plates were calculated in [1] by the method of integral equations, while the same method was used in [2, 3] to calculate the sound fields that are excited when sound waves Coming from outside or vorticity inhomogeneities of the oncoming flow act on the cascade. The problem of a two-dimensional cascade of vibrating plates in a supersonic flow was solved in [4, 5]. In [4] the solution was constructed on the basis of the well-known solution of the problem of vibrations of a single plate, while in [5] a variant of the method of integral equations was used which differed slightly from the usual formulation of this method [1–3]. The approach proposed in [5] is used below to calculate the unsteady flow past a two-dimensional cascade of plates in the case when vorticity inhomogeneities of a supersonic oncoming flow act on it. Equations are obtained for the strength of the unsteady pressure jumps arising in such a flow and the vortex wakes shed from the trailing edges of the plates. Examples of the calculations illustrating the accuracy of the method and its possibilities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp, 152–160, May–June, 1986.     

Variational model of a turbulent rotating flow

A model of effectively viscous turbulent flows satisfying the Navier-Stokes equations and certain slip conditions at the walls is analyzed. The turbulent viscosity is determined on the basis of the principle of minimum energy dissipation rate, whose significance and conditions of applicability are discussed in detail. A new separated turbulent flow model is outlined. The problem of turbulent flow in a porous rotating tube is solved. The existence of two metastable flow regimes is predicted: one with an axial circulation zone, the other straight-through. In the case of a strongly swirled flow the first of these has a greater probability of realization; however, as the rotation weakens, in a certain critical situation the circulation zone collapses, after which the flow can only be straight-through. Despite the absence of empirical content, every aspect of the proposed theory is in good agreement with the experimental research on vortex chamber flows.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 22–32, May–June, 1985.     

Estimates of nonequilibrium radiative heating of blunt bodies in a hypersonic flow of carbon dioxide gas

The solution to the problem of hypersonic nonequilibrium flow over cones with rounded noses [1, 2] is used to estimate the radiative heating of the surfaces of the bodies.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 158–160, November–December, 1979.     

Normal interaction of shock waves in the presence of vibrational relaxation

The effect of nonequilibrium physicochemical processes on the flow resulting from the normal collision and reflection of shock waves is studied by the example of nonequilibrium excitation of molecular oscillations in nitrogen. It is shown that the thermal effect of vibrational relaxation is small and the problem can be linearized around a known solution [1]. A similar approach to the solution of the problem of flow around a wedge and certain one-dimensional non-steady-state problems was used earlier in [2–4]. The solution of these problems was constructed in an angular domain, bounded by the shock wave and a solid wall (or the contact surface) and was reduced to a well-known functional equation [6]. The solution of this problem, because of the presence of two angular domains divided by a tangential discontinuity, reduces to a functional equation of more general form than in [6]. The results are obtained in finite form. In the special case of shocks of equal intensity, the normal reflection parameters are obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 90–96, July–August, 1976.     

Similarity autooscillations in a plane jet

Experimental investigations into the stability of a plane jet [1, 2] show that after the stationary flow has lost its stability a stable autooscillatory regime arises. In the present paper, an autooscillatory flow in a jet is studied theoretically on the basis of a plane-parallel flow in a fairly wide channel in the presence of a field of external forces. The external forces are such that at zero amplitude of the autooscillations they produce a Bickley—Schlichting velocity profile. The excitation of the secondary regimes is studied by the methods of bifurcation theory [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 26–32, May–June, 1979.We thank M. A. Gol'dshtik and V. N. Shtern for discussing the formulation of the problem and the results.     

Distribution of radiant thermal flux over the surface of a sphere for hypersonic flow of a nonviscous radiating gas past the sphere

In the present study using the Newtonian approximation [1] we obtain an analytical solution to the problem of flow of a steady, uniform, hypersonic, nonviscous, radiating gas past a sphere. The three-dimensional radiative-loss approximation is used. A distribution is found for the gasdynamic parameters in the shock layer, the withdrawal of the shock wave and the radiant thermal flux to the surface of the sphere. The Newtonian approximation was used earlier in [2, 3] to analyze a gas flow with radiation near the critical line. In [2] the radiation field was considered in the differential approximation, with the optical absorption coefficient being assumed constant. In [3] the integrodifferential energy equation with account of radiation was solved numerically for a gray gas. In [4–7] the problem of the flow of a nonviscous, nonheat-conducting gas behind a shock wave with account of radiation was solved numerically. To calculate the radiation field in [4, 7] the three-dimensional radiative-loss approximation was used; in [5, 6] the self-absorption of the gas was taken into account. A comparison of the equations obtained in the present study for radiant flow from radiating air to a sphere with the numerical calculations [4–7] shows them to have satisfactory accuracy.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 44–49, November–December, 1972.In conclusion the author thanks G. A. Tirskii and É. A. Gershbein for discussion and valuable remarks.     

Influence of nonequilibrium radiation on rarefied gas flow around blunt bodies

The hypersonic nonequilibrium rarefied gas flow is investigated in the neighborhood of the stagnation streamline ahead of a blunt body by taking into account nonequilibrium radiation due to electron excitation for air and carbon dioxide. The analysis is on the basis of a numerical solution of the Navier—Stokes equations simplified under the assumption of local self-similarity of the flow with the Shockwave structure taken into account. It is shown that at low densities, when the shock wave and shock layer thicknesses are of the same order of magnitudes, the two-layer Cheng model becomes inapplicable in the presence of radiation. In this case, the governing process is diffusion of the electronically excited molecules from the shock layer into the forward part of the shock front. The mechanism of the formation of a second luminous plateau on the diagram of the nonequilibrium radiation intensity density is discussed. The combined influence of the limit in collisions and the diffusion transport processes on the intensity of molecular band radiation is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 81–87, November–December, 1974.The author is grateful to V. V. Lunev for discussion and remarks during performance of the research.     

Unsteady aerodynamic interaction of two annular blade rows of thin lightly loaded blades rotating one relative to the other in a subsonic flow

The problem of subsonic unsteady ideal-gas flow over two annular blade rows of thin lightly loaded blades rotating one relative to the other is solved within the framework of linear small perturbation theory. As in the case of the interaction of two-dimensional cascades [1], the problem reduces to an infinite system of singular integral equations for the harmonic components of the oscillations in the distribution of the unknown aerodynamic load on one blade of each row. The system of integral equations for a finite number of harmonics is solved numerically by the collocation method. The kernels of the integral equations are regularized on the basis of the method proposed in [2].Translated from Izvestiya Akademii Nauk.SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 165–174, May–June, 1991.The author is grateful to A. A. Osipov and K. K. Butenko for their considerable assistance in the preparation of this paper.