High dimensional polynomial interpolation on sparse grids

We study polynomial interpolation on a d-dimensional cube, where d is large. We suggest to use the least solution at sparse grids with the extrema of the Chebyshev polynomials. The polynomial exactness of this method is almost optimal. Our error bounds show that the method is universal, i.e., almost optimal for many different function spaces. We report on numerical experiments for d = 10 using up to 652 065 interpolation points. This revised version was published online in June 2006 with corrections to the Cover Date.     

A convergent adaptive finite element method for an optimal design problem

The optimal design problem for maximal torsion stiffness of an infinite bar of given geometry and unknown distribution of two materials of prescribed amounts is one model example in topology optimisation. It eventually leads to a degenerate convex minimisation problem. The numerical analysis is therefore delicate for possibly multiple primal variables u but unique derivatives σ : = DW(D u). Even fine a posteriori error estimates still suffer from the reliability-efficiency gap. However, it motivates a simple edge-based adaptive mesh-refining algorithm (AFEM) that is not a priori guaranteed to refine everywhere. Its convergence proof is therefore based on energy estimates and some refined convexity control. Numerical experiments illustrate even nearly optimal convergence rates of the proposed AFEM. Supported by the DFG Research Center MATHEON “Mathematics for key technologies” in Berlin.     

A posteriori error estimates for elliptic problems with Dirac delta source terms

The aim of this paper is to introduce residual type a posteriori error estimators for a Poisson problem with a Dirac delta source term, in L ^p norm and W^1,p seminorm. The estimators are proved to yield global upper and local lower bounds for the corresponding norms of the error. They are used to guide adaptive procedures, which are experimentally shown to lead to optimal orders of convergence.     

Quadratic spline collocation methods for elliptic partial differential equations

We consider Quadratic Spline Collocation (QSC) methods for linear second order elliptic Partial Differential Equations (PDEs). The standard formulation of these methods leads to non-optimal approximations. In order to derive optimal QSC approximations, high order perturbations of the PDE problem are generated. These perturbations can be applied either to the PDE problem operators or to the right sides, thus leading to two different formulations of optimal QSC methods. The convergence properties of the QSC methods are studied. OptimalO(h ^3–j ) global error estimates for thejth partial derivative are obtained for a certain class of problems. Moreover,O(h ^4–j ) error bounds for thejth partial derivative are obtained at certain sets of points. Results from numerical experiments verify the theoretical behaviour of the QSC methods. Performance results also show that the QSC methods are very effective from the computational point of view. They have been implemented efficiently on parallel machines.This research was supported in part by David Ross Foundation (U.S.A) and NSERC (Natural Sciences and Engineering Research Council of Canada).     

Error estimates for the numerical solution of a time-periodic linear parabolic problem

In this paper we consider the numerical solution of a time-periodic linear parabolic problem. We derive optimal order error estimates inL ₂() for approximate solutions obtained by discretizing in space by a Galerkin finite-element method and in time by single-step finite difference methods, using known estimates for the associated initial value problem. We generalize this approach and obtain error estimates for more general discretization methods in the norm of a Banach spaceB L ₂(), e.g.,B=H ₀ ¹ () orL (). Finally, we consider some computational aspects and give a numerical example.     

Cubature formulas for function spaces with moderate smoothness

We construct simple algorithms for high-dimensional numerical integration of function classes with moderate smoothness. These classes consist of square-integrable functions over the d-dimensional unit cube whose coefficients with respect to certain multiwavelet expansions decay rapidly. Such a class contains discontinuous functions on the one hand and, for the right choice of parameters, the quite natural d-fold tensor product of a Sobolev space H^s[0,1] on the other hand.The algorithms are based on one-dimensional quadrature rules appropriate for the integration of the particular wavelets under consideration and on Smolyak's construction. We provide upper bounds for the worst-case error of our cubature rule in terms of the number of function calls. We additionally prove lower bounds showing that our method is optimal in dimension d=1 and almost optimal (up to logarithmic factors) in higher dimensions. We perform numerical tests which allow the comparison with other cubature methods.     

Interior penalty method for the indefinite time-harmonic Maxwell equations

Summary In this paper, we introduce and analyze the interior penalty discontinuous Galerkin method for the numerical discretization of the indefinite time-harmonic Maxwell equations in the high-frequency regime. Based on suitable duality arguments, we derive a-priori error bounds in the energy norm and the L²-norm. In particular, the error in the energy norm is shown to converge with the optimal order (h^min{s,}) with respect to the mesh size h, the polynomial degree , and the regularity exponent s of the analytical solution. Under additional regularity assumptions, the L²-error is shown to converge with the optimal order (h⁺¹). The theoretical results are confirmed in a series of numerical experiments.Supported by the EPSRC (Grant GR/R76615).Supported by the Swiss National Science Foundation under project 21-068126.02.Supported in part by the Natural Sciences and Engineering Council of Canada.     

Efficient computable error bounds for discontinuous Galerkin approximations of elliptic problems

We present guaranteed and computable both sided error bounds for the discontinuous Galerkin (DG) approximations of elliptic problems. These estimates are derived in the full DG-norm on purely functional grounds by the analysis of the respective differential problem, and thus, are applicable to any qualified DG approximation. Based on the triangle inequality, the underlying approach has the following steps for a given DG approximation: (1) computing a conforming approximation in the energy space using the Oswald interpolation operator, and (2) application of the existing functional a posteriori error estimates to the conforming approximation. Various numerical examples with varying difficulty in computing the error bounds, from simple problems of polynomial-type analytic solution to problems with analytic solution having sharp peaks, or problems with jumps in the coefficients of the partial differential equation operator, are presented which confirm the efficiency and the robustness of the estimates.     

Interpolation error estimates in W 1,p for degenerate Q 1 isoparametric elements

Optimal order error estimates in H ¹, for the Q ₁ isoparametric interpolation were obtained in Acosta and Durán (SIAM J Numer Anal37, 18–36, 1999) for a very general class of degenerate convex quadrilateral elements. In this work we show that the same conlusions are valid in W ^1,p for 1≤ p < 3 and we give a counterexample for the case p ≥ 3, showing that the result cannot be generalized for more regular functions. Despite this fact, we show that optimal order error estimates are valid for any p ≥ 1, keeping the interior angles of the element bounded away from 0 and π, independently of the aspect ratio. We also show that the restriction on the maximum angle is sharp for p ≥ 3.     

On conjugate gradient type methods and polynomial preconditioners for a class of complex non-hermitian matrices

Summary We consider conjugate gradient type methods for the solution of large linear systemsA x=b with complex coefficient matrices of the typeA=T+iI whereT is Hermitian and a real scalar. Three different conjugate gradient type approaches with iterates defined by a minimal residual property, a Galerkin type condition, and an Euclidean error minimization, respectively, are investigated. In particular, we propose numerically stable implementations based on the ideas behind Paige and Saunder's SYMMLQ and MINRES for real symmetric matrices and derive error bounds for all three methods. It is shown how the special shift structure ofA can be preserved by using polynomial preconditioning, and results on the optimal choice of the polynomial preconditioner are given. Also, we report on some numerical experiments for matrices arising from finite difference approximations to the complex Helmholtz equation.This work was supported in part by Cooperatives Agreement NCC 2-387 between the National Aeronautics and Space Administration (NASA) and the Universities Space Research Association (USRA) and by National Science Foundation Grant DCR-8412314     

<Emphasis Type="Italic">hp</Emphasis>-Discontinuous Galerkin methods for strongly nonlinear elliptic boundary value problems

In this paper, we have analyzed a one parameter family of hp-discontinuous Galerkin methods for strongly nonlinear elliptic boundary value problems with Dirichlet boundary conditions. These methods depend on the values of the parameter , where θ = + 1 corresponds to the nonsymmetric and θ = −1 corresponds to the symmetric interior penalty methods when and f(u,∇u) = −f, that is, for the Poisson problem. The error estimate in the broken H ¹ norm, which is optimal in h (mesh size) and suboptimal in p (degree of approximation) is derived using piecewise polynomials of degree p ≥ 2, when the solution . In the case of linear elliptic problems also, this estimate is optimal in h and suboptimal in p. Further, optimal error estimate in the L ² norm when θ = −1 is derived. Numerical experiments are presented to illustrate the theoretical results. Supported by DST-DAAD (PPP-05) project.     

A unified approach to a posteriori error estimation using element residual methods

Summary This paper deals with the problem of obtaining numerical estimates of the accuracy of approximations to solutions of elliptic partial differential equations. It is shown that, by solving appropriate local residual type problems, one can obtain upper bounds on the error in the energy norm. Moreover, in the special case of adaptiveh-p finite element analysis, the estimator will also give a realistic estimate of the error. A key feature of this is the development of a systematic approach to the determination of boundary conditions for the local problems. The work extends and combines several existing methods to the case of fullh-p finite element approximation on possibly irregular meshes with, elements of non-uniform degree. As a special case, the analysis proves a conjecture made by Bank and Weiser [Some A Posteriori Error Estimators for Elliptic Partial Differential Equations, Math. Comput.44, 283–301 (1985)].     

On the construction of stable curvilinear pp version elements for mixed formulations of elasticity and Stokes flow

Summary. The use of mixed finite element methods is well-established in the numerical approximation of the problem of nearly incompressible elasticity, and its limit, Stokes flow. The question of stability over curved elements for such methods is of particular significance in the p version, where, since the element size remains fixed, exact representation of the curved boundary by (large) elements is often used. We identify a mixed element which we show to be optimally stable in both p and h refinement over curvilinear meshes. We prove optimal p version (up to ) and h version (p = 2, 3) convergence for our element, and illustrate its optimality through numerical experiments. Received August 25, 1998 / Revised version received February 16, 1999 / Published online April 20, 2000 –? Springer-Verlag 2000     

Approximation of a bending plate problem with a boundary unilateral constraint

Summary We study a variational formulation of the unilaterally supported bent plate problem and we analyze the approximation of the problem by a mixed finite element method. We proveO(h) andO(h|lnh|^1/2) error bounds respectively for the moments and the displacement.Work partially supported by M.P.I., by G.N.I.M. of C.N.R. and by I.A.N. of C.N.R. of Pavia     

On additive Schwarz preconditioners for sparse grid discretizations

Summary Based on the framework of subspace splitting and the additive Schwarz scheme, we give bounds for the condition number of multilevel preconditioners for sparse grid discretizations of elliptic model problems. For a BXP-like preconditioner we derive an estimate of the optimal orderO(1) and for a HB-like variant we obtain an estimate of the orderO(k ² ·2 ^k/2 ), wherek denotes the number of levels employed. Furthermore, we confirm these results by numerically computed condition numbers.     

A posteriori error estimation of residual type for anisotropic diffusion-convection-reaction problems

This paper presents a robust a posteriori residual error estimator for diffusion-convection-reaction problems with anisotropic diffusion, approximated by a SUPG finite element method on isotropic or anisotropic meshes in R^d, d=2 or 3. The equivalence between the energy norm of the error and the residual error estimator is proved. Numerical tests confirm the theoretical results.     

Biquadratic finite volume element methods based on optimal stress points for parabolic problems

In this paper, the semi-discrete and full discrete biquadratic finite volume element schemes based on optimal stress points for a class of parabolic problems are presented. Optimal order error estimates in H¹ and L² norms are derived. In addition, the superconvergences of numerical gradients at optimal stress points are also discussed. A numerical experiment confirms some results of theoretical analysis.     

Robustness in a posteriori error analysis for FEM flow models

Summary. We derive a residual-based a posteriori error estimator for a stabilized finite element discretization of certain incompressible Oseen-like equations. We focus our attention on the behaviour of the effectivity index and we carry on a numerical study of its sensitiveness to the problem and mesh parameters. We also consider a scalar reaction-convection-diffusion problem and a divergence-free projection problem in order to investigate the effects on the robustness of our a posteriori error estimator of the reaction-convection-diffusion phenomena and, separately, of the incompressibility constraint. Received March 21, 2001 / Revised version received July 30, 2001 / Published online October 17, 2001     

A <Emphasis Type="Italic">V</Emphasis>-cycle Multigrid for multilevel matrix algebras: proof of optimality

We analyze the convergence rate of a multigrid method for multilevel linear systems whose coefficient matrices are generated by a real and nonnegative multivariate polynomial f and belong to multilevel matrix algebras like circulant, tau, Hartley, or are of Toeplitz type. In the case of matrix algebra linear systems, we prove that the convergence rate is independent of the system dimension even in presence of asymptotical ill-conditioning (this happens iff f takes the zero value). More precisely, if the d-level coefficient matrix has partial dimension n _r at level r, with , then the size of the system is , , and O(N(n)) operations are required by the considered V-cycle Multigrid in order to compute the solution within a fixed accuracy. Since the total arithmetic cost is asymptotically equivalent to the one of a matrix-vector product, the proposed method is optimal. Some numerical experiments concerning linear systems arising in 2D and 3D applications are considered and discussed.     

An asymptotically optimal Schur complement reduction for the Stokes equation

Summary. In this paper we develop an efficient Schur complement method for solving the 2D Stokes equation. As a basic algorithm, we apply a decomposition approach with respect to the trace of the pressure. The alternative stream function-vorticity reduction is also discussed. The original problem is reduced to solving the equivalent boundary (interface) equation with symmetric and positive definite operator in the appropriate trace space. We apply a mixed finite element approximation to the interface operator by iso triangular elements and prove the optimal error estimates in the presence of stabilizing bubble functions. The norm equivalences for the corresponding discrete operators are established. Then we propose an asymptotically optimal compression technique for the related stiffness matrix (in the absence of bubble functions) providing a sparse factorized approximation to the Schur complement. In this case, the algorithm is shown to have an optimal complexity of the order , q = 2 or q = 3, depending on the geometry, where N is the number of degrees of freedom on the interface. In the presence of bubble functions, our method has the complexity arithmetical operations. The Schur complement interface equation is resolved by the PCG iterations with an optimal preconditioner. Received March 20, 1996 / Revised version received October 28, 1997