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.     

Graphite combustion in a chemical-equilibrium boundary layer

The mechanism of graphite decomposition in a dissociated air stream in the presence of heterogeneous reactions on the surface and nonequilibrium evaporation for the case of a frozen boundary layer is studied in [1, 2]. Examples of the calculation of graphite decomposition in an oxygen stream with equilibrium [3] and nonequilibrium [4] chemical reactions in the boundary layer are presented in [3, 4]. The effect of the chemical reactions in the boundary layer on the rate of mass transfer and on the surface temperature with variation of the external flow parameters (p, T_e) remains unexplained. The present paper is devoted to a study of the mechanism of graphite ablation over a wide range of temperatures and pressures with air flow about the body in the case of an equilibrium boundary layer. The effect of the individual components on the heat and mass transfer processes is investigated.In conclusion, the author wishes to thank N. A. Anfimov for this constant interest and valuable advice, and also I. S. Epifanovskii for fruitful discussion of the results.     

Formation of a region of interaction and of a wake with the instantaneous start of a plate in an incompressible liquid

The flow arising in an incompressible liquid if, at the initial moment of time, a plate of finite length starts to move with a constant velocity in its plane, is discussed. For the case of an infinite plate, there is a simple exact solution of the Navier—Stokes equations, obtained by Rayleigh. The case of the motion of a semiinfinite plate has also been discussed by a number of authors. Approximate solutions have been obtained in a number of statements; for the complete unsteadystate equations of the boundary layer the statement was investigated by Stewartson (for example, [1–3]); a numerical solution of the problem by an unsteady-state method is given in [4]. The main stress in the present work is laid on investigation of the region of the interaction between a nonviscous flow and the boundary layer near the end of a plate. In passing, a solution of the problem is obtained for a wake, and a new numerical solution is also given for the boundary layer at the plate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–8, March–April, 1977.     

Allowance for longitudinal diffusion in the neighborhood of a discontinuity of catalytic surface properties

In the investigation of flow near surfaces with discontinuous changes in the catalytic properties the question arises of the applicability of parabolic boundary and viscous shock layer equations in the neighborhood of the discontinuity. In the present paper, three types of problem are solved in which longitudinal diffusion is taken into account. In the first an insertion with different catalytic properties is placed in the neighborhood of the stagnation point, in the second the discontinuity lines of the catalytic properties are perpendicular to the oncoming flow, while in the third they are parallel. On the main surface and on the insertion surface the heterogeneous catalytic reactions are assumed firstorder reactions with various rate constants whose values vary in a wide range. The data of the solution are compared with the solution obtained using the boundary layer approximation and the regions of influence of the longitudinal diffusion are estimated. In [1–4] a problem similar to the second one was solved by the numerical method of [1] and the Wiener-Hopf method for the case of transition from a noncatalytic to a perfectly catalytic surface and the region of applicability of the boundary layer was estimated [5].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 99–105, July–August, 1986.     

Interaction of a turbulent boundary layer with a shock wave at transonic flow velocities

An asymptotic theory of the interaction of a turbulent boundary layer on a plate with a normal shock wave of low intensity has been constructed in various studies [1–4] under the assumption that the averaged velocity of the particles in the boundary layer in front of the interaction region satisfies a logarithmic law. In the present paper a different approach to this problem is proposed based on a power law of the velocity in the undisturbed boundary layer. The obtained results give different estimates for not only the sizes of the characteristic flow regions in the interaction region but also for the shock intensity leading to boundary layer separation.     

Turbulent boundary layer on a catalytic surface in a nonequilibrium dissociating gas

The majority of the studies which consider the flow of a dissociating gas in a turbulent boundary layer are devoted to the investigation of either frozen or equilibrium flows on a flat plate.The frozen turbulent boundary layer has been studied by Dorrance [1], Kutateladze and Leont'ev [2], and Lapin and Sergeev [3]. A study of the effect of catalytic recombination processes at the plate surface on the heat transfer in a frozen turbulent boundary layer was made by Lapin [4].Kosterin and Koshmarov [5], Ginzburg [6], Dorrance [7], and Lapin [8] have studied the turbulent boundary layer on a plate in equilibrium dissociating gas.The calculation of the heat transfer in a turbulent boundary layer on a catalytic plate surface with nonequilibrium dissociation was made by Kulgein [9]. In this study the nonequilibrium nature of the dissociation process was taken into account only in the laminar sublayer, while the flow in the turbulent core was considered frozen. The solution was found numerically using a computer by means of a laborious iteration process.The present paper reports a method for calculating the turbulent boundary layer on a flat catalytic plate with arbitrary dissociation rate. The method, constructed using the assumptions customary for turbulent boundary layer theory, is a successive approximation method. Good convergence of the method is assured by the fact that the effect of the nonequilibrium nature of the dissociation process on the parameter distribution in the boundary layer and, consequently, on the friction and heat transfer may be allowed for merely by finding corrections, usually relatively small, to the distribution of these parameters in the equilibrium or frozen flows. The basis of the study is the two-layer scheme of the turbulent boundary layer. The Prandtl and Schmidt numbers and also their turbulent analogs are taken equal to unity. As the model of the dissociating gas we use the Lighthill model of the ideal dissociating gas [10], extended by Freeman [11] to nonequilibrium flows.     

Investigation of supersonic flow around a blunt body with injection for high Reynolds numbers

The paper discusses the supersonic flow around a blunt smooth body by a stream of viscous gas with subsonic injection from the surface of the body. The effect of various injection cycles on the physical flow characteristics ahead of the body are studied in [1, 2]; the problem is considered in the approximation of a boundary layer. The nonuniform composition of the gas ahead of the body, chemical reactions between the various components, and the effect of radiation are taken into account. For a number of flow cycles, which are of practical importance, it will be of interest to consider higher approximations in powers of [=1/Re, see Eq. (1.1) below] in the shock layer ahead of the body and, in particular, to explain the action of the displacement effect and also the limits of applicability of the boundary-layer approximation assumed in [1, 2]. Extensive literature has been devoted to the asymptotics of the problem of flow around a blunt body of a viscous gas at high Reynolds numbers (see, for example, Van Dyke's book [3]). An investigation of the problem, based on the method of M. I. Vishik and L. A. Lyusternik, is contained in [4–6]. (The advantage of the use of Vishlik and Lyusternik's method in comparison with the method of internal and external expansion is discussed in [4].) The effect of injection on the flow has not been considered in the papers listed. In this paper, approximate solutions are constructed with an error of order and ² which take into account the effect of the injectionf on the flow . The approximate solutions are compiled from a more accurate nonviscous flow (external solution) and boundary-layer corrections. The boundary-layer corrections are constructed on a shock wave and a contact boundary in such a way that the solution would be continuous and quite smooth. For the external solution at the contact boundary, conditions are obtained which take into account the effect of viscosity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 69–77, January–February, 1974.     

Investigation of hypersonic flow of rarefied gas past wings of infinite span

In the framework of the locally self-similar approximation of the Navier-Stokes equations an investigation is made of the flow of homogeneous gas in a hypersonic viscous shock layer, including the transition region through the shock wave, on wings of infinite span with rounded leading edge. The neighborhood of the stagnation line is considered. The boundary conditions, which take into account blowing or suction of gas, are specified on the surface of the body and in the undisturbed flow. A method of numerical solution of the problem proposed by Gershbein and Kolesnikov [1] and generalized to the case of flow past wings at different angles of slip is used. A solution to the problem is found in a wide range of variation of the Reynolds numbers, the blowing (suction) parameter, and the angle of slip. Flow past wings with rounded leading edge at different angles of slip has been investigated earlier only in the framework of the boundary layer equations (see, for example, [2], which gives a brief review of early studies) or a hypersonic viscous shock layer [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 150–154, May–June, 1984.     

A spatial laminar boundary layer on a streamline in conical external flow of a uniform gas in the absence of sources and sinks

Theoretical study of a three-dimensional laminar boundary layer is a complex problem, but it can be substantially simplified in certain particular cases and even reduced to the solution of ordinary differential equations.One such particular case is the flow of a compressible gas on a streamline in conical external flow. The case is of considerable practical importance because the local heat fluxes may take extremal values on such lines.Such flow, except for the conical case, has been examined [1–4], and an approximate method has been given [1] on the basis of integral relationships and a special form for the approximating functions. A numerical solution has been given [2, 3] for such flow around an infinite cylinder. It was assumed in [1–3] that the Prandtl number and the specific heats were constant, and that the dynamic viscosity was proportional to temperature. Heat transfer has been examined [4] near a cylinder exposed to a flow of dissociated air.Here we give results from numerical solution of a system of ordinary differential equations for the flow of a compressible gas in a laminar boundary layer on streamlines in conical external flow, with or without influx or withdrawal of a homogeneous gas. It is assumed that the gas is perfect and that the dynamic viscosity has a power-law temperature dependence.     

Method of mean-mass parameters for a three-dimensional boundary layer in a flow with vorticity

When a gas flows with hypersonic velocity over a slender blunt body, the bow shock induces large entropy gradients and vorticity near the wall in the disturbed flow region (in the high-entropy layer) [1]. The boundary layer on the body develops in an essentially inhomogeneous inviscid flow, so that it is necessary to take into account the difference between the values of the gas parameters on the outer edge of the boundary layer and their values on the wall in the inviscid flow. This vortex interaction is usually accompanied by a growth in the frictional stress and heat flux at the wall [2, 3]. In three-dimensional flows in which the spreading of the gas on the windward sections of the body causes the high-entropy layer to become narrower, the vortex interaction can be expected to be particularly important. The first investigations in this direction [4–6] studied the attachment lines of a three-dimensional boundary layer. The method proposed in the present paper for calculating the heat transfer generalizes the approach realized in [5] for the attachment lines and makes it possible to take into account this effect on the complete surface of a blunt body for three-dimensional laminar, transition, or turbulent flow regime in the boundary layer.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 80–87, January–February, 1981.     

A nonequilibrium viscous shock-wave layer in the vicinity of the critical point with the associated heat exchange taken into account

Results are presented in this article of numerical calculations of a viscous shock layer with the associated heat exchange in the vicinity of the critical point of a spherical blunt body taken into account in the presence of nonequilibrium chemical processes in the shock layer and on the surface of the body about which the flow takes place. A number of papers [1–4], in which specification of the surface temperature of the obstacle was utilized, have been devoted to the numerical investigation of a nonequilibrium viscous shock layer. At the same time the surface temperature of a body varies in actual flight due to heating, and together with this there is catalytic activity of the material, which appreciably complicates the problem and necessitates the simultaneous treatment of the course of processes in the gaseous and solid phases. The use of a separate formulation is difficult in this case, since the formulas for the thermal flux from the gaseous phase are of an estimative nature [5] when a volume nonequilibrium chemical reaction is present for a surface having an arbitrary catalytic activity. Taking account of the associated heat exchange has been done before for a number of problems of boundary-layer theory [6, 7], and in this case it has permitted determining the characteristics which are most important from the practical standpoint under conditions of flight along a specified trajectory, as well as under specified time-independent conditions of flight at altitudes at which the approximation of a viscous shock layer is valid. The effect of catalytic activity is discussed for a number of surface materials, and it is shown that the use of the formulas of boundary layer theory can appreciably distort the behavior of the surface temperature as a function of time for a certain altitude range.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 108–114, May–June, 1979.     

Unsteady three-dimensional viscous shock layer in nonuniform gas flow in the neighborhood of a stagnation point

The combined influence of unsteady effects and free-stream nonuniformity on the variation of the flow structure near the stagnation line and the mechanical and thermal surface loads is investigated within the framework of the thin viscous shock layer model with reference to the example of the motion of blunt bodies at constant velocity through a plane temperature inhomogeneity. The dependence of the friction and heat transfer coefficients on the Reynolds number, the shape of the body and the parameters of the temperature inhomogeneity is analyzed. A number of properties of the flow are established on the basis of numerical solutions obtained over a broad range of variation of the governing parameters. By comparing the solutions obtained in the exact formulation with the calculations made in the quasisteady approximation the region of applicability of the latter is determined. In a number of cases of the motion of a body at supersonic speed in nonuniform media it is necessary to take into account the effect of the nonstationarity of the problem on the flow parameters. In particular, as the results of experiments [1] show, at Strouhal numbers of the order of unity the unsteady effects are important in the problem of the motion of a body through a temperature inhomogeneity. In a number of studies the nonstationary effect associated with supersonic motion in nonuniform media has already been investigated theoretically. In [2] the Euler equations were used, while in [3–5] the equations of a viscous shock layer were used; moreover, whereas in [3–4] the solution was limited to the neighborhood of the stagnation line, in [5] it was obtained for the entire forward surface of a sphere. The effect of free-stream nonuniformity on the structure of the viscous shock layer in steady flow past axisymmetric bodies was studied in [6, 7] and for certain particular cases of three-dimensional flow in [8–11].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 175–180, May–June, 1990.     

Multicomponent chemically-reacting turbulent viscous shock layer on a catalytic surface

Turbulent flows past blunt bodies at high supersonic speeds are mainly investigated within the framework of the boundary layer model. However, even at large Reynolds numbers owing to the strong entropy gradient on the lateral surface it becomes necessary to take boundary layer corrections into account in the higher approximations [1]. The use of viscous shock layer theory makes it possible to obtain fairly accurate results over a broad interval of variation of the Reynolds numbers without organizing iterations with respect to vorticity and displacement thickness. The nonequilibrium nature of both homogeneous and heterogeneous catalytic reactions is taken into account. The results obtained are compared with the experimental data [2, 3]. Previously, in [4, 5] turbulent flow was investigated within the framework of viscous shock layer theory in the case of equilibrium homogeneous reactions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 144–149, March–April, 1989.     

Approximate solution of the problem of the three-dimensional boundary layer in an incompressible fluid

A method of successive approximations is proposed for the solution of the equations of the three-dimensional incompressible boundary layer on bodies of arbitrary shape. A coordinate system connected with the streamlines of the external nonviscous flow is used. It is assumed that the velocity across the external streamlines is small. When the intensity of secondary flow is low the equations describing the boundary layer in an incompressible fluid are reduced to a form analogous to the equations for the boundary layer on axially symmetrical bodies. An approximate analytical solution is obtained for the velocity and for the friction in the form of equations which can be used for any problems of a three-dimensional incompressible boundary layer. The method developed was applied to the problem of the three-dimensional boundary layer at a plate with a cylindrical obstacle in the presence of a slip angle.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

Hypersonic viscous flow around extended bodies

The characteristic feature of flow around an extended body is the interaction of the thickened boundary layer with the external nonviscous flow. This phenomenon becomes more significant at low Reynolds numbers and high Mach numbers. Theoretical investigation of this interaction is difficult because of the presence of shock waves, which are characteristic of hypersonic velocities; the position and curvature of these shock waves depend on the state of the boundary layer developing in conditions of pronounced vorticity of the external flow. With increasing rarefaction of the flow, the problem begins to take on an elliptic character, and this necessitates the use of methods of investigation of more general form than the classical boundary-layer theory.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 164–166, March–April, 1976.The authors thank V. G. Farafonov and V. N. Arkhipov for guidance and assistance in the work.     

Effect of turbulence of external flow on the turbulent boundary layer at a plate

In [1, 2] turbulence of the external flow was taken into account by specifying the turbulent energy at the external boundary of the boundary layer on integrating the energy-balance equation for the turbulence. In [3] a special correction that allowed the turbulence of the external flow to be taken into account was introduced in determining the mixture path. In [4, 5] the turbulent energy calculated from the energy-balance equation of the turbulence was added to the energy induced by turbulence of the external flow, the energy distribution of the induced turbulence being specified using an empirically selected function. In [6, 7] a method of taking into account the effect of turbulence of the external flow on a layer of mixing and a jet was proposed. In the present work, this method is applied to the boundary layer at a plate.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 26–31, May–June, 1977.     

Supersonic motion of bodies in a gas with shock waves

The authors consider the problem of supersonic unsteady flow of an inviscid stream containing shock waves round blunt shaped bodies. Various approaches are possible for solving this problem. The parameters in the shock layer on the axis of symmetry have been determined in [1, 2] by using one-dimensional theory. The authors of [3, 4] studied shock wave diffraction on a moving end plane and wedge, respectively, by the through calculation method. This method for studying flow around a wedge with attached shock was also used in [5]. But that study, unlike [4], used self-similar variables, and so was able to obtain a clearer picture of the interaction. The present study gives results of research into the diffraction of a plane shock wave on a body in supersonic motion with the separation of a bow shock. The solution to the problem was based on the grid characteristic method [6], which has been used successfully to solve steady and unsteady problems [7–10]. However a modification of the method was developed in order to improve the calculation of flows with internal discontinuities; this consisted of adopting the velocity of sound and entropy in place of enthalpy and pressure as the unknown thermodynamic parameters. Numerical calculations have shown how effective this procedure is in solving the present problem. The results are given for flow round bodies with spherical and flat (end plane) ends for various different values of the velocities of the bodies and the shock waves intersected by them. The collision and overtaking interactions are considered, and there is a comparison with the experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 141–147, September–October, 1984.     

Numerical investigation of convective motion in a closed cavity

The investigation of thermal convection in a closed cavity is of considerable interest in connection with the problem of heat transfer. The problem may be solved comparatively simply in the case of small characteristic temperature difference with heating from the side, when equilibrium is not possible and when slow movement is initiated for an arbitrarily small horizontal temperature gradient. In this case the motion may be studied using the small parameter method, based on expanding the velocity, temperature, and pressure in series in powers of the Grashof number—the dimensionless parameter which characterizes the intensity of the convection [1–4]. In the problems considered it has been possible to find only two or three terms of these series. The solutions obtained in this approximation describe only weak nonlinear effects and the region of their applicability is limited, naturally, to small values of the Grashof number (no larger than 10³).With increase of the temperature difference the nature of the motion gradually changes—at the boundaries of the cavity a convective boundary layer is formed, in which the primary temperature and velocity gradients are concentrated; the remaining portion of the liquid forms the flow core. On the basis of an analysis of the equations of motion for the plane case, Batchelor [4] suggested that the core is isothermal and rotates with constant and uniform vorticity. The value of the vorticity in the core must be determined as the eigenvalue of the problem of a closed boundary layer. A closed convective boundary layer in a horizontal cylinder and in a plane vertical stratum was considered in [5, 6] using the Batchelor scheme. The boundary layer parameters and the vorticity in the core were determined with the aid of an integral method. An attempt to solve the boundary layer equations analytically for a horizontal cylinder using the Oseen linearization method was made in [7].However, the results of experiments in which a study was made of the structure of the convective motion of various liquids and gases in closed cavities of different shapes [8–13] definitely contradict the Batchelor hypothesis. The measurements show that the core is not isothermal; on the contrary, there is a constant vertical temperature gradient directed upward in the core. Further, the core is practically motionless. In the core there are found retrograde motions with velocities much smaller than the velocities in the boundary layer.The use of numerical methods may be of assistance in clarifying the laws governing the convective motion in a closed cavity with large temperature differences. In [14] the two-dimensional problem of steady air convection in a square cavity was solved by expansion in orthogonal polynomials. The author was able to progress in the calculation only to a value of the Grashof numberG=10⁴. At these values of the Grashof numberG the formation of the boundary layer and the core has really only started, therefore the author's conclusion on the agreement of the numerical results with the Batchelor hypothesis is not justified. In addition, the bifurcation of the central isotherm (Fig. 3 of [14]), on the basis of which the conclusion was drawn concerning the formation of the isothermal core, is apparently the result of a misunderstanding, since an isotherm of this form obviously contradicts the symmetry of the solution.In [5] the method of finite differences is used to obtain the solution of the problem of strong convection of a gas in a horizontal cylinder whose lateral sides have different temperatures. According to the results of the calculation and in accordance with the experimental data [9], in the cavity there is a practically stationary core. However, since the authors started from the convection equations in the boundary layer approximation they did not obtain any detailed information on the core structure, in particular on the distribution of the temperature in the core.In the following we present the results of a finite difference solution of the complete nonlinear problem of plane convective motion in a square cavity. The vertical boundaries of the cavity are held at constant temperatures; the temperature varies linearly on the horizontal boundaries. The velocity and temperature distributions are obtained for values of the Grashof number in the range 0<G4·10⁵ and for a value of the Prandtl number P=1. The results of the calculation permit following the formation of the closed boundary layer and the very slowly moving core with a constant vertical temperature gradient. The heat flux through the cavity is found as a function of the Grashof number.