NUMERICAL METHOD AND APPLICATION FOR THREEDIMENSIONAL NONLINEAR SYSTEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

For the system of multilayer dynamics of fluids in porous media, the second order upwind finite difference fractional steps schemes applicable to parallel arithmetic are put forward. Some techniques, such as calculus of variations, energy method, multiplicative commutation rule of difference operators, decomposition of high order difference operators and prior estimates are adopted. Optimal order estimates are derived to determine the error in the second order approximate solution. These methods have already been applied to the numerical simulation of migration-accumulation of oil resources.     

Numerical method for nonlinear two-phase displacement problem and its application

For the three-dimensional nonlinear two-phase displacement problem, the modified upwind finite difference fractional steps schenles were put forward. Some techniques, such as calculus of variations, induction hypothesis, decomposition of high order difference operators, the theory of prior estimates and techniques were used. Optimal order estimates were derived for the error in the approximation solution. These methods have been successfully used to predict the consequences of seawater intrusion and protection projects.     

Modified characteristic finite difference fractional step method for moving boundary value problem of nonlinear percolation system

A fractional step scheme with modified characteristic finite differences running in a parallel arithmetic is presented to simulate a nonlinear percolation system of multilayer dynamics of fluids in a porous medium with moving boundary values. With the help of theoretical techniques including the change of regions, piecewise threefold quadratic interpolation, calculus of variations, multiplicative commutation rule of difference operators, multiplicative commutation rule of difference operators, decomposition of high order difference operators, induction hypothesis, and prior estimates, an optimal order in l ² norm is displayed to complete the convergence analysis of the numerical algorithm. Some numerical results arising in the actual simulation of migration-accumulation of oil resources by this method are listed in the last section.     

Upwind finite difference method for miscible oil and water displacement problem with moving boundary values

The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.     

TWO-GRID METHOD FOR CHARACTERISTICS FINITE-ELEMENT SOLUTION OF 2D NONLINEAR CONVECTION-DOMINATED DIFFUSION PROBLEM

Introduction Convection diffusionequationgovernssuchphenomenaastheflowofheatwithina movingfluid,thetransportofdissolvednutrientsorcontaminantswithinthegroundwater,andthetransportofasurfactantortracerwithinanincompressibleoilinapetroleum reservoir.Weconsid…     

Mathematical model of two-phase fluid nonlinear flow in low-permeability porous media with applications

IntroductionItisasuccessfulexampleinadevelopmentstoryofscienceandtechnologyformechanicsoffluidsinporousmediatocombinewithengineeringtechnology .Fieldsinfluencedbythemechanicsinvolveddevelopmentofoil_gasandgroundwaterresources,controlonseawaterintrusionandsubsidenceandgeologichazards,geotechnicalengineeringandbioengineering ,andairlineindustry[1～ 7].Aproblemonnonlinearflowinlow_permeabilityporousmediaisbutonlyabasiconeindifferentkindsofengineeringfields,butalsooneoffrontlineresearchfieldsofmod…     

Characteristic finite difference method and application for moving boundary value problem of coupled system

The coupled system of multilayer dynamics of fluids in porous media is to describe the history of oil-gas transport and accumulation in basin evolution.It is of great value in rational evaluation of prospecting and exploiting oil-gas resources.The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values.A kind of characteristic finite difference schemes is put forward,from which optimal order estimates in l^2 norm are derived for the error in the approximate solutions.The research is important both theoretically and practically for the model analysis in the field,the model numerical method and software development.     

Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

Characteristic fractional step finite difference method for nonlinear section coupled system

For the section coupled system of multilayer dynamics of fluids in porous media, a parallel scheme modified by the characteristic finite difference fractional steps is proposed for a complete point set consisting of coarse and fine partitions. Some tech- niques, such as calculus of variations, energy method, twofold-quadratic interpolation of product type, multiplicative commutation law of difference operators, decomposition of high order difference operators, and prior estimates, are used in theoretical analysis. Optimal order estimates in 12 norm are derived to show accuracy of the second order approximation solutions. These methods have been used to simulate the problems of migration-accumulation of oil resources.     

旋转机械系统多自由度非线性动力学数值分析

首先将转子系统的动力响应问题归结为2n维未知向量v的一阶非线性动力学方程dv／dt=Hv＋f（v,t）,并给出了求解这一方程的一次近似式法和三次多项式迫近法。在非稳态、非线性油膜力等作用下,以刚性Jeffcott转子与112个自由度的汽轮发电机组低压转子系统为例,用上述求解方法分析了它们的动力响应及非线性动力学特性;其问,还将计算结果与Runge—Kutta法、Newmark法的相应结果进行了比较,并深入讨论了数值稳定性问题。汽轮发电机组的算例表明对一些具有较复杂的非线性右端项,、同时规模又较大的问题,如果采用四阶Runge—Kutta法,才算几步就因数值骤然增大而失控;但若用同样步长的一次近似式,由于它是一种显式的无条件稳定算式,则计算过程迅速且结果合理可靠。     

Numerical methods for theoretical analysis of foil bearing dynamics

The bearing under consideration is a self-acting, gas bearing with a flexibly supported foil. In order to conduct the theoretical analysis, three systems were identified: a rotor, a gas film and a flexible structure. Each of these elements has been mathematically analysed and analytical equations have been formulated. It was found, that the gas flow in the bearing can be described by the Reynolds equation, whereas a spring-damper model was selected for the structural analysis. The Reynolds equation is a differential equation the exact solution to which is unknown. The work describes the finite difference method in detail, where the partial derivatives in the Reynolds equation are replaced by a system of algebraic equations. In order to solve the resulting system, the Alternating Direction Implicit method (ADI) was used. Based on that and on the analysis of interactions between the three systems the equations have been computed. The accuracy of methods has been verified by means of series of numerical tests. Chosen results are described in this paper.     

Step-like contrast structure for class of nonlinear singularly perturbed optimal control problem

The step-like contrast structure for a class of nonlinear singularly perturbed optimal control problems is considered. The existence of the step-like contrast structure for the singularly perturbed optimal control problem is proved by equivalence, which is based on the necessary conditions. The authors not only give the conditions under which there exists a step-like contrast structure, but also determine where the internal transition time is. Meanwhile, the uniformly valid asymptotic expansion of the step-like contrast structure solution is constructed by the direct scheme method. Finally, an example is presented to show the result.     

Numerical solution of nonlinear ordinary differential equation for a turning point problem

By using the method in[3],several useful estimations of the derivatives of the solutionof the boundary value problem for a nonlinear ordinary differential equation with a turningpoint are obtained.With the help of the technique in[4],the uniform convergence on thesmall parameterεfor a difference scheme is proved.At the end of this paper,a numericalexample is given.The numerical result coincides with theoretical analysis.     

基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

饱和多孔介质大变形分析的一种有限元-有限体积混合方法

建立了饱和多孔介质大变形分析的一种有限元-有限体积混合计算方法.将饱和多孔介质视为由固体骨架和孔隙水组成的两相体,其基本方程包括动力平衡方程和渗流连续方程.基于u-p假定和更新的Lagrange方法,饱和多孔介质的动力平衡方程在空间域内采用有限元方法进行离散,而渗流连续方程在空阃域内则采用有限体积法进行离散.通过两个数值算例,一维有限弹性固结和动力荷载作用下堤坝动力响应的计算,验证了该方法的有效性.     

Numerical method of Rayleigh-Stokes problem for heated generalized second grade fluid with fractional derivative

In this paper,we consider the Rayleigh-Stokes problem for a heated generalized second grade fluid（RSP-HGSGF）with fractional derivative.An effective numerical method for approximating RSP-HGSGF in a bounded domain is presented.The stability and convergence of the method are analyzed.Numerical examples are presented to show the application of the present technique.     

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.     

Modified characteristic finite difference fractional step method for moving boundary value problem of percolation coupled system

For the coupled system with moving boundary values of multilayer dynamicsof fluids in porous media,a characteristic finite difference fractional step scheme appli-cable to the parallel arithmetic is put forward.Some techniques,such as the change ofregions,the calculus of variations,the piecewise threefold quadratic interpolation,themultiplicative commutation rule of difference operators,the decomposition of high orderdifference operators,and the prior estimates,are adopted.The optimal order estimatesin the l2norm are derived to determine the error in the approximate solution.This nu-merical method has been successfully used to simulate the flow of migration-accumulationof the multilayer percolation coupled system.Some numerical results are well illustratedin this paper.     

可压缩流体流动多尺度混合有限元数值方法研究

可压缩流体是天然油藏中广泛存在的一种流体,研究其在多孔介质中的渗流规律对于油藏开发具有重要意义。本文采用多尺度混合有限元方法,对可压缩流体渗流问题进行了研究。考虑流体的可压缩性以及介质形变,推导得到了可压缩流体渗流问题的多尺度计算格式。数值计算结果表明,多尺度混合有限元适于求解非均质性和可压缩流问题,具有节省计算量、计算精度高等优势,对于实际大规模油藏模拟具有重要意义。     

储液容器内液体晃荡的非线性动力学分析

基于非线性波动理论模型,求解储液容器内液体晃动的固有频率、模态及动力学响应问题。流体使用u_s-u_p状态方程,利用ABAQUS软件的自适应网格技术,建立储液容器液体晃动数学模型,通过施加水平简谐激励得到液体晃动的固有频率和模态,并与解析解对比,验证了该方法的准确性与可行性。然后,分析了矩形储液容器在多种激励作用下液体非线性晃动响应特性。