A MIXED FINITE ELEMENT AND CHARACTERISTIC MIXED FINITE ELEMENT FOR INCOMPRESSIBLE MISCIBLE DARCY-FORCHHEIMER DISPLACEMENT AND NUMERICAL ANALYSIS

In this paper a mixed finite element-characteristic mixed finite element method is discussed to simulate an incompressible miscible Darcy-Forchheimer problem.The flow equation is solved by a mixed finite element and the approximation accuracy of Darch-Forchheimer velocity is improved one order.The concentration equation is solved by the method of mixed finite element,where the convection is discretized along the characteristic direction and the diffusion is discretized by the zero-order mixed finite element method.The characteristics can confirm strong stability at sharp fronts and avoids numerical dispersion and nonphysical oscillation.In actual computations the characteristics adopts a large time step without any loss of accuracy.The scalar unknowns and its adjoint vector function are obtained simultaneously and the law of mass conservation holds in every element by the zero-order mixed finite element discretization of diffusion flux.In order to derive the optimal 3/2-order error estimate in L^(2) norm,a post-processing technique is included in the approximation to the scalar unknowns.Numerical experiments are illustrated finally to validate theoretical analysis and efficiency.This method can be used to solve such an important problem.

A BLOCK-CENTERED UPWIND APPROXIMATION OF THE SEMICONDUCTOR DEVICE PROBLEM ON A DYNAMICALLY CHANGING MESH

The numerical simulation of a three-dimensional semiconductor device is a fundamental problem in information science. The mathematical model is defined by an initialboundary nonlinear system of four partial differential equations: an elliptic equation for electric potential, two convection-diffusion equations for electron concentration and hole concentration, and a heat conduction equation for temperature. The first equation is solved by the conservative block-centered method. The concentrations and temperature are computed by the block-centered upwind difference method on a changing mesh, where the block-centered method and upwind approximation are used to discretize the diffusion and convection, respectively. The computations on a changing mesh show very well the local special properties nearby the P-N junction. The upwind scheme is applied to approximate the convection, and numerical dispersion and nonphysical oscillation are avoided. The block-centered difference computes concentrations, temperature, and their adjoint vector functions simultaneously.The local conservation of mass, an important rule in the numerical simulation of a semiconductor device, is preserved during the computations. An optimal order convergence is obtained. Numerical examples are provided to show efficiency and application.

Theory and application of numerical simulation method of capillary force enhanced oil production

A kind of second-order implicit upwind fractional step finite difference methods are presented for the numerical simulation of coupled systems for enhanced （chemical） oil production with capillary force in the porous media. Some techniques, e.g., the calculus of variations, the energy analysis method, the commutativity of the products of difference operators, the decomposition of high-order difference operators, and the theory of a priori estimate, are introduced. An optimal order error estimate in the l2 norm is derived. The method is successfully used in the numerical simulation of the enhanced oil production in actual oilfields. The simulation results are satisfactory and interesting.

Theory and Application of Numerical Simulation of Chemical Flooding in High Temperature and High Salt Reservoirs

Applications, theoretical analysis and numerical methods are introduced for the simulation of mechanical models and principles of the porous flow in high temperature, high salt, complicated geology and large-scale reservoirs in this paper. Considering petroleum geology, geochemistry, computational permeation fluid mechanics and computer technology, we state the models of permeation fluid mechanics and put forward a sequence of implicit upwind difference iteration schemes based on refined fractional steps of the upstream, which can compute the pressures, the saturation and the concentrations of different chemistry components. A type of software applicable in major industries has been completed and carried out in numerical analysis and simulations of oil extraction in Shengli Oil-field, which brings huge economic benefits and social benefits. This software gives many characters: spatial steps are taken as ten meters, the number of nodes is up to hundreds of thousands and simulation time period can be tens of years and the high-order accuracy can be promised in numerical data. Precise analysis is present for simplified models of this type and that provides a tool to solve the international famous problem.

The Theory and Application of Upwind Finite Difference Fractional Steps Procedure for Seawater Intrusion

Numerical simulation and theoretical analysis of seawater intrusion is the mathematical basis for modern environmental science. Its mathematical model is the nonlinear coupled system of partial differential equations with initial-boundary problems. For a generic case of a three-dimensional bounded region, two kinds of finite difference fractional steps pro- cedures are put forward. Optimal order estimates in norm are derived for the error in the approximation solution. The present method has been successfully used in predicting the consequences of seawater intrusion and protection projects.

An Upwind Finite Volume Element Method for Nonlinear Convection Diffusion Problem

A class of upwind finite volume element method based on tetrahedron partition is put forward for a nonlinear convection diffusion problem. Some techniques, such as calculus of variations, commutating operators and the a priori estimate, are adopted. The a priori error estimate in L2-norm and H1-norm is derived to determine the error between the approximate solution and the true solution.

An Upwind Mixed Volume Element-Fractional Step Method on a Changing Mesh for Compressible Contamination Treatment from Nuclear Waste

In this paper the authors discuss the numerical simulation problem of threedimensional compressible contamination treatment from nuclear waste.The mathematical model is defined by an initial-boundary nonlinear convection-diffusion system of four partial differential equations:a parabolic equation for the pressure,two convection-diffusion equations for the concentrations of brine and radionuclide and a heat conduction equation for the temperature.The pressure appears within the concentration equations and heat conduction equation,and the Darcy velocity controls the concentrations and the temperature.The pressure is solved by the conservative mixed volume element method,and the order of the accuracy is improved by the Darcy velocity.The concentration of brine and temperature are computed by the upwind mixed volume element method on a changing mesh,where the diffusion is discretized by a mixed volume element and the convection is treated by an upwind scheme.The composite method can solve the convection-dominated diffusion problems well because it eliminates numerical dispersion and nonphysical oscillation and has high order computational accuracy.The mixed volume element has the local conservation of mass and energy,and it can obtain the brine and temperature and their adjoint vector functions simultaneously.The conservation nature plays an important role in numerical simulation of underground fluid.The concentrations of radionuclide factors are solved by the method of upwind fractional step difference and the computational work is decreased by decomposing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.By the theory and technique of a priori estimates of differential equations,we derive an optimal order result in L^(2) norm.Numerical examples are given to show the effectiveness and practicability and the composite method is testified as a powerful tool to solve the well-known actual problem.

THE UPWIND FINITE DIFFERENCE FRACTIONAL STEPS METHOD FOR NONLINEAR COUPLED SYSTEM OF DYNAMICS OF FLUIDS IN POROUS MEDIA

为在多孔的媒介的液体的多层的动力学的非线性的联合系统，第二订；首先订迎风的有限差别对平行算术适用的部分步计划被提出，；二维；三维的计划被用来形成一个完全的集合。一些技术例如变化的演算，差别操作员的趋于增加的交换统治，高顺序差别操作员的分解；优先的估计，被采用。在 L 的最佳的顺序估计[2 ] 标准被导出在第二份订单决定错误近似答案。这个方法已经被用于油资源的迁居累积的数字模拟。

AConservative Upwind Approximation on Block-Centered Difference for Chemical Oil Recovery Displacement Problem

A kind of conservative upwind method is discussed for chemical oil recovery displacement in porous media.The mathematical model is formulated by a nonlinear convection-diffusion system dependent on the pressure,Darcy velocity,concentration and saturations.The flow equation is solved by a conservative block-centered method,and the pressure and Darcy velocity are obtained at the same time.The concentration and saturations are determined by convection-dominated diffusion equations,so an upwind approximation is adopted to eliminate numerical dispersion and nonphysical oscillation.Block-centered method is conservative locally.An upwind method with block-centered difference is used for computing the concentration.The saturations of different components are solved by the method of upwind fractional step difference,and the computational work is shortened significantly by dividing a three-dimensional problem into three successive one-dimensional problems and using the method of speedup.Using the variation discussion,energy estimates,the method of duality,and the theory of a priori estimates,we complete numerical analysis.Finally,numerical tests are given for showing the computational accuracy,efficiency and practicability of our approach.