Nonisothermal multiphase flow of brine and gas through saline media

We propose a general formulation for nonisothermal multiphase flow of brine and gas through saline media. The balance equations include mass balance (three species), equilibrium of stresses and energy balance (total internal energy). Salt, water and air mass balance equations are established. The balance of salt allows the establishment of the equation for porosity evolution due to solid skeleton deformation, dissolution/precipitation of salt and migration of brine inclusions. Water and air mass balance equations are also obtained. Two equations are required for water: total water in the medium and water present in solid phase brine inclusions. The mechanical problem is formulated through the equation of stress equilibrium. Finally, the balance of internal energy is established assuming thermal equilibrium between phases. Some general aspects of the constitutive theory are also presented.     

Mechanics Applied to Skeletal Ontogeny and Phylogeny

The musculo-skeletal system serves the mechanical function of creating motion and transmitting loads. It is made up mainly of four components: bone, cartilage, muscle and fibrous connective tissue. These have evolved over millions of years into the complex and diverse shapes of the animal skeleton. The skeleton, however, is not built to a static plan: it can adapt to mechanical forces during growth, it can remodel if the forces change, and it can regenerate if it is damaged. In this paper, the regulation of skeletal construction by mechanical forces is analyzed from both ontogenetic and phylogenetic standpoints. In the first part, models of biomechanical processes that act during skeletal ontogenesis – tissue differentiation and bone remodeling – are presented and, in the second, the evolution of the middle ear is used as an example of biomechanical change in skeletal phylogenesis. Because the constitutive laws for skeletal tissues are relatively well understood, and because the skeleton is preserved in the fossil record, application of mechanics to skeletal evolution seems to present a good opportunity to explore the relationships governing ontogenetic adaptations and phylogenetic change.     

A simple model of the mechanical behavior of ceramic-like materials

This paper presents a macroscopic mechanical theory for ceramic-like materials undergoing isothermal deformations. The proposed model describes an elastic brittle material which is damageable only under tensile loading. The damage lowers the elastic stiffness in traction simulating hence the softening and the fracture (zero stillness) of the material. The basic idea is to consider the continuum as a mixture of two phases—a linear elastic phase and a masonry phase (which shows a linear elastic behavior under compression but cannot hold tractive loads at all). The damage is then related to the volume fraction β of the clastic constituent. The constitutive relations are derived from macroscopic thermodynamics with the volume fraction β and its gradient β taken as state variables.     

Linear waves in saturated porous media

The peculiarity of linear and nonlinear wave propagation in porous media saturated with liquid or gas has been investigated by the methods of multiphase media mechanics. It has been shown that for the analysis and interpretation of experimental data, it is expedient to build models taking into account the nonstationary powers of interaction between the solid and liquid phases and the viscouselastic behaviour of the porous media skeleton. Inertia and inertia-viscous powers principally influence wave attentuation in porous media. Two interphase mechanisms of momentum transfer (two stress tensors — in the solid phase and liquid) lead to two types of waves. Attenuation is determined not only by interphase friction, but also by dissipation resulting from intergrain friction in the solid phase, the influence of which multiplexly exceeds the liquid viscosity influence. The real decrement of attenuation may exceed the sphere restricted by the limiting curves corresponding to the frozen and equilibrium schemes of intergrain deformation. The attenuation of momentum perturbation has been studied. The method of discrete Fourier transform has been used. The analysis of experimental data contained in the literature and their comparison with the results of calculations has been carried out.     

Construction of two-phase equilibria in a non-elliptic hyperelastic material

This work focuses on the construction of equilibrated two-phase antiplane shear deformations of a non-elliptic isotropic and incompressible hyperelastic material. It is shown that this material can sustain metastable, two-phase equilibria which are neither piecewise homogeneous nor axisymmetric, but, rather, involve non-planar interfaces which completely segregate inhomogeneously deformed material in distinct elliptic phases. These results are obtained by studying a constrained boundary value problem involving an interface across which the deformation gradient jumps. The boundary value problem is recast as an integral equation and conditions on the interface sufficient to guarantee the existence of a solution to this equation are obtained. The contraints, which enforce the segregation of material in the two elliptic phases, are then studied. Sufficient conditions for their satisfaction are also secured. These involve additional restrictions on the interface across which the deformation gradient jumps — which, with all restrictions satisfied, constitutes a phase boundary. An uncountably infinite number of such phase boundaries are shown to exist. It is demonstrated that, for each of these, there exists a solution — unique up to an additive constant — for the constrained boundary value problem. As an illustration, approximate solutions which correspond to a particular class of phase boundaries are then constructed. Finally, the kinetics and stability of an arbitrary element within this class of phase boundaries are analyzed in the context of a quasistatic motion.     

Structure of weak shock waves with phase transitions of the first kind

An analytic solution is found to the problem of the structure of weak shock waves with phase transitions of the first kind; it is based on the use of the Mandel ' shtam—Leontovich—D ' yakov concept of second viscosity. The kinetic and thermodynamic characteristics of the problem are investigated, and also the condition of applicability of the second viscosity concept. It is shown that the obtained solution simplifies considerably if the vapor phase is regarded as a perfect gas. The results of calculations for the water—water vapor system are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 55–62, July–August, 1980.     

Nonlinear interaction of longitudinal—Transverse acoustic waves in a circular cylinder

The limiting amplitudes of acoustic oscillations in a cylindrical volume of a heat releasing medium in which one or several modes are unstable in the linear approximation are determined. One of the mechanisms limiting the amplitudes of unstable acoustic modes is the transfer of energy from them to damped modes by nonlinear interaction. The nonlinear interactions of plane acoustic waves in a long channel have been considered by Artamonov and Vorob'ev [1]; in the present paper, the interaction of mixed longitudinal—transverse acoustic modes in a closed cylindrical volume is considered. The equations describing the interaction of two and three longitudinal—transverse modes are derived and investigated in the quadratic approximation by the method of slowly varying amplitudes and phases of the oscillations [2]. The treatment is applicable to a high-temperature gas, for which general stability conditions in the linear approximation have been formulated by Artamonov [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 3–9, September–October, 1982.I should like to express my thanks to K. I. Artamonov (deceased) for suggesting the problem and for scientific supervision and A. P. Vorob'ev for constant interest in the work and helpful advice.     

Continuum thermo-mechanical model for homogeneous liquid–gas flows with internal surface tension effects

This work presents a rational thermodynamic model which consistently incorporates internal surface tension effects in homogeneous and isothermal two-component two-phase flows, without appealing to additional kinematical variables to describe its microstructure. The liquid and gaseous phases are regarded as a pseudo-fluid whose constitutive behavior is obtained from two thermodynamic potentials—the Helmholtz free energy and a pseudo-potential of dissipation. The structure of the flow is assumed to be described solely by the void fraction, which is treated as an internal variable. The constraint associated with it is regarded as a material property, being an effective part of the constitutive equations. Surface tension is taken by the model into account by adding a suitable extra term, with dependence on the void fraction, in the Helmholtz free energy. The capability of the proposed model in properly describing pertinent physical phenomenon is demonstrated by comparing the analytical predictions of sound speeds in the medium with available experimental data.     

Grouping of suspended particles in a standing acoustic wave field

Assuning a snall difference in phase densities averaged equations have been obtained for the motion of a two-phase medium in a standing acoustic wave field. The force causing the relative motion of the phases has been determined. On the basis of the averaged equations obtained the dynamics of the redistribution of particles in the suspension in a plane standing acoustic wave has been studied.Translated from Izvestiya Akadenii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 81–88, September–October, 1984.The authors wish to thank S. A. Regirer for his interest in this work and for useful advice.     

Nonlinear waves and localized structures in two-phase flow through porous media

Three problems of two-phase flow through porous media are considered. In the first two flows in the region of near-critical saturations are investigated. Since under these conditions the active saturation and hence the phase permeability of one of the phases are small, it is important to take into account the delay in phase redistribution — such types of flow as nonlinear waves and localized structures become important. In the third problem it is assumed that the capillary jump is insignificant as compared with the phase pressures. It is shown that in this case localized structures may also occur.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 71–76, January–February, 1988.     

Phase field dependent viscoplastic behaviour of solder alloys

A phase field model is presented to describe the mechanical behaviour of microstructure dependent materials, which is demonstrated by applying it to eutectic tin–lead solder. This solder material is known to be heavily influenced by its continuously evolving microstructure. For the constitutive behaviour of the different phases the elasto-viscoplastic Perzyna model has been used. It is coupled to the phase field model through the constitutive parameters which are taken dependent on the mass fraction field resulting from the solution of the phase field equations.     

Influence of the structure of a porous medium on the residual gas saturation in the case of capillary displacement by a liquid

A model of a porous medium consisting of randomly branching conical pores is used to investigate the quasistatic displacement of gas by a wetting liquid without application of an external pressure. Allowance is made for the circumstance that in the capillary process all the pores have at least one-sided permeability for the liquid phase. An expression is obtained that relates the residual gas saturation to the parameters which characterize the structure of the pores and the wetting properties of the system. Two new characteristics of the pore space are introduced — the branching parameter and the opening angle of the pores — and the influence of these parameters on the residual saturation is investigated. It is shown that for individual classes of natural media the residual gas saturation depends only on the porosity and the contact angle of wetting.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 128–133, September–October, 1981.     

Dynamical properties of heterogeneous surface layers. capillary wave scattering

The dynamics of a heterogeneous liquid surface constituting a two-dimensional disperse system is considered. One of the surface phases (the dispersed phase) forms circular regions of diameter comparable with the characteristic length of the mechanical disturbances within the continuous disperse medium. Inhomogeneous boundary conditions for the Navier-Stokes equations with a discontinuity on the surface phase contact line are formulated. Special attention is paid to the conditions on this line. An approximate method of solving the surface wave diffraction problem and the results for the case of transverse surface wave scattering are described. It is shown that if the wavelength is close to the dimensions of the two-dimensional dispersed particles and their concentration is sufficiently high, the energy of the scattered waves may exceed that dissipated in the vorticity layer. Thus, a new nonclassical mechanism of surfactant action on capillary wave damping is possible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 129–137, January–February, 1991.     

Displacement of oil by a slug of an active additive forced by water through a stratum

The paper gives a solution to the problem of the displacement of oil by a slug for different forms of the sorption isotherm and the distribution function of the additive between the phases and for different values of the initial flooding of the stratum. The process is considered under conditions of reversible sorption and also under conditions of partial retention of the additive by the skeleton of the porous medium. The behavior of slugs in the case of cyclic pumping of a solution of an active additive is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 3, pp. 102–111, May–June, 1982.I thank M. V. Lur'e and M. V. Filinov for suggesting the problem and constant interest in the work, and also V. M. Entov for helpful discussions.     

The constitutive equations of finite strain poroelasticity in the light of a micro-macro approach

After recalling the constitutive equations of finite strain poroelasticity formulated at the macroscopic level, we adopt a microscopic point of view which consists of describing the fluid-saturated porous medium at a space scale on which the fluid and solid phases are geometrically distinct. The constitutive equations of poroelasticity are recovered from the analysis conducted on a representative elementary volume of porous material open to fluid mass exchange. The procedure relies upon the solution of a boundary value problem defined on the solid domain of the representative volume undergoing large elastic strains. The macroscopic potential, computed as the integral of the free energy density over the solid domain, is shown to depend on the macroscopic deformation gradient and the porous space volume as relevant variables. The corresponding stress-type variables obtained through the differentiation of this potential turn out to be the macroscopic Boussinesq stress tensor and the pore pressure. Furthermore, such a procedure makes it possible to establish the necessary and sufficient conditions to ensure the validity of an ‘effective stress’ formulation of the constitutive equations of finite strain poroelasticity. Such conditions are notably satisfied in the important case of an incompressible solid matrix.     

Mechanics of a porous two-phase strongly viscous medium and some geophysical applications

According to modern geophysical ideas [1], the oceanic crust is formed when the lithospheric plates move apart in the axial regions of the rift zones of the Earth due to the penetration of molten matter from the asthenosphere into cracks and pores of the lithosphere and the subsequent solidification of this matter (the lithosphere and asthenosphere are the upper layers of the Earth). In these regions, which have a width of the order of several tens of kilometers, the liquid magmatic material penetrates the crystalline skeleton of the crust, which means that one must regard this region as a saturated fissured and porous deformable medium that permits phase transitions between the skeleton and the liquid filler. Geophysical data indicate that there is a similar two-phase fissured and porous structure of the medium in the asthenosphere (at least in certain regions) and in regions of high tectonic and volcanic activity. Because geological processes take place slowly, it can be assumed that the skeleton of the two-phase medium is deformed as a strongly viscous incompressible Newtonian liquid with a viscosity much greater than the viscosity of the filler.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 53–63, November–December, 1979.We thank V. P. Myasnikov and V. N. Nikolaevskii for helpful discussion of the present work.     

Hydraulic fracture in poro-hydro-elastic media

The prediction of fluid-driven crack propagation in deforming porous media has achieved increasing interest in recent years, in particular with regard to the modeling of hydraulic fracturing, the so-called “fracking”. Here, the challenge is to link at least three modeling ingredients for (i) the behavior of the solid skeleton and fluid bulk phases and their interaction, (ii) the crack propagation on not a priori known paths and (iii) the extra fluid flow within developing cracks. To this end, a macroscopic framework is proposed for a continuum phase field modeling of fracture in porous media that provides a rigorous approach to a diffusive crack modeling based on the introduction of a regularized crack surface. The approach overcomes difficulties associated with the computational realization of sharp crack discontinuities, in particular when it comes to complex crack topologies including branching. It shows that the quasi-static problem of elastically deforming, fluid-saturated porous media at fracture is related to a minimization principle for the evolution problem. The existence of this minimization principle for the coupled problem is advantageous with regard to a new unconstrained stable finite element design, while previous space discretizations of the saddle point principles are constrained by the LBB condition. This proposed formulation includes a generalization of crack driving forces from energetic definitions towards threshold-based criteria in terms of the effective stress related to the solid skeleton of a fluid-saturated porous medium. Furthermore, a Poiseuille-type constitutive continuum modeling of the extra fluid flow in developed cracks is suggested based on a deformation-dependent permeability, that is scaled by a characteristic length.     

Conditional Stochastic Analysis for Multiphase Trsnsport in Randomly Heterogeneous,variably Saturated Media

Groundwater contamination of organics has recently become a problem of growing concern over the resulting health and environmental problems. In general, the multiphase system of nonaqueous phase liquid (NAPL), water and air has to be studied in order to realistically describe the movement of such materials in the subsurface. Numerous models have been developed to study multiphase flow and/or multispecies transport in porous media. However, using models to study the influence of medium heterogeneity on such flow and transport is only a recent event. It has been demonstrated for single-phase flow and transport in saturated and unsaturated media that the study of medium heterogeneity is amenable to stochastic analysis. In this paper, we extend our Eulerian–Lagrangian stochastic theory for single-phase transport to the problem of multiphase–multispecies transport in randomly heterogeneous media under the conditions that the flow is steady-state and the phases are in local chemical equilibrium. We present theoretical expressions to describe the first two conditional moments of the random concentration of any species in any phase. Though they reveal some of the fundamental properties and help gaining insight into the nature of the problem, these expressions cannot be evaluated without either high resolution Monte Carlo simulation or approximation (closure). Therefore, we propose two sets of workable approximations, one being a weak approximation and the other being a linearized pseudo-Fickian approximation. The former yields a nonlinear integro-differential equation for the first conditional moment and the latter yields a linear differential equation. Then the second moments can be computed from explicit expressions from either the weak or pseudo-Fickian approximation.     

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.     

Nonstationary model of immiscible fluid flow through porous media

Unsteady two-phase flow through a microinhomogeneous porous medium is considered. A forest growth model — a percolation model that enables nonequilibrium effects to be taken into account — is proposed for describing the dynamics of the process. In the context of the plane problem expressions are obtained for determining the saturation and the characteristic dimensions of the stagnation zones of trapped phase behind the displacement front.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 73–80, November–December, 1993.