New concepts in relative permeabilities at high capillary numbers for surfactant flooding

Flooding oil reservoirs with surfactant solutions can increase the amount of oil that can be recovered. Macroscopic modelling of the process requires relative permeabilities to be functions of saturation and capillary number. With only limited experimental data, relative permeabilities have usually been assumed to be linear functions of saturation at high capillary numbers. The experimental data is reviewed, some of which suggest that this assumption is not necessarily correct. The basis for the assumption is therefore reviewed and it is concluded that the linear model corresponds to microscopically segregated flow in the porous medium. Based on new but equally plausible complementary assumptions about the flow pattern, a mixed flow model is derived. These models are then shown to be limiting cases of a droplet model which represents the mixing scale within the porous medium and gives a physical basis for interpolating between the models. The models are based on physical concepts of flow in a porous medium and so the approach described here represents a significant improvement in the understanding of high capillary number flow. This is shown by the fact that fewer parameters are needed to describe experimental data.Notation A total cross-sectional area assigned to capillary bundle - A ⁽ⁱ⁾ physical cross-sectional area of tube i - c ⁽ⁱ⁾ ordered configurational label for droplets in tube i - c configuration label for tube i (order not considered) - D defined by Equation (26) - E(...) expectation value with respect to the trinomial distribution - S _r ⁽⁾ fractional flow of phase - k absolute permeability - k _r relative permeability of phase - k _r ⁰ endpoint relative permeability of phase - L capillary tube length in bundle model - m ⁽ⁱ⁾ number of droplets of phase a occupying tube i - n exponent for phase a in Equation (2) - N number of droplets in bundle model - N _c capillary number - p pressure - p(c') probability of configuration c - Q ⁽ⁱ⁾ total volume flow rate in tube i - S saturation of phase - S flowing saturation of phase - S _r residual saturation of phase - S _r ⁽⁾ saturations when fractional flow of phase is 1 in the case of varying residual saturations for three-phase flow ( ) - t _c residence time for droplet configuration c - v ⁽ⁱ⁾ total fluid velocity in bundle tube i - , phase label - p pressure differential across capillary bundle - ⁽ⁱ⁾ tube conductivity defined by Equation (7) - viscosity of phase - interfacial tension - gradient operator - ... average over tube droplet configurations     

Reservoir transport equations by compositional approach

The objective of this paper is to present an overview of the fundamental equations governing transport phenomena in compressible reservoirs. A general mathematical model is presented for important thermo-mechanical processes operative in a reservoir. Such a formulation includes equations governing multiphase fluid (gas-water-hydrocarbon) flow, energy transport, and reservoir skeleton deformation. The model allows phase changes due to gas solubility. Furthermore, Terzaghi's concept of effective stress and stress-strain relations are incorporated into the general model. The functional relations among various model parameters which cause the nonlinearity of the system of equations are explained within the context of reservoir engineering principles. Simplified equations and appropriate boundary conditions have also been presented for various cases. It has been demonstrated that various well-known equations such as Jacob, Terzaghi, Buckley-Leverett, Richards, solute transport, black-oil, and Biot equations are simplifications of the compositional model.Notation List B reservoir thickness - B formation volume factor of phase - Cⁱ mass of component i dissolved per total volume of solution - C ⁱ mass fraction of component i in phase - C heat capacity of phase at constant volume - C_p heat capacity of phase at constant pressure - D ⁱ hydrodynamic dispersion coefficient of component i in phase - D_MTf thermal liquid diffusivity for fluid f - F = F(x, y, z, t) defines the boundary surface - f_p fractional flow of phase - g gravitational acceleration - H_p enthalpy per unit mass of phase - J_p volumetric flux of phase - k_rf relative permeability to fluid f - k₀ absolute permeability of the medium - M_p ⁱ mass of component i in phase - n porosity - N rate of accretion - P_f pressure in fluid f - p_ca capillary pressure between phases and =p-p - Rⁱ rate of mass transfer of component i from phase to phase - Rⁱ _source source rate of component i within phase - S saturation of phase - s gas solubility - T temperature - t time - U displacement vector - u velocity in the x-direction - v velocity in the y-direction - V volume of phase - V_s velocity of soil solids - W_i body force in coordinate direction i - x horizontal coordinate - z vertical coordinate Greek Letters _p volumetric coefficient of compressibility - _T volumetric coefficient of thermal expansion - _ij Kronecker delta - volumetric strain - _m thermal conductivity of the whole matrix - internal energy per unit mass of phase - gf suction head - density of phase - _ij tensor of total stresses - _ij tensor of effective stresses - volumetric content of phase - f viscosity of fluid f     

Determination of permeability tensors for two-phase flow in homogeneous porous media: Theory

In this paper we continue previous studies of the closure problem for two-phase flow in homogeneous porous media, and we show how the closure problem can be transformed to a pair of Stokes-like boundary-value problems in terms of pressures that have units of length and velocities that have units of length squared. These are essentially geometrical boundary value problems that are used to calculate the four permeability tensors that appear in the volume averaged Stokes' equations. To determine the geometry associated with the closure problem, one needs to solve the physical problem; however, the closure problem can be solved using the same algorithm used to solve the physical problem, thus the entire procedure can be accomplished with a single numerical code.Nomenclature a a vector that maps V onto , m^-1. - A a tensor that maps V onto . - A area of the - interface contained within the macroscopic region, m². - A area of the -phase entrances and exits contained within the macroscopic region, m². - A area of the - interface contained within the averaging volume, m². - A area of the -phase entrances and exits contained within the averaging volume, m². - Bo Bond number (= (=(–)g²/). - Ca capillary number (= v/). - g gravitational acceleration, m/s². - H mean curvature, m^-1. - I unit tensor. - permeability tensor for the -phase, m². - viscous drag tensor that maps V onto V. - ^* dominant permeability tensor that maps onto v , m². - ^* coupling permeability tensor that maps onto v , m². - characteristic length scale for the -phase, m. - l characteristic length scale representing both and , m. - L characteristic length scale for volume averaged quantities, m. - n unit normal vector directed from the -phase toward the -phase. - n unit normal vector representing both n and n . - n unit normal vector representing both n and n . - P pressure in the -phase, N/m². - p superficial average pressure in the -phase, N/m². - p intrinsic average pressure in the -phase, N/m². - p –p , spatial deviation pressure for the -phase, N/m². - r ₀ radius of the averaging volume, m. - r position vector, m. - t time, s. - v fluid velocity in the -phase, m/s. - v superficial average velocity in the -phase, m/s. - v intrinsic average velocity in the -phase, m/s. - v –v , spatial deviation velocity in the -phase, m/s. - V volume of the -phase contained within the averaging volmue, m³. - averaging volume, m³. Greek Symbols V /, volume fraction of the -phase. - viscosity of the -phase, Ns/m². - density of the -phase, kg/m³. - surface tension, N/m. - (v +v ^T ), viscous stress tensor for the -phase, N/m².     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

Stability of stationary traveling waves on the surface of a vertical film of viscous fluid

The stability of stationary traveling waves of the first and second families with respect to infinitesimal perturbations of arbitrary wavelength is subjected to a detailed numerical investigation. The existence of a unique region of stability of the first family is established for wave numbers (₁, ₁) corresponding to the optimal wave regime. There are several regions of stability of the second family ( _k , _k),k=2,3,..., lying close to the local flow rate maxima. In the regions of instability the growth rates of perturbations of the first family are several times greater than for the second family. This difference increases with increase in the Reynolds number. The calculations make it possible to explain a number of experimental observations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 33–41, May–June, 1989.The authors are grateful to V. Ya. Shkadov for his constant interest, and to A. G. Kulikovskii, A. A. Barmin and their seminar participants for useful discussions and suggestions.     

Chaos and integrability in a nonlinear wave equation

We consider the parametrized family of equations _tt ,u- _xx u-au+u ₂ ² u=O,x(0,L), with Dirichlet boundary conditions. This equation has finite-dimensional invariant manifolds of solutions. Studying the reduced equation to a four-dimensional manifold, we prove the existence of transversal homoclinic orbits to periodic solutions and of invariant sets with chaotic dynamics, provided that =2, 3, 4,.... For =1 we prove the existence of infinitely many first integrals pairwise in involution.     

Mixed convection flow in narrow vertical ducts

The mixed convection flow in a vertical duct is analysed under the assumption that , the ratio of the duct width to the length over which the wall is heated, is small. It is assumed that a fully developed Poiseuille flow has already been set up in the duct before heat from the wall causes this to be changed by the action of the buoyancy forces, as measured by a buoyancy parameter . An analytical solution is derived for the case when the Reynolds numberRe, based on the duct width, is of 0 (1). This is extended to the case whenRe is 0 (^–1) by numerical integrations of the governing equations for a range of values of representing both aiding and opposing flows. The limiting cases, || 1 andR=Re of 0 (1), andR and both large, with of 0 (R ^1/3) are considered further. Finally, the free convection limit, large with R of 0 (1), is discussed.

---

Mischkonvektion in engen senkrechten Rohren

Zusammenfassung Mischkonvektion in einem senkrechten Rohr wird unter der Voraussetzung untersucht, daß das Verhältnis der Rohrbreite zur Länge, über welche die Wand beheizt wird, klein ist. Es wird angenommen, daß sich bereits eine voll entwickelte Poiseuille-Strömung in dem Rohr eingestellt hat, bevor Antriebskräfte, gemessen mit dem Auftriebsparameter , aufgrund der Wandbeheizung die Strömung verändern. Es wird eine analytische Lösung für den Fall erhalten, daß die mit der Rohrbreite als charakteristische Länge gebildete Reynolds-ZahlRe konstant ist. Dies wird mittels einer numerischen Integration der wichtigsten Gleichungen auf den FallRe =f (^–1) sowohl für Gleich- als auch für Gegenstrom ausgedehnt. Weiterhin werden die beiden Grenzfälle betrachtet, wenn || 1 undR=Re konstant ist, sowieR und beide groß mit proportionalR ^1/3. Schließlich wird der Grenzfall der freien Konvektion, großes mit konstantem R, diskutiert.

---

Nomenclature g acceleration due to gravity - Gr Grashof number - G modified Grashof number - h duct width - l length of the heated section of the duct wall - p pressure - Pr Prandtl number - Q flow rate through the duct - Q ₀ heat transfer on the wally=0 - Q ₁ heat transfer on the wally=1 - Re Reynolds number - R modified Reynolds number - T temperature of the fluid - T ₀ ambient temperature - T applied temperature difference - u, velocity component in thex-direction - v, velocity component in they-direction - x, co-ordinate measuring distance along the duct - y, co-ordinate measuring distance across the duct - buoyancy parameter - ₀ modified buoyancy parameter, ₀=R ^–1/3 - coefficient of thermal expansion - ratio of duct width to heated length, =h/l - (non-dimensional) temperature - _w applied temperature on the wally=0 - kinematic viscosity - density of the fluid - ₀ shear stress on the wally=0 - ₁ shear stress on the wally=1 - stream function     

Effect of the Nature of the Interaction between a Gas and a Surface on the Relations for the Knudsen Layer in the Case of Strong Evaporation

The effect of the temperature accommodation coefficient _T on the relations at the Knudsen layer edge is investigated for strong evaporation using the moment method. An explicit expression for the dimensionless density as a function of the temperature and the Mach number M is obtained for 0 < _T < 1. For _T = 0 the entire solution is obtained in explicit form. It is shown that for = 0 and a condensation coefficient << 1 the temperature outside the Knudsen layer changes sharply as M varies from 0 to a certain value much less than unity after which the temperature ceases to depend on . For the model of specular reflection of the molecules from the surface the density and the temperature outside the Knudsen layer are found in explicit form as functions of the Mach number.     

The wedge subjected to tractions: a paradox resolved

The classical two-dimensional solution provided by Lévy for the stress distribution in an elastic wedge, loaded by a uniform pressure on one face, becomes infinite when the opening angle 2 of the wedge satisfies the equation tan 2 = 2. Such pathological behavior prompted the investigation in this paper of the stresses and displacements that are induced by tractions of O(r ^–) as r0. The key point is to choose an Airy stress function which generates stresses capable of accommodating unrestricted loading. Fortunately conditions can be derived which pre-determine the form of the necessary Airy stress function. The results show that inhomogeneous boundary conditions can induce stresses of O(r ^–), O(r ^– ln r), or O(r ^– ln² r) as r0, depending on which conditions are satisfied. The stress function used by Williams is sufficient only if the induced stress and displacement behavior is of the power type. The wedge loaded by uniform antisymmetric shear tractions is shown in this paper to exhibit stresses of O(ln r) as r0 for the half-plane or crack geometry. At the critical opening angle 2, uniform antisymmetric normal and symmetric shear tractions also induce the above type of stress singularity. No anticipating such stresses, Lévy used an insufficiently general Airy stress function that led to the observed pathological behavior at 2.     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Solutions for the anisotropic elastic wedge at critical wedge angles

The classical solution for an isotropic elastic wedge loaded by uniform tractions on the sides of the wedge becomes infinite everywhere in the wedge when the wedge angle 2 equals , 2 or 2_* where tan 2_* = 2_*. When the wedge is loaded by a concentrated couple at the wedge apex the solution also becomes infinite at 2 = 2_*. A similar situation occurs when the wedge is anisotropic except that 2_* is governed by a different equation and depends on material properties. Solutions which do not become infinite everywhere in the wedge are available for isotropic elastic wedges. In this paper we present solutions for the anisotropic elastic wedge at critical wedge angles. The main feature of the solutions obtained here is that they are in a real form even though Stroh's complex formalism is employed.     

Loss of stability and bifurcation of auto-oscillatory regimes close to Poiseuille flow

Bifurcation of Poiseuille flow in a flat channel is used as an example to analyze the problem of determining variables that permit study of bifurcation of a main steady flow of a viscous incompressible liquid for parameters close to the values of the coordinates of a point on the curve of neutral stability at which the first Lyapunov exponent d₀ vanishes and there is a changeover from subcritical to supercritical bifurcation. For Poiseuille flow, such a point (R₂,₂, where R₂ is the Reynolds number, and ₂ is the wave number, occurs on the lower branch of the neutral curve. In this paper, it is shown by the Lyapunov-Schmidt method that for < ₂ the stable time-periodic solution that bifurcates into the subcritical region loses stability in the case of slight supercriticality, and a fold singularity is formed in the amplitude surface. The nature of this additional bifurcation is determined by the sign of the second Lyapunov exponent d₁. For its calculation, the value of ₂ is fixed, and the bifurcation that occurs when the Reynolds number is changed is considered. A solution is sought in the form of a convergent series in powers of = ((R – R₀)^1/4, = ±1. The condition of solvability, which serves to determine the coefficient of ⁴, makes it possible to determine the value of d₂. This procedure is entirely general and makes it possible to study bifurcation in the neighborhood of a point of degeneracy on the neutral curve in other hydrodynamic problems too.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 41–48, March–April, 1991.     

Strömungstechnische Gesichtspunkte der Flußbelastung durch Abwasser aus Kläranlagen

Zusammenfassung Für ein im Durchlauf betriebenes System bestehend aus einem Fluß (Vorfluter) und den angeschlossenen Kläranlagen wird eine Methode zur Bestimmung der Vorfluterbelastung durch die eingeleiteten Klärwässer angegeben. Die Methode erfaßt mit Rücksicht auf die Anwendung des Verursacherprinzips im Gewässerschutz die Belastung durch jede Kläranlage für sich, und zwar in Abhängigkeit von der Wasserführung, den Emissionsraten der betreffenden Kläranlage und dem Selbstreinigungsvermögen von den organischen Stoffen aus der betreffenden Kläranlage. Die abhängigen Veränderlichen sind mit der Fließgeschwindigkeit gewichtete Mittelwerte von Schmutzstoffdichten über den Vorfluterquerschnitt. Im Falle konstanter Vorflutertemperatur und zeitunabhängiger Struktur der Klärwässer ergeben sich beispielsweise für die abhängigen Veränderlichen einfache analytische Darstellungen, welche sich als spezielle Formen des -Theorems erweisen. Es wird gezeigt, bei einem unendlich langen Vorfluter mit konstantem Volumenstrom stromabwärts der Klärwassereinleitungen stimmen die erwähnten gewichteten Mittelwerte mit den entsprechenden ungewichteten stromabwärts der Klärwassereinleitungen überein. Die entwickelte Methode kann leicht erweitert werden, um den Sauerstoffschwund im Vorfluter durch jede Kläranlage für sich zu bestimmen.

---

Fluid mechanical aspects of river pollution by effluents from waste treatment plants

The pollution of a river by effluent inflows from waste treatment plants is modeled under steady-state conditions. With respect to modern policies of environmental protection the method describes the river pollution by each plant separately, depending on the flow conditions, the emission rates of the plant and the microbiological decomposition of the biodegradable matter from the plant. Each dependent variable is a weighted cross-sectional mean of a density of organic matter. If the water temperature is constant and the composition of each effluent is independent of time the method gives simple analytic expressions for the dependent variables, which prove to be special versions of the -theorem. It is shown for an infinitely long river of constant volume rate of flow downstream of the effluent inflows: the weighted means mentioned agree with the corresponding nonweighted downstream of the effluent inflows. The present paper can easily the extended to determine the oxygen deficit in the river due to each plant.

---

Bezeichnungen a Anzahl der Kläranlagen - D(t_b) Kennzahl, Einführung in 4.3 - ^eA Emissionsrate der abbaubaren or ganischen Verschmutzung aus der -ten Kläranlage - ^eU Emissionsrate der nichtabbaubaren organischen Verschmutzung aus der -ten Kläranlage - Vorfluterquerschnitt, Einführung in Gl. (4) - F Flächeninhalt von - dF Betrag eines Flächenelements, Einführung in Gl. (6) - J_A Diffusionsstromdichten, Einführung in Gl. (2) bzw. Gl. (3) - L Anzahl der Stromstrecken - M Gesamtmasse der abbaubaren or- ganischen Verschmutzung in den N Teilchen, Einführung in Gl. (17) - N Anzahl der verschmutzten Flußwasserteilchen, welche die -te Nahfeldvermischungszone während des Zeitintervalles t_a t_b für immer verlassen - P(x, t, x, t_c) Teilchendichte, Einführung in Gl. (11) und Gl.(12) - Q Selbstreinigungsvermögen, Einführung in Gl.(26) - t Zeitpunkt, Einführung in Gl.(11) - t, t_b Intervallgrenzen, Einführung in 4.1 - t_c Zeitpunkt, Einführung in Gl.(11) - t Zeitdifferenz, Einführung im Anschluß an Gl.(10) - t^* charakteristische Zeit, Einführung in 4.3 - Strömungsgeschwindigkeit Komponente von ¯b in Richtung der zu Tal weisenden Oberflächennormalen eines Vorfluterquerschnitts, Einführung in Gl. (5) und Gl. (6) - Volumenstrom, Einführung in Gl. (7) - x Ortsvektor - x Ortsvektor eines bestimmten markierten Teilchens zur Zeit t_c, Einführung in Gl.(11) - x längs der Stromachse gemessene Längenkoordinate - x x-Koordinate des Vorfluterquerschnitts durch x - x,x₊₁ x-Koordinaten der Vorfluterquerschnitte, welche die -te Stromstrecke stromaufwärts bzw. stromabwärts begrenzen. Einführung in 4.2. - transformierte Variable, Einführung in Gl.(65) - Zeitvariable - (t_b) Kennzahl, Einführung in 4.3. - Masse der abbaubaren organischen Verschmutzung in dem markierten Teilchen, Einführung in Gl.(14) - , Integrationsvariablen, Einführung in Gl.(38) bzw. Gl.(28) - _A durch die -te Kläranlage bedingte Dichte der abbaubaren organischen Verschmutzung - _U durch die -te Kläranlage bedingte Dichte der nichtabbaubaren organischen Verschmutzung - Mittelwerte von bzw· , Einführung in Gl.(31) bzw. Gl.(8) - _m -Wert zu einem Maximum, Einführung in Gl.(31) - Verhältnis zweier Mittelwerte, Einführung in Gl.(64) - stochastischer Mittelwert einer Zufallsgröße Y - Y Schwankung einer Zufallsgröße Y um den stochastisehen Mittelwert - Mittlung über den Vorfluterquerschnitt Der saubere Vorfluter sei definiert durch Standardwerte für Mindestanforderungen an die Flußwasserqualität. Vorschläge für solche Standardwerte werden in jüngster Zeit unter Berücksichtigung des Umweltschutzes ausführlich diskutiert ([1]; [2], S.- K 13 -).     

Hypersonic flow over a narrow delta plate at large angles of attack

The problem of hypersonic flow over a flat delta plate with a high sweepback anglex at angles of attack close to /2 is solved using a numerical algorithm based on transition to the conical solution. The existence of conical flow at /2 with the velocity vector directed towards the apex of the plate is established. Values ofC _p/sin² and the thickness of the shock layer in the plane of symmetry of the plate are given as functions of the hypersonic similarity parameterk=tan tanx. A comparison of the calculated and experimental data shows that they are in good agreement.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.5, pp. 183–185, September–October, 1992.     

Stabilization of the uniform rotations of a rigid body around a principal axis

Conclusion We sum up the results regarding the stabilization of the investigated motion. The system (1.2) is stabilizable in the linear approximation if the conditions (2.2) on the parameter are not satisfied, and also in the case 4 if a₂₂₃₀>0. The system (1.2) is stabilizable in the nonlinear approximation for ₀=0, ₁0. The system (1.2) is nonstabilizable in the case 4 if a₂₂₃₀<0. In the remaining cases the investigated motion is nonasymptotically stable.Comparing the results regarding stabilizability and controllability, we note that the relation between these properties is, possibly, more complex than for linear systems, since in the cases 2, 3, in which the necessary conditions for the nonlinear system (1.2) are satisfied, while the sufficient ones are not, we do have a nonasymptotic stability. The further investigation of these cases requires the determination of more general sufficient, possibly necessary and sufficient, conditions of controllability of nonlinear systems; this seems to be possible in any case for systems that are linear with respect to control.In conclusion we note that the critical cases of stabilization as well as the problem of the control of the motion of a rigid body by a reactive force have been investigated long ago (we mention [2, 5]) and, as shown by this paper, they have not been definitively solved and continue to present interest for both theory and practice.Institute of Applied Mathematics and Mechanics, Academy of Sciences of Ukraine, Kiev. University of Science, Sebha, Libya. Translated from Prikladnaya Mekhanika, Vol. 28, No. 9, pp. 73–79, September, 1992.     

Calculation of the pressure on the side surface of bodies consisting of a spherical segment joined to an inverted cone in a supersonic stream

The results of investigations of inviscid flow over inverted cones with nose consisting of a spherical segment were published for the first time in Soviet literature in [1–4]. In the present paper, a numerical solution to this problem is obtained using the improved algorithms of [5, 6], which have proved themselves well in problems of exterior flow over surfaces with positive angles of inclination to the oncoming flow. It is shown that the Mach number 2 M , equilibrium and nonequilibrium physicochemical transformations in air (H = 60 km, V = 7.4 km/sec, R₀ = 1 m), and the angle of attack 0 40° influence the investigated pressure distributions. A comparison of the results of the calculations with drainage experiments for M = 6, = 0-25° confirms the extended region of applicability of the developed numerical methods. Also proposed is a simple correlation of the dependence on the Mach number in the range 1.5 M of the shape of the shock wave near a sphere in a stream of ideal gas with adiabatic exponent = 1.4.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 178–183, January–February, 1981.     

Optical co-ordinates in general relativity 2 — the problem of distance

Summary Let denote the congruence of null geodesics associated with a given optical observer inV ₄. We prove that determines a unique collection of vector fieldsM₍₎ ( =1, 2, 3) and ₍₀₎ overV ₄, satisfying a weak version of Killing's conditions.This allows a natural interpretation of these fields as the infinitesimal generators of spatial rotations and temporal translation relative to the given observer. We prove also that the definition of the fieldsM₍₎ and ₍₀₎ is mathematically equivalent to the choice of a distinguished affine parameter f along the curves of, playing the role of a retarded distance from the observer.The relation between f and other possible definitions of distance is discussed.

---

Sommario Sia la congruenza di geodetiche nulle associata ad un osservatore ottico assegnato nello spazio-tempoV ₄. Dimostriamo che determina un'unica collezione di campi vettorialiM₍₎ ( =1, 2, 3) e ₍₀₎ inV ₄ che soddisfano una versione in forma debole delle equazioni di Killing. Ciò suggerisce una naturale interpretazione di questi campi come generatori infinitesimi di rotazioni spaziali e traslazioni temporali relative all'osservatore assegnato. Dimostriamo anche che la definizione dei campiM₍₎, ₍₀₎ è matematicamente equivalente alla scelta di un parametro affine privilegiato f lungo le curve di, che gioca il ruolo di distanza ritardata dall'osservatore. Successivamente si esaminano i legami tra f ed altre possibili definizioni di distanza in grande.

---

---

Work performed in the sphere of activity of: Gruppo Nazionale per la Fisica Matematica del CNR.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

Hydrodynamic models as applied to the investigation of magnetic plasma stability

In the present paper magnetohydrodynamic models are employed to investigate the stability of an inhomogeneous magnetic plasma with respect to perturbations in which the electric field may be regarded as a potential field (rot E 0). A hydrodynamic model, actually an extension of the well-known Chew-Goldberg er-Low model [1], is used to investigate motions transverse to a strong magnetic field in a collisionless plasma. The total viscous stress tensor is given; this includes, together with magnetic viscosity, the so-called inertial viscosity.Ordinary two-fluid hydrodynamics is used in the case of strong collisions=. It is shown that the collisional viscosity leads to flute-type instability in the case when, collisions being neglected, the flute mode is stabilized by a finite Larmor radius. A treatment is also given of the case when epithermal high-frequency oscillations (not leading immediately to anomalous diffusion) cause instability in the low-frequency (drift) oscillations in a manner similar to the collisional electron viscosity, leading to anomalous diffusion.Notation f particle distribution function - E electric field component - H₀ magnetic field - density - V particle velocity - e charge - m, M electron and ion mass - _i, _e ion and electron cyclotron frequencies - viscous stress tensor - P pressure - r_i Larmor radius - P pressure tensor - t time - frequency - T temperature - collision frequency - collision time - j current density - _i, _e ion and electron drift frequencies - k_x, k_y, k_z wave-vector components - n₀ particle density - g acceleration due to gravity. The authors are grateful to A. A. Galeev for valuable discussion.     

Saint-Venant's principle in nonlinear plane elasticity with sufficiently small strains

A homogeneous, isotropic cylinder in an equilibrium state of plane strain, whose cross-section is a rectangle R : [0 < y ₁ < 2L; 0 < y ₂ < h] with h/L 1, is considered. There are no body forces and the long sides are stress free. At y ₁ = 0 and y ₁ = 2L, there are arbitrary loadings, each statically equivalent to a uniformly distributed tensile or compressive stress c. Within the theory of nonlinear elasticity and with the strains and strain gradients assumed to be sufficiently small (but with no such assumptions on the displacement gradients), it is proved that if (,=1,2) represents the Cauchy stress tensor and the Kronecker delta, then |–c₁₁| decays exponentially to zero in R with distance from the nearer end, and the decay constant depends only upon the material but is independent of L.