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

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

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.     

压力梯度对边界层两相流动稳定性的影响

分析了有压力梯度的边界层两相流动稳定性,推导出类似于Saffman理论的修正的稳定性方程,数值计算采用高精度的谱方法。结果说明,压力梯度对边界层两相流动稳定性有显著的影响,顺压梯度增强流动稳定性,而逆压梯度则促进流动失稳。在不同的压力梯度和浓度下,Stokes数对流动稳定性的影响是一致的,存在一个临界Stokes数,小Stokes数促进流动失稳,而大Stokes数则提高临界雷诺数,抑制流动失稳的最佳Stokes数为10的量级。     

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.     

Reynolds Stress Budgets in Couette and Boundary Layer Flows

Reynolds stress budgets for both Couette and boundary layer flows are evaluated and presented. Data are taken from direct numerical simulations of rotating and non-rotating plane turbulent Couette flow and turbulent boundary layer with and without adverse pressure gradient. Comparison of the total shear stress for the two types of flows suggests that the Couette case may be regarded as the high Reynolds number limit for the boundary layer flow close to the wall. The limit values of turbulence statistics close to the wall for the boundary layer for increasing Reynolds number approach the corresponding Couette flow values. The direction of rotation is chosen so that it has a stabilizing effect, whereas the adverse pressure gradient is destabilizing. The pressure-strain rate tensor in the Couette flow case is presented for a split into slow, rapid and Stokes terms. Most of the influence from rotation is located to the region close to the wall, and both the slow and rapid parts are affected. The anisotropy for the boundary layer decreases for higher Reynolds number, reflecting the larger separation of scales, and becomes close to that for Couette flow. The adverse pressure gradient has a strong weakening effect on the anisotropy. All of the data presented here are available on the web [36]. This revised version was published online in July 2006 with corrections to the Cover Date.     

Bulk flow parameters in turbulent wall boundary layers near separation

Experiments were conducted in a turbulent boundary layer near separation along a flat plate. The pressure gradient in flow direction was varied such that three significant boundary layer configurations could be maintained. The flow in the test section thus had simultaneously a region of favourable pressure gradient, a region of strong adverse pressure gradient with boundary layer separation and a region of reattached boundary layer. Specially designed fine probes facilitated the measurements of skin friction and velocity distribution very close to the wall. Bulk flow parameters such as skin friction coefficient C _f, Reynold's number Re_δ2 and shape factors H and G, which are significant characteristics of wall boundary layers were evaluated. The dependence of these parameters on the Reynolds number and along the test section was explored and the values were compared with other empirical and analytical formulae known in the literature.     

Application of the Dorodnitsyn finite element method to swirling boundary layer flow

The Dorodnitsyn finite element method for turbulent boundary layer flow with surface mass transfer is extended to include axisymmetric swirling internal boundary layer flow. Turbulence effects are represented by the two-layer eddy viscosity model of Cebeci and Smith¹ with extensions to allow for the effect of swirl. The method is applied to duct entry flow and a 10 degree included-angle conical diffuser, and produces results in close agreement with experimental measurements with only 11 grid points across the boundary layer. The introduction of swirl (w_e/u_e = 0.4) is found to have little effect on the axial skin friction in either a slightly favourable or adverse pressure gradient, but does cause an increase in the displacement area for an adverse pressure gradient. Surface mass transfer (blowing or suction) causes a substantial reduction (blowing) in axial skin friction and an increase in the displacement area. Both suction and the adverse pressure gradient have little influence on the circumferential velocity and shear stress components. Consequently in an adverse pressure gradient the flow direction adjacent to the wall is expected to approach the circumferential direction at some downstream location.     

Wake layer in a thermal turbulent boundary layer with pressure gradient

The open equations of thermal turbulent boundary layer subjected to pressure gradient have been analysed by method of matched asymptotic expansions at large Reynolds number. The flow is divided into outer wake layer and inner wall layer. The asymptotic expansions are matched by Millikan-Kolmogorov hypothesis. The temperature profile in overlap region yields composite law which reduce to log. law for moderate pressure gradient and inverse half power law for strong adverse pressure gradient. In case of a shallow thermal wake, the matching result of outer wake layer reduces to composite temperature defect law, which is more general than the classical log. law. The comparison of data for thermal boundary layer with strong adverse pressure gradient is also considered. Received on 26 May 1998     

Global pressure relaxation for laminar two-dimensional internal flow

This study extends the reduced Navier–Stokes (RNS) global pressure relaxation procedure developed by Rubin and co-workers for external flow to internal flow applications. The streamwise pressure gradient is split into a backward-differenced or initial value component, as in boundary layer marching, and a forward-differenced or boundary value component that represents the elliptic downstream effects. The streamwise convection terms are upwind-differenced and all other streamwise derivatives are backward-differenced. We thus obtain a standard boundary layer marching technique imbedded in a conventional line relaxation technique. For compressible flow the pressure iteration determines the interior flow interation as well as the inlet mass flux that is consistent with the outflow pressure boundary condition. Results have been computed for incompressible flow in both rectangular and curved channels, and for subsonic compressible flow in the simulation of an aerofoil in a wind tunnel. Converged solutions were obtained over a range of Reynolds numbers generating small to moderately large separation bubbles.     

动边界双重介质油藏低速非达西渗流试井模型

裂缝性油藏中基质岩块的渗透率一般很低,大量岩心测试实验证实在基质岩块内的液体渗流和在一定含水饱和度下的气体渗流将偏离达西渗流,往往出现低速非达西渗流,表现出启动压力梯度以及流体流动边界不断向外扩展等特殊现象。本文充分考虑启动压力梯度与动边界的影响,建立了微可压缩双重介质油藏低速非达西渗流的试井数学模型,对时间和空间变量...     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

Analytical solutions of the boundary-value problem of nonstationary flow of viscoplastic medium between two plates

Summary Nonstationary flow of a viscoplastic medium between two parallel plates is considered for the case of a varying pressure gradient. The problem is reduced to the Stephan problem, with the condition on the boundary separating the flow domain from the quasi-rigid domain. Four multiparameter families of exact solutions are found. The first family describes the flow decelerations up to a full stop. The second family determines the development of the flow from the state of rest as the pressure gradient increases. The third family describes the development of the flow for the case where (1) the pressure gradient is constant and exceeds the threshold value related to the yield stress, (2) the upper plate does not move, and (3) the lower plate moves with a constant acceleration. Finally, the fourth family determines the flow retardation, when the pressure gradient is constant and is less than the threshold value. The decrease in the flow of the viscoplastic medium can be achieved for certain values of parameters by increasing the quasi-rigid domain, whereas the viscoplastic flow remains unchanged. Received 7 October 1998; accepted for publication 8 April 1999     

Analytical investigation of boundary layer growth and swirl intensity decay rate in a pipe

In this research, the developing turbulent swirling flow in the entrance region of a pipe is investigated analytically by using the boundary layer integral method. The governing equations are integrated through the boundary layer and obtained differential equations are solved with forth-order Adams predictor-corrector method. The general tangential velocity is applied at the inlet region to consider both free and forced vortex velocity profiles. The comparison between present model and available experimental data demonstrates the capability of the model in predicting boundary layer parameters (e.g. boundary layer growth, shear rate and swirl intensity decay rate). Analytical results showed that the free vortex velocity profile can better predict the boundary layer parameters in the entrance region than in the forced one. Also, effects of pressure gradient inside the boundary layer is investigated and showed that if pressure gradient is ignored inside the boundary layer, results deviate greatly from the experimental data.     

Permissible Height of the Surface Roughness in a Turbulent Boundary Layer with a Streamwise Pressure Gradient

The results of an experimental investigation of the effect of the streamwise pressure gradient in a turbulent boundary layer on the permissible height of the surface roughness of bodies in an incompressible fluid flow are presented. The permissible roughness Reynolds number for which the characteristics of the turbulent boundary layer remain the same as in the case of flow past a smooth surface is determined.     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

Computing axisymmetric,laminar internal flows

A simple computational scheme is developed to compute laminar flows inside axisymmetric ducts. It is based on the Keller box method where the equations are approximated at the centre of the downstream face of each computational box. The coupling between the pressure gradient and the velocities for internal flow has been observed to introduce stability problems for the Keller box method that are not present for external, boundary layer flow problems. The difference scheme for the velocities is coupled to an iterative scheme to solve for the pressure gradient at each axial step. Example results for developing flow in a pipe and in a 2° conical diffuser are presented.     

Response of a skewed turbulent boundary layer to favourable pressure gradient

Experimental results are reported for the response to a favourable pressure gradient of an initially turbulent boundary layer (Re _θ?≈?1600) developing on a flat plate with its leading edge skewed at 60° to the approach flow. The pressure gradient orthogonal to the leading edge is nominally the same as that which was shown by Escudier et?al. [(1998) Exp Fluids 25: 491–502] to cause extreme thinning of a two-dimensional (2D) (i.e. unskewed) turbulent boundary layer and the intermittency in the immediate vicinity of the surface to fall to zero, i.e. an apparent laminarisation of the boundary layer. In the case of the skewed boundary layer, the responses of the turbulence and mean-flow structures are qualitatively similar to those for the 2D situation. However, the streamwise pressure gradient is much weaker than for the 2D experiment and the extent of the changes it produces is much reduced. Even so, the changes are considerably greater than would be expected from the magnitude of the streamwise pressure gradient.     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

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.     

An experimental study of turbulent boundary layers approaching separation

The present paper deals with the experimental analysis of a strong decelerated turbulent boundary layer developed on a flat plate. The aim of the study was to examine the effects of pressure gradient on a non-equilibrium boundary layer while indicating local areas of equilibrium flow. The effect of the Reynolds number on a turbulent boundary layer developed with matching the external pressure gradient conditions was also analysed. The emphasis was on the analysis of mean flow statistics i.e. mean velocity profiles, streamwise Reynolds stress and the effect of large- and small-scale interactions by analysing the skewness factor and energy isocontours maps. The comparative analysis of the external data indicated that the structure of the turbulent boundary layer depends not only on local effects of pressure gradient but also on the upstream history of the flow. For the same condition of pressure gradient, the increased momentum is observed near the wall with the increase of the Reynolds number at the Incipient Detachment, where increased turbulence production is also observed, leading to the failure of the outer scaling methods. Surprisingly, the effect of the Reynolds number decays at the intermittent transitory detachment where similar profiles were observed. The upper inflection point in the mean profile corresponded well with the outer maximum of the Reynolds stress and zero crossing of skewness factor. Position of this point occurs at different locations, depending on the flow history effects. The last observation demonstrates that the inflection points results from large- and small-scale interactions, which led to the increased convection velocity of small scales near the wall.