EXPRESSIONS FOR PRESSURE-VELOCITY-GRADIENT CORRELATIONS

The term for pressure-velocity-gradient correlation was initiated by Rott’srewriting the correlation between the pressure fluctuation gradient and velocity fluctuation.However,it is very difficult to consider the effect of this term.Since Rotta’s work,Launderet al.has made some estimates of this term.In this paper according to the equations forvelocity fluctuation,the pressure fluctuation is solved so that the average value of theproduct of the pressure fluctuation and the velocity fluctuation gradient is obtained.Thus,the whole expressions for the pressure-velocity-gradient correlation are derived.The resultexplains that the limited expressions by Rotta and Launder are reasonable to a certaindegree.The whole expressions in this paper are discussed respectively in two situations:oneis without a separate consideration of large and small vortexes;the other is with a separateconsideration of three kinds of vortexes.Therefore,the paper gives the whole expressions forpressure-velocity-gradient correlation     

Peristaltic flow of Johnson-Segalman fluid in asymmetric channel with convective boundary conditions

This work is concerned with the peristaltic transport of the Johnson-Segalman fluid in an asymmetric channel with convective boundary conditions. The mathematical modeling is based upon the conservation laws of mass, linear momentum, and energy. The resulting equations are solved after long wavelength and low Reynolds number are used. The results for the axial pressure gradient, velocity, and temperature profiles are obtained for small Weissenberg number. The expressions of the pressure gra-dient, velocity, and temperature are analyzed for various embedded parameters. Pumping and trapping phenomena are also explored.     

Transient electro-osmotic and pressure driven flows of two-layer fluids through a slit microchannel

By method of the Laplace transform, this article presents semi-analytical solutions for transient electroosmotic and pressure-driven flows (EOF/PDF) of two-layer fluids between microparallel plates. The linearized Poisson-Boltzmann equation and the Cauchy momentum equation have been solved in this article. At the interface, the Maxwell stress is included as the boundary condition. By numerical computations of the inverse Laplace transform, the effects of dielectric constant ratio ε , density ratio ρ , pressure ratio p, viscosity ratio μ of layer II to layer I, interface zeta potential difference △ψ, interface charge density jump Q, the ratios of maximum electro-osmotic velocity to pressure velocity α , and the normalized pressure gradient B on transient velocity amplitude are presented.We find the velocity amplitude becomes large with the interface zeta potential difference and becomes small with the increase of the viscosity. The velocity will be large with the increases of dielectric constant ratio; the density ratio almost does not influence the EOF velocity. Larger interface charge density jump leads to a strong jump of velocity at the interface. Additionally, the effects of the thickness of fluid layers (h1 and h2 ) and pressure gradient on the velocity are also investigated.     

Lattice Boltzmann simulation of flows in bifurcate channel at rotating inflow boundary conditions and resulted different outflow fluxes

The Lattice Boltzmann method(LBM) is used to simulate the flow field in a bifurcate channel which is a simplified model of the draft tube of hydraulic turbine machine.According to the simulation results,some qualitative conclusions can be deduced.The reason of uneven flux in different branches of draft tube is given.Not only the vortex rope itself,but also the attenuation of the rotation strength is important in bringing on the uneven flux.The later leads to adverse pressure gradient,and changes the velocity profile.If the outlet contains more than one exit,the one that contains the vortex rope will lose flux because of this adverse pressure gradient.Several possible methods can be used to minimize the adverse pressure gradient domain in order to improve the efficiency of turbine machine.     

PRESSURE AND PRESSURE GRADIENT IN AN AXISYMMETRIC RIGID VESSEL WITH STENOSIS

Based on an improvement of the Karman-Pohlhausen's method, using nonlinear polynomial fitting and numerical integral, the axial distributions of pressure and its gradient in an axisymmetric rigid vessel with stenosis were obtained, and the distributions related to Reynolds number and the geometry of stenotic vessel were discussed. It shows that with the increasing of stenotic degree or Reynolds number, the fluctuation of pressure and its gradient in stenotic area is intense rapidly, and negative pressure occurs subsequently in the diverging part of stenotic area. Especially when the axial range of stenosis extends, the flow of blood in the diverging part will be more obviously changed. In higher Reynolds number or heavy stenosis, theoretical calculation is mainly in accordance with nast experiments.     

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.     

Effect of an adverse pressure gradient on the streamwise Reynolds stress profile maxima in a turbulent boundary layer

It is widely accepted that in a turbulent boundary layer （TBL） with adverse pressure gradient （APG） an outer peak usually appears in the profile of streamwise Reynolds stress. However, the effect of APG on this outer peak is not clearly understood. In this paper, the effect of APG is analysed using the numerical and experimental results in the literature. Because the effect of upstream flow is inherent in the TBL, we first analyse this effect in TBLs with zero pressure gradient on flat plates. Under the individual effect of upstream flow, an outer peak already appears in the profile of streamwise Reynolds stress when the TBL continues developing in the streamwise direction. The APG accelerates the appearance of the outer peak, instead of being a trigger.     

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.     

ON THE DYNAMICAL STRUCTURE OF THE WIND FIELD OF JUPITER''''S GREAT RED SPOT

Based on the geographic approximation the two-dimensional dynamical structure of the wind fields of Jupiterls Great Red Spot and White Oval BC is obtained. The results of calculation are in good agreement with the observations. Thus, an explanation of the observed dispersion of the velocities along the horizontal streamline is given. The major physical mechanism of this dispersion is as follwos: The distance between two adjacent elliptical streamlines varies along the elliptical streamline; leading to the variance of the normal pressure gradient. Thus, the horizontal velocity V_T has to vary correspondingly so that the Coriolis force can approximately balance the normal pressure gradient Another less important factor, i.e., the change of the Coriolis force parameter with the latitude, is also taken into account. The distributions of the vorticities of GRS and White Oval BC are also calculated.     

ANALYTICAL AND NUMERICAL INVESTIGATIONS OF TWO SPECIAL CLASSES OF GENERALIZED CONTINUUM MEDIA

In this paper,the micromorphic theory and the second gradient theory are proposedwhere the micromorphic model can be reduced to the second gradient model with the vanishing relative deformation between macrodeformation gradient and microdeformation.Analytical solutions for the simple shear problem in the case of a general small strain isotropic elasticity micromorphic model and the second gradient model are presented,respectively.Besides,uniaxial tension of a constrained layer with two different boundary conditions is also analytically solved.Finally,the micromorphic theory is implemented numerically within a two-dimensional plane strain finite element framework by developing two isoparametric elements.     

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed.Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model(DPM). The comparisons results show that the DPM overestimates the inter-porosity coefficient λ and the storage factor ω. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density.The pressure propagation is slower in the direction of larger fracture density.     

Time-Dependent Shape Factors for Uniform and Non-Uniform Pressure Boundary Conditions

Matrix–fracture transfer functions are the backbone of any dual-porosity or dual-permeability formulation. The chief feature within them is the accurate definition of shape factors. To date, there is no completely accepted formulation of a matrix–fracture transfer function. Many formulations of shape factors for instantly-filled fractures with uniform pressure distribution have been presented and used; however, they differ by up to five times in magnitude. Based on a recently presented transfer function, time-dependent shape factors for water imbibing from fracture to matrix under pressure driven flow are proposed. Also new matrix–fracture transfer pressure-based shape factors for instantly-filled fractures with non-uniform pressure distribution are presented in this article. These are the boundary conditions for a case for porous media with clusters of parallel and disconnected fractures, for instance. These new pressure-based shape factors were obtained by solving the pressure diffusivity equation for a single phase using non-uniform boundary conditions. This leads to time-dependent shape factors because of the transient part of the solution for pressure. However, approximating the solution with an exponential function, one obtains constant shape factors that can be easily implemented in current dual-porosity reservoir simulators. The approximate shape factors provide good results for systems where the transient behavior of pressure is short (a case commonly encountered in fractured reservoirs).     

A dual-porosity model for gas reservoir flow incorporating adsorption behaviour—part I. Theoretical development and asymptotic analyses

In this paper a rigorous dual-porosity model is formulated, which accurately represents the coupling between large-scale fractures and the micropores within dual porosity media. The overall structure of the porous medium is conceptualized as being blocks of diffusion dominated micropores separated by natural fractures (e.g. cleats for coal) through which Darcy’s flow occurs. In the developed model, diffusion in the matrix blocks is fully coupled to the pressure distribution within the fracture system. Specific assumptions on the pressure behaviour at the matrix boundary, such as step-time function employed in some earlier studies, are not invoked. The model involves introducing an analytical solution for diffusion within a matrix block, and the resultant combined flow equation is a nonlinear integro-(partial) differential equation. Analyses to the equation in this text, in addition to the theoretical development of the proposed model, include: (1) discussion on the “fading memory” of the model; (2); one-dimensional perturbation solution subject to a specific condition; and (3) asymptotic analyses of the “long-time” and “short-time” responses of the flow. Two previous models, the Warren-Root and the modified Vermeulen models, are compared with the proposed model. The advantages of the new model are demonstrated, particularly for early time prediction where the approximations of these other models can lead to significant error.     

Transient Pressure Behavior of Reservoirs with Discrete Conductive Faults and Fractures

Fractures and faults are common features of many well-known reservoirs. They create traps, serve as conduits to oil and gas migration, and can behave as barriers or baffles to fluid flow. Naturally fractured reservoirs consist of fractures in igneous, metamorphic, sedimentary rocks (matrix), and formations. In most sedimentary formations both fractures and matrix contribute to flow and storage, but in igneous and metamorphic rocks only fractures contribute to flow and storage, and the matrix has almost zero permeability and porosity. In this study, we present a mesh-free semianalytical solution for pressure transient behavior in a 2D infinite reservoir containing a network of discrete and/or connected finite- and infinite-conductivity fractures. The proposed solution methodology is based on an analytical-element method and thus can be easily extended to incorporate other reservoir features such as sealing or leaky faults, domains with altered petrophysical properties (for example, fluid permeability or reservoir porosity), and complicated reservoir boundaries. It is shown that the pressure behavior of discretely fractured reservoirs is considerably different from the well-known Warren and Root dual-porosity reservoir model behavior. The pressure behavior of discretely fractured reservoirs shows many different flow regimes depending on fracture distribution, its intensity and conductivity. In some cases, they also exhibit a dual-porosity reservoir model behavior.     

Flow processes with change of regime in fractured reservoirs

Experimental and industrial observations indicate a strong nonlinear dependence of the parameters of the flow processes in a fractured reservoir on its state of stress. Two problems with change of boundary condition at the well — pressure recovery and transition from constant flow to fixed bottom pressure — are analyzed for such a reservoir. The latter problem may be formulated, for example, so as not to permit closure of the fractures in the bottom zone. For comparison, the cases of linear [1] and nonlinear [2] fractured porous media and a fractured medium [3] are considered, and solutions are obtained in a unified manner using the integral method described in [1]. Nonlinear elastic flow regimes were previously considered in [3–6], where the pressure recovery process was investigated in the linearized formulation. Problems involving a change of well operating regime were examined for a porous reservoir in [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 67–73, May–June, 1991.     

Solution of some problems of flow through fractured porous media

Problems of fluid flow through a fractured porous medium consisting of fractures and blocks with different filter characteristics are solved. The mass exchange between fractures and blocks is assumed to be proportional to the pressure difference between them. The porosity in the fractures is assumed to be negligibly small. Under these assumptions the determination of the pressure fields reduces to the integration of a system of linear differential equations. The solution is found by the operational method using the Efros theorem. The cases of oil reservoir operation by means of both galleries and wells are considered. The solutions are obtained in an analytical form convenient for calculations.Kazan'. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, pp. 94–102, January–February, 1995.     

New Correlative Models to Improve Prediction of Fracture Permeability and Inertial Resistance Coefficient

Presence of fracture roughness and occurrence of nonlinear flow complicate fluid flow through rock fractures. This paper presents a qualitative and quantitative study on the effects of fracture wall surface roughness on flow behavior using direct flow simulation on artificial fractures. Previous studies have highlighted the importance of roughness on linear and nonlinear flow through rock fractures. Therefore, considering fracture roughness to propose models for the linear and nonlinear flow parameters seems to be necessary. In the current report, lattice Boltzmann method is used to numerically simulate fluid flow through different fracture realizations. Flow simulations are conducted over a wide range of pressure gradients through each fracture. It is observed that creeping flow at lower pressure gradients can be described using Darcy’s law, while transition to inertial flow occurs at higher pressure gradients. By detecting the onset of inertial flow and regression analysis on the simulation results with Forchheimer equation, inertial resistance coefficients are determined for each fracture. Fracture permeability values are also determined from Darcy flow as well. According to simulation results through different fractures, two parametric expressions are proposed for permeability and inertial resistance coefficient. The proposed models are validated using 3D numerical simulations and experimental results. The results obtained from these two proposed models are further compared with those obtained from the conventional models. The calculated average absolute relative errors and correlation coefficients indicate that the proposed models, despite their simplicity, present acceptable outcomes; the models are also more accurate compared to the available methods in the literature.     

Large-scale averaging analysis of multiphase flow in fractured reservoirs

The method of large-scale averaging is applied to derive and analyze a dual-porosity model of multiphase flow in naturally fractured reservoirs. The dual-porosity model contains the usual equations based on Darcy's law, and the coupling terms representing the fluid transfer between the matrix and the fractures. Both quasisteady and transient closure schemes are considered to obtain and analyze the fracture and matrix permeability tensors and the fluid transfer terms. The techniques developed here are not restricted to regular geometric fractures. Computational work aimed at showing the implications of the theory behind the derivation of the present dual-porosity model is also described. In particular, comparisons among the dual-porosity model, the single porosity model, and other dual-porosity models are presented through numerical experiments.