数值积分对特征值有限元外推的影响

§1.引言在分片一致三角形剖分下,用线性有限元法解特征值问题求得近似特征值λ、λ~(h/2).[1]证明了对λ~h.λ~(h/2)作外推可提高收敛阶:     

关于特征值问题有限元方法研究的一个说法

受林群的微积分哲学公式((相对真理)/(绝对真理)=0.9)的启示,总结了近些年来关于特征值问题有限元方法的研究,并发现其背后同样隐藏着该哲学公式,换句话说,所追求的是特征值问题的有限元数值解和真解的零距离,其实就是追求真解的过程,要经多道(即0.9,0.99,0.999,…),再将比例,即(数值解)/(真解),提到1.     

特征值问题有限元逼近位移超收敛的一个定理

本文建立了二阶椭圆特征值问题 m 次有限元逼近位移超收敛定理:‖Pu-u_h‖_(0,∞)≤ch~(m+2),m≥2‖Pu-u_h‖_(0,∞)≤ch~(m+1+(2/ξ),m≥3其中,1<ξ<2,u_h 是 m 次有限元特征函数,u 是某一精确特征函数,P 是 Ritz 投影算子。利用这个定理就可以把二阶线性椭圆方程有限元逼近位移超收敛的结果移植到特征值问题,     

特征值问题的Lagrange型各向异性有限元方法

在各向异性网格下首先研究了二阶椭圆特征值问题算子谱逼近的若干抽象结果.然后将这些结果具体应用于线性和双线性Lagrange型协调有限元,得到了与传统有限元网格剖分下相同的最优误差估计,从而拓宽了已有的成果.     

用Fourier Galerkin方法求解间断特征值问题

本文我们考虑一个带间断系数的特征值问题,使用Fourier G a lerk in方法求解.数值试验表明特征值收敛速度达到三阶,而特征函数收敛速度达到2.5阶.表明此方法对间断系数问题非常有效.     

杆有限元模型的特征值反问题

1引言非均匀杆的轴向自由振动方程为其中L为杆的长度,A(x)为杆的横截面积,E(x)为弹性模量,ρ(x)是杆的密度．杆的主振动的一般形式为u(x,t)=u(x)sin(wt φ), (2)其中u(x)满足d/dx[E(x)A(x)du/dx] λρ(x)A(x)u=0,λ=ω2．(3)     

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

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

抛物型界面问题的变网格有限元方法

本文针对抛物型界面问题,提出了一种线性三角形变网格有限元方法.其主要思路是针对空间变量采用有限元离散,对时间变量采用差分离散,但是不同时刻的有限元剖分网格可以不同.在不引入Ritz投影这一传统分析工具的情况下,得到了最优误差估计结果,使得证明过程更加简洁.给出的数值算例验证了理论分析的正确性.     

各向异性网格下特征值问题的有限元分析

主要目的是在各向异性网格下研究二阶椭圆特征值问题的两类非协调有限元—类Wilson矩形元和Carey三角形元—的收敛性分析.通过新的技巧和方法,得到了与传统有限元网格剖分下相同的特征对的最优误差估计.推广了已有的结果.     

Finite Volume Methods for Eigenvalue Problems

Based on a linear finite element space, in this paper, two symmetric finite volume schemes are proposed for self-adjoint elliptic boundary eigenvalue problems. Both convergence and superconvergence are discussed.This revised version was published online in October 2005 with corrections to the Cover Date.     

Immersed Interface Finite Element Methods for Elasticity Interface Problems with Non-Homogeneous Jump Conditions

<正>In this paper,a class of new immersed interface finite element methods (IIFEM) is developed to solve elasticity interface problems with homogeneous and non-homogeneous jump conditions in two dimensions.Simple non-body-fitted meshes are used.For homogeneous jump conditions,both non-conforming and conforming basis functions are constructed in such a way that they satisfy the natural jump conditions. For non-homogeneous jump conditions,a pair of functions that satisfy the same non-homogeneous jump conditions are constructed using a level-set representation of the interface.With such a pair of functions,the discontinuities across the interface in the solution and flux are removed;and an equivalent elasticity interface problem with homogeneous jump conditions is formulated.Numerical examples are presented to demonstrate that such methods have second order convergence.     

非协调元特征值渐近下界

利用有限元收敛速度下界的结果获得某些非协调元方法新的Aubin-Nitsche估计形式,然后再结合非协调元特征值的展开式获得不需要额外条件下非协调元特征值渐近下界的结果.     

Local and Parallel Finite Element Algorithms for Eigenvalue Problems

Abstract Some new local and parallel finite element algorithms are proposed and analyzed in this paper foreigenvalue problems.With these algorithms, the solution of an eigenvalue problem on a fine grid is reduced tothe solution of an eigenvalue problem on a relatively coarse grid together with solutions of some linear algebraicsystems on fine grid by using some local and parallel procedure.A theoretical tool for analyzing these algorithmsis some local error estimate that is also obtained in this paper for finite element approximations of eigenvectorson general shape-regular grids.     

Model Reduction Methods for Solving Symmetric Rational Eigenvalue Problems

We consider large and sparse rational eigenproblems where the rational term is of low rank compared to the dimension of the problem. Exploiting model‐order reduction techniques the dimension can be reduced considerably, and problems of this type can be solved very efficiently. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Mixed Finite Element Method for Fourth Order Eigenvalue Problems

A new mixed finite element method for approximating eigenpairs of IV order elliptic eigenvalue problems with Dirichlet boundary conditions has been given. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)