Modeling of retrograde condensation in the process of steady radial flow through a porous medium

The problem of the steady axisymmetric two-phase flow of a multicomponent mixture through a porous medium with phase transitions is considered. It is shown that the system of equations for the two-phase multicomponent flow process, together with the equations of phase equilibrium, reduces to a system of two ordinary differential equations for the pressures in the gas and liquid phases. A family of numerical solutions is found under certain assumptions concerning the pressure dependence of the molar fraction of the liquid phase.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 92–97, November–December, 1994.     

Two-phase three-component flow through a porous medium with variable total flow

In analyzing the processes of the displacement of oil, in which intensive interphase mass transfer takes place, it is normally assumed that the partial volumes of the components as they mix are additive (Amagat's Law) [1, 2]. Then the equations of motion have an integral, which is the total volume flow rate through the porous medium, and the basic problems of frontal displacement, if there are not too many components in the system, permit an exact analytical study to be made [3–5]. If this assumption is rejected, the total flow becomes variable [3, 6, 7]. It appears that the consequences of this as applied to the processes of the displacement of oil by high pressure gases have not previously been considered. The results of such a study, developing the approach outlined in [4], are given below. The initial multicomponent system is simulated by a three-component one which contains oil (the component being displaced), gas (the neutral or main displacing component), and intermediate hydrocarbon fractions or solvent (the active component). It is shown that instead of the triangular phase diagram (TPD) normally used where the partial volumes of the components are additive, in this case it is convenient to use a special spatial phase diagram (SPD) of the apparent volume concentrations of the components to construct the solutions and to interpret them graphically. The method of constructing the SPD and its main properties are explained. A corresponding graphoanalytical technique is developed for constructing the solutions of the basic problems of frontal displacement which correspond to motions with variable total flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 113–120, May–June, 1985.     

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.     

Theory of diffraction of weak shock waves

A numerical solution is considered to the universal nonlinear boundary-value diffraction problem which occurs in various problems of weak interaction [1, 2] in the asymptotic analysis of the flow in a region with large gradients of the parameters near the point of intersection of the incident, diffracted, and reflected waves. The analytical solutions to this type of problem usually approximately satisfy the conditions on the diffracted front, the position of which is not known beforehand, but is found along with the solution. In the present paper, the problem is solved by the numerical method of [3], which reduces the initial boundary-value problem for the system of short-wave equations with an unknown boundary to the solution of a series of boundary-value problems with a fixed boundary. The problem of the diffraction of a weak shock wave on a wedge with a finite apex angle is considered as an application of the solution. The data calculated by the asymptotic theory agree significantly better with the experimental data [5] than the theoretical data of [4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 176–178, November–December, 1984.     

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.     

Flow of a viscous fluid in an annular gap with intensive blowing from the walls

Self-similar solutions are obtained in [1, 2] to the Navier-Stokes equations in gaps with completely porous boundaries and with Reynolds number tending to infinity. Approximate asymptotic solutions are also known for the Navier-Stokes equations for plane and annular gaps in the neighborhood of the line of spreading of the flow [3, 4]. A number of authors [5–8] have discovered and studied the effect of increase in the stability of a laminar flow regime in channels of the type considered and a significant increase in the Reynolds number of the transition from the laminar regime to the turbulent in comparison with the flow in a pipe with impermeable walls. In the present study a numerical solution is given to the system of Navier-Stokes equations for plane and annular gaps with a single porous boundary in the neighborhood of the line of spreading of the flow on a section in which the values of the local Reynolds number definitely do not exceed the critical values [5–8]. Generalized dependences are obtained for the coefficients of friction and heat transfer on the impermeable boundary. A comparison is made between the solutions so obtained and the exact solutions to the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 21–24, January–February, 1987.     

Numerical study of flows of a viscous compressible gas and of an equilibrium gas mixture in a flat nozzle in the presence of strong blowing

The problem of the interaction of a viscous supersonic stream in a flat nozzle with a transverse gas jet of the same composition blown through a slot in one wall of the nozzle is examined. The complete Navier-Stokes equations are used as the initial equations. The statement of the problem in the case of the absence of blowing coincides with [1]. The conditions at the blowing cut are obtained on the assumption that the flow of the blown jet up to the blowing cut is described by one-dimensional equations of ideal gasdynamics. The proposed model of the interaction is generalized to the case of flow of a multicomponent gas mixture in chemical equilibrium. The exact solutions found in [2] are used as the boundary conditions at the entrance to the section of the nozzle under consideration. The results of numerical calculations of the flows of a homogeneous nonreacting gas and of an equilibrium mixture of gases consisting of four components (H₂, H₂O, CO, CO₂) are given for different values of the parameters of the main stream and of the blown jet. In the latter case it is assumed that the effect of thermo- and barodiffusion can be neglected.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 55–63, July–August, 1974.     

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.     

Theory of nonequilibrium inductive high-frequency discharge in a gas flow

The global nonequilibrium flow in the discharge chamber of an induction plasma generator is modeled. The problem for an equilibrium discharge was considered in [2, 3]. Here, on the basis of a numerical solution of the combined system of Navier-Stokes, Maxwell, energy, ionization kinetics and electron-gas energy balance equations, the structure of the nonequilibrium discharge is analyzed and the results obtained within the framework of the local one-dimensional approach [1] and on the basis of global numerical modeling of the flow are compared. As distinct from [2, 3], in finding the electromagnetic field distribution in the discharge chamber the boundary-value problem for the two-dimensional Maxwell equations is solved.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 153–160, March–April, 1991.The author is grateful to V. V. Lunev and G. N. Zalogin for their constant interest and useful discussions.     

The uniqueness of the solution of boundary-value problems of a thermal model of discharge

The paper is devoted to a nonlinear analysis of superheating [1, 2] instability of an electric discharge stabilized by electrodes [3] in the framework of a thermal model [4] where the stability of the discharge relative to the long-wave and short-wave perturbations is proved in a linear approximation. Similar boundary-value problems arise in the theories of chemically and biologically reacting mixtures [5–7], thermal breakdown of dielectrics [8], thermal explosion [9], in the investigation of nonlinear waves in semiconductors and superconductors [10, 11], and in the investigation of Couette flow with variable viscosity [12]. The uniqueness of the one-dimensional steady solutions of the thermal model of discharge and the stability relative to the small spatial perturbations, respectively, for the exponential and step dependence of the electrical conductivity on the temperature are proved in [3, 13]. The uniqueness of the solutions in the one-dimensional case for the same electrode temperature and arbitrary dependences of the electrical and thermal conductivity on the temperature is established in paper [14]. In the present paper, the existence and uniqueness of steady solutions of the thermal model of discharge in a three-dimensional formulation for arbitrary fairly smooth electrical and thermal conductivity functions of the temperature in the case of isothermal isopotential electrodes are proved analytically.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 140–145, January–February, 1986.The author expresses his gratitude to A. G. Kulikovskii and A. A. Barmin for the formulation of the problem and their discussions.     

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.     

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.     

Flow processes with change of regime in fractured reservoirs

Experimental and industrial observations indicate a strong nonlinear dependence of the parameters of the flow processes in a fractured reservoir on its state of stress. Two problems with change of boundary condition at the well — pressure recovery and transition from constant flow to fixed bottom pressure — are analyzed for such a reservoir. The latter problem may be formulated, for example, so as not to permit closure of the fractures in the bottom zone. For comparison, the cases of linear [1] and nonlinear [2] fractured porous media and a fractured medium [3] are considered, and solutions are obtained in a unified manner using the integral method described in [1]. Nonlinear elastic flow regimes were previously considered in [3–6], where the pressure recovery process was investigated in the linearized formulation. Problems involving a change of well operating regime were examined for a porous reservoir in [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 67–73, May–June, 1991.     

Region of applicability and the main features of the generalized Chapman-Enskog method

The present paper is made necessary by the publication of the foregoing paper in this issue by Kolesnichenko [1]. It considers the basic propositions of the generalized Chapman-Enskog method and analyzes the arguments put forward by Kolesnichenko [1] and the validity of the method. The position of the results obtained by Kolesnichenko [14–17] is indicated. Nonequilibrium flows of multiatomic gases in which there occur processes of exchange of internal energy of the molecules in collisions between them and chemical reactions (such processes are called inelastic) are encountered frequently in nature and technology. It is therefore naturally of interest to derive gas-dynamic equations for such flows. The methods of the kinetic theory of gases were first used to obtain equations describing the limiting cases of very fast inelastic processes that take place in times of the order of the molecule-molecule collision times (equilibrium case) and very slow inelastic processes that take place over times of the order of the characteristic flow time (relaxation case). In [2–5], an algorithm was proposed for deriving gas-dynamic equations valid for arbitrary ratios of the rates of the elastic and inelastic processes and reducing to the well-known equations for the limiting cases already mentioned. The algorithm is called the generalized Chapman-Enskog method (abbreviated to the generalized method). The development, modification, and analysis of its properties can be found in [4, 6–13]. In [1], Kolesnichenko has questioned the validity of this algorithm.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 126–136, May–June, 1984.We thank V. A. Rykov for helpful and constructive discussions of the work.     

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.     

Studies of quality of geometrically nonlinear elasticity theory for small strains and arbitrary displacements

Earlier it was shown in [1, 2] that the equations of classical nonlinear elasticity constructed for the case of small strains and arbitrary displacements are ill posed, because their use in specific problems may result in the appearance of “spurious” bifurcation points. A detailed analysis of these equations and the construction, in their stead, of consistent equations of geometrically nonlinear theory of elasticity can be found in [3]. Certain steps in this direction were also made in [4, 5]. In [3], it was also stated that the methods and applied program packages (APPs) based on the use of the classical relations of nonlinear elasticity require some revision and correction. In the present paper, this conclusion is justified and confirmed by numerical finite-element solutions of several three-dimensional geometrically nonlinear deformation problems and linearized problems on the stability of equilibrium of a rectilinear beam. These solutions were obtained by using two APPs developed by the authors and the well-known APP “ANSYS.” It is shown that the classical equations of the geometrically nonlinear theory of elasticity, which underly the first of the developed APP and the well-known APP “ANSYS,” often lead to overestimated buckling loads for structural members as compared with the consistent equations proposed in [1–3].     

Wave motions in viscoelastic tubes. Forced oscillations

A characteristic of small blood and lymphatic vessels is the capacity of the wall to change its rheological properties and lumen by active contraction of the annular muscle cells contained in it [1–3]. A model of flow in the vessels taking this feature into account has been proposed in [4, 5], where a linear stability analysis is also given. A consequence of wall activity is the existence of auto-oscillatory flow conditions [6–8], which have also been discovered in the numerical solutions of the corresponding problems [9, 10]. Up to the present time flows have only been studied under steady conditions at the ends of the vessel and in the environment. The wall of an actual blood vessel is subject to various actions, frequently of a periodic nature: pressure pulsations at entry and rhythmically changing external forces applied from the surrounding tissues. Data exist on the sensitivity of vessels to transient actions [11–13], in particular on the relationship of their hydraulic resistance to frequency and amplitude of the action. There has been frequent discussion of the hypothesis that bv contraction of muscles in its walls or by external compression the vessel can act as a valveless pump [14, 15]. Within the framework of the quasione-dimensional approximation given below [4] the movement of liquid along a viscoelastic tube in the presence of small amplitude periodic external actions has been studied. A general solution of the problem has been constructed and concrete examples are given illustrating the features of forced wave motions in a tube having passive and active properties.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 4, pp. 94–99, July–August, 1984.     

Asymptotic laws of degeneracy of the shape of thin cavities

Differential equations are obtained in the approximation of a slender body for the shape of thin nonstationary axisymmetric cavities. In contrast to the equations obtained earlier by Grigoryan [1] and the author [2], the equations derived in the present paper contain the well-known asymptotic behavior for a flow of Kirchhoff type in the stationary case. For the case of compressible fluid, the dependence found for the asymptotic behavior of the cavity on the Mach number differs from the result obtained by Gurevich [3]. When compressibility is taken into account, the cavity area grows less rapidly. It is shown that the solutions of the obtained differential equations satisfy the well-known separation conditions [3, 4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 3–10, May–June, 1981.I thank L. I. Sedov, S. S. Grigoryan, V. P. Karlikov, V. A. Eroshin and others who participated in a discussion of the work, and also L. A. Barmina for assisting in the preparation of the paper.     

Nonisothermal viscoplastic fluid flow through porous media

The Alishaev model [1] is extended to the case of nonisothermal flow. Neglecting conductive heat transfer, it is shown that for the model in question in the plane of the complex potential not only are the problems linear but the decoupling of the thermal and hydrodynamic problems is also allowed. The latter is reduced to a mixed problem for an analytic function. This makes it possible to use the wellknown methods and results of the theory of limiting equilibrium pillars for isothermal flow [2–5]. It is also established that the solutions of the unsteady problems tend asymptotically to the solutions of the corresponding steady-state problems and can be obtained from the latter by simpler conversion. The effectiveness of the approach proposed is illustrated with reference to the problem of a source-sink system [1–4].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 117–122, July–August, 1990.     

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.