求解两类Maxwell方程组棱元离散系统的快速算法和自适应方法

Maxwell方程组棱元离散系统的快速算法和自适应方法是当前计算电磁场中的研究热点和难点. 首先, 针对H(curl)椭圆方程组的棱元离散系统, 通过建立棱元空间的稳定性分解, 设计了相应的快速迭代法和高效预条件子, 并且证明了迭代算法的收敛率和预条件子的条件数均不依赖于模型参数和网格规模. 其次, 针对时谐Maxwell方程组的棱有限元方法, 利用离散的Helmholtz分解, 连续散度为零函数对离散散度为零函数的逼近性和对偶论证, 获得了在L²和H(curl)范数下的拟最优误差估计. 进而设计和分析了相应的两网格法. 最后, 分别针对变系数H(curl)椭圆方程组和不定时谐Maxwell方程组, 考虑了一种不需要标记振荡项和加密单元不需要满足“内节点” 性质的自适应棱有限元法(AEFEM), 并证明了AEFEM的收敛性. 进一步, 当初始网格和Dörfler标记策略参数满足一定的假设条件时, 利用AEFEM的收敛性、误差的整体下界和局部上界估计, 证明了AEFEM的拟最优复杂性.     

定常Navier-Stokes方程的一个二阶非协调三角型混合元格式（英文）

In this paper, a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure. The convergence analysis is presented and optimal error estimates of both broken H~1-norm and L~2-norm for velocity as well as the L~2-norm for the pressure are derived.     

THE FULL DISCRETE DISCONTINUOUS FINITE ELEMENTANALYSIS FOR FIRST-ORDER LINEAR HYPERBOLICEQUATION

1.IntroductionLetfibeaboundeddomaininRZwithpiecewisesmoothboundaryOff,[0,T]beatimeinterval.Considerthefirst-orderhyperbolicproblemasfollowingwhereac~(%,%),p(x,t)~(gi(x,t),pZ(x,t)),Off--(t)~{xEOff:fi(x,t)'ac<0},700istheoutwardunitnormaltoOff;fi(t)~fi\Ofl--(t).Asusual,Off--(t)isreferedtoasinflowboundaryattimet,andOff (t)~Off\Ofl--(t)iscalledoutflowboundaryattimet.FOrsimplicityinfiniteelementanalysis,supposethatboundaryOff--(t)isindependentoft.ThusforalltE(0,T]wecanwriteandproblem(1.0)can…     

基于Crouzeix-Raviart元的界面浸入有限元方法及其收敛性分析

本文对具间断系数的二阶椭圆界面问题提出一种浸入有限元方法(theimmersed finite element method), 即在界面单元上采用依赖于界面的线性多项式空间离散, 而在非界面单元上采用Crouzeix-Raviart非协调元离散. 论证表明, 该方法具有对界面问题解的最优L²-模和H¹-模收敛精度.     

非线性抛物方程的质量集中非协调元分析

将一个低阶Crouzeix-Raviart型非协调三角形元应用到一类非线性抛物方程,并建立了质量集中的半离散和向后Euler全离散逼近格式,在一般各向异性网格上利用插值算子导出了L²-模的最优误差估计.     

拟线性抛物型积分微分方程的一个新最低阶混合元格式

对一类拟线性抛物型积分微分方程构造了一个新的最低阶三角形协调混合元格式,并直接利用单元插值的性质,给出了相应的收敛性分析和H~1-模及L~2-模意义下的最优误差估计.     

电报方程H~1-Galerkin非协调混合有限元分析

主要研究一类电报方程的H~1-Galerkin非协调混合有限元方法,在任意四边形网格剖分下,其逼近空间分别取为类Wilson元与双线性Q_1元,在不需要满足LBB相容性条件及不采用传统的Ritz投影的情况下,得到了与常规有限元方法相同的L~2-模和H~1-模的误差估计,进一步拓展了H~1-Galerkin混合有限元和类Wilson元的应用范围.     

非定常Navier-Stokes方程的稳定化特征有限元法

1引言特征线有限元法是求解对流扩散问题的有效方法。在处理对流占优问题时,表现出了很好的稳定性[8]。对于求解Navier-Stokes方程,文[9]建立了特征有限元格式,并进行了详细分析,但得到的收敛阶O(h~m △t (h~(m 1)/△t))只是拟丰满的。文[10]对此作了非线性稳定性的进一步分析,给出了关于速度和压力的最优误差估计。但目前所有的特征有限元法都要求有限元空间满足inf-sup条件,这就排除了工程实际应用计算方便的低阶有     

Stokes方程的压力梯度局部投影间断有限元法

本文对定常的Stokes方程提出了一种新的间断有限元法,通过将通常的间断Galerkin有限元法与压力梯度局部投影相结合,建立了一个稳定的间断有限元格式,对速度和压力的任意分片多项式空间P_l(K),P_m(K)的间断有限元逼近证明了解的存在唯一性,给出了关于速度和压力的L~2范数的最优误差估计.     

热传导方程一个新的H~1-Galerkin非协调混合有限元格式

对热传导方程提出了一个新的H~1-Galerkin非协调混合有限元格式,其逼近空间不需满足LBB相容性条件,且在不引进传统的Rutz投影的情况下,得到了与以往协调有限元方法相同的L~2-模和H~1-模的误差估计.     

ERROR ESTIMATES FOR SPARSE OPTIMAL CONTROL PROBLEMS BY PIECEWISE LINEAR FINITE ELEMENT APPROXIMATION

Optimization problems with L¹-control cost functional subject to an elliptic partial differential equation(PDE)are considered.However,different from the finite dimensiona l¹-regularization optimization,the resulting discretized L¹norm does not have a decoupled form when the standard piecewise linear finite element is employed to discretize the continuous problem.A common approach to overcome this difficulty is employing a nodal quadrature formula to approximately discretize the L¹-norm.In this paper,a new discretized scheme for the L¹-norm is presented.Compared to the new discretized scheme for L¹-norm with the nodal quadrature formula,the advantages of our new discretized scheme can be demonstrated in terms of the order of approximation.Moreover,finite element error estimates results for the primal problem with the new discretized scheme for the L¹-norm are provided,which confirms that this approximation scheme will not change the order of error estimates.To solve the new discretized problem,a symmetric Gauss-Seidel based majorized accelerated block coordinate descent(sGS-mABCD)method is introduced to solve it via its dual.The proposed sGS-mABCD algorithm is illustrated at two numerical examples.Numerical results not only confirm the finite element error estimates,but also show that our proposed algorithm is efficient.     

一类半线性反应对流扩散模型的特征差分方法和分析

1.引 言如下形式的半线性反应对流扩散方程组分别在生命科学、化学和环境科学中,有大量的应用模型[1-3].其中文献[2-6]分别讨论了方程组(1.1)的各种特殊模型的定性性质.文献[6]讨论了一类线性模型的流线扩散有限元分析.作者在文[7]中,分别利用标准有限元方法和交替方向有限元方法,对(1.1)的一些特殊情形作了数值分析.     

带有记忆项和阻尼项的非线性双曲型方程的有限元方法

This paper studies the finite element method for some nonlinear hyperbolic partial differential equations with memory and dampling terms.A Crank-Nicolson approximation for this kind of equations is presented.By using the elliptic Ritz-Volterra projection,the analysis of the error estimates for the finite element numerical solutions and the optimal H1-norm error estimate are demonstrated.     

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 NEW NONCONFORMING MIXED FINITE ELEMENT SCHEME FOR THE STATIONARY NAVIER-STOKES EQUATIONS

In this article, a new stable nonconforming mixed finite element scheme is proposed for the stationary Navier-Stokes equations, in which a new low order CrouzeixRaviart type nonconforming rectangular element is taken for approximating space for the velocity and the piecewise constant element for the pressure. The optimal order error estimates for the approximation of both the velocity and the pressure in L2-norm are established, as well as one in broken H1-norm for the velocity. Numerical experiments are given which are consistent with our theoretical analysis.     

Weak Galerkin finite element methods for electric interface model with nonhomogeneous jump conditions

In this paper, the weak Galerkin finite element method (WG-FEM) is applied to a pulsed electric model arising in biological tissue when a biological cell is exposed to an electric field. A fitted WG-FEM is proposed to approximate the voltage of the pulsed electric model across the physical media involving an electric interface (surface membrane), and heterogeneous permittivity and a heterogeneous conductivity. This method uses totally discontinuous functions in approximation space and allows the usage of finite element partitions consisting of general polygonal meshes. Optimal pointwise-in-time error estimates in L²-norm and H¹-norm are shown to hold for the semidiscrete scheme even if the regularity of the solution is low on the whole domain. Furthermore, a fully discrete approximation based on backward Euler scheme is analyzed and related optimal error estimates are derived.     

一类非线性双曲型方程有限元方法的误差分析

关于二阶双曲型方程有限元方法的理论研究,已有不少工作,如[1]—[5]。[5]对具Dirichlet边界条件且初边值均取0值的一类非线性双曲方程定解问题的有限元方法,导出了H~1-逼近阶估计,其中,对有关辅助函数u([5],p,151)施加了||?u||_(L~∞(Ω×[0,T]))< ∞的假定。 本文对[5]中研究过的方程,就Dirichlet边界及第三类边界两种情况,给出了半离散Galerkin方法H~1及L~2误差估计。得到的逼近阶都是最佳的,而且,在建立H~1估计的     

鱼骨形网格上二阶方程混合元的超收敛

In this paper, superconvergence of the lowest order Raviart-Thomas mixed finite element approximation for second order Neumann boundary value problem on fishbone shape meshes is analyzed. The main term of the error between the exact solution and the finite element interpolating function is determined by Bramble-Hilbert lemma on the individual finite element. A part of the main term of the error on two adjacent finite elements can be cancelled along the special direction, and thus the higher order error estimate is obtained on the whole domain by summation. Compared with the general finite element error estimate,the convergence rate can be increased from order one to order two in L2-norm by postprocessing superconvergence technique.     

非线性双相滞热传导方程的新混合元超收敛分析

针对非线性双相滞热传导方程,建立了一种自由度少且自然满足B-B条件的新混合元逼近格式.在半离散格式下,基于双线性元的高精度结果,分别导出了原始变量的H~1模及中间变量的L~2模的超逼近性质,进而,借助于插值后处理算子,得到了原始及中间变量比传统误差高一阶的整体超收敛结果.     

A Second Order Nonconforming Triangular Mixed Finite Element Scheme for the Stationary Navier-Stokes Equations

In this paper,a nonconforming triangular mixed finite element scheme with second order convergence behavior is proposed for the stationary Navier-Stokes equations.The new nonconforming triangular element is taken as approximation space for the velocity and the linear element for the pressure.The convergence analysis is presented and optimal error estimates of both broken H1-norm and L2-norm for velocity as well as the L2-norm for the pressure are derived.