Modeling of cementitious materials exposed to isothermal calcium leaching,considering process kinetics and advective water flow. Part 1: Theoretical model

Part 1 of the paper presents a new mathematical model of hydro-chemo-mechanical behavior of cementitious materials exposed to contact with the deionized water. For analyzing the calcium leaching process, an equation describing kinetics of the process is proposed instead of an equilibrium curve commonly used in previous models. This allows taking directly into account the characteristic time of calcium leaching which affects the time evolution of the process, especially that with a relatively fast decrease of the concentration of calcium in pore water. Constitutive relationships describing the transport and strength properties of concrete during chemical degradation are discussed. They are expressed in terms of the leaching degree, an internal variable of the model describing the material decalcification. The governing equations of the model and boundary conditions are expressed in terms of the state variables: gas pressure, capillary pressure, calcium concentration and displacement vector. Numerical solution of the model equations with the finite element method, as well as several examples of application of the model for analysis of various problems, are presented in the part 2 paper.     

Boundary layer on the moving surface of a cylinder

A study is made of the problem of the boundary layer on a cylinder with a moving surface when the cylinder moves with constant velocity in an incompressible fluid. Expressions are obtained for the distributions of the frictional stress on the surface of the cylinder and the coordinate of the singular point in the solution of the boundary layer equations that indicates the appearance of a region of reverse flow for different values of the relative velocity of the motion of the surface of the cylinder. Numerical calculations have been made of the work of the force of friction associated with displacement of the cylinder, the work expended on the motion of its surface, and, in the case of flow separation, the work of the pressure forces (it being assumed here that the pressure and friction on the wall behind the singular point are constant and equal to the pressure and friction at the singular point).     

Effect of residual oil saturation on the flow through a porous medium in the neighborhood of an injection well

Models of the residual oil saturation and models of its effect on the flow in injection wells are proposed. The threshold nature of the dependence of the residual oil saturation on the capillary number determines a change in the flow regimes in the neighborhood of the injection well. The cases of pure, contaminated, and compressible reservoirs are considered. The dependences of the basic problem parameters on the displacement conditions and the state of the reservoir are obtained, together with formulas for the pressure distribution and well injectivity. The topicality of such a simulation for field calculations is demonstrated.     

Upscaling of Two-Phase Immiscible Flows in Communicating Stratified Reservoirs

A semi-analytical method for upscaling two-phase immiscible flows in heterogeneous porous media is described. This method is developed for stratified reservoirs with perfect communication between layers (the case of vertical equilibrium), in a viscous dominant regime, where the effects of capillary forces and gravity may be neglected. The method is discussed on the example of its basic application: waterflooding in petroleum reservoirs. We apply asymptotic analysis to a system of two-dimensional (2D) mass conservation equations for incompressible fluids. For high anisotropy ratios, the pressure gradient in vertical direction may be set zero, which is the only assumption of our derivation. In this way, the 2D Buckley–Leverett problem may be reduced to a one-dimensional problem for a system of quasi-linear hyperbolic equations, of a number equal to the number of layers in the reservoir. They are solved numerically, based on an upstream finite difference algorithm. Self-similarity of the solution makes it possible to compute pseudofractional flow functions depending on the average saturation. The computer partial differential equation solver COMSOL is used for comparison of the complete 2D solutions with averaged 1D simulations. Cases of both discrete and continuous (log-normal) permeability distribution are studied. Generally, saturation profiles of the 1D model are only slightly different from the 2D simulation results. Recovery curves and fractional flow curves fit well. Calculations show that at a favorable mobility ratio (displaced to displacing phase) crossflow increases the recovery, while at an unfavorable mobility ratio, the effect is the opposite. Compared with the classical Hearn method, our method is more general and more precise, since it does not assume universal relative permeabilities and piston-like displacement, and it presumes non-zero exchange between layers. The method generalizes also the study of Yortsos (Transp Porous Media 18:107–129, 1995), taking into account in a more consistent way the interactions between the layers.     

Investigation of the influence of blowing on the flow in a hypersonic viscous shock layer near the stagnation streamline on a blunt body

In the framework of the approximation of local similarity to the Navier-Stokes equations, an investigation is made of the axisymmetric flow of homogeneous gas in a hypersonic shock layer, this including the region of transition through the shock wave. Boundary conditions, which take into account blowing of gas, are specified on the surface of the body and in the undisturbed flow. A numerical solution to the problem is obtained in a wide range of variation of the Reynolds number and the blowing parameter. Expressions are found for the dependences on the blowing parameter usually employed in boundary layer theory of the coefficients of friction and heat transfer on the surface of the body, which are divided by their values obtained for blowing parameter equal to zero. It is shown that these dependences are universal and the same as the dependences obtained from the solution of the equations of a hypersonic viscous shock layer with modified Rankin-Hugoniot relations across the shock wave and from the solution of the boundary layer equations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 199–202, January–February, 1980.     

Theoretical and Numerical Investigation of Gas Lubrication Processes on the Basis of Aerodynamic Equations

The dimensionless parameters of the complete system of Navier-Stokes equations of a compressible gas are estimated with reference to a typical gas bearing. It is found that the three-dimensional compressible boundary layer equations should be used as the determining equations for describing gas lubrication processes. After introducing certain assumptions with respect to the dimensionless parameters in the determining equations, an equation for the pressure, the generalized Reynolds equation, is obtained.Use of the spectral method of analysis makes it possible to transform the generalized Reynolds equation into a system of ordinary differential equations. An analytic solution of the entire boundary value problem is obtained for a journal bearing with fairly small eccentricity. By comparing the numerical results obtained using both the solution of the generalized Reynolds equation and the traditional theory it is possible to estimate the effect of the inertia forces, dissipation processes, and heat transfer.     

The displacement solution of conical shell and its application

In this paper,the displacement solution method of the conical shell is presented.Fromthe differential equations in displacement form of conical shell and by introducing adisplacement function,U(s,θ),the differential equations are changed into an eight-ordersoluble partial differential equation about the displacement function U(s,θ)in which thecoefficients are variable.At the same time,the expressions of the displacement and internalforce components of the shell are also given by the displacement function.As special casesof this paper,the displacement function introduced by V.Z.Vlasov in circular cylindricalshell,the basic equation of the cylindrical shell and that of the circular plate are directlyderived.Under the arbitrary loads and boundary conditions,the general bending problem of theconical shell is reduced to finding the displacement function U(s,θ),and the generalsolution of the governing equation is obtained in generalized hypergeometric function,Forthe axisymmetric bending deformation of the     

New method for solving the one-dimensional problem of the displacement of petroleum by water taking account of capillary forces

The article describes a method for the solution of the one-dimensional problem of the displacement of petroleum by water taking account of capillary forces; on the basis of the method, in the equations of two-phase filtration it is proposed to use a new, soughtfor function, i.e., the fraction of water in the total flow of liquid. As a result, the problem is simplified, the labor consumption of the calculations is reduced, and it becomes possible to make a qualitative evaluation of the behavior of the solutions, including the analytical justification of a solution with a stabilized zone.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 83–91, May–June, 1975.     

Problem of filling of a spherical cavity in a Kelvin-Voigt medium

This paper is concerned with mathematical modeling and solution of the problem of the collapse of a spherical cavity in a viscoelastic medium under the action of constant pressure at infinity. A differential equation of motion for the cavity boundary is constructed and solved numerically. The existence of three modes of motion of the boundary is established, and a map of these modes in the plane of the determining parameters is constructed. Asymptotic forms of the solutions of the problem for all modes are constructed. The problem of cavity collapse with capillary forces taken into account is formulated and solved. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 93–101, September–October, 2008.     

Percolation model of two-phase flow through porous media

A physical model of the process of two-phase flow of immiscible fluids through a porous medium is developed and used to make an analytical calculation of the dependence of the relative phase permeabilities on the saturation of the medium by one of the phases. The theory is compared qualitatively with experiment for a model capillary radius frequency function and quantitatively with numerical calculations made on a computer. In both cases good agreement is obtained. The pressure dependences of the phase permeabilities are analyzed. The question of residual saturation with the wetting fluid after completion of the displacement process is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 88–95, January–February, 1987.     

AN ALGEBRAIC MULTIGRID METHOD FOR COUPLED THERMO-HYDRO-MECHANICAL PROBLEMS

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓ,ａｇｒｅａｔｄｅａｌｏｆａｔｔｅｎｔｉｏｎｈａｓｂｅｅｎｄｉｒｅｃｔｅｄｔｏｗａｒｄｓｆｕｌｌｃｏｕｐｌｅｄｔｈｅｒｍｏ＿ｈｙｄｒｏ＿ｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｉｎｄｅｆｏｒｍｉｎｇｐｏｒｏｕｓｍｅｄｉａｄｕｅｔｏｅｘｔｒａｃｔｉｏｎｏｆｕｎｄｅｒｇｒｏｕｎｄｆｌｕｉｄｓ(ｗａｔｅｒ,ｐｅｔｒｏｌｅｕｍ ,ｎａｔｕｒａｌｇａｓ)ｉｎｒｅｓｅｒｖｏｉｒ,ｅｎｖｉｒｏｎｍｅｎｔａｎｄｃｏｎｓｔｒｕｃｔｉｏｎｅｎｇｉｎｅｅｒｉｎｇ (Ｌｅｗｉｓａｎ…     

Translationally nonequilibrium free jet boundary layer in the wake behind a wedge

Macroscopic equations obtained as a thin-layer version of the 13-moment Grad equations derived from kinetic considerations are used for describing the translationally nonequilibrium monatomic gas flow in a hypersonic free jet boundary layer formed in the wake behind a wedge. This model makes it possible to investigate flows with strong violations of equilibrium with respect to the translational degrees of freedom. A method of constructing the solution of this kinetically justified problem based on the solution of an analogous problem in the Navier-Stokes interpretation is proposed. It is established that for the kinetic variant of the problem considered the gas flow velocity distribution along the separating streamline in a plane orthogonal to the wedge generator coincides with the distribution obtained in solving the Navier-Stokes variant. It is found that taking into account the nonequilibrium nature of the flow with respect to the translational degrees of freedom of the gas particles has no effect on the base pressure and the wake angle.     

Microscale Visual Study of End Effects at Permeability Discontinuities

The physical effect of multiphase fluid distribution and flow at permeability boundaries has not been fully investigated, particularly at the pore scale (1–100 μm), although such behaviour can significantly affect the overall scaled-up reservoir trapping capacity and production performance. In this article, microscale physical models have been used to qualitatively study the pore scale flow events at permeability boundaries, both high to low and vice versa, to gain a better understanding of the role of these boundaries and water saturation on multiphase displacement behaviour at the pore scale. We have used etched glass models of stripes of large and small (a factor of two) pores with circular matrix. Capillary pressure, which is the controlling parameter is itself dependant on pore size and its spatial distribution, the magnitude of the interfacial tensions and the wettability between the fluids and the solid surface of the models. Sometimes, the only way the non-wetting fluid can penetrate the boundary is through a fortuitous leakage, whereby the presence of an initial saturation reduces the controlling capillary pressure. Examples are demonstrated including mechanisms of end-effects and how capillary boundary resistance (due to capillary forces) can be broken down and fluid movement across the boundary can develop. These micromodel experiments show vividly that connate water can assist in these processes, particularly oil trapping and leakage of water across a permeability boundary.     

Godunov–mixed methods for immiscible displacement

The immiscible displacement problem in reservoir engineering can be formulated as a system of partial differential equations which includes an elliptic pressure–velocity equation and a degenerate parabolic saturation equation. We apply a sequential numerical scheme to this problem where time splitting is used to solve the saturation equation. In this procedure one approximates advection by a higher-order Godunov method and diffusion by a mixed finite element method. Numerical results for this scheme applied to gas–oil centrifuge experiments are given.     

Modeling and Numerical Simulations of Immiscible Compressible Two-Phase Flow in Porous Media by the Concept of Global Pressure

A new formulation is presented for the modeling of immiscible compressible two-phase flow in porous media taking into account gravity, capillary effects, and heterogeneity. The formulation is intended for the numerical simulation of multidimensional flows and is fully equivalent to the original equations, contrary to the one introduced in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986). 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. A finite volume method is used to solve the global pressure equation and the saturation equation for the water and gas phase in the context of gas migration through engineered and geological barriers for a deep repository for radioactive waste. Numerical results for the one-dimensional problem are presented. The accuracy of the fully equivalent fractional flow model is demonstrated through comparison with the simplified model already developed in Chavent and Jaffré (Mathematical Models and Finite Elements for Reservoir Simulation, 1986).     

一种全耦合多相流分析的并行计算方法

研究了孔隙介质中热、水和汽流全耦合分析的并行计算方法．模型中采用了考虑毛细压力关系的修正有效应力概念,并考虑了相变和潜热传递．基本变量为位移、毛细压力、汽压和温度．并行程序是在国家高性能计算中心（北京）的曙光１０００Ａ上借助ＰＶＭ（ＰａｒａｌｅｌＶｉｒｔｕａｌＭａｃｈｉｎｅ）软件系统实现的,考题显示出较高的并行加速比和效率     

Refined differential equations of deflections in axial symmetrical bending problems of spherical shell and their singular perturbation solutions

This paper deals with the research of accuracy of differential equations of deflections.The basic idea is as follows.Firstly,considering the boundary effect the meridianmidsurface displacement u=0,thus we derive the deflection differential equations;secondly we accurately prove that by use of the deflection differential equations or theoriginal differential equations the same inner forces solutions are obtained;finally,weaccurately prove that considering the boundary effect the meridian surface displacementu=0 is an exact solution.In this paper we give the singular perturbation solution of thedeflection differential equations.Finally we check the equilibrium condition and prove theinner forces solved by perturbation method and the outer load are fully equilibrated.Itshows that perturbation solution is accurate.On the other hand,it shows again that thedeflection differential equation is an exact equation.The features of the new differential equations are as follows:1.The accuracies of the new differentia     

An Analytical Solution for Capillary Gravity Drainage with Dominant Viscous Forces

The development of the capillary fringe during gravity drainage has a significant influence on saturation and pressure distributions in porous formations (Sarkarfarshi et al. in Int J Greenh Gas Control 23:61–71, 2014). This paper introduces an analytical solution for gravity drainage in an axisymmetric geometry with significant capillary pressure. The drainage process results from the injection of a lighter and less viscous injectant into a porous medium saturated with a heavier and more viscous pore fluid. If the viscous force dominates the capillary and the buoyancy forces, then the flow regime is approximated by differential equations and the admissible solution comprises a front shock wave and a trailing simple wave. In contrast to existing analytical solutions for capillary gravity drainage problems (e.g., Nordbotten and Dahle in 47(2) 2011; Golding et al. in J Fluid Mech 678:248–270 2011), this solution targets the saturation distribution during injection at an earlier point in time. Another contribution of this analytical solution is the incorporation of a completely drained flow regime close to the injection well. The analytical solution demonstrates the strong dependency of the saturation distribution upon relative permeability functions, gas entry capillary pressure, and residual saturation. The analytical results are compared to results from a commercial reservoir engineering software package (\(\hbox {CMG } \hbox {STARS}^{\mathrm{TM}}\)).     

Calculations on the flow of heterogeneous mixtures through porous media

Summary Muskat and Meres¹) have formulated basic differential equations governing the motion of heterogeneous fluid mixtures through porous media. They obtained a solution of these equations by numerical integration for the problem of a column initially filled with liquid saturated with gas, which is closed at one end and kept at a constant low pressure at the other¹).Buckley and Leverett²) obtained an analytical solution for the problem of a column initially filled with liquid which is flooded with a second immiscible liquid. They found a solution in which the saturation is a three-valued function of the coordinate along the column.In our paper in the first place a discussion is given of Buckley and Leverett's solution. It appeared that the true solution which contains a discontinuity may be derived from the three-valued solution by a discussion of the integral relation which represents the total liquid recovery from the column. This discussion bears a formal resemblance to that occurring in the theory of van der Waal's equation of state.The second problem treated in our paper is that of a vertical column initially filled with liquid saturated with gas under a high pressure, which is opened at its lower end. For this problem we succeeded in finding an analytical solution for low values of the pressure gradient. Here again the paradox of a three-valued solution occurred and led to a discontinuity in the saturation.A general argument is given to the effect that three-valued solutions are unavoidable for a theoretical treatment based on Muskat's equations of problems which lead to discontinuities in the saturation.Paper presented at the VIIth Congres of Applied Mechanics, London 1948.     

非牛顿流体两相渗流问题的摄动解

本文用奇异摄动法结合正则摄动法求解了考虑毛管力因素时多孔介质中弱非牛顿流体的两相驱替问题,得到了分流函数和湿相饱和度的渐近解析解。所得结果同数值解和经典的牛顿流体两相渗流结果进行了比较,并着重讨论了非牛顿因素的影响。