特征值问题协调/非协调有限元后验误差分析

本文研究对称椭圆特征值问题的有限元后验误差估计,包括协调元和非协调元,具有下列特色:(1)对协调/非协调元建立了有限元特征函数uh的误差与相应的边值问题有限元解的误差在局部能量模意义下的恒等关系式,该边值问题的右端为有限元特征值λh与uh的乘积,有限元解恰好为uh.从而边值问题有限元解在能量模意义下的局部后验误差指示子,包括残差型和重构型后验误差指示子,成为有限元特征函数在能量模意义下的局部后验误差指示子.(2)讨论了协调有限元特征函数的基于插值后处理的梯度重构型后验误差估计,对有限元特征函数的导数得到了最大模意义下的渐近准确局部后验误差指示子.     

平面弹性的一个新的Locking-free非协调有限元

针对平面弹性问题构造了一个Locking-free的矩形非协调有限元,并证明该有限元格式关于λ有一致最优收敛阶,其离散误差的能量模为O（h2）,L2模为O（h3）.最后给出了数值实验对理论结果进行了验证.     

Sine-Gordon方程的混合有限体积元方法及数值模拟

研究了在Dirichlet边界条件和Neumann边界条件下一维sine-Gordon方程的混合有限体积元方法.通过引入将试探函数空间映射到检验函数空间的迁移算子γh,结合混合有限元方法和有限体积元方法,构造了半离散格式,时间显式和隐式全离散混合有限体积元格式.给出了显格式离散解的稳定性分析,并得到了三种格式的最优阶误差估计.最后,给出数值算例来验证理论分析结果和数值格式的有效性.     

广义Rayleigh商与有限元法2-网格离散方案

研究广义Rayleigh商和高效率有限元计算方案,做了下列工作：1）把Rayleigh商加速技巧推广到非自共轭问题,定义了算子型广义Rayleigh商和弱形式型广义Rayleigh商,并建立了近似特征向量及其广义Rayleigh商之间的基本关系式.2）在误差估计式中用有限元特征值的陡度取代准确特征值的陡度,得到新的误差估计式.3）在许进超和周爱辉工作的基础上建立了解非自共轭椭圆微分算子特征值问题的有限元2-网格离散方案,并用于协调有限元法和非协调有限元法.从理论分析和数值实验两个方面证明了2-网格方案的有效性.4）把解自共轭椭圆微分算子特征值问题的迭代Galerkin法、插值校正法和梯度重构法推广到非自共轭椭圆微分算子特征值问题.     

二阶椭圆问题的混合有限元法的泡函数稳定性及其后验误差估计

本文讨论二阶椭圆问题的混合有限元逼近的一种泡函数稳定性,并给出其基于简化的稳定化格式的先验误差估计和后验误差估计.该方法较通常的格式(例如,Raviaxt-Thomas方法的同阶格式)节省大量的自由度.     

伪双曲方程的新分裂式正定混合元方法

提出一类二阶伪双曲型方程的新的分裂正定式混合有限元方法.给出了半离散和全离散格式误差估计及其格式的稳定性.与传统的混合元相比,所提出的格式有几个优点:首先所提出的格式能够分裂成两个独立的积分微分子格式并且不需要求解匹配方程组系统;其次不必满足LBB相容性条件.     

一类抛物型方程组的特征修正区域分解有限元方法

在实际生产和科学研究中,有许多物理问题的数学模型为抛物型方程组问题,如可压缩核废料污染问题,地下水资源问题,杨青提出了差分格式和有限元格式,应用先验估计得到了最优的l^2和L^2模误差估计,江城顺等利用交替方向有限元方法得到了H^1模和L^2模误差估计．杨国强等采用显式可解的三层差分格式求解二维方程组得到了H^1模误差估计．     

热传导—对流问题的混合有限元分析注记

本文研究流体力学中的热传导－对流问题,分析其混合有限元解的存在性和误差估计,并给出一些混合有限元格式的例子。     

环形空腔内自然对流问题的混合有限元方法的误差估计

本文从研究环形空腔内自然对流问题的流体运动状态及温度分布规律出发,构造出Boussinesq方程组初、边值问题,给出了两种混合有限元计算格式及它们的有限元误差估计.     

三维热传导型半导体问题的交替方向特征有限元方法

提出交替方向特征有限元方法,对电场位势方程采用混合元格式,对电子,空穴浓度方程采用交替方向特征有限元格式,对温度方程提出交替方向格式.应用向量积计算及先验估计理论和技巧,得到最佳的L2误差估计.     

A posteriori estimates for eigenvalue/vector approximations

We combine abstract eigenvalue/eigenvector estimates (from our earlier work) with a saturation assumption for finite element solution of associated stationary problem to obtain a posteriori estimates of the accuracy of finite element Rayleigh–Ritz approximations. Attention will be payed to the interplay between the accuracy estimate for the finite element method and a strategy for generating an adapted mesh. The obtained results use a preconditioned residuum of Neymeyr and extend his study of eigenvalue approximations with eigenvector estimates. We also prove that this eigenvalue estimator is equivalent to the global error. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Multiscale finite element discretizations based on local defect correction for the biharmonic eigenvalue problem of plate buckling

In this paper, we study the multiscale finite element discretizations about the biharmonic eigenvalue problem of plate buckling. On the basis of the work of Dai and Zhou (SIAM J. Numer. Anal. 46[1] [2008] 295‐324), we establish a three‐scale scheme, a multiscale discretization scheme, and the associated parallel version based on local defect correction. We first prove a local priori error estimate of finite element approximations, then give the error estimates of multiscale discretization schemes. Theoretical analysis and numerical experiments indicate that our schemes are suitable and efficient for eigenfunctions with local low smoothness.     

Asymptotic expansion and extrapolation for the Eigenvalue approximation of the biharmonic eigenvalue problem by Ciarlet–Raviart scheme

We propose and analyze the Ciarlet–Raviart mixed scheme for solving the biharmonic eigenvalue problem with bilinear finite elements. We derive a higher order convergence rate for eigenvalue and eigenfunction approximations. Furthermore, we give an asymptotic expansion of the eigenvalue error, from which an efficient extrapolation and an a posteriori error estimate for the eigenvalue are given. Finally, numerical experiments illustrating the theoretical results are reported. This author was supported by China Postdoctoral Sciences Foundation.     

Error analysis for the computation of zeros of regular Coulomb wave function and its first derivative

In 1975 one of the coauthors, Ikebe, showed that the problem of computing the zeros of the regular Coulomb wave functions and their derivatives may be reformulated as the eigenvalue problem for infinite matrices. Approximation by truncation is justified but no error estimates are given there.

The class of eigenvalue problems studied there turns out to be subsumed in a more general problem studied by Ikebe et al. in 1993, where an extremely accurate asymptotic error estimate is shown.

In this paper, we apply this error formula to the former case to obtain error formulas in a closed, explicit form.

     

非稳态电导介质中Maxwell方程组的反问题

本文主要考虑非稳态电导介质的Maxwell 方程组. 本文考查通过有限组的边界区域观测值决定关于本构方程中系数ε, ζ, μ 和电导率系数σ 的反问题, 利用Carleman 估计证明该反问题的Lipschitz稳定性.     

An accurate a posteriori error estimator for the Steklov eigenvalue problem and its applications

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine mesh and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.

     

Error Estimation for Eigenvalues of Unbounded Linear Operators and an Application to Energy Levels in Graphene Quantum Dots

The eigenvalue problem for linear differential operators is important since eigenvalues correspond to the possible energy levels of a physical system. It is also important to have good estimates of the error in the computed eigenvalues. In this work, we use spline interpolation to construct approximate eigenfunctions of a linear operator using the corresponding eigenvectors of a discretized approximation of the operator. We show that an error estimate for the approximate eigenvalues can be obtained by evaluating the residual for an approximate eigenpair. The interpolation scheme is selected in such a way that the residual can be evaluated analytically. To demonstrate that the method gives useful error bounds, we apply it to a problem originating from the study of graphene quantum dots where the goal was to investigate the change in the spectrum from incorporating electron–electron interactions in the potential.     

Central polynomials and matrix invariants

LetK be a field, charK=0 andM _n (K) the algebra ofn×n matrices overK. If λ=(λ₁,…,λ _m ) andμ=(μ ₁,…,μ _m ) are partitions ofn ² let wherex ₁,…,x _n ²,y ₁,…,y _n ² are noncommuting indeterminates andS _n ² is the symmetric group of degreen ². The polynomialsF ^{λ, μ} , when evaluated inM _n (K), take central values and we study the problem of classifying those partitions λ,μ for whichF ^{λ, μ} is a central polynomial (not a polynomial identity) forM _n (K). We give a formula that allows us to evaluateF ^{λ, μ} inM(K) in general and we prove that if λ andμ are not both derived in a suitable way from the partition δ=(1, 3,…, 2n−3, 2n−1), thenF ^{λ, μ} is a polynomial identity forM _n (K). As an application, we exhibit a new class of central polynomials forM _n (K). In memory of Shimshon Amitsur Research supported by a grant from MURST of Italy.     

Acceleration of two-grid stabilized mixed finite element method for the Stokes eigenvalue problem

This paper provides an accelerated two-grid stabilized mixed finite element scheme for the Stokes eigenvalue problem based on the pressure projection. With the scheme, the solution of the Stokes eigenvalue problem on a fine grid is reduced to the solution of the Stokes eigenvalue problem on a much coarser grid and the solution of a linear algebraic system on the fine grid. By solving a slightly different linear problem on the fine grid, the new algorithm significantly improves the theoretical error estimate which allows a much coarser mesh to achieve the same asymptotic convergence rate. Finally, numerical experiments are shown to verify the high efficiency and the theoretical results of the new method.     

Steklov eigenvalue problem and its applications An accurate a posteriori error estimator for the

In this paper, a type of accurate a posteriori error estimator is proposed for the Steklov eigenvalue problem based on the complementary approach, which provides an asymptotic exact estimate for the approximate eigenpair. Besides, we design a type of cascadic adaptive finite element method for the Steklov eigenvalue problem based on the proposed a posteriori error estimator. In this new cascadic adaptive scheme, instead of solving the Steklov eigenvalue problem in each adaptive space directly, we only need to do some smoothing steps for linearized boundary value problems on a series of adaptive spaces and solve some Steklov eigenvalue problems on a low dimensional space. Furthermore, the proposed a posteriori error estimator provides the way to refine meshes and control the number of smoothing steps for the cascadic adaptive method. Some numerical examples are presented to validate the efficiency of the algorithm in this paper.