Exact Solutions of the Problem of Unsteady Flow of a Viscoplastic Medium in a Circular Pipe

Two one-parameter series of real solutions describing the process of deceleration and acceleration of a viscoplastic medium under the action of a time-varying pressure gradient are obtained. The problem of axisymmetric unsteady viscoplastic flow is reduced to the solution of the Stefan boundary-value problem for the heat conduction equation with a nonlinear condition on the boundary of the quasi-rigid core. By a self-similar change of variables the problem can be reduced to a second-order ordinary differential equation. The solutions of this equation are represented in terms of Bessel and elementary functions. As a result, two one-parameter series of solutions, the first of which describes the acceleration and the second the deceleration of a viscoplastic medium in a pipe under the action of a time-varying pressure gradient are obtained.     

低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy’s law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equations with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be stably numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy’s flow problem, the exact analytical solution for the case of Darcy’s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclusion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equations; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pressure gradient increasing for the one-dimensional problem.     

Straight flow in an elastoviscoplastic cylindrical layer with possible two-sided slip

The solution of the boundary value problem of the theory of large strains dealing with the viscoplastic flow of an elastoviscoplasticmaterial in the gap between two rigid coaxial cylindrical surfaces in the case of motion of the inner surface is obtained under the assumption that the material can slip on both surfaces. The reversible deformation, the development of a viscoplastic flow in the cases of uniformly accelerated and constant-velocity motions of the surface, the flow deceleration in the case of uniformly decelerated motion of the surface until full stop, and the unloading of the medium are considered.     

Exact solutions for nonlinear transient flow model including a quadratic gradient term

The models of the nonlinear radial flow for the infinite and finite reservoirs including a quadratic gradient term were presented. The exact solution was given in real space for flow equation including quadratic gradiet term for both constant-rate and constant pressure production cases in an infinite system by using generalized Weber transform.Analytical solutions for flow equation including quadratic gradient term were also obtained by using the Hankel transform for a finite circular reservoir case. Both closed and constant pressure outer boundary conditions are considered. Moreover, both constant rate and constant pressure inner boundary conditions are considered. The difference between the nonlinear pressure solution and linear pressure solution is analyzed. The difference may be reached about 8% in the long time. The effect of the quadratic gradient term in the large time well test is considered.     

On shear thinning fluid flow induced by continuous mass injection in porous media with variable conductivity

A new formulation is proposed to examine the propagation of the pressure disturbance induced by the injection of a time-variable mass of a weakly compressible shear thinning fluid in a porous domain with generalized geometry (plane, radial, or spherical). Medium heterogeneity along the flow direction is conceptualized as a monotonic power-law permeability variation. The resulting nonlinear differential problem admits a similarity solution in dimensionless form which provides the velocity of the pressure front and describes the pressure field within the domain as a function of geometry, fluid flow behavior index, injection rate, and exponent of the permeability variation. The problem has a closed-form solution for an instantaneous injection, generalizing earlier results for constant permeability. A parameter-dependent upper bound to the permeability increase in the flow direction needs to be imposed for the expression of the front velocity to retain its physical meaning. An example application to the radial injection of a remediation agent in a subsurface environment demonstrates the impact of permeability spatial variations and of their interplay with uncertainties in flow behavior index on model predictions.     

EXACT SOLUTION AND ITS BEHAVIOR CHARACTERISTIC OF NONLINEAR DUAL-POROSITY MODEL

IntroductionIn the computation of petroleum reservoir engineering design,the nonlinear quadraticgradient term is neglected by assuming small pressure gradient or small compressibility.Theassumption of small pressure gradient may cause significant errors i…     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

Helical Viscoplastic Flow in a Gap between Rigid Cylinders

In this paper, the problem of viscoplastic flow of an incompressible non-Newtonian material between two rigid coaxial cylindrical surfaces in the case of helical motion of each of the cylinders is solved. The deformation under an increasing, constant, and decreasing loading is considered. The problem is solved by using the model of large elastoviscoplastic strains. The parameters of the processes under study are calculated for the domains of reversible strain, for the domains of viscoplastic flow, and for the unloading domains.     

Estimate of the residual volume of viscoplastic oil after flooding

The problem of calculating regions of residual viscoplastic oil with a critical pressure gradient is considered. This problem has hitherto been solved mainly for uniform strata (see, for example, [1–3]) but in the present paper, which follows the work of Entov, Pankov, and Pan'ko [4], the problem is considered not only for uniform strata but also for layered strata and vertically nonuniform strata. It is shown that for sufficiently thin strata allowing free flow between individual seams the formulation of the problem is the same for the piston model and the Buckley-Leverett model and reduces to a two-dimensional boundary value problem with a continuous nonlinear effective flow law for the water. An approximate numerical method of calculating regions of residual oil that is suitable for the three listed types of strata is described and illustrated by examples.     

Spatial stability of the incompressible corner flow

The linear spatial stability of the incompressible corner flow under pressure gradient has been studied. A self-similar form has been used for the mean flow, which reduces the related problem to the solution of a two-dimensional problem. The stability problem was formulated using the parabolised stability equations (PSE) and results were obtained for the viscous modes at medium and high frequencies. The related N-factors indicate that the flow is stable at these frequencies, but probably unstable for small frequencies. Furthermore the inviscid mode for each mean flow was obtained and the results indicate that its importance increases considerably with an increase in the adverse pressure gradient. Finally the dependence of the stability characteristics on the extent of the domain is also considered.     

Viscoplastic boundary layer

Boundary layer (BL) solutions around a flat plate for viscoplastic fluids are re-examined, after the precursory study by Oldroyd of low inertia Bingham fluid mechanics. Due consideration is paid to the admissible stress fields far from the obstacle. Normalized Cauchy equations are introduced in the flowing regions. They initially contain inertia, pressure, and viscoplastic terms obtained by adding the yield value and a viscous stress excess in the flowing regions. New similarity solutions for the BL along a flat plate are derived, and the VPBL properties are given, in the limiting case of creeping flows with a dominant yield value. Improvements with respect to Oldroyd solution and relevant aspects are presented and discussed.     

Law of flow of a viscoplastic fluid through a porous medium with allowance for inertial losses

A macroscopic law of flow of a viscoplastic Schwedoff-Bingham fluid through a porous medium is obtained on the basis of percolation theory with allowance for viscous and inertial losses. The asymptotics of the flow law are estimated and expressions for determining the limiting pressure gradient as a function of the microinhomogeneity parameters are given. Satisfactory qualitative agreement between the theoretical and known experimental data is observed. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 68–73, January–February, 1999.     

考虑二次梯度项影响的三孔双渗模型渗流特征

传统的试井模型与物质平衡方程都是不一致的.在非线性偏微分方程中,根据弱可压缩液体的假设,忽略二次梯度项,在试井较长时间内将产生误差.由于二次梯度项的存在,描述多孔介质微可压缩液体流动的压力分布的偏微分方程是非线性的.本文对三重介质渗流特征系统,保留了非线性偏微分方程中的所有项,没有忽略二次梯度项,建立了由基岩系统、裂缝系统及溶洞系统组成的考虑井筒储集和表皮效应对压力影响的三孔双渗模型.采用一阶向前差分,二阶中心差分的方法获得了此类三重介质模型的差分方程,用解非线性方程组的Newton迭代法求得了方程组的数值解.分别讨论了三重介质的参数变化时的无因次压力的变化规律,做出了典型无因次压力曲线图.     

Flow instability in a curved porous channel formed by two concentric cylindrical surfaces

Stability of laminar flow in a curved channel formed by two concentric cylindrical surfaces is investigated. The channel is occupied by a fluid saturated porous medium; the flow in the channel is driven by a constant azimuthal pressure gradient. The momentum equation takes into account two drag terms: the Darcy term that describes friction between the fluid and the porous matrix, and the Brinkman term, which allows imposing the no-slip boundary condition at the channel walls. An analytical solution for the basic flow velocity is obtained. Numerical analysis is carried out using the collocation method to investigate the onset of instability leading to the development of a secondary motion in the form of toroidal vortices. The dependence of the critical Dean number on porosity and the channel width is analyzed.     

Forced Convection Flow of Power-Law Fluids Over a Flat Plate Embedded in a Darcy-Brinkman Porous Medium

The steady forced convection flow of a power-law fluid over a horizontal plate embedded in a saturated Darcy-Brinkman porous medium is considered. The flow is driven by a constant pressure gradient. In addition to the convective inertia, also the “porous Forchheimer inertia” effects are taken into account. The pertinent boundary value problem is investigated analytically, as well as numerically by a finite difference method. It is found that far away from the leading edge, the velocity boundary layer always approaches an asymptotic state with identically vanishing transverse component. This holds for pseudoplastic (0 < n < 1), Newtonian (n = 1), and dilatant (n > 1) fluids as well. The asymptotic solution is given for several particular values of the power-law index n in an exact analytical form. The main flow characteristics of physical and engineering interest are discussed in the paper in some detail.     

Effect of a periodic action on average flow in an inhomogeneous liquid-saturated porous medium

The problem of the average flow of a viscous incompressible fluid saturating a stationary porous incompressible matrix under a periodic action is considered. It is shown that a spatial inhomogeneity of the medium porosity leads to an average fluid flow, quadratically dependent on the action amplitude, in the direction of increase in porosity. In particular, this means that wave action on an oil reservoir could lead to fluid flow on the interfaces from low-porosity,weakly permeable collector regions into high-porosity regions, for example, to flow from blocks to fractures in fractured porous reservoirs, which makes it possible to enhance oil production. It is shown that in the presence of a constant pressure gradient the flow component generated by a periodic action can provide a substantial proportion of the total flow, especially on the boundaries with low-porosity strata or blocks. With increase in amplitude this may significantly exceed the component associated with the constant pressure gradient.     

Cavitation bubble collapse in nonlinear viscous and viscoplastic media

The problem of the collapse of a bubble in nonlinear viscous and viscoplastic media under the influence of pressure at infinity is solved numerically. The pressure at which the radius of the bubble begins to decrease, the limiting radius of the bubble in the viscoplastic case and the critical collapse pressure where there is no plastic component are found. The critical pressure is found to be an order smaller than the corresponding value for a Newtonian viscous fluid. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.2, pp. 181–184, March–April, 1994.     

LES of the adverse-pressure gradient turbulent boundary layer

We describe large-eddy simulations (LES) of the flat-plate turbulent boundary layer in the presence of an adverse pressure gradient. The stretched-vortex subgrid-scale model is used in the domain of the flow coupled to a wall model that explicitly accounts for the presence of a finite pressure gradient. The LES are designed to match recent experiments conducted at the University of Melbourne wind tunnel where a plate section with zero pressure gradient is followed by section with constant adverse pressure gradient. First, LES are described at Reynolds numbers based on the local free-stream velocity and the local momentum thickness in the range 6560–13,900 chosen to match the experimental conditions. This is followed by a discussion of further LES at Reynolds numbers at approximately 10 times and 100 times these values, which are well out of range of present day direct numerical simulation and wall-resolved LES. For the lower Reynolds number runs, mean velocity profiles, one-point turbulent statistics of the velocity fluctuations, skin friction and the Clauser and acceleration parameters along the streamwise, adverse pressure-gradient domain are compared to the experimental measurements. For the full range of LES, the relationship of the skin-friction coefficient, in the form of the ratio of the local free-stream velocity to the local friction velocity, to both Reynolds number and the Clauser parameter is explored. At large Reynolds numbers, a region of collapse is found that is well described by a simple log-like empirical relationship over two orders of magnitude. This is expected to be useful for constant adverse-pressure gradient flows. It is concluded that the present adverse pressure gradient boundary layers are far from an equilibrium state.     

Lower bound for LCF lifetime and its application to safe design of elastic viscoplastic structures

The theory of shakedown is applied to obtain an upper estimation of LCF lifetime of structures. A model of elastic viscoplastic material similar to the Perzyna one with isotropic strain hardening and isotropic damage is adopted. Assumptions: viscoplastic strain rate is proportional to the access of the yield function over zero; the rate of damage evolution is equal to a function of hardening and damage parameters with the coefficient of fluidity, as a factor of proportionality; damage process is coupled with the viscoplastic deformation process; the hardening parameter is equal to accumulated viscoplastic deformation. The yield surfaces form a family of self-similar surfaces with the diameter as the parameter. The shakedown condition of the Melan type is formulated relatively to the initial yield surface. Features of the stress path lead to an equation with min–max problem of the mathematical programming in the left side, which determines a safe value of the virtual residual stress. The equation provides an opportunity to compute the maximal value of the strain hardening parameter possible under the prescribed loading program. This value allows to obtain an upper estimate to safe work time of the structure, which results in a sufficient condition of the structure integrity during the prescribed time period. An example of the developed theory application to resolve various problems arising from designing of structures is considered.