A Note on the Energy Norm for a Singularly Perturbed Model Problem

This paper considers a singularly perturbed reaction diffusion problem. It is investigated whether adaptive approaches are successful to design robust solution procedures. A key ingredient is the a posteriori error estimator. Since robust and mathematically analysed error estimation is possible in the energy norm, the focus is on this choice of norm and its implications. The numerical performance for several model problems confirms that the proposed adaptive algorithm (in conjunction with an energy norm error estimator) produces optimal results. Hence the energy norm is suitable for the purpose considered here. The investigations also provide valuable justification for forthcoming research. Received October 25, 2001; revised July 12, 2002 Published online: October 24, 2002     

ℋ-matrices for Convection-diffusion Problems with Constant Convection

L ^∞-coefficients. This paper analyses the application of hierarchical matrices to the convection-dominant convection-diffusion equation with constant convection. In the case of increasing convection, the convergence of a standard ℋ-matrix approximant towards the original matrix will deteriorate. We derive a modified partitioning and admissibility condition that ensures good convergence also for the singularly perturbed case. Received January 1, 2003; revised March 4, 2003 Published online: May 2, 2003     

A priori and a posteriori error analysis of an augmented mixed finite element method for incompressible fluid flows

In this paper we extend recent results on the a priori and a posteriori error analysis of an augmented mixed finite element method for the linear elasticity problem, to the case of incompressible fluid flows with symmetric stress tensor. Similarly as before, the present approach is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and from the relations defining the pressure in terms of the stress tensor and the rotation in terms of the displacement, all of them multiplied by stabilization parameters. We show that these parameters can be suitably chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well-posed for any choice of finite element subspaces. Next, we present a reliable and efficient residual-based a posteriori error estimator for the augmented mixed finite element scheme. Finally, several numerical results confirming the theoretical properties of this estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are reported.     

Reliable Approximation of Weight Factors Entering Residual-based Error Bounds for FE-discretisations

In this note we refine strategies of the so called dual-weighted-residual (DWR) approach to a posteriori error control for FE-schemes. We derive rigorous error bounds, especially we control the approximation process of the (unknown) dual solution entering the proposed estimate.     

Two-level methods for the single layer potential in ℝ3

We consider weakly singular integral equations of the first kind on open surface pieces Γ in ℝ³. To obtain approximate solutions we use theh-version Galerkin boundary element method. Furthermore we introduce two-level additive Schwarz operators for non-overlapping domain decompositions of Γ and we estimate the conditions numbers of these operators with respect to the mesh size. Based on these operators we derive an a posteriori error estimate for the difference between the exact solution and the Galerkin solution. The estimate also involves the error which comes from an approximate solution of the Galerkin equations. For uniform meshes and under the assumption of a saturation condition we show reliability and efficiency of our estimate. Based on this estimate we introduce an adaptive multilevel algorithm with easily computable local error indicators which allows direction control of the local refinements. The theoretical results are illustrated by numerical examples for plane and curved surfaces. Supported by the German Research Foundation (DFG) under grant Ste 238/25-9.     

Bounds for theN lowest eigenvalues of fourth-order boundary value problems

We describe a method for the calculation of theN lowest eigenvalues of fourth-order problems inH ₀ ² (Ω). In order to obtain small error bounds, we compute the defects inH ⁻²(Ω) and, to obtain a bound for the rest of the spectrum, we use a boundary homotopy method. As an example, we compute strict error bounds (using interval arithmetic to control rounding errors) for the 100 lowest eigenvalues of the clamped plate problem in the unit square. Applying symmetry properties, we prove the existence of double eigenvalues.     

Weighted Error Estimates for Finite Element Solutions of Variational Inequalities

In this note the studies begun in Blum and Suttmeier (1999) on adaptive finite element discretisations for nonlinear problems described by variational inequalities are continued. Similar to the concept proposed, e.g., in Becker and Rannacher (1996) for variational equalities, weighted a posteriori estimates for controlling arbitrary functionals of the discretisation error are constructed by using a duality argument. Numerical results for the obstacle problem demonstrate the derived error bounds to be reliable and, used for an adaptive grid refinement strategy, to produce economical meshes. Received September 6, 1999; revised February 8, 2000     

Error Estimates of a FEM with Lumping for Parabolic PDEs

L ²-norm. This well-known FEM is given by the use of the vertical line method and conforming linear finite elements on a triangulation. The main result of the paper are new estimates in the L ²-norm for the additional error term originated by lumping. Using these ones, for the FEM with lumping we can apply directly the proof technique of error estimates known for conforming FEMs. Received May 17, 2001; revised November 2, 2001 Published online February 18, 2002     

Superconvergence of a Mixed Covolume Method for Elliptic Problems

We consider a mixed covolume method for a system of first order partial differential equations resulting from the mixed formulation of a general self-adjoint elliptic problem with a variable full diffusion tensor. The system can be used to model the transport of a contaminant carried by a flow. We use the lowest order Raviart-Thomas mixed finite element space. We show the first order convergence in L ² norm and the superconvergence in certain discrete norms both for the pressure and velocity. Finally some numerical examples illustrating the error behavior of the scheme are provided. Supported by the National Natural Science Foundation of China under grant No. 10071044 and the Research Fund of Doctoral Program of High Education by State Education Ministry of China.     

An a Posteriori Error Estimator for Two-Body Contact Problems on Non-Matching Meshes

A posteriori error estimates for two-body contact problems are established. The discretization is based on mortar finite elements with dual Lagrange multipliers. To define locally the error estimator, Arnold–Winther elements for the stress and equilibrated fluxes for the surface traction are used. Using the Lagrange multiplier on the contact zone as Neumann boundary conditions, equilibrated fluxes can be locally computed. In terms of these fluxes, we define on each element a symmetric and globally H(div)-conforming approximation for the stress. Upper and lower bounds for the discretization error in the energy norm are provided. In contrast to many other approaches, the constant in the upper bound is, up to higher order terms, equal to one. Numerical examples illustrate the reliability and efficiency of the estimator. This work was supported in part by the Deutsche Forschungsgemeinschaft, SFB 404, B8.     

Multiresolution Visualization of Higher Order Adaptive Finite Element Simulations

h –p–adaptive projection with respect to any prescribed threshold value for the visual error. This projection can then be processed by various local rendering methods, e.g. color coding of data or isosurface extraction. Especially for color coding purposes modern texture capabilities are used to directly render higher polynomial data by superposition of polynomial basis function textures and final color look-up tables. Numerical experiments from CFD clearly demonstrate the applicability and efficiency of our approach. Received September 25, 2001; revised March 31, 2003 Published online: May 26, 2003 The authors acknowledge the valuable hints of the anonymous referees.     

A Posteriori Error Estimates for Weak Galerkin Finite Element Methods for Second Order Elliptic Problems

A residual type a posteriori error estimator is presented and analyzed for Weak Galerkin finite element methods for second order elliptic problems. The error estimator is proved to be efficient and reliable through two estimates, one from below and the other from above, in terms of an $H^1$ -equivalent norm for the exact error. Two numerical experiments are conducted to demonstrate the effectiveness of adaptive mesh refinement guided by this estimator.     

Verification of Reduced Convergence Rates

In this short article, we recalculate the numerical example in Kíek and Neittaanmäki (1987) for the Poisson solution u=x(1–x)siny in the unit square S as . By the finite difference method, an error analysis for such a problem is given from our previous study by where h is the meshspacing of the uniform square grids used, and C₁ and C₂ are two positive constants. Let =u–u_h, where u_h is the finite difference solution, and is the discrete H¹ norm. Several techniques are employed to confirm the reduced rate of convergence, and to give the constants, C₁=0.09034 and C₂=0.002275 for a stripe domain. The better performance for arises from the fact that the constant C₁ is much large than C₂, and the h in computation is not small enough.     

A Posteriori Error Estimates and an Adaptive Finite Element Method for the Allen–Cahn Equation and the Mean Curvature Flow

This paper develops an a posteriori error estimate of residual type for finite element approximations of the Allen–Cahn equation u_t − Δu+ ε⁻² f(u)=0. It is shown that the error depends on ε⁻¹ only in some low polynomial order, instead of exponential order. Based on the proposed a posteriori error estimator, we construct an adaptive algorithm for computing the Allen–Cahn equation and its sharp interface limit, the mean curvature flow. Numerical experiments are also presented to show the robustness and effectiveness of the proposed error estimator and the adaptive algorithm.     

Iterative Stokes Solvers in the Harmonic Velte Subspace

We explore the prospects of utilizing the decomposition of the function space (H ¹ ₀) ⁿ (where n=2,3) into three orthogonal subspaces (as introduced by Velte) for the iterative solution of the Stokes problem. It is shown that Uzawa and Arrow-Hurwitz iterations – after at most two initial steps – can proceed fully in the third, smallest subspace. For both methods, we also compute optimal iteration parameters. Here, for two-dimensional problems, the lower estimate of the inf-sup constant by Horgan and Payne proves useful and provides an inclusion of the spectrum of the Schur complement operator of the Stokes problem. We further consider the conjugate gradient method in the third Velte subspace and derive a corresponding convergence estimate. Computational results show the effectiveness of this approach for discretizations which admit a discrete Velte decomposition. Received June 11, 1999; revised October 27, 2000     

Error trends in Quasi-Monte Carlo integration

Several test functions, whose variation could be calculated, were integrated with up to 10¹⁰ trials using different low-discrepancy sequences in dimensions 3, 6, 12, and 24. The integration errors divided by the variation of the functions were compared with exact and asymptotic discrepancies. These errors follow an approximate power law, whose constant is essentially given by the variance of the integrand, and whose power depends on its effective dimension. Included were also some calculations with scrambled Niederreiter sequences, and with Niederreiter-Xing sequences. A notable result is that the pre-factors of the asymptotic discrepancy function D∗(N), which are so often used as a quality measure of the sequences, have little relevance in practical ranges of N.     

Finite Volume Element Approximations of Nonlocal in Time One-Dimensional Flows in Porous Media

Various finite volume element schemes for parabolic integro-differential equations in 1-D are derived and studied. These types of equations arise in modeling reactive flows or material with memory effects. Our main goal is to develop a general framework for obtaining finite volume element approximations and to study the error analysis. We consider the lowest-order (linear and L-splines) finite volume elements, although higher-order volume elements can be considered as well under this framework. It is proved that finite volume element approximations are convergent with optimal order in H ¹-norms, suboptimal order in the L ²-norm and super-convergent order in a discrete H ¹-norm. Received August 3, 1998; revised October 11, 1999     

An A Posteriori Error Estimate for the Local Discontinuous Galerkin Method Applied to Linear and Nonlinear Diffusion Problems

In this paper we present a new residual-based reliable a posteriori error estimator for the local discontinuous Galerkin approximations of linear and nonlinear diffusion problems in polygonal regions of R ². Our analysis, which applies to convex and nonconvex domains, is based on Helmholtz decompositions of the error and a suitable auxiliary polynomial function interpolating the Dirichlet datum. Several examples confirming the reliability of the estimator and providing numerical evidences for its efficiency are given. Furthermore, the associated adaptive method, which considers meshes with and without hanging nodes, is shown to be much more efficient than a uniform refinement to compute the discrete solutions. In particular, the experiments illustrate the ability of the adaptive algorithm to localize the singularities of each problem.Mathematics Subject Classifications (1991). 65N30This revised version was published online in July 2005 with corrected volume and issue numbers.     

On the Convergence of a Multigrid Method for Linear Reaction-Diffusion Problems

In this note we consider discrete linear reaction-diffusion problems. For the discretization a standard conforming finite element method is used. For the approximate solution of the resulting discrete problem a multigrid method with a damped Jacobi or symmetric Gauss-Seidel smoother is applied. We analyze the convergence of the multigrid V- and W-cycle in the framework of the approximation- and smoothing property. The multigrid method is shown to be robust in the sense that the contraction number can be bounded by a constant smaller than one which does not depend on the mesh size or on the diffusion-reaction ratio. Received June 15, 2000     

<Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-interpolation on quadrilateral and hexahedral meshes with hanging nodes

We propose a Scott-Zhang type finite element interpolation operator of first order for the approximation of H ¹-functions by means of continuous piecewise mapped bilinear or trilinear polynomials. The novelty of the proposed interpolation operator is that it is defined for general non-affine equivalent quadrilateral and hexahedral elements and so-called 1-irregular meshes with hanging nodes. We prove optimal local approximation properties of this interpolation operator for functions in H ¹. As necessary ingredients we provide a definition of a hanging node and a rigorous analysis of the issue of constrained approximation which cover both the two- and three-dimensional case in a unified fashion.