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.     

Optimization of body shape at small reynolds numbers

The problem of the optimization of the shape of a body in a stream of viscous liquid or gas was treated in [1–5]. The necessary conditions for a body to offer minimum resistance to the flow of a viscous gas past it were derived in [1], The necessary optimality conditions when the motion of the fluid is described by the approximate Stokes equations were derived in [2], The shape of a body of minimum resistance was found numerically in [3] in the Stokes approximation. The optimality conditions when the motion of the fluid is described by the Navier—Stokes equations were derived in [4, 5], and in [4] these conditions were extended to the case of a fluid whose motion is described in the boundary-layer approximation. The necessary optimality conditions when the motion of the fluid is described by the approximate Oseen equations were derived in [5] and an asymptotic analysis of the behavior of the optimum shape near the critical points was performed for arbitrary Reynolds numbers.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp, 87–93, January–February, 1978.     

Choice of a self-similar solution in boundary-layer theory

If the speed of the outer flow at the edge of the boundary layer does not depend on the time and is specified in the form of a power-law function of the longitudinal coordinate, then a self-similar solution of the boundary-layer equations can be found by integrating a third-order ordinary differential equation (see [1–3]). When the exponent of the power in the outerflow velocity distribution is negative, a self-similar solution satisfying the equations and the usually posed boundary conditions is not uniquely determinable [4], A similar result was obtained in [5] for flows of a conducting fluid in a magnetic field. In the present paper we study the behavior of non-self-similar perturbations of a self-similar solution, enabling us to provide a basis for the choice of a self-similar solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 42–46, July–August, 1974.     

Aerodynamics of end-wall boundary layer in a vortex chamber

The need for the inclusion of end-wall boundary layers in the study of the aerodynamics of vortex chambers has been frequently mentioned in the literature. However, owing to limited experimental data [1–3] with reliable information on the wall layers, the existing computational methods for end-wall boundary layers are not well-founded. The question of which parameters determine the formation of end-wall flow remains debatable. In some studies [4, 5], the vortex chambers are conditionally divided into short and long chambers. However, there is no unique opinion on the role of end-wall flows in vortex chambers of different lengths. It has also not been established for what geometric and flow parameters the chamber could be considered long or short. In the present study, as in [1, 5–8], solution is obtained for the end-wall boundary-layer equations using integral methods, considering the boundary layer in the radial direction in the form of a submerged wall jet. Such an approach made it possible to use the laws for the development of wall jets [9], and obtain fairly simple relations for integral parameters, skin friction, mass flow in the boundary layer, and other characteristics. Results are compared with available experimental data and computations of others authors; turbulent flow is considered; results for laminar boundary layer are given in [10].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 117–126, September–October, 1986.     

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.     

Small vibrations of a sphere in a spherical volume of viscous fluid

Problems of the vibration of bodies in confined viscous fluids have been solved to determine the added masses and damping coefficients of rods [1–3] and floats [4–5]. The solutions of these problems, based on the use of simplifications of the boundary-layer method [4–6], are obtained analytically in general form and are in good agreement with the experimental data. However, in each specific case the possibility of using such solutions for given values of the fluid viscosity and vibration frequency must be justified either experimentally [2, 4, 5] or theoretically as, for example, in [1], where an analytic solution was obtained for concentric cylinders. The present paper offers a general solution of the problem of the small vibrations of a sphere in a spherical volume of fluid valid over a broad range of variation of the dimensionless kinematic viscosity. The limiting cases of this solution for both high and low viscosity are considered. The asymptotic expressions obtained are compared with calculations based on the analytic solution.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 29–34, March–April, 1986.     

Stability of heat transfer regimes at the front stagnation point of a body in a stream of dissociated gas

The formulation and solution of the stationary problem of heat transfer in the neighborhood of the front point of a body at constant temperature in a stream of dissociated air are given in [1]. In [2], the results are given of numerical solution of this problem in the nonstationary formulation; the establishment of a stationary heat transfer regime was established for all the calculated variants. In the present paper, we investigate the stability of stationary heat transfer regimes at the front stagnation point of a body in a stream of dissociated air using the Lyapunov functional method [3, 4] and the method of [2, 5], which is based on the use of Meksyn's method in boundary-layer theory [6, 7]. It is established that an arbitrarily strong growth of the Damköhler number does not lead to instability and multiplicity of the stationary regimes, in contrast to the case when a hot mixture of gases flows over the front point of a thermostat [2, 5, 8]. Numerical solution of the boundary-layer equations for a wide range of Damköhler numbers confirms the results of the approximate qualitative analysis and shows that in a number of cases the time of establishment of the stationary state is a nonmonotonic function of the Damköhler number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 97–106, September–October, 1979.     

Self-similar solutions of three-dimensional laminar magnetohydrodynamic boundary-layer equations

Self-similar solutions of three-dimensional boundary-layer equations of an incompressible fluid in ordinary hydrodynamics were considered in [1–3] et al. The present work looks for self-similar solutions of three-dimensional magnetohydrodynamic boundary-layer equations.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 9, No. 4, pp. 10–17, July–August, 1968.     

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.     

Rayleigh stability of shear flow in relaxing media

The stability of steady-state flow is considered in a medium with a nonlocal coupling between pressure and density. The equations for perturbations in such a medium are derived in the linear approximation. The results of numerical integration are given for shear motion. The stability of parallel layered flow in an inviscid homogeneous fluid has been studied for a hundred years. The mathematics for investigating an inviscid instability has been developed, and it has been given a physical interpretation. The first important results in flow stability of an incompressible fluid were obtained in the papers of Helmholtz, Rayleigh, and Kelvin [1] in the last century. Heisenberg [2] worked on this problem in the 1920's, and a series of interesting papers by Tollmien [3] appeared subsequently. Apparently one of the first problems in the stability of a compressible fluid was solved by Landau [4]. The first investigations on the boundary-layer stability of an ideal gas were carried out by Lees and Lin [5], and Dunn and Lin [6]. Mention should be made of a series of papers which have appeared quite recently [7–9]. In all the papers mentioned flow stability is investigated in the framework of classical single-phase hydrodynamics. Meanwhile, in recent years, the processes by which perturbations propagate in media with relaxation have been intensively studied [10–12].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 87–93, May–June, 1976.     

Laminar boundary layer on a partly moving surface in the presence of blowing or suction

Planar and axisymmetric flows of a multicomponent compressible gas in a laminar boundary layer with nonzero tangential component of the velocity on a permeable surface are considered. The asymptotic solutions of the boundary-layer equations obtained earlier [1–4] for large values of the blowing and suction parameters are generalized to the case when the velocity vector of the blown or extracted gas makes an acute angle with the surface of the body, this angle depending on the longitudinal coordinate. The region of applicability of the asymptotic formulas is estimated on the basis of the results of numerical solution of the boundary-layer equations. The results are given of some calculations of the boundary layer on a partly moving surface.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 28–36, September–October, 1979.We thank G. A. Tirskii and G. G. Chernyi for a helpful discussion of the results.     

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-layer separation on a cooled body and interaction with a hypersonic flow

In previous papers, e.g., [1, 2], boundary-layer separation was investigated by analyzing solutions of the boundary-layer equations with a given external pressure distribution. In general, this kind of solution cannot be continued after the separation point. Study of the asymptotic behavior of solutions of the Navier-Stokes equations [3–5] shows that, in boundarylayer separation in supersonic flow over a smooth surface, the main effect on the flow in the immediate vicinity of the separation point is a large local pressure gradient induced by interaction with the external flow. The solution can be continued beyond the separation point and linked to the solutions in the other regions, located downstream [5]. Similar results for incompressible flow were recently obtained in [6]. We can assume that in general there is always a small region near the separation point in which separation is self-induced, and where the limiting solution of the Navier-Stokes equations does not contain unattainable singular points. However, this limiting slope picture can be more complex and can contain more regions where the behavior of the functions differed from that found in [3–6]. The present paper investigates separation on a body moving at hypersonic speed, where the ratio of the stagnation temperature to the body temperature is large. It is shown that both. for hypersonic and supersonic speeds the flow near the separation point is appreciably affected by the distribution of parameters over the entire unperturbed boundary layer, and not only in a narrow layer near the body, as was true in the flows studied earlier [3–6]. Regions may appear with appreciable transverse pressure drops within the zone occupied by layers of the unperturbed boundary layer. Similarity parameters are given, the boundary problems are formulated, and the results of computer calculation are presented. The concept of subcritical and supercritical boundary layers is refined, and the dependence of pressure coefficients responsible for separation on the temperature factor is established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 99–109, November–December 1973.     

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.     

Boundary layer on a blunt body in a flow of dusty gas

The problem of interaction of gas-dust flows with solid surfaces arose in connection with the study of the motion of aircraft in a dusty atmosphere [1–2], the motion of a gas suspension in power generators, and in a number of other applications [3]. The presence of a disperse admixture may lead to a significant increase in the heat fluxes [4] and to erosion of the surface [5]. These phenomena are due to the joint influence of several factors — the change in the structure of the carrier-phase boundary layer due to the presence of the particles, collisions of the particles with the surface, roughness of the ablating surface, and so forth. This paper continues an investigation begun earlier [6–7] into the influence of particles on the structure of the dynamical and thermal two-phase boundary layer formed around a blunt body in a flow. The model of the dusty gas [8] has an incompressible carrier phase. The method of matched asymptotic expansions [9] is used to obtain the equations of the two-phase boundary layer. In the frame-work of the refined classification made by Stulov [6], it is shown that the form of the boundary layer equations is different in the presence and absence of inertial precipitation of the particles. The equations are solved numerically in the neighborhood of the stagnation point of the blunt body. The temperature and phase velocity distributions in the boundary layer, and also the friction coefficients and the heat transfer of the carrier phase are found for a wide range of the determining parameters. In the case of an admixture of low-inertia particles that are not precipitated on the body, it is shown that even when the mass concentration of the particles in the undisturbed flow is small their accumulation in the boundary layer can lead to a sharp increase in the thermal fluxes at the stagnation point.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 99–107, September–October, 1985.I thank V. P. Strulov for a discussion.     

Relaxation of anharmonic molecules

The inclusion of anharmonic effects is important in vibrational population inversions in CO-lasers [1, 2], in relaxation processes in jets [3], in thermal dissociation [1], in the kinetics of chemical reactions with high thresholds [4], etc. Usually these effects are studied by including anharmonic corrections to the kinetic constants in the discrete model of single-quantum transitions or in the diffusion approximation [1, 2]. In [5] a method was given of solving the relaxation equations fro arbitrary forms of the rate constants and the spectrum of the molecule. The method is valid when the ratio of the population densities of neighboring levels varies smoothly with quantum number. It was shown in [5] that this approximation can be used to construct analytical solutions for a wide class of problems. In the present paper the method of [5] is extended to the case of equations with variable coefficients. The properties of the solutions for VT-relaxation of anharmonic molecules are analyzed, and the inclusion of sources is considered. A simple method of taking into account multiple-quantum transitions is given, as well as an extension of the method to an arbitrary mixture of gases. The population densities are calculated and the possibility of using our solutions in relaxation gas dynamics is discussed.Translated from Zhurnal Prikladnoi Mekhaniki Tekhnicheskoi Fiziki, No. 3, pp. 22–31, May–June, 1986.     

Computation of base flow characteristics with uniform blowing

Blowing at bluff body base was considered under different conditions and for small amount of blowing this problem was solved using dividing streamline model [1]. The effect of supersonic blowing on the flow characteristics of the external supersonic stream was studied in [2–4]. The procedure and results of the solution to the problem of subsonic blowing of a homogeneous fluid at the base of a body in supersonic flow are discussed in this paper. Analysis of experimental results (see, e.g., [5]) shows that within a certain range of blowing rate the pressure distribution along the viscous region differs very little from the pressure in the free stream ahead of the base section. In this range the flow in the blown subsonic jet and in the mixing zones can be described approximately by slender channel flow. This approximation is used in the computation of nozzle flows with smooth wall inclination [6, 7]. On the other hand, boundary layer equations are used to compute separated stationary flows with developed recirculation regions [8] in order to describe the flow at the throat of the wake. The presence of blowing has significant effect on the flow structure in the base region. An increasing blowing rate reduces the size of the recirculation region [9] and increases base pressure. This leads to a widening of the flow region at the throat, usually described by boundary-layer approximations. At a certain blowing rate the recirculation region completely disappears which makes it possible to use boundary-layer equations to describe the flow in the entire viscous region in the immediate neighborhood of the base section.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 76–81, January–February, 1984.     

Investigation of the flow and heat transfer of viscous reacting fluids in long tubes

Heat transfer and resistance in the case of laminar flow of inert gases and liquids in a circular tube were considered in [1–4], the justification of the use of boundary-layer type equations for investigating two-dimensional flows in tubes being provided in [4]. The flow of strongly viscous, chemically reacting fluids in an infinite tube has been investigated analytically and numerically in the case of a constant pressure gradient or constant flow rate of the fluid [5–8]. An analytic analysis of the flow of viscous reacting fluids in tubes of finite length was made in [9, 10]. However, by virtue of the averaging of the unknown functions over the volume of the tube in these investigations, the allowance for the finite length of the tube reduced to an analysis of the influence of the time the fluid remains in the tube on the thermal regime of the flow, and the details of the flow and the heat transfer in the initial section of the tube were not taken into account. In [11], the development of chemical reactions in displacement reactors were studied under the condition that a Poiseuille velocity profile is realized and the viscosity does not depend on the temperature or the concentration of the reactant; in [12], a study was made of the regimes of an adiabatic reactor of finite length, and in [13] of the flow regimes of reacting fluids in long tubes in the case of a constant flow rate. The aim of the present paper is to analyze analytically and numerically in the two-dimensional formulation the approach to the regimes of thermal and hydrodynamic stabilization in the case of the flow of viscous inert fluids and details of the flow of strongly viscous reacting fluids.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 17–25, January–February, 1930.     

System of equations for describing dynamical problems associated with the mechanics of lung parenchyma

In [1], which is an expanded version of the paper [2], the equations of conservation of mass and momentum are shown to be valid for dynamical problems of lung parenchyma. This system of equations is now closed by means of the heat flux equation. As in [1, 2], the heat flux equation is obtained on the basis of the methods of the mechanics of heterogeneous media [3] and anatomical data on the structure of lung parenchyma.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 21–29, May–June, 1988.     

Character of the boundary layer near its origin on a body rotating in a fluid at rest

The evolution of the boundary layer on bodies of revolution rotating about the symmetry axis in a fluid at rest is largely determined by the position of its origin with respect to the axis of rotation. If the origin of the boundary layer coincides with a pole of the rotating body, then under fairly general assumptions as to the shape of the body the boundary layer has a nonzero thickness in the vicinity of the pole, and the flow in it is described by a particular self-similar solution of the boundary-layer equations [1, 2]. The applicability of existing approximate methods for calculation of the boundary layer [2, 3] is restricted to this case. The results of the present article refer to the case in which the boundary originates at the leading edge at a finite distance from the rotation axis. The behavior of the solution of the boundary-layer equations near the edge is determined. A transformation of variables that reduces the system of boundary-layer equations to a form suitable for analysis and solution is derived. This transformation is used to obtain universal equations determining the local behavior of the boundary layer in the vicinity of its origin in both of the cases indicated above.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 137–140, July–August, 1976.