Asymptotic behavior of solutions of the Orr-Sommerfeld equations in regions adjacent to two branches of the neutral curve

The asymptotic behavior of the upper and lower branches of the neutral stability curve of a boundary layer found by Lin [1] was determined more accurately by various authors [2–4], who, on the basis of the linearized Navien-Stokes equations, analyzed the higher approximations in the Reynolds number R. In the limit R , neutral perturbations have wavelengths that exceed in order of magnitude the boundary layer thickness. The long-wavelength asymptotic behavior of the Orr-Sommerfeld equation is, in particular, of interest because the characteristic solutions of the linearized equations of free interaction (triple-deck theory) [5–7] are a limiting form of Tollmierr-Schlichting waves in an incompressible fluid with critical layers next to the wall [8–9]. At the same time, the dispersion relation, which is identical to the secular equation of the Orr-Sommerfeld problem, contains an entire spectrum of solutions not considered in the earlier studies [2–4]. The first oscillation mode in the spectrum may be either stable or unstable. In the present paper, solutions are constructed for each of the subregions (including the critical layer) into which the perturbed velocity field in the linear stability problem is divided at large Reynolds numbers. Dispersion relations describing the neighborhood of the upper and lower branches of the neutral curve for the boundary layer are derived. These relations, which contain neutral solutions as a special case, go over asymptotically into each other in the unstable region between the two branches.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 3–11, July–August, 1984.     

On the stability of certain nonparallel flows of a viscous incompressible liquid in a channel

The stability of very simple nonparallel flows of a viscous incompressible liquid in an infinite plane channel described by the exact solutions of the Navier-Stokes equations is studied. Such solutions are realized between two parallel porous plates when the liquid (or gas) is forced in at one wall and drawn out at the same velocity at the other, with a steady flow of liquid along the channel. In this case the transverse velocity component is constant, and the profile of the longitudinal velocity component is independent of the longitudinal con-ordinate x, being an asymmetric function of the transverse coordinate y. A study of the hydrodynamic stability then reduces to the solution of an equation differing from the Orr-Sommerfeld equation by virtue of the presence of additional terms containing the transverse velocity component of the main flow. By numerically solving both this equation and the ordinary Orr-Sommerfeld equation and comparing the corresponding results for various inflowing Reynolds numbers R₀=v₀h/ (v₀ is the inflow velocity, h is the width of the channel), the effect of the nonparallel and asymmetrical nature of such flows on their stability is discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 125–129, July–August, 1970.     

Linear stability in the flow of a liquid film concurrent with a gas flow

The two-phase flow of liquid films are often encountered in practice, but the number of theoretical papers devoted to this problem is limited. The problem of the linear stability of a viscous liquid film subjected to a gas flow has been formulated in [1] and, in somewhat different form, in [2]. The linear stability of plane-parallel motion in films has been studied analytically in [1–8] for some limiting cases. The range of validity of the analytic approaches remains an open question. Therefore, an exact numerical analysis of flow stability over a fairly broad range is required. In the present paper a separate solution of the problem for the gas and the liquid is shown to be possible. The Orr-Sommerfeld equation has been integrated numerically, and the results are compared to the results of analytic calculations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 143–146, January–February, 1976.The author is grateful to É. É. Markovich for directing the work and to V. Ya. Shkadov for his interest in the work and many useful comments.     

Normal interaction of shock waves in the presence of vibrational relaxation

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

Korteweg-de vries-type equations in the theory of surface waves in electrically conducting liquids

Equations are obtained which describe the propagation of long waves of small, but finite amplitude in an ideal weakly conducting liquid and on the basis of these equations the influence of MHD interaction effects on the characteristics of the solitary waves is investigated. The wave equations are derived under less rigorous constraints on the external magnetic field and the MHD interaction parameter than in [1–3]. It is shown that the evolution of the free surface is described by the KdV-Burgers or KdV equations with a dissipative perturbation, and that the propagation velocity of the solitary waves depends on the strength of the external magnetic field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 177–180, November–December, 1989.     

A theory of hydroshock in a two-phase gas-liquid mixture

A hydroshock in a two-phase gas-liquid mixture is usually calculated by an analog of the Zhukovskii formula for a mixture (see, for example, [1, 2]) which establishes the relation of the hydro-shock intensity to the velocity of sound in the mixture. However, as the experimental data [1] show, the hydroshock intensity in the mixture can significantly exceed the calculated values, a fact which is explained by the increase in the propagation velocity of the perturbation wave in comparison to the velocity of sound in the mixture [3, 4]. In the present paper, an equilibrium model of shock transition, similar to [5, 6], is used to calculate the attenuation of the hydroshock as the gas content of the mixture increases in a bubble flow regime. It is shown that owing to the high compressibility of the mixture the effect of the elasticity of the pipe-line walls is small, and the dependence of the propagation velocity of the perturbation wave on the intensity and gas content becomes the main effect. A simple dependence of the hydroshock intensity on the gas content and two similarity parameters is obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 202–204, September–October, 1984.     

Steady two-phase flow in a vertical tube subject to combined free and forced convection

Using the two-velocity, two-temperature model of a continuous medium, the viscousgravitational flow of a mixture of incompressible liquid and solid particles in a vertical round tube is considered. The free-convection equations are written down on the basis of the general equation of motion and the energy equation of a two-phase medium [1, 2]. Using a finite Hankel integral transformation, a solution is constructed for the case of a linear wall-temperature distribution along the tube. The results of some practical calculations of the velocity and temperature fields over the cross section of the tube are presented, together with the dimensionless heat-transfer coefficient expressed as a function of the Rayleigh number and phase concentration. Here it is assumed that the dynamic and thermal-interaction coefficients between the phases correspond to the Stokes mode of flow for each particle, as a result of which the velocity and thermal phase lag is very small [3].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 132–136, July–August, 1975.     

Propagation of acoustic perturbations along a supersonic jet flowing out into the atmosphere

The propagation of acoustic perturbations (specified in the outlet cross section of a particular channel) along a supersonic jet flowing out of the channel is considered; also considered is acoustic emission from the surface of the jet into the atmosphere. The solution of these problems is obtained by a numerical method on the linear approximation. The laws governing the propagation of the perturbations as a function of the perturbation frequency and other determining parameters are investigated; these parameters include the velocity and temperature of the jet, the velocity of the subsonic accompanying flow in the external medium, and the character of the perturbation in the initial cross section of the jet.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 92–99, March–April, 1977.     

Stability of a laminar boundary layer of a power-law no n-newtonian liquid

The stability of a laminar boundary layer of a power-law non-Newtonian fluid is studied. The validity of the Squire theorem on the possibility of reducing the flow stability problem for a power-law fluid relative to three-dimensional disturbances to a problem with two-dimensional disturbances is demonstrated. A numerical method of integrating the generalized Orr-Sommerfeld equation is constructed on the basis of previously proposed [1] transformations. Stability characteristics of the boundary layer on a longitudinally streamlined semiinfinite plate are considered.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 101–106, January–February, 1976.     

Strong recondensation in a one-component and a two-component rarefied gas for an arbitrary value of Knudsen number

As is known, surface phenomena such as evaporation, absorption, and reflection of molecules from the surface of a body depend strongly on its temperature [1–5]. This leads to the establishment of a flow of a substance between two surfaces maintained at different temperatures (recondensation). The phenomenon of recondensation was studied in kinetic theory comparatively long ago. However, up to the present, only the case of small mass flows in a onecomponent gas has been investigated completely [3,4]. Meanwhile it is clear that by the creation of appropriate conditions we can obtain considerable flows of the recondensing substance, so that the mass-transfer rate will be of the order of the molecular thermal velocity. Such a numerical solution of the problem with strong mass flows along the normal to the surface for small Knudsen numbers for a model Boltzmann kinetic equation was obtained in [7]. In this study we numerically solve the problem of strong recondensation between two infinite parallel plates over a wide range of Knudsen numbers for a one-component and a two-component gas, on the basis of the model Boltzmann kinetic equation [6] for a one-component gas and the model Boltzmann kinetic equation for a binary mixture in the form assumed by Hamel [8], for a ratio of the plate temperatures equal to ten. We also investigate the effect of the relative plate motion on the recondensation flow.Moscow. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 130–138, September–October, 1972.     

Self-similar problems of propagation of an extended hydrofracture in a permeable medium

The propagation of an extended hydrofracture in a permeable elastic medium under the influence of an injected viscous fluid is considered within the framework of the model proposed in [1, 2]. It is assumed that the motion of the fluid in the fracture is turbulent. The flow of the fluid in the porous medium is described by the filtration equation. In the quasisteady approximation and for locally one-dimensional leakage [3] new self-similarity solutions of the problem of the hydraulic fracture of a permeable reservoir with an exponential self-similar variable are obtained for plane and axial symmetry. The solution of this two-dimensional evolution problem is reduced to the integration of a one-dimensional integral equation. The asymptotic behavior of the solution near the well and the tip of the fracture is analyzed. The difficulties of using the quasisteady approximation for solving problems of the hydraulic fracture of permeable reservoirs are discussed. Other similarity solutions of the problem of the propagation of plane hydrofractures in the locally one-dimensional leakage approximation were considered in [3, 4] and for leakage constant along the surface of the fracture in [5–7].Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 91–101, March–April, 1992.     

Stability of quasistatic equilibrium state of an inhomogeneous viscous layer on an inclined plane

A study is made of the stability against small perturbations [1] of a slow flow of an incompressible inhomogeneous linearly viscous liquid under the influence of a force of gravity on an unbounded inclined plane. Problems of such kind arise in glaciology when one estimates the stability of snow on mountain slopes or determines the catastrophic movement of a glacier; the results can also be applied to solifluction phenomena [2, 3]. Equations for perturbations of parallel flows of linearly viscous fluids in the case of a continuous variation of the viscosity and density across the flow were derived in [4]. An attempt to solve the hydrodynamic problem with allowance for a perturbation of the viscosity was made in [5]; however, in the equations for the perturbations, simplifications resulted in the neglect of terms that take into account perturbations of the viscosity. In the quasistatic formulation considered here in the case when allowance is made for perturbation of the density and viscosity, the equation for the perturbation amplitudes is an ordinary differential equation with variable coefficients; analytic solution of the equation is very difficult, even for long-wave perturbations. In this connection a study is made of the stability of a laminar model; the viscosity and density are constant within each layer. A similar hydrodynamic problem in the long-wave approximation was considered in [6]. In the present paper an exact solution is constructed in the quasistatic formulation for a two-layer model; the solution shows that the viscosity of the lower layer has an important influence on the stability. Essentially, instability is observed when the lower layer acts as a lubricant.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 20–24, November–December, 1973.     

Interaction between the boundary layer and a supersonic flow when part of the wall is rapidly heated

There have been many theoretical studies of aspects of the unsteady interaction of an exterior inviscid flow with a boundary layer [1–9]. The mathematical flow models obtained in these studies by the method of matched asymptotic expansions describe a wide range of phenomena observed experimentally. These include boundary layer separation near the hinge of a flap, the flow in the neighborhood of the trailing edge of an oscillating airfoil [1–2], and the development and propagation of perturbations in a boundary layer excited by an oscillating wall or some other way [3–5]. The present paper studies the interaction of an unsteady boundary layer with a supersonic flow when a small part of the surface of a body in the flow is rapidly heated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 66–70, January–February, 1984.     

Wave flow regimes of a thin layer of viscous fluid subject to gravity

The studies of Kapitsa initiated the detailed experimental and theoretical study of the flow of a thin layer of viscous liquid (liquid film) over a solid surface [1–2]. Extensive experimental data on this question have now been accumulated. As a rule, the existing theories are based on linearization of the problem and diverge considerably from the experimental results. The present paper is also addressed to the theoretical solution of this problem. The solution method used enables consideration of the wave flow of the liquid as a nonlinear problem and on this basis permits determining all the parameters of the wave regime-amplitude, wavelength, wave propagation speed, frequency.     

Velocity field near a wing—Fuselage combination at an angle of attack in a supersonic stream

To investigate interference between the wing and fuselage at supersonic flight velocities, one can, besides numerical methods based on the exact equations of motion, make effective use of the theory of small perturbations [1]. This is the direction adopted, in particular, in [2–4], in which the problem is solved in the framework of linear theory. In [5], the results obtained in the first approximations are corrected by taking into account the following term in the expansion of the potential function in a series in a small parameter. The present paper considers the velocity field near an arbitrarily profiled wing with supersonic edges and the features due to the presence of the fuselage. A general expression is found for the singular term of the asymptotic expansion of the solution of the linear equation in the neighborhood of the Mach cone with apex at the point of intersection of the leading edge of the wing with the surface of the fuselage. A uniformly exact solution for the linear differential equation for the additional velocity potential is constructed. The position and intensity of the shock wave on the upper surface of the wing are determined. Analytic dependences and quantitiative estimates are obtained for the local downwashes below the wing in the region of the flow where the linear theory leads to the largest errors. The obtained results are important for the correct determination of the aerodynamic characteristics of aircraft in the three-dimensional velocity field produced by the wing-fuselage combination.Translated from Izvestlya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 136–148, November–December, 1980.I am grateful to M. F. Pritulo for discussing the results of the work.     

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.     

Example of parametric resonance in a rotating flow

Parametric resonance is one of the common types of instability of mechanical systems [1]. A standard example of the equations describing parametric oscillations is the Mathieu equation and its generalizations. In hydrodynamics these oscillations have been closely studied in connection with the problem of the vertical oscillations of a vessel containing an incompressible fluid in a uniform gravity field [1–5]. In this paper a new example of a flow whose stability problem reduces to the Mathieu equation is given. This is a flow of special type in a rotating cylindrical channel. The direction of the angular velocity is perpendicular to the channel axis, and its magnitude varies periodically with time. Flows with this geometry are of potential interest in technical applications [6, 7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 175–177, March–April, 1987.     

Qualitative mathematical analysis of the Richards equation

The Richards equation is widely used as a model for the flow of water in unsaturated soils. For modelling one-dimensional flow in a homogeneous soil, this equation can be cast in the form of a specific nonlinear partial differential equation with a time derivative and one spatial derivative. This paper is a survey of recent progress in the pure mathematical analysis of this last equation. The emphasis is on the interpretation of the results of the analysis. These are explained in terms of the qualitative behaviour of the flow of water in an unsaturated soil which is described by the Richards equation.Nomenclature a coefficient in second-order diffusion term of equation - b coefficient in first-order advection term of equation - D soil-moisture diffusivity [L²T^-1] - h pressure head [L] - H quarter-plane domain for Cauchy-Dirichlet problem [L] x [T] - K hydraulic conductivity scalar [LT^–1] - K hydraulic conductivity tensor [LT^–1] - q soil-moisture flux scalar [LT^–1] - q soil-moisture flux vector [LT^–1] - r dummy variable - R rectangle [L] x [T] - s dummy variable - s* representative value of dummy variable - S half-plane domain for Cauchy problem [L] x [T] - t time [T] - u unknown solution of partial differential equation - u₀ initial-value function - v soil-moisture velocity scalar [LT^–1] - v soil-moisture velocity vector [LT^–1]     

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.