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

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

可压可溶两相驱动问题的迎风混合元方法

研究迎风混合元方法求解多孔介质中含弥散的可压缩可混溶两相渗流驱动问题, 利用变分形式和先验估计的理论技巧,得到饱和度近似的L²模和压力近似的H¹模最优阶误差估计,数值实验证实该方法在克服数值扩散和非物理振荡方面是很有效的.     

半正定两相驱动问题的多步有限体积方法及其理论分析

考虑多维半正定两相驱动方程的初边值问题,在非结构网格上构造多步的迎风有限体积格式,利用微分方程先验估计理论证明了格式的离散模形式的误差估计为D(△t~2 h),其中△t和h分别表示时空步长．数值算例进一步验证了格式的有效性．     

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

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

拟抛物方程的两类分裂正定混合元方法

本文对拟抛物方程构造两种分裂对称正定混合元方法.通过适当选取变分形式,格式分裂成两个独立对称正定子格式,并且方法不需要验证LBB条件.收敛性分析表明方法关于变量u和引进的变量σ分别具有L 2(Ω)和H(div;Ω)范数意义下的最优收敛阶.最后,通过数值实验验证了方法的有效性.     

一四元数矩阵方程组的实部半正定解

设Ａ是一个ｎ阶实四元数方阵，若对任意的非零ｎ元列向量ｘ，有ｘＡｘ的实部非负，则称Ａ是一个实部半正定阵，本文给出了实四元数矩阵方程组。     

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

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

各向异性网格下二阶椭圆问题的一个混合元形式

提出了二阶椭圆问题的一个混合变分形式,同时证明了Rariart-Thomas元的各向异性插值性质,并给出了单元的对二阶问题的最优误差估计。     

半正定矩阵及矩阵方程AX=B的反问题

文从研究一类控制系统的实际背景提出对已知实向量x,b求满足Ax=b的对称正定阵A的一类反问题。文[2]与[3]研究了上述反问题在对称正定类、正定类中有解的充要条件及解的一般形式。本文讨论复矩阵方程 AX=B(1)(X,B为m×n阵,A为m×m阵)在半正定、正定、H半正定、H正定类中反问题有解的充要条件及其解集的一般形式。如无特别申明,本文总考虑复矩阵和复向量,其共轭转置用“*”表     

两个四元数自共轭半正定矩阵乘积的特征估计

设A和B均非0的n阶实四元数自共轭矩阵,λ_i及μ_i分别为共特征值(i=1,…,n),且规定|λ₁|≥|λ₂|≥…≥|λ_n|,|μ₁|≥|μ₂|≥…≥|μ_n|,又λ为AB之任意特征值,则λ为实数,且(1)若A≥0,A(?)GL_n(Q),B≥0,B GL_n(Q),则λ≤λ₁μ₁;(2)若A>0或B>0,则|λ|≤|λ₁μ₁|,特别当A>0且B>0时有λ≤λ₁μ₁;(3)若A>0,B∈GL_n(Q),或B>0,A∈GL_n(Q)则|λ|≥|λ_nμ_n|,特别当A>0且B>0时有λ≥λ_nμ_n。     

Mixed Method for Compressible Miscible Displacement with Dispersion in Porous Media

Compressible miscible displacement of one fluid by another in porous media is modelled by a nonlinear parabolic system. A finite element procedure is introduced to approximate the concentration of one fluid and the pressure of the mixture. The concentration is treated by a Galerkin method while the pressure is treated by a parabolic mixed finite element method. The effect of dispersion, which is neglected in [1], is considered. Optimal order estimates in L2 are derived for the errors in the approximate solutions.     

A Characteristic Mixed Finite Element Two-Grid Method for Compressible Miscible Displacement Problem

A nonlinear parabolic system is derived to describe compressible miscible displacement in a porous medium. The concentration equation is treated by a mixed finite element method with characteristics (CMFEM) and the pressure equation is treated by a parabolic mixed finite element method (PMFEM). Two-grid algorithm is considered to linearize nonlinear coupled system of two parabolic partial differential equations. Moreover, the $L^q$ error estimates are conducted for the pressure, Darcy velocity and concentration variables in the two-grid solutions. Both theoretical analysis and numerical experiments are presented to show that the two-grid algorithm is very effective.     

Stokes问题的各向异性误差分析

对二维空间Stokes问题提出了各向异性平行四边形混合有限元逼近格式,证明了其在不要求满足正则性和拟一致条件下的收敛性以及在各向异性条件下对相容误差部分的超收敛估计.     

Stabilized Nonconforming Mixed Finite Element Method for Linear Elasticity on Rectangular or Cubic Meshes

Based on the primal mixed variational formulation, a stabilized nonconforming mixed finite element method is proposed for the linear elasticity on rectangular and cubic meshes. Two kinds of penalty terms are introduced in the stabilized mixed formulation, which are the jump penalty term for the displacement and the divergence penalty term for the stress. We use the classical nonconforming rectangular and cubic elements for the displacement and the discontinuous piecewise polynomial space for the stress, where the discrete space for stress are carefully chosen to guarantee the well-posedness of discrete formulation. The stabilized mixed method is locking-free. The optimal convergence order is derived in the $L^2$-norm for stress and in the broken $H^1$-norm and $L^2$-norm for displacement. A numerical test is carried out to verify the optimal convergence of the stabilized method.     

Design of Finite Element Tools for Coupled Surface and Volume Meshes

Many problems with underlying variational structure involve a coupling of volume with surface effects.A straight-forward approach in a finite element discretiza- tion is to make use of the surface triangulation that is naturally induced by the volume triangulation.In an adaptive method one wants to facilitate"matching"local mesh modifications,i.e.,local refinement and/or coarsening,of volume and surface mesh with standard tools such that the surface grid is always induced by the volume grid. We describe the concepts behind this approach for bisectional refinement and describe new tools incorporated in the finite element toolbox ALBERTA.We also present several important applications of the mesh coupling.     

关于Ciarlet-Raviart混合有限元法的一个注记

本文推广解双调合方程的Ciarlet-Raviart混合有限元方案：用二次元逼近流函数φ．一次元逼近涡度-Δφ．在拟一致三角形剖分的条件下,证明了推广方案具有φ和-Δφ都用二次元逼近的标准Ciarlet－Raviart方案同样的精度阶．     

基于非结构自适应网格的复合有限体积法

利用文献[1]中将Lax-Wendroff格式和Lax-Friedrichs格式整体复合作用构成二维无结构网格上的复合型有限体积法,同时利用Delaunay方法,根据流场流动特性变化的梯度值为指示器对网格进行加密和粗化,实现自适应,并将此方法应用到二维浅水波方程的求解上,进行了二维部分溃坝,倾斜水跃的数值实验.结果表明,该方法是一个计算稳定、能适应复杂的求解域、能很好地捕捉激波、且计算速度快的算法.