双参数任意四边形非协调板元的构造及收敛性分析

本文利用双参数有限元方法[1]构造出十二参和十三参任意四边形板元，并对其收敛性进行分析证明．     

A new robust C-type nonconforming triangular element for singular perturbation problems

In this paper, a new robust C⁰ triangular element is proposed for the fourth order elliptic singular perturbation problem with double set parameter method and bubble function technique, and a general convergence theorem for C⁰ nonconforming elements is presented. The convergence of the new element is proved in the energy norm uniformly with respect to the perturbation parameter. Numerical experiments are also carried out to demonstrate the efficiency of the new element.     

Convergence and stability of finite element nonlinear Galerkin method for the Navier-Stokes equations

In this article, we present a new fully discrete finite element nonlinear Galerkin method, which are well suited to the long time integration of the Navier-Stokes equations. Spatial discretization is based on two-grid finite element technique; time discretization is based on Euler explicit scheme with variable time step size. Moreover, we analyse the boundedness, convergence and stability condition of the finite element nonlinear Galerkin method. Our discussion shows that the time step constraints of the method depend only on the coarse grid parameter and the time step constraints of the finite element Galerkin method depend on the fine grid parameter under the same convergence accuracy. Received February 2, 1994 / Revised version received December 6, 1996     

Composite penalty method of a low order anisotropic nonconforming finite element for the Stokes problem

Composite penalty method of a low order anisotropic nonconforming quadrilateral finite element for the Stokes problem is presented. This method with a large penalty parameter can achieve the same accuracy as the stand method with a small penalty parameter and the convergence rate of this method is two times as that of the standard method under the condition of the same order penalty parameter. The superconvergence for velocity is established as well. The results of this paper are also valid to the most of the known nonconforming finite element methods.     

双参数Morley有限元多重网格法

1、引言 多重网格方法是求解偏微分方程的高效快速算法,在实际中得到广泛应用．[2][6]中考察了Morley元的多重网格方法,并用于双调和方程问题。     

九参数广义协调元的收敛性

众知周所,解板弯曲问题的Zienkiewicz不协调三次元只对特殊的单元剖分才收敛.但由于这种元采用单元顶点的函数值及二个一阶导数值作为节点参数,计算简单,总体自由度少,所以相继出现一些对Zienkiewicz元的改进形式,使之对任意剖分均收敛,如拟协调元,TRUNC元,Specht元.对这些元的分析见[7—9].最近龙     

Error analysis of finite element approximations of the inverse mean curvature flow arising from the general relativity

This paper proposes and analyzes a finite element method for a nonlinear singular elliptic equation arising from the black hole theory in the general relativity. The nonlinear equation, which was derived and analyzed by Huisken and Ilmanen in (J Diff Geom 59:353–437), represents a level set formulation for the inverse mean curvature flow describing the evolution of a hypersurface whose normal velocity equals the reciprocal of its mean curvature. We first propose a finite element method for a regularized flow which involves a small parameter ɛ; a rigorous analysis is presented to study well-posedness and convergence of the scheme under certain mesh-constraints, and optimal rates of convergence are verified. We then prove uniform convergence of the finite element solution to the unique weak solution of the nonlinear singular elliptic equation as the mesh size h and the regularization parameter ɛ both tend to zero. Computational results are provided to show the efficiency of the proposed finite element method and to numerically validate the “jumping out” phenomenon of the weak solution of the inverse mean curvature flow. Numerical studies are presented to evidence the existence of a polynomial scaling law between the mesh size h and the regularization parameter ɛ for optimal convergence of the proposed scheme. Finally, a numerical convergence study for another approach recently proposed by R. Moser (The inverse mean curvature flow and p-harmonic functions. preprint U Bath, 2005) for approximating the inverse mean curvature flow via p-harmonic functions is also included.     

Weak Galerkin finite element methods for the biharmonic equation on polytopal meshes

A new weak Galerkin (WG) finite element method is introduced and analyzed in this article for the biharmonic equation in its primary form. This method is highly robust and flexible in the element construction by using discontinuous piecewise polynomials on general finite element partitions consisting of polygons or polyhedra of arbitrary shape. The resulting WG finite element formulation is symmetric, positive definite, and parameter‐free. Optimal order error estimates in a discrete H² norm is established for the corresponding WG finite element solutions. Error estimates in the usual L² norm are also derived, yielding a suboptimal order of convergence for the lowest order element and an optimal order of convergence for all high order of elements. Numerical results are presented to confirm the theory of convergence under suitable regularity assumptions. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1003–1029, 2014     

NONCONFORMING FINITE ELEMENT PENALTY METHOD FOR STOKES EQUATION

a special penalty method is presented to improve the accuracy of the standard penaltymethod (or solving Stokes equation with nonconforming finite element, It is shown that thismethod with a larger penalty parameter can achieve the same accuracy as the staodaxd methodwith a smaller penalty parameter. The convergence rate of the standard method is just hall order of this penalty method when using the same penalty parameter, while the extrapolationmethod proposed by Faik et al can not yield so high accuracy of convergence. At last, we alsoget the super-convergence estimates for total flux.     

双参数十二参矩形板元的对称列式

1 引言 在位移有限元中,九参数三角形板元的研究取得了丰硕成果,根据不同方法已构造出众多收敛性能很好的单元（见[1]、[2]、[3]）。相比之下,矩形板元的研究却较少报道,ACM元及广义协调元RGC—12是其中比较成功的单元.但是ACM是C~0元。其位移形函数的外法向导数平均值在单元间不连续。广义协调元是基于势能原理建立单元协调的,其自由度(协调条件)不对称是其本身的一个弱点,陈万吉研究表明。这种不对称性会破坏单元的几何不变性。     

SOR-like Methods for Augmented Systems

Several SOR-like methods are proposed for solving augmented systems. These have many different applications in scientific computing, for example, constrained optimization and the finite element method for solving the Stokes equation. The convergence and the choice of optimal parameter for these algorithms are studied. The convergence and divergence regions for some algorithms are given, and the new algorithms are applied to solve the Stokes equations as well.     

一阶双曲问题的间断流线扩散法

1．引言众所周知，求解一阶双曲问题的Galerkin有限元法，仅具有次最优LZ一收敛阶估计，难于建立H＼误差估计【‘1，且Gderkill有限元解常呈现伪数值振荡．为改善计算精度与稳定性，诸多非标准有限元解法相继提出，其中，间断（Discontinuous）Galerkin有限元法（以下简称DG方法）与流线扩散（StreamlineDiffusion）有限元法（以下简称SD法）是两种具有鲜明特点，较为成功的算法．具体地，DG方法是一种迎风型显式算法，它从入流边界开始，沿流场方向，自上游往下游，逐个单元进行解算，计算十分简便且可局部并行化．SD方法则是一种P…     

TRAPEZOIDAL PLATE BENDING ELEMENT WITH DOUBLE SET PARAMETERS

Using double set parameter method, a 12-parameter trapezoidal plate bending element is presented. The first set of degrees of freedom, which make the element convergent, are the values at the four vertices and the middle points of the four sides together with the mean values of the outer normal derivatives along four sides. The second set of degree of freedom, which make the number of unknowns in the resulting discrete system small and computation convenient are values and the first derivatives at the four vertices of the element. The convergence of the element is proved.     

Continuous-discontinuous finite element method for convection-diffusion problems with characteristic layers

We study convergence properties of a numerical method for convection-diffusion problems with characteristic layers on a layer-adapted mesh. The method couples standard Galerkin with an h-version of the nonsymmetric discontinuous Galerkin finite element method with bilinear elements. In an associated norm, we derive the error estimate as well as the supercloseness result that are uniform in the perturbation parameter. Applying a post-processing operator for the discontinuous Galerkin method, we construct a new numerical solution with enhanced convergence properties.     

Numerical modelling of semi-coercive beam problem with unilateral elastic subsoil of Winkler’s type

A non-linear semi-coercive beam problem is solved in this article. Suitable numerical methods are presented and their uniform convergence properties with respect to the finite element discretization parameter are proved here. The methods are based on the minimization of the total energy functional, where the descent directions of the functional are searched by solving the linear problems with a beam on bilateral elastic “springs”. The influence of external loads on the convergence properties is also investigated. The effectiveness of the algorithms is illustrated on numerical examples.     

Sufficient conditions for uniform convergence on layer-adapted meshes for one-dimensional reaction–diffusion problems

We study convergence properties of a finite element method with lumping for the solution of linear one-dimensional reaction–diffusion problems on arbitrary meshes. We derive conditions that are sufficient for convergence in the L_∞ norm, uniformly in the diffusion parameter, of the method. These conditions are easy to check and enable one to immediately deduce the rate of convergence. The key ingredients of our analysis are sharp estimates for the discrete Green function associated with the discretization. AMS subject classification 65L10, 65L12, 65L15     

A relaxation procedure for domain decomposition methods using finite elements

Summary We present the convergence analysis of a new domain decomposition technique for finite element approximations. This technique was introduced in [11] and is based on an iterative procedure among subdomains in which transmission conditions at interfaces are taken into account partly in one subdomain and partly in its adjacent. No global preconditioner is needed in the practice, but simply single-domain finite element solvers are required. An optimal strategy for an automatic selection of a relaxation parameter to be used at interface subdomains is indicated. Applications are given to both elliptic equations and incompressible Stokes equations.     

A LOCKING-FREE ANISOTROPIC NONCONFORMING FINITE ELEMENT FOR PLANAR LINEAR ELASTICITY PROBLEM

The main aim of this article is to study the approximation of a locking-free anisotropic nonconforming finite element for the pure displacement boundary value prob-lem of planar linear elasticity. The optimal error estimates are obtained by using some novel approaches and techniques. The method proposed in this article is robust in the sense that the convergence estimates in the energy and L2-norms are independent of the Lame parameterλ.     

Local projection stabilization of equal order interpolation applied to the Stokes problem

The local projection stabilization allows us to circumvent the Babuška-Brezzi condition and to use equal order interpolation for discretizing the Stokes problem. The projection is usually done in a two-level approach by projecting the pressure gradient onto a discontinuous finite element space living on a patch of elements. We propose a new local projection stabilization method based on (possibly) enriched finite element spaces and discontinuous projection spaces defined on the same mesh. Optimal order of convergence is shown for pairs of approximation and projection spaces satisfying a certain inf-sup condition. Examples are enriched simplicial finite elements and standard quadrilateral/hexahedral elements. The new approach overcomes the problem of an increasing discretization stencil and, thus, is simple to implement in existing computer codes. Numerical tests confirm the theoretical convergence results which are robust with respect to the user-chosen stabilization parameter.

     

A direct discontinuous Galerkin finite element method for convection-dominated two-point boundary value problems

The direct discontinuous Galerkin (DDG) finite element method, using piecewise polynomials of degree k ≥ 1 on a Shishkin mesh, is applied to convection-dominated singularly perturbed two-point boundary value problems. Consistency, stability and convergence of order k (up to a logarithmic factor) are proved in an energy-type norm appropriate to the method and problem. The results are robust, i.e., they hold uniformly for all values of the singular perturbation parameter. Numerical experiments confirm the theoretical convergence rate.