油水驱动半正定问题迎风混合元格式及其数值分析

研究了不可压缩油水两相渗透流驱动问题.在扩散矩阵仅是半正定的假设条件下,提出了迎风混合元方法.混合元方法近似压力方程,饱和度方程的对流项用Godunov迎风格式来处理,扩散项则用推广的混合元来逼进,并推导出格式的误差估计.此种格式的优越性表现在两个方面:首先是饱和度方程的扩散矩阵仅是半正定的;二是摒弃了特征格式所限制的周期性条件,更适用于实际问题.     

半导体瞬态问题的混合迎风有限体积格式

半导体器件的瞬时状态由包含3个拟线性偏微分方程所组成的方程组的初边值问题来描述.在三角剖分的基础上,对椭圆型的电子位势方程采用混合有限体积元法来逼近,对对流扩散型的电子浓度和空穴浓度方程采用迎风有限体积元方法来逼近,并进行了详细的理论分析,得到了最优阶的误差估计结果.最后,针对混合有限体积元法和迎风有限体积元法分别单独使用以及两种方法结合使用的情形给出了不同的数值算例.     

半导体器件瞬态模拟的对称正定混合元方法

提出具有对称正定特性的混合元格式求解非稳态半导体器件瞬态模拟问题。提出一个最小二乘混合元方法、一个新的具有分裂和对称正定性质的混合元格式和一个解经典混合元方程的对称正定失窃工格式求解电场位势和电场强度方程;提出一个最小二乘混合元格式求解关于电子与空穴浓度的非稳态对流扩散方程,浓度函数和流函数被同时求解;采用标准的有限元方法求解热传导方程。建立了误差分析理论。     

土壤水流与溶质耦合运移问题的混合元-迎风广义差分法研究

本文研究了一维非饱和土壤水流与溶质耦合运移问题的数学模型, 建立了求其数值解的守恒混合元-迎风广义差分格式. 对非线性土壤水分入渗方程, 采用守恒混合元法进行离散模拟, 同时得到了土壤含水量和水分通量; 而对对流-扩散形式的溶质运移方程, 利用迎风的广义差分法离散求解. 且分析了解的存在唯一性, 并讨论了误差估计. 最后给出数值算例, 模拟结果表明利用本文格式来求解非饱和土壤水流与溶质耦合运移问题是可靠的, 且该格式具有稳定性和可实用性.     

多孔介质中两相不可压混溶驱替的差分流线扩散法及其误差分析

多孔介质中两相不可压混熔驱替问题可描述为椭圆和抛物耦合的非线性偏微分方程组，对椭圆方程采用混合元方法，而对抛物方程采用差分流线扩散法，本文构造了求解该问题的差分流线扩散-混合元格式，最后，给出所构造格式按L^∞(L^2）模的拟最优误差阶估计。     

不可压混溶驱替问题的流线扩散──混合元数值模拟

采用标准元模拟不可压混溶流问题，当扩散系数矩阵小过剖分参数时，有限元格式仅能给出比最优精度低一阶的逼近解，格式稳定性差并伴有强烈的数值弥散现象．为了克服上述缺陷，本文对压力方程采用混合元，而对浓度方程采用流线扩散格式，在扩散矩阵为线性的假定下，证明了该格式具有较标准元更高的逼近精度（比最优阶低1／2）和更好的稳定性．     

求解对流扩散方程的低耗散中心迎风格式

以四阶CWENO重构为基础,通过将对流项采用低耗散中心迎风格式离散,扩散项采用四阶中心差分格式离散,对得到的半离散格式采用四阶龙格库塔方法在时间方向上推进,得到一种求解对流扩散方程的高阶有限差分格式.数值结果验证了该格式的四阶精度和基本无振荡特性.     

三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

迎风紧致格式与驱动方腔流动问题的直接数值模拟

本文给出了一种求解不可压缩流动问题的高精度差分格式,即迎风紧致格式.出发方程采用二维非定常原始变量Naiver-Stokes方程组.在差分方程中,对流项采用三阶精度的迎风紧致差分,其余空间导数项采用四阶紧致差分.本文利用该差分格式在等距网格上数值模拟了驱动方腔流动中的分离涡运动.在257×257的细网格上,Re数最高计算到10000.Re≤5000时的计算结果与前人结果符合得很好.当Re≥7500时发现流动不存在定常层流解而为非定常周期性解,并首次给出了非定常解的结果。     

一维定常型对流占优扩散方程的一类迎风有限体积格式

本文针对一维定常型对流占优扩散方程提出了一类迎风有限体积格式 .该格式对对流项具有二阶精度 ,对扩散项保持一阶精度 ,符合对流占优扩散问题强对流、弱扩散的特点 .     

An upwind-mixed method on changing meshes for two-phase miscible flow in porous media

Two-phase miscible flow in porous media is governed by a system of nonlinear partial differential equations. In this paper, the upwind-mixed method on dynamically changing meshes is presented for the problem in two dimensions. The pressure is approximated by a mixed finite element method and the concentration by a method which upwinds the convection and incorporates diffusion using an expanded mixed finite element method. The method developed is shown to obtain almost optimal rate error estimate. When the method is modified we can obtain the optimal rate error estimate that is well known for static meshes. The modification of the scheme is the construction of a linear approximation to the solution, which is used in projecting the solution from one mesh to another. Finally, numerical experiments are given.     

Superconvergence analysis of splitting positive definite nonconforming mixed finite element method for pseudo-hyperbolic equations

In this paper, a new splitting positive definite nonconforming mixed finite element method is proposed for pseudo-hyperbolic equations, in which a quasi-Wilson quadrilateral element is used for the flux p, and the bilinear element is used for u. Superconvergence results in ｜｜·｜｜div,h norm for p and optimal error estimates in L2 norm for u are derived for both semi-discrete and fully discrete schemes under almost uniform meshes.     

双线形元的各向异性后验误差估计

The main aim of this paper is to give an anisotropic posteriori error estimator. We firstly study the convergence of bilinear finite element for the second order problem under anisotropic meshes.By using some novel approaches and techniques,the optimal error estimates and some superconvergence results are obtained without the regularity assumption and quasi-uniform assumption requirements on the meshes.Then,based on these results, we give an anisotropic posteriori error estimate for the second problem.     

A Posteriori Error Estimators of Gradient Recovery Type for Fem of a Model Optimal Control Problem

In this paper, we derive an a posteriori error estimator of gradient recovery type for a model optimal control problem. We show that the a posteriori error estimator is equivalent to the discretization error in a norm of energy type on general meshes. Furthermore, when the solution of the control problem is smooth and the meshes are uniform, it is shown to be asymptotically exact.     

Raviart-Thomas混合元的超收敛

考虑二阶椭圆方程Dirichlet边值问题在正则矩形网格上k阶RaviartThomas混合有限元的超收敛.对有限元解经插值处理后,与通常的有限元最优误差估计相比,收敛速度提高了两阶.     

CONVERGENCE ANALYSIS FOR A NONCONFORMING MEMBRANE ELEMENT ON ANISOTROPIC MESHES

Regular assumption of finite element meshes is a basic condition of most analysis offinite element approximations both for conventional conforming elements and nonconform-ing elements.The aim of this paper is to present a novel approach of dealing with theapproximation of a four-degree nonconforming finite element for the second order ellipticproblems on the anisotropic meshes.The optimal error estimates of energy norm and L~2-norm without the regular assumption or quasi-uniform assumption are obtained based onsome new special features of this element discovered herein.Numerical results are givento demonstrate validity of our theoretical analysis.     

A posteriori error estimates for optimal control problems constrained by convection-diffusion equations

We propose a characteristic finite element discretization of evolutionary type convection-diffusion optimal control problems. Nondivergence-free velocity fields and bilateral inequality control constraints are handled. Then some residual type a posteriori error estimates are analyzed for the approximations of the control, the state, and the adjoint state. Based on the derived error estimators, we use them as error indicators in developing efficient multi-set adaptive meshes characteristic finite element algorithm for such optimal control problems. Finally, one numerical example is given to check the feasibility and validity of multi-set adaptive meshes refinements.     

非线性抛物型方程组的二次有限体积元方法及其误差估计

讨论基于三角形网格的二维非线性抛物型方程组的有限体积元方法,其中试探函数空间为二次Lagrange元,检验函数空间为分片常数函数空间,对问题的全离散格式证明了最优的能量模误差估计。最后给出一个相关数值算例以验证格式的有效性。     

ON THE ANISOTROPIC ACCURACY ANALYSIS OF ACM''''S NONCONFORMING FINITE ELEMENT

The main aim of this paper is to study the superconvergence accuracy analysis of thefamous ACM's nonconforming finite element for biharmonic equation under anisotropicmeshes. By using some novel approaches and techniques, the optimal anisotropic inter-polation error and consistency error estimates are obtained. The global error is of orderO(h~2). Lastly, some numerical tests are presented to verify the theoretical analysis.     

多孔介质驱动问题的混合元最小二乘特征有限元方法及其收敛分析

１．引言多孔介质二相驱动问题的数学模型是偶合的非线性偏微分方程组的初边值问题．该问题可转化为压力方程和浓度方程［１－４］．浓度方程一般是对流占优的对流扩散方程,它的对流速度依赖于比浓度方程的扩散系数大得多的Ｆａｒｃｙ速度．因此Ｄａｒｃｙ速度的求解精度直接影响着浓度的求解精度．为了提高速度的求解精度,７０年代Ｐ．Ａ．Ｒａｖｉａｔ和Ｊ．Ｍ．Ｔｈｏｍａｓ提出混合有限元方法［５］．Ｊ．ＤｏｕｇｌａｓＪｒ,Ｔ．Ｆ．Ｒｕｓｓｅｌｌ,Ｒ．Ｅ．Ｅｗｉｎｇ,Ｍ．Ｆ．Ｗｈｅｅｌｅｒ［１］－［４］,［９］,［１２］袁…