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

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

Hamilton系统的连续有限元法

利用常微分方程的连续有限元法,对非线性Hamilton系统证明了连续一次、二次有限元法分别是2阶和3阶的拟辛格式,且保持能量守恒;连续有限元法是辛算法对线性Hamilton系统,且保持能量守恒．在数值计算上探讨了辛性质和能量守恒性,与已有的辛算法进行对比,结果与理论相吻合．     

Analysis of finite element approximations of a phase field model for two-phase fluids

This paper studies a phase field model for the mixture of two immiscible and incompressible fluids. The model is described by a nonlinear parabolic system consisting of the nonstationary Stokes equations coupled with the Allen-Cahn equation through an extra phase induced stress term in the Stokes equations and a fluid induced transport term in the Allen-Cahn equation. Both semi-discrete and fully discrete finite element methods are developed for approximating the parabolic system. It is shown that the proposed numerical methods satisfy a discrete energy law which mimics the basic energy law for the phase field model. Error estimates are derived for the semi-discrete method, and the convergence to the phase field model and to its sharp interface limiting model are established for the fully discrete finite element method by making use of the discrete energy law. Numerical experiments are also presented to validate the theory and to show the effectiveness of the combined phase field and finite element approach.

     

Interior estimates for a low order finite element method for the Reissner–Mindlin plate model

Interior error estimates are obtained for a low order finite element introduced by Arnold and Falk for the Reissner–Mindlin plates. It is proved that the approximation error of the finite element solution in the interior domain is bounded above by two parts: one measures the local approximability of the exact solution by the finite element space and the other the global approximability of the finite element method. As an application, we show that for the soft simply supported plate, the Arnold–Falk element still achieves an almost optimal convergence rate in the energy norm away from the boundary layer, even though optimal order convergence cannot hold globally due to the boundary layer. Numerical results are given which support our conclusion. This revised version was published online in June 2006 with corrections to the Cover Date.     

Convection-diffusion-reaction problems on a B-type mesh

One-dimensional singularly perturbed problems with two small parameters are considered. Numerical methods for such problems are discussed in several papers, but on a Shishkin-type mesh. The first optimal result of convergence in an energy norm on a Bakhvalov-type mesh for one-dimensional convection-diffusion problem was given by Roos in [9]. In this paper we analyze Galerkin finite element method on a Bakhvalov-type mesh for two-parameter convection-diffusion-reaction problems. In the interpolation error analysis, instead of the usual interpolation operator in the finite element space we use a quasi-interpolant with improved stability properties. We prove that the finite element method for these problems is uniformly convergent in the energy norm. Numerical results confirm our theoretical analysis and show first-order convergence rate. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the calculation of energy release rate for curved cracks of rubbery materials

The energy release rate G for 2-D rubbery material problems with curved cracks is calculated with finite element solutions in this paper. Two approaches, a generalized domain integral method and the virtual crack extension method, are investigated. The generalized domain integral method is demonstrated to be “patch independent” and therefore a complicated finite element mesh adjacent to the crack tip area is not required.     

The modified method of characteristics with mixed finite element domain decomposition procedures for the transient behavior of a semiconductor device

For the transient behavior of a semiconductor device, the modified method of characteristics with mixed finite element domain decomposition procedures applicable to parallel arithmetic is put forward. The electric potential equation is described by the mixed finite element method, and the electric, hole concentration and heat conduction equations are treated by the modified method of characteristics finite element domain decomposition methods. Some techniques, such as calculus of variations, domain decomposition, characteristic method, energy method, negative norm estimate and prior estimates and techniques are employed. Optimal order estimates in L² norm are derived for the error in the approximation solution. Thus the well‐known theoretical problem has been thoroughly and completely solved.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 353–368 2012     

A constitutive model for granular materials with microstructures using the concept of energy relaxation

In this paper, we present a constitutive model for granular materials exhibiting microstructures using the concept of energy relaxation. Within the framework of Cosserat continuum theory the free energy of the material is enriched with an interaction energy potential taking into account the counter rotations of the particles. The enhanced energy potential fails to be quasiconvex. Energy relaxation theory is employed to compute the relaxed energy which yields all possible displacement and micro-rotations field fluctuations as minimizers. Based on a two-field variational principle the constitutive response of the material is derived. The developed constitutive model is then implemented in a finite element analysis program using the finite element method. Numerical simulations are presented to observe the localized deformation phenomenon in a granular medium. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Analysis of some low-order finite element schemes for Mindlin-Reissner and Kirchhoff plates

Summary We set up a framework for analyzing mixed finite element methods for the plate problem using a mesh dependent energy norm which applies both to the Kirchhoff and to the Mindlin-Reissner formulation of the problem. The analysis techniques are applied to some low order finite element schemes where three degrees of freedom are associated to each vertex of a triangulation of the domain. The schemes proceed from the Mindlin-Reissner formulation with modified shear energy.Dedicated to Professor Ivo Babuka on the occasion of his 60th birthday     

An object-oriented hybrid symbolic/numerical approach for the development of finite element codes

The association of object-oriented programming and symbolic computation techniques introduces certain changes in finite element code organization. The purpose of this approach is to speed up the design of new formulations. Previous papers have described the basic concepts of the method. In this paper, the focus is placed on functional aspects of symbolic tools for the development of finite element formulations. Two practical examples are used to illustrate this point. The first is a space-time formulation for an incompressible flow driven by the Navier–Stokes equations, and the second is a finite element derivation of the total potential energy for linear elasticity.     

An Error Analysis of Discontinuous Finite Element Methods for the Optimal Control Problems Governed by Stokes Equation

In this article, an abstract framework for the error analysis of discontinuous finite element method is developed for the distributed and Neumann boundary control problems governed by the stationary Stokes equation with control constraints. A priori error estimates of optimal order are derived for velocity and pressure in the energy norm and the L²-norm, respectively. Moreover, a reliable and efficient a posteriori error estimator is derived. The results are applicable to a variety of problems just under the minimal regularity possessed by the well-posedness of the problem. In particular, we consider the abstract results with suitable stable pairs of velocity and pressure spaces like as the lowest-order Crouzeix–Raviart finite element and piecewise constant spaces, piecewise linear and constant finite element spaces. The theoretical results are illustrated by the numerical experiments.     

抛物方程时间连续有限元全离散格式的超收敛性

利用张量积分解和时间方向单元正交分解,证明了线性抛物型方程的时间连续全离散有限元在单元节点和内部的特征点的超收敛性.并用连续有限元计算了非线性Schrodinger方程,验证了能量的守恒性.计算结果与理论相吻合.     

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 posteriori error estimation with finite element methods of lines for one-dimensional parabolic systems

Summary Consider the solution of one-dimensional linear initial-boundary value problems by a finite element method of lines using a piecewiseP ^th -degree polynomial basis. A posteriori estimates of the discretization error are obtained as the solutions of either local parabolic or local elliptic finite element problems using piecewise polynomial corrections of degreep+1 that vanish at element ends. Error estimates computed in this manner are shown to converge in energy under mesh refinement to the exact finite element discretization error. Computational results indicate that the error estimates are robust over a wide range of mesh spacings and polynomial degrees and are, furthermore, applicable in situations that are not supported by the analysis.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR 90-0194; by the U.S. Army Research Office under Contract Number DAAL03-91-G-0215; and by the National Science Foundation under Institutional Infrastructure Grant Number CDA-8805910     

非协调有限元的全塑性分析

本文采用满足相容条件的非协调有限元模型以解决全塑性分析中有限元解的数值精度问题.文中讨论了该模型适用于全塑性分析的机理和判据,还设计了一个确定塑性极限载荷的算法.     

On an Exact Boundary Technique

In Marshall & Mitchell (1973) a finite element was introducedwhich matches Dirichlet boundary data exactly and was comparedwith the standard bilinear finite element which matches suchdata at only a finite number of points on the boundary. Ourpurpose here is to compare these same two elements with respectto the energy norm and the process of minimizing the energyfunctional.     

Adaptive numerical solution of thick plates using first‐order shear deformation theory. Part II: Adaptivity

This is the second in a pair of articles concerned with the adaptive finite element solution of Riessner‐Mindlin thick plates modeled using first‐order shear deformation theory. This article is concerned with enhancing the a posteriori energy‐error estimators developed in Part I in order to accomodate transition elements in the finite element mesh. The resulting estimators are then used in an adaptive finite element model employing transition elements and the subsequent results discussed and compared with those in Part I. A major part of the article is devoted to identifying a novel patch assembly node algorithm for using the ZZ recovery‐type estimator with transition elements. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 227–253, 2003.     

Numerical analysis of relaxed micromagnetics by penalised finite elements

Summary. Some micromagnetic phenomena in rigid (ferro-)magnetic materials can be modelled by a non-convex minimisation problem. Typically, minimising sequences develop finer and finer oscillations and their weak limits do not attain the infimal energy. Solutions exist in a generalised sense and the observed microstructure can be described in terms of Young measures. A relaxation by convexifying the energy density resolves the essential macroscopic information. The numerical analysis of the relaxed problem faces convex but degenerated energy functionals in a setting similar to mixed finite element formulations. The lowest order conforming finite element schemes appear instable and nonconforming finite element methods are proposed. An a priori and a posteriori error analysis is presented for a penalised version of the side-restriction that the modulus of the magnetic field is bounded pointwise. Residual-based adaptive algorithms are proposed and experimentally shown to be efficient. Received June 24, 1999 / Revised version received August 24, 2000 / Published online May 4, 2001     

A two-grid finite element method for the Allen-Cahn equation with the logarithmic potential

In this paper, we present a two-grid finite element method for the Allen-Cahn equation with the logarithmic potential. This method consists of two steps. In the first step, based on a fully implicit finite element method, the Allen-Cahn equation is solved on a coarse grid with mesh size H. In the second step, a linearized system whose nonlinear term is replaced by the value of the first step is solved on a fine grid with mesh size h. We give the energy stabilities of the traditional finite element method and the two-grid finite element method. The optimal convergence order of the two-grid finite element method in H¹ norm is achieved when the mesh sizes satisfy h = O(H²). Numerical examples are given to demonstrate the validity of the proposed scheme. The results show that the two-grid method can save the CPU time while keeping the same convergence rate.