Homogenization of Wall-Slip Gas Flow Through Porous Media

The permeability of reservoir rocks is most commonly measured with an atmospheric gas. Permeability is greater for a gas than for a liquid. The Klinkenberg equation gives a semi-empirical relation between the liquid and gas permeabilities. In this paper, the wall-slip gas flow problem is homogenized. This problem is described by the steady state, low velocity Navier–Stokes equations for a compressible gas with a small Knudsen number. Darcy's law with a permeability tensor equal to that of liquid flow is shown to be valid to the lowest order. The lowest order wall-slip correction is a local tensorial form of the Klinkenberg equation. The Klinkenberg permeability is a positive tensor. It is in general not symmetric, but may under some conditions, which we specify, be symmetric. Our result reduces to the Klinkenberg equation for constant viscosity gas flow in isotropic media.     

Dynamic model of liquid motion in a porous medium

A theory which takes account of the role of inertial effects in liquid motion in a porous medium is developed. For a compressible liquid, not only the hydrodynamic equations but also the thermodynamic transfer equation are formulated. The initiation, propagation, and dissipation characteristics of eddy motion are considered. Matching conditions at the interface between the media and boundary conditions are obtained. An approximate formulation of the problem is given, isolating in the porous medium a basic flow region in which the classical Darcy law is valid. As an illustration, weakly perturbed liquid flow in a plane channel with an insert of porous material is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 89–95, November–December, 1978.     

Simulation of Single-Phase Multicomponent Flow Problems in Gas Reservoirs by Eulerian–Lagrangian Techniques

Over the past two decades most discussions of the simulation of miscible displacement in porous media were related to incompressible flow problems; recently, however, attention has shifted to compressible problems. The first goal of this paper is the derivation of the governing equations (mathematical models) for a hierarchy of miscible isothermal displacements in porous media, starting from a very general single-phase, multicomponent, compressible flow problem; these models are then compared with previously proposed models. Next, we formulate an extension of the modified method of characteristics with adjusted advection to treat the transport and dispersion of the components of the miscible fluid; the fluid displacement must be coupled in a two-stage operator-splitting procedure with a pressure equation to define the Darcy velocity field required for transport and dispersion, with the outer stage incorporating an implicit solution of the nonlinear parabolic pressure equation and an inner stage for transport and diffussion in which the mass fraction equations are solved sequentially by first applying a globally conservative Eulerian–Lagrangian scheme to solve for transport, followed by a standard implicit procedure for including the diffusive effects. The third objective is a careful investigation of the underlying physics in compressible displacements in porous media through several high resolution numerical experiments. We consider real binary gas mixtures, with realistic thermodynamic correlations, in homogeneous and heterogeneous formations.     

Nonlinear waves in a viscous compressible fluid with relaxation

The propagation and stability of nonlinear waves in a viscous compressible fluid with relaxation that satisfies a Theological equation of state of Oldroyd type are investigated. An equation that describes the structure of the wave perturbations and its evolution is derived subject to the condition of balance of the nonlinear dissipative and relaxation effects, and its solutions of the solitary wave type are analyzed.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 31–35, May–June, 1993.     

A Class of Physically Stable Non-linear Models of Flow Through Anisotropic Porous Media

A linear stability analysis of the single-phase conservation equation in multidimensional porous media is performed, for both weakly compressible and compressible fluids. Non-Newtonian and non-Darcy effects are accounted for using a non-linear Darcy-like form for the superficial velocity, where the mobility tensor is velocity-dependent and proportional to the permeability. It is found that under this hypothesis, flows at an angle with respect to the principal axes of the permeability tensor can be unstable, unless the mobility is a function of the velocity magnitude in terms of the inverse permeability norm. As shown by previous authors, for steady-state incompressible flows this is also the condition ensuring that the governing equation derives from the minimization of a dissipation potential.     

Development of taylor instability in compressible conducting media in a magnetic field

Article [1] is devoted to the investigation of the interface of viscous incompressible conducting media in the presence of a current and with a magnetic field with a small magnetic Reynolds number. Article [2] discusses the stability of a contact discontinuity in compressible media. In this article, in an analysis of the case of long-wave vibrations in the region z<0, the boundary condition is unsatisfactorily met. Therefore, in this part of it the author actually considered a problem with a mass force f=(0, 0, –g sign z) instead of f=(0, 0, –g). The present, article considers the stability of the interface of compressible conducting media in a magnetic field. It is postulated that the magnetic intensity can undergo a discontinuity at this boundary. The article gives the dependence of the maximal increment of the rise in the instability on the determining parameters. An analysis is made of the stability of a contact discontinuity as a function of the angles formed by the wave vector and the intensity of the magnetic field. The stabilizing effect of the walls on the stability is demonstrated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, 15–19, September–October, 1975.The author thanks A. G. Kulikovskii for his aid in stating the problem and for his continuing interest in the work.     

Linear problem of water-ice phase transitions in unsaturated soils

The approach proposed by Menot is developed. The problem of the freezing of an unsaturated soil is formulated for a compressible gas and equal component pressures. An analytic solution is obtained in the linear approximation. The results indicate the strong dependence of the ice saturation on the permeability of the soil and the applied pressure gradient.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 68–73, May–June, 1990.     

Derivation of the Forchheimer Law via Homogenization

In this paper we derive the Forchheimer law via the theory of homogenization. In particular, we study the nonlinear correction to Darcy's law due to inertial effects on the flow of a Newtonian fluid in rigid porous media. A general formula for this correction term is derived directly from the Navier–Stokes equation via homogenization. Unlike other studies based on the same approach that concluded for the nonlinear correction to be cubic in velocity for isotropic media, the present work shows that the nonlinear correction is quadratic. An example is constructed to illustrate our theory. In this example, the analytic solution to the Navier–Stokes equation is obtained and is utilized to show the validity of the quadratic correction. Both incompressible and compressible fluids are considered.     

Influence of sound waves on steady duct flow

The article discusses the problem of determining the secondary steady flow in a plane duct when a sound field is superimposed on an undisturbed compressible laminar flow. It is shown that under certain simplifying conditions the velocity distribution of the secondary flow in the wall region is given by a simple analytical expression. In the rest of the duct the problem is reduced to the solution of a linear fourth-order ordinary differential equation (in complex variables); this problem is solved numerically. The indicated equation is transformed to an Airy equation for large Reynolds numbers Re of the undisturbed flow. The results are presented graphically.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 2, pp. 57–64, March–April, 1976.The author is indebted to V. E. Nakoryakov for valuable comments and interest.     

Constitutive Relations for Finite Elastic-Inelastic Strains

The evolutionary constitutive elastic-inelastic relation with its compatible objective derivative is derived in general form using the kinematics of superposition of small elastic and inelastic strains on finite elastic-inelastic strains. The equation is rendered concrete using the elastic law for a slightly compressible material.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 5, pp. 138–149, September–October, 2005.     

Waves from an underground explosion

The problem of the propagation of a spherical detonation wave in water-saturated soil was solved in [1, 2] by using a model of a liquid porous multicomponent medium with bulk viscosity. Experiments show that soils which are not water saturated are solid porous multicomponent media having a viscosity, nonlinear bulk compression limit diagrams, and irreversible deformations. Taking account of these properties, and using the model in [2], we have solved the problem of the propagation of a spherical detonation wave from an underground explosion. The solution was obtained by computer, using the finite difference method [3]. The basic wave parameters were determined at various distances from the site of the explosion. The values obtained are in good agreement with experiment. Models of soils as viscous media which take account of the dependence of deformations on the rate of loading were proposed in [4–7] also. In [8] a model was proposed corresponding to a liquid multicomponent medium with a variable viscosity.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 34–41, May–June, 1984.     

Variational method for solving problems of the plasticity of compressible media

Variational methods used in the theory of plastic flow are formulated on the assumption of the incompressibility of the deformable medium. In solving problems of the mechanics of soils and friable media and technological problems of the plastic shaping of uncompacted materials it is very important to take account of irreversible volumetric change. Extremum and variational theorems are proved in [1, 2] for rigid-plastic and viscoplastic expanding bodies. A variational equation equivalent to a complete system of differential equations is derived for a compressible plastic body.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 153–155, September–October, 1977.     

Inertia Effects in High-Rate Flow Through Heterogeneous Porous Media

The paper deals with the effects of large scale permeability–heterogeneity on flows at high velocities through porous media. The media is made of a large number of homogeneous blocks where the flow is assumed to be governed by the Forchheimer equation with a constant inertial coefficient. By assuming the validity of the Forchheimer equation at the large scale, an effective inertial coefficient is deduced from numerical simulations. Different media are investigated: serial-layers, parallel-layers and correlated media. The numerical results show that: (i) for the serial-layers, the effective inertial coefficient is independent of the Reynolds number and decreases when the variance and the mean permeability ratio increases; (ii) for the parallel-layers and the correlated media, the effective inertial coefficient is function of the Reynolds number and increases when the variance and the mean permeability ratio increases. Theoretical relationships are proposed for the inertial coefficient as function of the Reynolds number and the characteristics of the media.     

Shock Compression of a Plate on a Wedged–Shaped Target

Problems of compression of a plate on a wedge–shaped target by a strong shock wave and plate acceleration are studied using the equations of dissipationless hydrodynamics of compressible media. The state of an aluminum plate accelerated or compressed by an aluminum impactor with a velocity of 5—15 km/sec is studied numerically. For a compression regime in which a shaped–charge jet forms, critical values of the wedge angle are obtained beginning with which the shaped–charge jet is in the liquid or solid state and does not contain the boiling liquid. For the jetless regime of shock–wave compression, an approximate solution with an attached shock wave is constructed that takes into account the phase composition of the plate material in the rarefaction wave. The constructed solution is compared with the solution of the original problem. The temperature behind the front of the attached shock wave was found to be considerably (severalfold) higher than the temperature behind the front of the compression wave. The fundamental possibility of initiating a thermonuclear reaction is shown for jetless compression of a plate of deuterium ice by a strong shock wave.     

Ventilated entry of a slender profile into a compressible liquid at subsonic velocity

The problem of vertical entry at subsonic velocity into an ideal compressible fluid is solved in the linear formulation for a slender profile with open attached cavity. An integral equation is obtained for the potential of the accelerations. Expressions are given for the calculation of the drag of a thin wedge and also some results of calculations which show that in the limiting case of infinite depth of penetration, which corresponds to stationary flow past the thin wedge with separation of a jet, an analog of the Prandtl—Glauert theorem holds.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 9–17, July–August, 1980.     

A new formulation of immiscible compressible two-phase flow in porous media

A new formulation is proposed to describe immiscible compressible two-phase flow in porous media. The main feature of this formulation is the introduction of a global pressure. The resulting equations are written in a fractional flow formulation and lead to a coupled system which consists of a nonlinear parabolic (the global pressure equation) and a nonlinear diffusion–convection one (the saturation equation) which can be efficiently solved numerically. To cite this article: B. Amaziane, M. Jurak, C. R. Mecanique 336 (2008).     

Exact solutions of linearized equations of convection of a weakly compressible fluid

A mathematical model of fluid convection under microgravity conditions is considered. The equation of state is used in a form that allows considering the fluid as a weakly compressible medium. Based on the previously proposed mathematical model of convection of a weakly compressible fluid, unsteady convective motion in a vertical band, with a heat flux periodic in time set on the solid boundaries of this band, is considered. This model of convection allows one to study the problem with the boundary thermal model oscillating in an antiphase rather than in-phase mode, while the latter was required for the model of microconvection of an isothermally incompressible fluid. Exact solutions for velocity components and temperature are derived, and the trajectories of fluid particles are constructed. For comparison, the trajectories predicted by the classical Oberbeck-Boussinesq model of convection and by the microconvection model are presented.Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 2, pp. 52–63, March–April, 2005.     

Wave Dynamics of Saturated Porous Media and Evolutionary Equations

Nonlinear wave dynamics of an elastically deformed saturated porous media is investigated following the Biot approach. Mathematical models under research are the Biot model and its generalization by consideration of viscous stresses inside liquids. Using two-scales and linear WKB methods, the classical Biot system is transformed to a first-order wave equation. To construct the solution of the other system, an asymptotic modified two-scales method is developed. Initial system of equations is transformed to a nonlinear generalized Korteweg–de Vries–Burgers equation for quick elastic wave. Distinctions of wave propagation in the context of the Biot model and its generalization are shown.     

Nonlinear long wave modulation in symmetric waves of a plane magnetic layer

The special features of the distribution of the magnetic field in the photosphere of the Sun and the experimental discovery of waves which propagate along magnetic tubes in the solar atmosphere have brought about the publication recently of a large number of articles which study the wave-conducting properties of media with a magnetic structure. One of the simplest cases was that of a plane magnetic layer, which was studied in detail in the linear approximation [1–3]. Starting from the dispersion properties of such a structure, [4] indicates the possibility of the existence in it of solitons in the approximation of waves of low amplitude which are long in relation to the layer. The present study has used the method of different-scale expansions to obtain the Schrödinger equation describing the propagation of nonlinear modulations of a symmetric harmonic mode over a plane magnetic layer in an incompressible fluid. A similar equation has been deduced, for example, for waves in water [5–9].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, 164–168, March–April, 1985.The author wishes to thank M. S. Ruderman for formulating the problem and for useful discussions, and V. B. Baranov for his attention to the study.     

Hydrodynamic features of the exploitation of highly nonuniform source-type oil formations

The flow of a weakly compressible fluid in a highly nonuniform formation with block structure, classified as source-type, is considered. An analytic solution of the problem of fluid flow into a well in a bounded circular reservoir is obtained. On the basis of this solution the effect of the fluid offtake rate on the depletion of the reservoir is investigated. It is shown that in highly nonuniform media a number of unsteady effects, which cannot be described by the classical model with steady mass transfer, occur.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 113–120, September–October, 1993.