Wilson砖特征值的渐进展开式

在有限元方法中,特征值的渐进展开式对于研究数值特征值的逼近性质有着非常重要的应用,使用林群等发展的广义Bramble-Hilbert引理,证明了Wilson砖特征值的渐进展开式.     

用有限元方法求解界面特征值问题

用有限元方法计算椭圆型界面特征值问题,实验数据显示近似特征值的变化规律:界面特征值问题中系数的间断性对协调和非协调Crouzeix-Raviart有限元特征值的收敛性并无影响,而且对协调有限元特征值外推以后得到高精度的解,相应的外推值还提供特征值下界;Crouzeix-Raviart元特征值提供特征值下界,这对一般有界区域如"镂空"型区域也成立.另外,还展示近似特征函数的图形.     

非协调元特征值渐近下界

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

混合网格下的有限元特征值重构

利用有限元后处理技术在混合网格上重构了线性有限元解,使其梯度具有超收敛性,在此基础上利用Rayleigh商重构特征值,获得了线元特征值的四阶超收敛结果.     

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

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

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

本文建立了二阶椭圆特征值问题 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 投影算子。利用这个定理就可以把二阶线性椭圆方程有限元逼近位移超收敛的结果移植到特征值问题,     

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

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

浅谈微积分所蕴含的特征值计算思想方法

数学思想方法是数学的灵魂.主要阐述微积分中泰勒公式的应用,泰勒公式不但能用来近似计算数学常数π,而且能用来理论证明两个粗糙的近似值通过简单的组合加工技术(现代人称之为"外推法")产生更准确的近似值,也把外推法推广用于求解Steklov特征值问题,数值算例表明非协调Crouzeix-Raviart有限元外推法既能提高解的精度又能提供准确特征值的下界.     

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

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

关于有限元离散方程特征值的界

各种有限元离散方程(包括协调元、非协调元及混合元等)的特征值的上、下界以及条件数的估计,早已引起了注意(见[1],[2],[5]).这里我们用一种简明的方法来处理这类问题.本文用到的关于Sobolev空间中的标准符号见[3].下面出现的常数c,c_1,c_2等在不同地方可能取不同的值.     

二阶椭圆问题的一类广义有限元法

本文提出了求解二阶椭圆问题的一类广义有限元方法,分析了广义有限元方法的优越性,证明了二阶椭圆问题的广义有限元方法具有比标准的Galerkin有限元方法更高阶的收敛速度,根据插值算子的性质,进一步证明了有限元解的亏量迭代校正收敛到广义有限元解,并用数值例子说明广义有限元方法是有效的.     

AN ADAPTIVE INVERSE ITERATION FEM FOR THE INHOMOGENEOUS DIELECTRIC WAVEGUIDES

We introduce an adaptive finite element method for computing electromagnetic guided waves in a closed, inhomogeneous, pillared three-dimensional waveguide at a given frequency based on the inverse iteration method. The problem is formulated as a generalized eigenvalue problems. By modifying the exact inverse iteration algorithm for the eigenvalue problem, we design a new adaptive inverse iteration finite element algorithm. Adaptive finite element methods based on a posteriori error estimate are known to be successful in resolving singularities of eigenfunctions which deteriorate the finite element convergence. We construct a posteriori error estimator for the electromagnetic guided waves problem. Numerical results are reported to illustrate the quasi-optimal performance of our adaptive inverse iteration finite element method.     

三维Poisson方程特征值的四种有限元解及比较

应用三维EQ1rot元、三维Crouzeix-Raviart元、八节点等参数元、四面体线性元计算三维Poisson方程的近似特征值.计算结果表明:三维EQ1rot元和三维Crouzeix-Raviart元特征值下逼近准确特征值,八节点等参数元、四面体线性元特征值上逼近准确特征值,三维EQr1ot元和三维Crouzeix-Raviart元外推特征值下逼近准确特征值.计算结果还表明三维Crouzeix-Raviart元是一种计算效率较高的非协调元.     

Symmetric modified finite volume element methods for self-adjoint elliptic and parabolic problems

The finite volume element method is a discretization technique for partial differential equations, but in general case the coefficient matrix of its linear system is not symmetric, even for the self-adjoint continuous problem. In this paper we develop a kind of symmetric modified finite volume element methods both for general self-adjoint elliptic and for parabolic problems on general discretization, their coefficient matrix are symmetric. We give the optimal order energy norm error estimates. We also prove that the difference between the solutions of the finite volume element method and symmetric modified finite volume element method is a high order term.     

A new mixed finite element method for the Stokes problem

We propose a new mixed formulation of the Stokes problem where the extra stress tensor is considered. Based on such a formulation, a mixed finite element is constructed and analyzed. This new finite element has properties analogous to the finite volume methods, namely, the local conservation of the momentum and the mass. Optimal error estimates are derived. For the numerical implementation of this finite element, a hybrid form is presented. This work is a first step towards the treatment of viscoelastic fluid flows by mixed finite element methods.     

Finite-element method for time-dependent euler equation

This paper presents finite element methods to approximate inviscid incompressible flow problems. First we emphasize the conservation properties of these problems, and we show that finite element methods appear as a very natural way to find conservative schemes such as Arakawa's scheme. We give convergence theorems and an error analysis of finite element discretization schemes. We turn then to the time differencing problem. We derive stability and convergence results for a second-order semi-implicit scheme and for the leap-frog scheme.     

Weighted quadrature rules for finite element methods

We discuss the numerical integration of polynomials times non-polynomial weighting functions in two dimensions arising from multiscale finite element computations. The proposed quadrature rules are significantly more accurate than standard quadratures and are better suited to existing finite element codes than formulas computed by symbolic integration. We validate this approach by introducing the new quadrature formulas into a multiscale finite element method for the two-dimensional reaction–diffusion equation.     

Some remarks on the flux-free finite element method for immiscible two-fluid flows

We give some theoretical considerations on the the flux-free finite element method for the generalized Stokes interface problem arising from the immiscible two-fluid flow problems. In the flux-free finite element method, the flux constraint is posed as another Lagrange multiplier to keep the zero-flux on the interface. As a result, the mass of each fluid is expected to be preserved at every time step. We first study the effect of discontinuous coefficients (viscosity and density) on the error of the standard finite element approximations very carefully. Then, the analysis is extended to the flux-free finite element method.     

A sparse finite element method with high accuracy Part I

Summary. In this paper, we develop and analyze a new finite element method called the sparse finite element method for second order elliptic problems. This method involves much fewer degrees of freedom than the standard finite element method. We show nevertheless that such a sparse finite element method still possesses the superconvergence and other high accuracy properties same as those of the standard finite element method. The main technique in our analysis is the use of some integral identities. Received October 1, 1995 / Revised version received August 23, 1999 / Published online February 5, 2001