Discontinuous Galerkin finite element methods for variational inequalities of first and second kinds

We develop the error analysis for the h‐version of the discontinuous Galerkin finite element discretization for variational inequalities of first and second kinds. We establish an a priori error estimate for the method which is of optimal order in a mesh dependant as well as L²‐norm.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Convergence analysis of the LDG method applied to singularly perturbed problems

Considering a two‐dimensional singularly perturbed convection–diffusion problem with exponential boundary layers, we analyze the local discontinuous Galerkin (DG) method that uses piecewise bilinear polynomials on Shishkin mesh. A convergence rate O(N^‐1 lnN) in a DG‐norm is established under the regularity assumptions, while the total number of mesh points is O(N²). The rate of convergence is uniformly valid with respect to the singular perturbation parameter ε. Numerical experiments indicate that the theoretical error estimate is sharp. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Superconvergence analysis of the Galerkin FEM for a singularly perturbed convection–diffusion problem with characteristic layers

We analyze the superconvergence property of the Galerkin finite element method (FEM) for elliptic convection–diffusion problems with characteristic layers. This method on Shishkin meshes is known to be almost first‐order accurate (up to a logarithmic factor) in the energy norm induced by the bilinear form of the weak formulation, uniformly in the perturbation parameter. In the present paper the method is shown to be almost second‐order superconvergent in this energy norm for the difference between the FEM solution and the bilinear interpolant of the exact solution. This supercloseness property is used to improve the accuracy to almost second order by means of a postprocessing procedure. Numerical experiments confirm these results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

A posteriori error estimation for convection dominated problems on anisotropic meshes

A singularly perturbed convection–diffusion problem in two and three space dimensions is discretized using the streamline upwind Petrov Galerkin (SUPG) variant of the finite element method. The dominant convection frequently gives rise to solutions with layers; hence anisotropic finite elements can be applied advantageously. The main focus is on a posteriori energy norm error estimation that is robust in the perturbation parameter and with respect to the mesh anisotropy. A residual error estimator and a local problem error estimator are proposed and investigated. The analysis reveals that the upper error bound depends on the alignment of the anisotropies of the mesh and of the solution. Hence reliable error estimation is possible for suitable anisotropic meshes. The lower error bound depends on the problem data via a local mesh Peclet number. Thus efficient error estimation is achieved for small mesh Peclet numbers. Altogether, error estimation approaches for isotropic meshes are successfully extended to anisotropic elements. Several numerical experiments support the analysis. Copyright © 2003 John Wiley & Sons, Ltd.     

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)     

Global error bounds for the Petrov–Galerkin discretization of the neutron transport equation

In this paper, we prove that the piecewise bilinear Petrov‐Galerkin discretization for the mono‐directional neutron transport equation described in (J. Comput. Phys. 1986; 64 :96–111) is convergent and second‐order accurate, provided that the true solution to the problem has continuous partial derivatives of all orders up through three. We do this by giving a bound on the 2‐norm of the inverse of the system matrix that is independent of the mesh size. This shows that the global error is of the same order as the local truncation error. Copyright © 2010 John Wiley & Sons, Ltd.     

The combination technique for a two-dimensional convection-diffusion problem with exponential layers

Convection-diffusion problems posed on the unit square and with solutions displaying exponential layers are solved using a sparse grid Galerkin finite element method with Shishkin meshes. Writing N for the maximum number of mesh intervals in each coordinate direction, our “combination” method simply adds or subtracts solutions that have been computed by the Galerkin FEM on and meshes. It is shown that the combination FEM yields (up to a factor ln N) the same order of accuracy in the associated energy norm as the Galerkin FEM on an N × N mesh, but it requires only (N ^3/2) degrees of freedom compared with the (N ²) used by the Galerkin FEM. An analogous result is also proved for the streamline diffusion finite element method. This work was supported by the National Natural Science Foundation of China (10701083 and 10425105), the Chinese National Basic Research Program (2005CB321704) and the Boole Centre for Research in Informatics at National University of Ireland Cork.     

Convergence of adaptive FEM for some elliptic obstacle problem

In this work, we treat the convergence of adaptive lowest-order FEM for some elliptic obstacle problem with affine obstacle. For error estimation, we use a residual error estimator from [D. Braess, C. Carstensen, and R. Hoppe, Convergence analysis of a conforming adaptive finite element method for an obstacle problem, Numer. Math. 107 (2007), pp. 455–471]. We extend recent ideas from [J. Cascon, C. Kreuzer, R. Nochetto, and K. Siebert, Quasi-optimal convergence rate for an adaptive finite element method, SIAM J. Numer. Anal. 46 (2008), pp. 2524–2550] for the unrestricted variational problem to overcome the lack of Galerkin orthogonality. The main result states that an appropriately weighted sum of energy error, edge residuals and data oscillations satisfies a contraction property within each step of the adaptive feedback loop. This result is superior to a prior result from Braess et al. (2007) in two ways: first, it is unnecessary to control the decay of the data oscillations explicitly; second, our analysis avoids the use of some discrete local efficiency estimate so that the local mesh-refinement is fairly arbitrary.     

Coupling of FEM and BEM for a nonlinear interface problem: The h–p version

This article presents some numerical examples for coupling the finite element method (FEM) and the boundary element method (BEM) as analyzed in [11]. This coupling procedure combines the advantages of boundary elements (problems in unbounded regions) and of finite elements (nonlinear problems with inhomogeneous data). In [28], experimental rates of convergence for the h version are presented, where the accuracy of the Galerkin approximation is achieved by refining the mesh. In this article we treat the h–p version, combining an increase of the degree of the piecewise polynomials with a certain mesh refinement. In our model examples, we obtain theoretically and numerically exponential convergence, which indicates a great efficiency in particular if singularities appear. © 1995 John Wiley & Sons, Inc.     

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.     

Sensitivities for topological graph changes in finite element meshes and application to adaptive refinement

We propose a novel approach to adaptivity in FEM based on local sensitivities for topological mesh changes. To this end, we consider refinement as a continuous operation on the edge graph of the finite element discretization, for instance by splitting nodes along edges and expanding edges to elements. Thereby, we introduce the concept of a topological mesh derivative for a given objective function that depends on the discrete solution of the underlying PDE. These sensitivities may in turn be used as refinement indicators within an adaptive algorithm. For their calculation, we rely on the first-order asymptotic expansion of the Galerkin solution with respect to the topological mesh change. As a proof of concept, we consider the total potential energy of a linear symmetric second-order elliptic PDE, minimization of which is known to decrease the approximation error in the energy norm. In this case, our approach yields local sensitivities that are closely related to the reduction of the energy error upon refinement and may therefore be used as refinement indicators in an adaptive algorithm. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

拟线性拋物方程半离散变网格有限元方法的误差估计

§1.导言近年来,变网格方法正日益为人们所重视与应用,但理论性分析文献仍不多见。文献[1]讨论了某些发展型方程变网格方法的误差估计,但未给出收敛阶估计;文献[2,3]仅对全离散方法讨论了收敛阶问题。本文对一类拟线性抛物问题,于第二节中给出了半离散Galerkin变网格计算格式及其可解性定理;第三节中建立了对称误差估计;第四节给     

Superconvergence Analysis of Galerkin FEM and SDFEM for Elliptic Problems with Characteristic Layers

We analyse the superconvergence properties of the Galerkin FEM and of the streamline-diffusion finite element method (SDFEM) using bilinear functions in the case of elliptic problems with characteristic layers. To resolve the layers we use appropriate Shishkin meshes. For the SDFEM we give an optimal choice for the streamline-diffusion parameter δ for maximal stability in the induced streamline-diffusion norm. In the characteristic layers we are able to show that δ can be chosen of order δ = Cε- N^–2 which is con.rmed by numerical results. (© 2006 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Discontinuous Stable Elements for the Incompressible Flow

In this paper, we derive a discontinuous Galerkin finite element formulation for the Stokes equations and a group of stable elements associated with the formulation. We prove that these elements satisfy the new inf–sup condition and can be used to solve incompressible flow problems. Associated with these stable elements, optimal error estimates for the approximation of both velocity and pressure in L ² norm are obtained for the Stokes problems, as well as an optimal error estimate for the approximation of velocity in a mesh dependent norm.     

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.

     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

Uniform superconvergence of a Galerkin finite element method on Shishkin‐type meshes

We consider a Galerkin finite element method that uses piecewise bilinears on a class of Shishkin‐type meshes for a model singularly perturbed convection‐diffusion problem on the unit square. The method is shown to be convergent, uniformly in the diffusion parameter ϵ, of almost second order in a discrete weighted energy norm. As a corollary, we derive global L₂‐norm error estimates and local L_∞‐norm estimates. Numerical experiments support our theoretical results. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16:426–440, 2000     

A TWO-LEVEL FINITE ELEMENT GALERKIN METHOD FOR THE NONSTATIONARY NAVIER-STOKES EQUATIONS I： SPATIAL DISCRETIZATION

In this article we consider a two-level finite element Galerkin method using mixed finite elements for the two-dimensional nonstationary incompressible Navier-Stokes equations. The method yields a $H^1$-optimal velocity approximation and a $L_2$-optimal pressure approximation. The two-level finite element Galerkin method involves solving one small, nonlinear Navier-Stokes problem on the coarse mesh with mesh size $H$, one linear Stokes problem on the fine mesh with mesh size $h << H$. The algorithm we study produces an approximate solution with the optimal, asymptotic in $h$, accuracy.     

Fully discrete interior penalty discontinuous Galerkin methods for nonlinear parabolic equations

In this article, we investigate interior penalty discontinuous Galerkin (IPDG) methods for solving a class of two‐dimensional nonlinear parabolic equations. For semi‐discrete IPDG schemes on a quasi‐uniform family of meshes, we obtain a priori bounds on solutions measured in the L² norm and in the broken Sobolev norm. The fully discrete IPDG schemes considered are based on the approximation by forward Euler difference in time and broken Sobolev space. Under a restriction related to the mesh size and time step, an hp ‐version of an a priori l^∞(L²) and l²(H¹) error estimate is derived and numerical experiments are presented.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 288–311, 2012     

Local discontinuous galerkin approximation of convection‐dominated diffusion optimal control problems with control constraints

In this article, we investigate local discontinuous Galerkin approximation of stationary convection‐dominated diffusion optimal control problems with distributed control constraints. The state variable and adjoint state variable are approximated by piecewise linear polynomials without continuity requirement, whereas the control variable is discretized by variational discretization concept. The discrete first‐order optimality condition is derived. We show that optimization and discretization are commutative for the local discontinuous Galerkin approximation. Because the solutions to convection‐dominated diffusion equations often admit interior or boundary layers, residual type a posteriori error estimate in L² norm is proved, which can be used to guide mesh refinement. Finally, numerical examples are presented to illustrate the theoretical findings. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 339–360, 2014