Analysis of a FEM for a coupled system of singularly perturbed reaction–diffusion equations

We consider a system of ℓ ≥ 2 one-dimensional singularly perturbed reaction–diffusion equations coupled at the zero-order term. The second derivative of each equation is multiplied by a distinct small parameter. We present a convergence theory for conforming linear finite elements on arbitrary meshes. As a result convergence independently of the perturbation parameters on a wide class of layer-adapted meshes is established.      

ROBUSTNESS OF AN UPWIND FINITE DIFFERENCE SCHEME FOR SEMILINEAR CONVECTION-DIFFUSION PROBLEMS WITH BOUNDARY TURNING POINTS

We consider a singularly perturbed semilinear convection-diffusion problem with a boundary layer of attractive turning-point type. It is shown that its solution can be decomposed into a regular solution component and a layer component. This decomposi-tion is used to analyse the convergence of an upwinded finite difference scheme on Shishkin meshes.     

MULTIGRID METHOD FOR FLUID DYNAMIC PROBLEMS

This paper covers the dynamics problems. The review and some aspects of main development stages of using Multigrid method for fluid multigrid technics are presented. Some approaches for solving Navier-Stokes equations and convection- diffusion problems are considered.     

DISCONTINUOUS FINITE ELEMENT METHOD FOR CONVECTION-DIFFUSION EQUATIONS

1. IlltroductionThe finite element approximation of the convection--diffusin equations has been investigated using several different approaches (see e.g. [3] [4] and the references therein).Previous analysis in primal formulation of these problems was done for two types ofapproximation schemes: one which produces a continuous piecewise polynomial approximation and one which produces a piecewise polynomial approximation which arecontinuous for certain number of moments accross interelement edge…     

Graded Meshes for Higher Order FEM

A singularly perturbed one-dimensional convection-diffusion problem is solved numerically by the finite element method based on higher order polynomials. Numerical solutions are obtained using S-type meshes with special emphasis on meshes which are graded (based on a mesh generating function) in the fine mesh region. Error estimates in the ε-weighted energy norm are proved. We derive an 'optimal' mesh generating function in order to minimize the constant in the error estimate. Two layer-adapted meshes defined by a recursive formulae in the fine mesh region are also considered and a new technique for proving error estimates for these meshes is presented. The aim of the paper is to emphasize the importance of using optimal meshes for higher order finite element methods. Numerical experiments support all theoretical results.     

分层网格上奇异摄动问题的一致NIPG分析

本文采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在分层网格上至多相差一个关于摄动参数对数因子的拟最优阶的一致收敛性,即在能量范数度量下其误差估计为O((log~2(1/e))/N),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.     

UNIFORM SUPERCONVERGENCE OF A FINITE ELEMENT METHOD WITH EDGE STABILIZATION FOR CONVECTION-DIFFUSION PROBLEMS

In the present paper the edge stabilization technique is applied to a convection-diffusion problem with exponential boundary layers on the unit square, using a Shishkin mesh with bilinear finite elements in the layer regions and linear elements on the coarse part of the mesh. An error bound is proved for ‖πu-u^h‖Е, where πu is some interpolant of the solution u and uh the discrete solution. This supercloseness result implies an optimal error estimate with respect to the L2 norm and opens the door to the application of postprocessing for improving the discrete solution.     

对流扩散方程一类改进的特征线修正有限元方法

１引言在地下水污染,地下渗流驱动,核污染,半导体等问题的数值模拟中,均涉及抛物型对流扩散方程（或方程组）的数值求解问题．这些对流扩散型偏微分方程（或方程组）具有共同的特点：对流的影响远大于扩散的影响,即对流占优性,对流占优性给问题的数值求解带来许多困难,因此对流占优问题的有效数值解法一直是计算数学中重要的研究内容．用通常的差分法或有限元法进行数值求解将出现数值振荡．为了克服数值振荡,提出各种迎风方法和修正的特征方法并在这些问题上得到成功的实际应用、８０年代,Ｄｏｕｇｌａｓ和Ｒｕｓｓｅｌｌ［２］等…     

Hamilton-Jacobi方程的单调有限元格式

A monotone finite element scheme is obtained by applying the finite element method to the viscosity equation of the Hamilton-Jacobi equation on unstructured meshes. Under some constraints, we show that this scheme is monotone and its numerical solution converges to the viscosity solution of the Hamilton-Jacobi equa-tion. Numerical examples test the stability and the convergence of this scheme.     

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 modified Bakhvalov mesh

We consider a few numerical methods for solving a one-dimensional convection–diffusion singularly perturbed problem. More precisely, we introduce a modified Bakvalov mesh generated using some implicitly defined functions. Properties of this mesh and convergence results for several methods on it are given. Numerical results are presented in support of the theoretical considerations.     

SUPERCONVERGENCE OF GRADIENT RECOVERY SCHEMES ON GRADED MESHES FOR CORNER SINGULARITIES

For the linear finite element solution to the Poisson equation, we show that supercon- vergence exists for a type of graded meshes for corner singularities in polygonal domains. In particular, we prove that the L^2-projection from the piecewise constant field △↓UN to the continuous and piecewise linear finite element space gives a better approximation of △↓U in the Hi-norm. In contrast to the existing superconvergence results, we do not assume high regularity of the exact solution.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

INFINITE ELEMENT METHOD FOR PROBLEMS ON UNBOUNDED AND MULTIPLY CONNECTED DOMAINS

1 IntroductionDifferent kinds of numerical n1etl1ods llavc been apPlied to so1lle proble1l1s o11 exteriordolllai11s successf1lly, e.g. tlle bou11dary element metl1od, the absorbillg boundary condi-tion lnethod, the spectrunl 11letllod, a11d the i11fillite eleIue1lt 1nethod. Tlle il1finite elementllletllod llas beell applied to tl1e Laplace equatioll['], the Stokes equation[']['], the plane elastic-ity systeln['1, alld the Heln1ho1tz equatioll[n. We study tl1e infinite ele111ent llletl1od tbr…     

RATE OF CONVERGENCE OF SCHWARZ ALTERNATING METHOD FOR TIME-DEPENDENT CONVECTION-DIFFUSION PROBLEM

This paper provides a theoretical justification to a overlapping domain decomposition method applied to the solution of time-dependent convection-diffusion problems. The method is based on the partial upwind finite element scheme and the discrete strong maximum principle for steady problem. An error estimate in L∞(0, T; L∞ (Ω)) is obtained and the fact that convergence factor ρ(τ, h) →0 exponentially as τ, h→0 is also proved under some usual conditions.     

MULTIGRID METHODS FOR MORLEY ELEMENT ON NONNESTED MESHES

1.IntroductionWeconsidersomemultigridalgorithmsforthebiharmonicequationdiscretizedbyMoneyelementonnonnestedmeshes.TOdefineamultigridalgorithm,certainintergridtransferoperatorhastobeconstructed.Throughtakingtheaveragesofthenodalvariables,weconstructanintergridtransferoperatorforMoneyelementonnonnestedmeshesthatsatisfiesacertainstableapproximationpropertywhichplaysakeyroleinmultigridmethodsfornonconformingplateelementsonnonnestedmeshes.Theso--calledregularity-approximaticnassurnptionisestablis…     

A SUPERONVERGENCE ANALYISIS FOR FINITE ELEMENT SOLUTION BY THE INTERPOLANT POSTPROCESSING ON IRREGULAR MESHES FOR SMOOTH PROBLEM

1. IntroductionThe results in this paper are based on the idea of interpolation postprocessing in [11 andthe techniques of LZ projection processing in [2]For simlicity) we consider the model problem: Finds e Hi(fl),such thatSuppose that Jh and JH are irregular triangulations (or quadrilateral partitions). Theirsizes satisfy h << H, (H - 0). Construct piecewise k-order and r-order finite element spaceSh and SH respectively. Let ah E Sh be the Galerkin approximation of u E HJ(fl), andbe …