Nonlinear problem of flow past a source with the formation of a stagnant zone given different Bernoulli constants

The problem of free flow past a point source is considered for two streams with different Bernoulli constants whose encounter creates a bounded region of constant pressure. The theory and method of solving problems of plane ideal jet flows with different Bernoulli constants in the jets were developed in [1]. Here, in conformity with [1], a nonlinear system of equations is derived, the question of the construction of a high-accuracy numerical solution is considered, and certain calculation results are presented for various values of the Bernoulli and cavitation numbers, which are dimensionless parameters of the problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 55–60, May–June, 1986.     

Stability of nonplane-parallel three-dimensional jet flows

The problem of the stability of nonplane-parallel flows is one of the most difficult and least studied problems in the theory of hydrodynamic stability [1]. In contrast to the Heisenberg approximation [1], the basic state whose stability is investigated depends on several variables, and the stability problem reduces to the solution of an eigenvalue problem for partial differential equations in which the coefficients depend on several variables [2–7]. In the case of a periodic dependence of these coefficients on the time [2] or the spatial coordinates [3, 4], the analog of Floquet theory for the partial differential equations is constructed. With rare exceptions, the case of a nonperiodic dependence has usually been considered under the assumption of weak nonplane-parallelism, i.e., a fairly small deviation from the plane-parallel case has been assumed and the corresponding asymptotic expansions in the linear [6] and nonlinear [7] stability analyses considered. The present paper considers the case of an arbitrary dependence of the velocity profile of the basic flow on two spatial variables. The deviation from the plane-parallel case is not assumed to be small, and the corresponding eigenvalue problem for the partial differential equations is solved by means of the direct methods of [5], which were introduced for the first time and justified in the theory of hydrodynamic stability by Petrov [8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 21–28, May–June, 1987.     

Computation of the pressure on the surface of a fuselage with a wing in an ideal fluid

The flows around complex three-dimensional bodies by an ideal fluid are computed by methods [1–3] using approximation of the surface by a set of plane elements. A layer of surface singularities, whose intensity is found by solving a system of linear algebraic equations of very high order, is distributed continuously over each element. Evaluation of the system coefficients and its solution require significant machine time expenditures on powerful electronic computers. If in the method of [2] the total system of equations is separated successfully into several subsystems by simplifications and an approximate solution of the problem is obtained more rapidly than by the method in [1], then the same author practically used the method in [1] to design specific fuselages in [3]. A method [4] developed for a fuselage is expanded in this paper to design a wing-fuselage combination. This method turns out to be less tedious, without being inferior in accuracy, by being different from the method in [1] in the means of solving the fundamental integral equation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 110–115, May–June, 1977.     

The zonal linearization method for multidimensional mass transfer problems in porous media

A number of methods have been proposed in recent years for calculating the combined flows of immiscible and miscible liquids in strata to systems of boreholes. We propose a method which can naturally be called the zonal linearization method [1]. It is more compact than the usual finite-difference method and has high accuracy, in particular, in the neighborhood of a borehole, since it is closely similar to the method of characteristics. The method can be applied to both continuous and discontinuous flows and in principle makes it possible to investigate the formation and breakdown of discontinuities. As distinct from the method of characteristics, it is well suited to programming and implementation on a computer, and it also makes it possible to obtain an approximate analytic solution of the problem in many cases and to estimate the accuracy of the solution. The method is based on the zonal linearization of the equation for mass conservation in the total flow between chosen surfaces or contour lines (lines of equal saturation or concentration). Determination of the dynamics of the contour surfaces leads to a Cauchy problem for a system of integrodifferential equations involving partial derivatives. The zonal linearization method is a development of the scheme described in [2–4], and the method of solving the Cauchy problem is a generalization of the methods described in [4–13]. The essence of the method and its convergence are illustrated by means of two-dimensional problems in two-phase filtration.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 66–80, July–August, 1973.     

On radiative heat transfer in a plane layer of an absorbing medium

Recently there has arisen increased interest in the study of radiative heat transfer between geometrically simple systems, both as autonomous problems and as elements of more complex problems.Problems of this kind have been treated by many authors [1–111 who have considered gray, diffusely emitting and absorbing boundaries and gray nonscattering media. In most cases these investigations were restricted either to the derivation of approximate formulas for the net radiative flux, without an exact analysis of the temperature distribution in the layer [5–7], or to numerical computation [1–4], In the latter case, with the exception of [8], which contains a numerical analysis for the case of optical symmetry, no attempt was made to analyze the effect of the optical properties of the boundaries on the temperature field in the layer.These papers can be divided into two groups according to the method of analysis used. The first group includes papers based on the integral equations of radiative transfer, with the corresponding integral analytical methods [1, 2], Similar in nature are [3, 4] which use the slab method, applicable to electrical-analog computation, as well as a recent paper [8] based on probability methods.The second group of papers [5–7] is based on the so-called differential methods. Of particular interest is [7], which develops these methods to an advanced degree. In several papers the problem of radiative transfer is analyzed in conjunction with more complex problems (cf., e.g. [10, 11]).In the present work we shall attempt to carry out an approximate analytical study of problems connected with radiative heat transfer in a plane layer of an absorbing, emitting, nonscattering gray medium with temperature-independent optical properties. The layer is bounded by two parallel, diffusely emitting and diffusely reflecting, isothermal, gray planes.The paper presents the fundamental formulation of the problem, which consists in: (a) the determination of the net heat flux on the basis of given temperature distribution (direct formulation), and (b) the determination of the temperature distribution on the basis of given distribution of the net radiative heat source per unit volume and boundary temperatures (inverse formulation). The analysis is based on integral methods appropriate to the integral equations which represent the net total and hemispherical radiation flux densities [12].The author would like to thank S. S. Kutateladze for his interest in this work.     

Boundary-value problems with external matching conditions and their application to flow through porous media

A method of solving plane problems of flow through soils with curved level lines and lines of discontinuity of the permeability function, fissures and curtain walls is proposed with reference to the example of m inhomogeneous zones separated by ellipses. The method is based on the solution of boundary-value problems with external matching conditions and is more efficient than the method of constructing flows on Riemann surfaces for two homogeneous zones [1], the method of reducing problems for homogeneous zones to the solution of a system of integral equations [2] and the circle theorem method for four homogeneous zones [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 181–183, July–August, 1991.     

Some problems in the theory of plane jets of conducting fluid

Using the boundary-layer equations as a basis, the author considers the propagation of plane jets of conducting fluid in a transverse magnetic field (noninductive approximation).The propagation of plane jets of conducting fluid is considered in several studies [1–12]. In the first few studies jet flow in a nonuniform magnetic field is considered; here the field strength distribution along the jet axis was chosen in order to obtain self-similar solutions. The solution to such a problem given a constant conductivity of the medium is given in [1–3] for a free jet and in [4] for a semibounded jet; reference [5] contains a solution to the problem of a free jet allowing for the dependence of conductivity on temperature. References [6–8] attempt an exact solution to the problem of jet propagation in any magnetic field. An approximate solution to problems of this type can be obtained by using the integral method. References [9–10] contain the solution obtained by this method for a free jet propagating in a uniform magnetic field.The last study [10] also gives a comparison of the exact solution obtained in [3] with the solution obtained by the integral method using as an example the propagation of a jet in a nonuniform magnetic field. It is shown that for scale values of the jet velocity and thickness the integral method yields almost-exact values. In this study [10], the propagation of a free jet is considered allowing for conduction anisotropy. The solution to the problem of a free jet within the asymptotic boundary layer is obtained in [1] by applying the expansion method to the small magnetic-interaction parameter. With this method, the problem of a turbulent jet is considered in terms of the Prandtl scheme. The Boussinesq formula for the turbulent-viscosity coefficient is used in [12].This study considers the dynamic and thermal problems involved with a laminar free and semibounded jet within the asymptotic boundary layer, propagating in a magnetic field with any distribution. A system of ordinary differential equations and the integral condition are obtained from the initial partial differential equations. The solution of the derived equations is illustrated by the example of jet propagation in a uniform magnetic field. A similar solution is obtained for a turbulent free jet with the turbulent-exchange coefficient defined by the Prandtl scheme.     

Generalized variables in boundary-layer theory

Independent variables are widely used in boundary-layer theory to construct efficient methods of solving problems. The Dorodnitsyn variables in Lees' form [1] are the most common and general. This form combines the transformations proposed by Dorodnitsyn [2], Blasius [3], and Mangler-Stepanov [4, 5]. As is well known, transformation of the boundary-layer equations to Dorodnitsyn variables in Lees' form leads to a generalized single system of equations describing plane and axisymmetric gas flows. An analogous generalization of the Mises [6] and Crocco [7] variables is carried out below.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 166–168, September–October, 1976.     

Flow of a viscous fluid in a cylindrical vessel with a rotating cover

The problem considered arises in solving various technical problems associated with flows of a viscous fluid in a closed space near rotating plane surfaces, turbomachine disks, thrust bearings, rotational viscosimeters, etc. The approximate solution of the problem on the basis of a simplified flow scheme was first obtained by Schultz-Grunow [1], The most complete investigation has been made recently by Grohne [2], who outlined a program for solving the problem by joining several partial solutions on the basis of definite hypotheses concerning the flow core.With the development of electronic digital computers and the necessary numerical methods, the most effective means of solving the considered problem is the use of the grid methods for solving partial differential equations. The present paper is devoted to presenting the results of the solution of the problem using the grid method on a digital computer.     

Solution to the problem of a swirling jet from an annular orifice

The problem of the propagation of a laminar immersed fan jet with swirling was considered in [1–3]. In [1], the jet source scheme was used to find a self-similar solution for a weakly swirling jet. An attempt to solve by an integral method the analogous problem for a jet emanating from a slit of finite size was made in [2]. In [3], the equations of motion for a jet with arbitrary swirling were reduced under a number of assumptions to the equations that describe the flow of a flat immersed jet. This paper gives the numerical solution to the problem of the propagation of a radial jet emanating with arbitrary swirling from a slit of finite size and an analytic solution for the main section of the jet.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 49–54, March–April, 1991.     

Dynamic problem for a two-layer medium with a vibrational source at the boundary interface

The present study is concerned with an analysis of gravitational and acoustic waves which are excited by a vibrational source deeply placed in a liquid covered by ice. An analysis of the rigidity characteristics of ice modeled by an elastic layer or by a Kirchhoff plate is done by factorization of the solution to the integral equation equivalent to an initially combined boundary value problem. The uncombined boundary condition is used to solve problems for unrestricted ice fields in [1–3], whereas combined conditions with vibrational sources positioned at the boundary of the medium are used in [4].Translated from Zhurnal Prikladnoi Mekhaniki, No. 3, pp. 125–129, May–June, 1986.     

Asymptotic method for analysis of a conical shell under local loading

Differential equations of the general theory describing the stress-strain state of conical shells are very complicated, and when computing the exact solution of the problem by analytic methods one encounters severe or even so far insurmountable difficulties. Therefore, in the present paper we develop an approach based on the method of asymptotic synthesis of the stressed state, which has already proved efficient when solving similar problems for cylindrical shells [1, 2]. We essentially use the fourth-order differential equations obtained by Kan [3], which describe the ground state and the boundary effect. Earlier, such equations have already been used to solve problems concerning force and thermal actions on weakly conical shells [4–6]. By applying the asymptotic synthesis method to these equations, we manage to obtain sufficiently accurate closed-form solutions.     

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.     

Method of generalized similitude in free convection problems with arbitrary distribution of the temperature or heat flux on a vertical wall

The self-similar problem of free convection near a heated vertical plate was solved for the first time in [1] for the simplest case of a constant wall temperature. In [2], Yang proved the existence of a self-similar solution to the problem of free convection for vertical plates and cylinders on the surfaces of which the temperature has a power-law distribution. In [3], Yang's solution was generalized to the case of free convection near a slender figure of revolution, but also only in the self-similar case of a power-law distribution of the temperature on the wall. In [4], this problem was solved in an extended nonsimilar formulation but by an artificial and not general method similar to Gertler's, the convergence of the approximations being slow. The present paper contains the solution to the problem of free convection near a vertical plate with arbitrary distribution of the temperature or heat flux on its surface. Rigorous application of the method of generalized similitude [5] leads in this case to universal equations that present insuperable computational difficulties, which forces one to use a simplified but fairly general method to solve this class of problems.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 167–170, May–June, 1980.I thank L. G. Loitsyanskii and E. M. Smirnov for discussing the results and for valuable comments.     

Calculation of the unsteady-state aerodynamic characteristics of three-dimensional bearing systems

At the present times effective numerical methods have been developed for solving both steady-state and unsteady-state linear problems on the determination of the aerodynamic characteristics of wings of complex form, in a plan view [1–3]. The transition to three-dimensional systems requires the development of new methods. There are a number of known investigations in this direction, for simple three-dimensional configurations [4–7]. The present article proposes a method for solving supersonic problems of the determination of the steady-state and unsteady-state three-dimensional bearing systems. It is a development of known methods for calculating wings of complex form in a plan view. The most effective route is one based on a solution of the problem for stepwise dependences on the time. After this, the transition to any other laws is effected using a packet integral.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 173–176, January–February, 1976.     

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.     

Nonlinear problems of the thermal conductivity equation

The asymptotic behavior of solutions of parabolic equations at infinite times has been investigated for various cases [1–6]. Two initial boundary-value problems are considered in this paper. The solution of the thermal conductivity equation with a nonlinear right-hand side is found, including also nonlinear boundary conditions. It is shown that the solution of the corresponding problem tends either to a stable, steady-state solution, or to a periodic solution, depending on the initial values of the functions and constants appearing in the conditions of the problem. Other papers [7, 8] are devoted to finding the periodic solutions of these two problems encountered in hydrodynamics (diffusion, underground hydrodynamics), and to studying the asymptotic behavior of the corresponding initial boundary problems.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 123–128, May–June, 1972.     

Calculation of aerodynamic characteristics of arrays of profiles of arbitrary shape

It is well known that the problem on nonseparating potential flow of an incompressible fluid about an array of profiles reduces to an integral equation for a certain real function, determined on the contours of the profiles of the array. As such a function one can take, as was done, for instance, in [1–5], the relative velocity of the fluid on the profiles of the array. For arrays of profiles of arbitrary shape it is necessary to solve the corresponding integral equation numerically. In the particular examples of the calculation of aerodynamic arrays that are available [1–3] the numerical methods used were based on the approximate evaluation of contour integrals by rectangle formulas. As investigations showed, sizeable errors arose thereby in the approximate solution obtained, these being especially significant in the case of curved profiles of relatively small bulk. In the present paper a method for the numerical solution of the integral equation obtained in [5] is proposed. The method is based on the replacement of a profile of the array with an inscribed N polygon, the length of whose sides is of the order N^–1 and whose internal angles are close to . Convergence with increasing N of the numerical solution to an exact solution of the integral equations at the reference points is demonstrated. Examples of the calculation are given.Novosibirsk. Translated from Izvestiya Akademii Nauk SSSR. Mekhanika Zhidkosti i Gaza, No. 2, pp. 105–112, March–April, 1972.     

Nonlinear dispersion of long internal waves

The propagation of long weakly nonlinear waves in an atmospheric waveguide is considered. A model system of Kadomtsev-Petviashvili equations [1], which describes the propagation of such waves, is derived. In the case of one excited wave mode the system of model equations goes over into the Kadomtsev-Petviashvili equation, in which, however, the variables x and t are interchanged. The reasons for this are clarified. In the two-dimensional case an approximate solution of the model equations is constructed, and steady nonlinear waves and their interaction in a collision are considered. The results of a numerical verification of the stability of the approximate steady solutions and of the solution to the problem of decay of the wave into quasisolitons are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 151–157, May–June, 1988.     

Cavitation streamline flow of lamina in transverse gravitational field

The problem of cavitation streamline flow located on the linear base of a lamina in a gravity solution current is solved by the systems of Ryabushinskii and Zhukovskii-Roshko. The method of fragment-continuum approximation of the boundary condition at the free boundary was used, in which this condition is exactly satisfied at a finite number of points. In this way the original problem comes down to a solution of a system of nonlinear equations whose solvability can be shown by the method of V. N. Monakhov [1]. The main consideration in the present work was given to a numerical solution of this system of equations on a computer. The problem is similar to the type for large Froude numbers, when the effect of weight on the flow is small, studied in [2-5]. In [6, 7] the flow problems were solved by the method of finite differences. The approximations of the boundary condition at the free boundary used earlier are based on the use of the smallness of these or other characteristics of flow. Thus, for example, the linearization of Levi-Chivit [8] is rightly used in the assumption of smallness of the change in the modulus and angle of inclination of the velocity at the free flow line; a stronger linearization is based on the requirement of smallness of additional velocities caused by an obstacle in comparison with the velocity of the undisturbed current [9]. In the given work the problems studied lead to a range of cavitation and Froude numbers when the gravitational force exerts a considerable effect on the main characteristics of the flow. As an example of one of the possible applications of the calculation, the solution of the problem of choice of the form of a body of zero buoyancy with a zone of constant pressure is given.Translated from Zhurnal Prikladnoi Mekhanik i Tekhnicheskoi Fiziki, No. 5, pp. 132–136, September–October, 1971.