Computation of the viscous shock layer on blunted cones

A numerical method is described for solving the equations of the compressible viscous shock layer on smooth spherically blunted axisymmetric cones at zero angle of attack and flow of a perfect gas. Effective use is made of the scheme of separating the original system of equations into parabolic (second order) and inviscid (first order) subsystems, which are solved by intrinsic methods. The results of the computations are presented. The method is capable of natural generalization to the case of nonequilibrium physical and chemical processes and diffusion. In most published papers dealing with computation of the compressible shock layer, the authors examine either the vicinity of the stagnation point or a certain region of spherical blunting [1–5]. In all the papers except [4, 5], a number of simplified assumptions have been made regarding the flow picture. Very few papers [6–8] have calculated the viscous shock layer on the forward surface of blunted bodies. In [6, 7] an approximate examination was made only of hyperboloids and paraboloids of revolution, which have very favorable geometry. Reference [8] used a approximate Karman—Polhausen integral method for a very simple system of equations. The method proposed here is essentially an accurate numerical method for solution of the viscous shock layer equations.     

Approximate symmetry group classification for a nonlinear fractional filtration equation of diffusion-wave type

A one-dimensional nonlinear fractional filtration equation with the Riemann–Liouville time-fractional derivative is proposed for modeling fluid flow through a porous medium. This equation is derived under an assumption that the fluid has a fractional equation of state in which the fluid density depends on the time-fractional derivative of pressure. The obtained equation belongs to the diffusion-wave type of equations. A case when the order of fractional differentiation is close to an integer number is considered, and a small parameter is introduced into the fractional filtration equation under consideration. An expansion of the Riemann–Liouville time-fractional derivative into the series with respect to this small parameter is obtained. Using this expansion, a first-order approximation of the derived fractional filtration equation is performed. Next, the problem of approximate Lie point symmetry group classification for this approximate nonlinear filtration equation with a small parameter is studied. It is shown that approximate symmetry groups admitted by different realizations of the approximate filtration equation have much more dimensions than the corresponding exact Lie point symmetry groups admitted by unperturbed fractional diffusion-wave equations. Obtained classification results permit to construct approximate invariant solutions for the considered nonlinear time-fractional filtration equations.     

Construction of approximate solutions of boundary value impact strain problems

The method for constructing approximate solutions of boundary value problems of impact strain dynamics in the form of ray expansions behind the strain discontinuity fronts is generalized to the case of curvilinear and diverging rays. This proposed generalization is illustrated by an example of dynamics of an antiplane motion of an elastic medium. The ray method is one of the methods for constructing approximate solutions of nonstationary boundary value problems of strain dynamics. It was proposed in [1, 2] and then widely used in nonstationary problems of mathematical physics involving surfaces on which the desired function or its derivatives have discontinuities [3–7]. A complete, qualified survey of papers in this direction can be found in [8]. This method is based on the expansion of the solution in a Taylor-type series behind the moving discontinuity surface rather than in a neighborhood of a stationary point. The coefficients of this series are the jumps of the derivatives of the unknown functions, for which, as a consequence of the compatibility conditions, one can obtain ordinary differential equations, i.e., discontinuity damping equations. In the case where the problem with velocity discontinuity surfaces is considered in a nonlinear medium, this method cannot be used directly, because one cannot obtain the damping equation. A modification of this method for the purpose of using it to solve problems of that type was proposed in [9–11], where, as an example, the solutions of several one-dimensional problems were considered. In the present paper, we show how this method can be transferred to the case of multidimensional impact strain problems in which the geometry of the ray is not known in advance and the rays become curvilinear and diverging. By way of example, we consider a simple problem on the antiplane motion of a nonlinearly elastic incompressible medium.     

Exponential filtration of a liquid in the presence of sources

A study is made of the plane exponential filtration of an incompressible liquid under the action of two sources (sinks). The solution is based on an S. A. Chaplygin transformation, the possibility of whose use in the investigation of nonlinear filtration was first noted in [1]. In [2–5] this transformation was used in a consideration of filtration with a limiting gradient. In the present article, another nonlinear law of resistance, an exponential law, is used to construct an exact solution. The use of S. A. Chaplygin variables makes it possible to transform the starting system of equations to a Helmholtz equation, which then reduces to a functional relation which is solvable by the Wiener-Hopf method. The results obtained point to the possibility of using the proposed method to solve other problems of plane exponential filtration, generated by sources or sinks, particularly when they are arranged symmetrically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 91–96, September–October, 1973.     

The spherical problem of the filtration of a liquid and a gas with a change in the permeability of rocks in accordance with an exponential law

The article considers the problem of the filtration of liquids (or gases), pumped through a borehole at a constant rate with elastic filtration conditions. The permeability of the stratum is assumed to be an exponential function of the coordinates. The viscosities of the injected and displaced liquids are assumed to be different. To increase the capacity of strata, i.e., of collectors used for the burial of industrial waste flows and gases, various methods are employed to increase the fracturing and the permeability of the rocks (hydro-pulse techniques, explosions, and other methods). As a result of this, a spherical region is formed in the rocks, in which the permeability varies along the radius. The character of this change is well described by an exponential function. The pumping of waste flows or industrial gases into such a cavity leads to the displacement of the stratum liquid (or gas). The problem of the displacement of one liquid by another liquid not miscible with it under rigid filtration conditions was first discussed in [1–5]. Here a study was made of a region of finite dimensions, bounded by two boundaries, with given pressures or mass flow rates (the linear and axisymmetric flow problems). The permeability of the stratum was assumed to be independent of the coordinates. A special characteristic of these problems is the fact that it is impossible to consider unbounded or semi-bounded filtration conditions in them since, under rigid filtration conditions, the condition of bounded character of the pressure (the head) is not satisfied at infinity. Elastic filtration conditions for two immiscible liquids were first discussed in [6], and later in [7, 8] and other reports. Here an investigation was made of the linear and axisymmetric problems for an unbounded region. In [9, 10] solutions are given to some problems with spherical symmetry for an unbounded region, with rigid filtration conditions and a jumpwise change of the permeability along the radius. In the problems of [6–10] the condition of the bounded character of the pressure is satisfied. In [11] the case of a hyperbolic change in the permeability of the rocks is considered.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 42–51, November–December, 1974.     

Method of outer and inner asymptotic expansions in the theory of Brownian motion of aerosol particles

We examine the Brownian motion of particles in a gaseous medium, complicated by the influence of inertial forces. The equation for the distribution function in phase space describing motion of this type was obtained in [1]. Also presented in [1] are the solutions of this equation for certain simple particular cases. The approximate equations of motion of aerosol particles in coordinate space were first obtained in [2] and solved for certain concrete problems in [3,4]. More exact equations of motion in coordinate space, and also the limits of applicability of the equations of [2], are presented in [5].     

Calculation of principal characteristics of turbulent streams in a state of structural equilibrium

A review of articles on the study of turbulent streams having transverse displacement, in which a turbulent energy balance equation is used, is contained in [1]. Levin [2] proposed a certain development of Rotta's method [3] making it possible to determine the characteristics of the average flow and the radial distribution of pulsation magnitudes. However, in this article the scale of the turbulence (the quantityl) was given as an empirical function of the coordinates. At the same time it is clear that the distribution of the turbulence scale depends on the conditions of the problem. A special differential equation proposed in [4,5] describing the variation in time and space of the quantityl has the drawback that in deriving this equation it is necessary to invoke additional hypotheses which are difficult to test experimentally. In the present article, along with the velocity of the average flow, the pressure, and the pulsation magnitudes, the scale of the turbulence is considered as an important characteristic of the stream, determined by the reference system which consists of the Reynolds equations, continuity equations, and equations for the component of the Reynolds stress tensor. Rotta's approximate semiempirical relations and an experimental relation for the single-point correlation coefficient between the turbulent pulsations in velocity are used for closure of the system obtained. An approximate calculation is given for the principal average and pulsation characteristics of the flow for the region of the stream where the turbulence is in a state of structural equilibrium [6]. A comparison of the calculated and experimental data is presented.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 95–99, January–February, 1973.     

Steady seepage flow under aprons in anisotropic ground

Seepage under a flat apron in a two-layer ground consisting of homogeneous isotropic layers of equal thickness was investigated by means of the analytic theory of linear differential equations by Polubarinova-Kochina [1]. An approximate solution of the same problem for layers of different thickness by the methods of the variational calculus was given in [2]. The solution to such a problem for homogeneous anisotropic layers of different thickness by means of a special complex variable representation and linear conjugation was given in [3] by Prusov and the author of the present paper. The method of [3] is now used to investigate two-dimensional steady seepage flow under two aprons in a multilayer ground consisting of n homogeneous anisotropic layers of different thickness, the principal directions of anistropy of each layer having an arbitrary position relative to the horizontal.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 155–158, October–December, 1981.     

Free Fluid Convection in a Closed Contour in a Heat-Conducting Half-Space

The steady-state convection of a fluid in a thin porous vertical ring located in a heat-conducting half-plane is considered. For this problem approximate equations are derived. For a circular ring an analytic solution is obtained. For an elliptic ring a numerical-analytic solution is found. The Nusselt number and the fluid flow rate as functions of the Rayleigh number, the aspect ratio, and the contour depth are investigated.Many studies have been devoted to fluid convection in a porous ring [1–3]. In [1] two-dimensional convection with an isothermal internal boundary was considered when a temperature stratification is given on the outer boundary. A feature of this problem is the fact that the ring is located inside an impermeable heat-conducting medium in which a thermal gradient directed vertically downward is specified at a large distance from the ring. In [2, 3] two-dimensional convection in an annular ring occupied by a porous medium was investigated. From the results obtained in these studies it follows that in the formulation considered the hydraulic approximation can be used with satisfactory accuracy. In the present study this question is discussed more concretely and the necessary estimates are found. The results obtained could be useful for investigating hydrothermal convection in the Earth's crust, which has important geophysical applications [4–6].     

Dynamics of a flat cavity in a low-viscosity liquid

The problem of the motion of a cavity in a plane-parallel flow of an ideal liquid, taking account of surface tension, was first discussed in [1], in which an exact equation was obtained describing the equilibrium form of the cavity. In [2] an analysis was made of this equation, and, in a particular case, the existence of an analytical solution was demonstrated. Articles [3, 4] give the results of numerical solutions. In the present article, the cavity is defined by an infinite set of generalized coordinates, and Lagrange equations determining the dynamics of the cavity are given in explicit form. The problem discussed in [1–4] is reduced to the problem of seeking a minimum of a function of an infinite number of variables. The explicit form of this function is found. In distinction from [1–4], on the basis of the Lagrauge equations, a study is also made of the unsteady-state motion of the cavity. The dynamic equations are generalized for the case of a cavity moving in a heavy viscous liquid with surface tension at large Reynolds numbers. Under these circumstances, the steady-state motion of the cavity is determined from an infinite system of algebraic equations written in explicit form. An exact solution of the dynamic equations is obtained for an elliptical cavity in the case of an ideal liquid. An approximation of the cavity by an ellipse is used to find the approximate analytical dependence of the Weber number on the deformation, and a comparison is made with numerical calculations [3, 4]. The problem of the motion of an elliptical cavity is considered in a manner analogous to the problem of an ellipsoidal cavity for an axisymmetric flow [5, 6]. In distinction from [6], the equilibrium form of a flat cavity in a heavy viscous liquid becomes unstable if the ratio of the axes of the cavity is greater than 2.06.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 15–23, September–October, 1973.The author thanks G. Yu. Stepanov for his useful observations.     

Fluid flow in an arbitrary cascade on an axisymmetric stream surface in a variable thickness layer

A solution is given for the problem of flow past a cascade on an axisymmetric stream surface in a layer of variable thickness, which is a component part of the approximate solution of the three-dimensional problem for a three-dimensional cascade. Generalized analytic functions are used to obtain the integral equation for the potential function, which is solved via iteration method by reduction to a system of linear algebraic equations. An algorithm and a program for the Minsk-2 computer are formulated. The precision of the algorithm is evaluated and results are presented of the calculation of an example cascade.In the formulation of [1, 3] the problem of flow past a three-dimensional turbomachine cascade is reduced approximately to the joint solution of two-dimensional problems of the averaged axisymmetric flow and the flow on an axisymmetric stream surface in an elementary layer of variable thickness.In the following we solve the second problem for an arbitrary cascade with finite thickness rotating with constant angular velocity in ideal fluid flow: the solution is carried out on a Minsk-2 computer.Many studies have been devoted to this problem. A method for solving the direct problem for a cascade of flat plates in a hyperbolic layer was presented in [2]. Methods were developed in [1, 3] for constructing the flow for the case of a channel with variable thickness; these methods are approximately applicable for dense cascades but yield considerable error for small-load turbomachine cascades. The solution developed in [4], somewhat reminiscent of that of [2], is applicable for thin, slightly curved profiles in a layer with monotonically varying thickness. A solution has been given for a circular cascade for layers varying logarithmically [5] and linearly [6]. Approximate methods for slightly curved profiles in a monotonically varying layer with account for layer variability only in the discharge component were examined in [7–9]. A solution is given in [10] for an arbitrary layer by means of the relaxation method, which yields a roughly approximate flow pattern. The general solution of the problem by means of potential theory and the method of singularities presented in [11] is in error because of neglect of the crossflow through the skeletal line. The computer solution of [12] contains an unassessed error for the calculations in an arbitrary layer. The finite difference method is used in [13] to solve the differential equation of flow, which is illustrated by numerical examples for monotonie layers of axial turbomachines. The numerical solution of [13] is very complex.The solution presented below is found in the general formulation with respect to the geometric parameters of the cascade and the axisymmetric surface and also in terms of the layer thickness variation law.The numerical solution requires about 15 minutes of machine time on the Minsk-2 computer.     

Three-dimensional fully established flow of an ideal liquid around the blade systems of turbomachines

The present-day methods for the quasi-three-dimensional calculation of flow around the working organs of turbomachines were developed in [1, 2]. The most widely used in practice is the method of [2] (see, for example, [3, 4]), which, in a complete statement, proposes the solution of three known two-dimensional problems at the coordinate surfaces of a triorthogonal natural system of coordinate surfaces in the flow-through part of a turbomachine. A method is proposed below for the hydrodynamic calculation of blade systems directly in the three-dimensional region, based on solutions of the integral equations of the three-dimensional theory of a field, and not connected with any kind, of metric schematization of the flow of the liquid. Comparative analyses show that, with the use of an electronic computer, the present method and the method of [2] (in the complete statement) are practically equivalent.     

Calculation of the flow of an ideal fluid past a wing of finite thickness

The problem of irrotational flow past a wing of finite thickness and finite span can be reduced by Green's formula to the solution of a system of Fredholm equations of the second kind on the surface of the wing [1]. The wake vortex sheet is represented by a free vortex surface. Besides panel methods (see, for example, [2]) there are also methods of approximate solution of this problem based on a preliminary discretization of the solution along the span of the wing in which the two-dimensional integral equations are reduced to a system of one-dimensional integral equations [1], for which numerical methods of solution have already been developed [3–6]. At the same time, a discretization is also realized for the wake vortex sheet along the span of the wing. In the present paper, this idea of numerical solution of the problem of irrotational flow past a wing of finite span is realized on the basis of an approximation of the unknown functions which is piecewise linear along the span. The wake vortex sheet is represented by vortex filaments [7] in the nonlinear problem. In the linear problem, the sheet is represented both by vortex filaments and by a vortex surface. Examples are given of an aerodynamic calculation for sweptback wings of finite thickness with a constriction, and the results of the calculation are also compared with experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 124–131, October–December, 1981.     

Asymptotic analysis of transonic flows

In this paper we derive the equations of the second and third approximations for the stream function of two-dimensional and axisymmetric potential transonic flow of an inviscid gas and find their particular solutions corresponding to certain transonic flows.A similar study concerning the second approximation of subsonic and supersonic flow was made by Van Dyke [1] and Hayes [2]. The second approximation for the velocity potential of transonic flow has been examined in detail by Hayes [3]. Euvrard [4, 5] has investigated the asymptotic behavior of transonic flow far from a body, while Fal'kovich, Chernov, and Gorskii [6] have studied the flow in a nozzle throat.The transonic asymptotic analysis for the stream function is presented in this paper.     

板的非线性热弹耦合振动（Ⅰ）：近似解析解

本文以文[2,3,4]为基础,导出了板的热弹耦合非线性振动控制方程,在采用Galerkin法离散化以后,按各个变量性质分别用多尺度法或正则摄动法求得近似解析解。籍此可揭示系统各参数对非线性热弹耦合振动影响的机理和作出必要的近似计算,对工程实际具有较大的参考价值。