hp–Adaptive composite discontinuous Galerkin methods for elliptic problems on complicated domains

In this article, we develop the a posteriori error estimation of hp–version discontinuous Galerkin composite finite element methods for the discretization of second‐order elliptic partial differential equations. This class of methods allows for the approximation of problems posed on computational domains which may contain a huge number of local geometrical features, or microstructures. Although standard numerical methods can be devised for such problems, the computational effort may be extremely high, as the minimal number of elements needed to represent the underlying domain can be very large. In contrast, the minimal dimension of the underlying composite finite element space is independent of the number of geometric features. Computable bounds on the error measured in terms of a natural (mesh‐dependent) energy norm are derived. Numerical experiments highlighting the practical application of the proposed estimators within an automatic hp–adaptive refinement procedure will be presented. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1342–1367, 2014     

Exponential convergence of mixed hp-DGFEM for Stokes flow in polygons

Summary. We analyze mixed hp-discontinuous Galerkin finite element methods (DGFEM) for Stokes flow in polygonal domains. In conjunction with geometrically refined quadrilateral meshes and linearly increasing approximation orders, we prove that the hp-DGFEM leads to exponential rates of convergence for piecewise analytic solutions exhibiting singularities near corners. Mathematics Subject Classification (2000):65N30     

Exponential convergence in a Galerkin least squares hp-FEM for Stokes flow

A stabilized hp-finite element method (FEM) of Galerkin leastsquares (GLS) type is analysed for the Stokes equations in polygonaldomains. Contrary to the standard Galerkin FEM, the GLSFEM admitsthe implementationally attractive equal-order interpolationin the velocity and the pressure. In conjunction with geometricallyrefined meshes and linearly increasing approximation ordersit is shown that the hp-GLSFEM leads to exponential rates ofconvergence for solutions exhibiting singularities near corners.To obtain this result a novel hp-interpolation result is provedthat allows the approximation of pressure functions in certainweighted Sobolev spaces in a conforming way and at exponentialrates of convergence on geometric meshes. Received 6 June 1999. Accepted 14 March 2000.     

Analysis of a semidiscrete discontinuous Galerkin scheme for compressible miscible displacement problem

Numerical approximation of the coupled system of compressible miscible displacement problem in porous media is considered in this paper. A continuous in time discontinuous Galerkin scheme is developed. The symmetric interior penalty discontinuous Galerkin method is used to solve both the flow and transport equations. Upwind technique is used to treat the convection term in the transport equation. The hp-a priori error bounds are derived.     

hp-Adaptive Finite Element Methods for Variational Inequalities

In this work, we combine an hp–adaptive strategy with a posteriori error estimates for variational inequalities, which are given by contact problems. The a posteriori error estimates are obtained using a general approach based on the saddle point formulation of contact problems and making use of a yposteriori error estimates for variational equations. Error estimates are presented for obstacle problems and Signorini problems with friction. Numerical experiments confirm the reliability of the error estimates for finite elements of higher order. The use of the hp–adaptive strategy leads to meshes with the same characteristics as geometric meshes and to exponential convergence. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two-Grid hp-Version DGFEMs for Strongly Monotone Second-Order Quasilinear Elliptic PDEs

In this article we develop the a priori error analysis of so-called two-grid hp-version discontinuous Galerkin finite element methods for the numerical approximation of strongly monotone second-order quasilinear partial differential equations. In this setting, the fully nonlinear problem is first approximated on a coarse finite element space V(𝒯_H, P ). The resulting ‘coarse’ numerical solution is then exploited to provide the necessary data needed to linearize the underlying discretization on the finer space V(𝒯_h, p ); thereby, only a linear system of equations is solved on the richer space V(𝒯_h, p ). Numerical experiments confirming the theoretical results are presented. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A hp-discontinuous Galerkin method for the time-dependent Maxwell’s equation: a priori error estimate

A discontinuous Galerkin method for the numerical approximation for the time-dependent Maxwell’s equations in “stable medium” with supraconductive boundary, is introduced and analysed. its hp-analysis is carried out and error estimates that are optimal in the meshsize h and slightly suboptimal in the approximation degree p are obtained.     

On the p-version approximation in the boundary element method for a variational inequality of the second kind modelling unilateral contact and given friction

In this paper we complement recent work of Maischak and Stephan on adaptive hp-versions of the BEM for unilateral Signorini problems, respectively on FEM-BEM coupling in its h-version for a nonlinear transmission problem modelling Coulomb friction contact. Here we focus on the boundary element method in its p-version to treat a scalar variational inequality of the second kind that models unilateral contact and Coulomb friction in elasticity together. This leads to a nonconforming discretization scheme. In contrast to the work cited above and to a related paper of Guediri on a boundary variational inequality of the second kind modelling friction we take the quadrature error of the friction functional into account of the error analysis. At first without any regularity assumptions, we prove convergence of the BEM Galerkin approximation in the energy norm. Then under mild regularity assumptions, we establish an a priori error estimate that is based on a novel Céa–Falk lemma for abstract variational inequalities of the second kind.     

A priori error analysis of a discontinuous Galerkin approximation for a kind of compressible miscible displacement problems

A kind of compressible miscible displacement problems which include molecular diffusion and dispersion in porous media are investigated.A symmetric interior penalty discontinuous Galerkin (SIPG) method is applied to the coupled system of flow and transport.Using the induction hypotheses instead of the cut-off operator and the interpolation projection properties,a priori hp error estimates are presented.The error bounds in L2(H1) norm for concentration and in L∞(L2) norm for velocity are optimal in h and suboptimal in p with a loss of power 1/2.     

Benchmark results for testing adaptive finite element eigenvalue procedures

A discontinuous Galerkin method, with hp-adaptivity based on the approximate solution of appropriate dual problems, is employed for highly-accurate eigenvalue computations on a collection of benchmark examples. After demonstrating the effectivity of our computed error estimates on a few well-studied examples, we present results for several examples in which the coefficients of the partial-differential operators are discontinuous. The problems considered here are put forward as benchmarks upon which other adaptive methods for computing eigenvalues may be tested, with results compared to our own.     

The hp-MITC finite element method for the Reissner–Mindlin plate problem

The popular MITC finite elements used for the approximation of the Reissner–Mindlin plate are extended to the case where elements of non-uniform degree p distribution are used on locally refined meshes. Such an extension is of particular interest to the hp-version and hp-adaptive finite element methods. A priori error bounds are provided showing that the method is locking-free. The analysis is based on new approximation theoretic results for non-uniform Brezzi–Douglas–Fortin–Marini spaces, and extends the results obtained in the case of uniform order approximation on globally quasi-uniform meshes presented by Stenberg and Suri (SIAM J. Numer. Anal. 34 (1997) 544). Numerical examples illustrating the theoretical results and comparing the performance with alternative standard Galerkin approaches are presented for two new benchmark problems with known analytic solution, including the case where the shear stress exhibits a boundary layer. The new method is observed to be locking-free and able to provide exponential rates of convergence even in the presence of boundary layers.     

Mixed FEM of higher-order for a frictional contact problem

This paper presents mixed finite element methods of higher-order for an idealized frictional contact problem in linear elasticity. The approach relies on a saddle point formulation where the frictional contact condition is captured by a Lagrange multiplier. The convergence of the mixed scheme is proven and some a priori estimates for the h- and p-method are derived. Furthermore, a posteriori error estimates are presented which rely on the estimation of the discretization error of an auxiliary problem and some further terms capturing the error in the friction and complementary conditions. Numerical results confirm the applicability of the a posteriori error estimates within h- and hp-adaptive schemes. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

DG and hp-DG for highly indefinite Helmholtz problems

We develop a stability and convergence theory for a Discontinuous Galerkin formulation (DG) of a highly indefinite Helmholtz problem in . We prove that the DG-method admits a unique solution under much weaker conditions than for conventional Galerkin methods. It is shown that for the case of hp-DGFEM the optimal convergence order estimate can be obtained under the conditions that is sufficiently small and the polynomial degree p is at least O(log k). (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Adaptive higher-order space-time discontinuous Galerkin method for the computer simulation of variably-saturated porous media flows

This paper is concerned with the numerical simulation of time-dependent variably-saturated Darcian flow problems described by the Richards equation. We present the adaptive higher-order space-time discontinuous Galerkin (hp-STDG) method which optimizes accuracy and efficiency by balancing the errors that arise from the space and time discretizations and from the resulting nonlinear algebraic system. Convergence problems related to the transition between unsaturated flow and saturated flow are eliminated by regularizing the constitutive formulas. We also present an hp-anisotropic mesh adaptation technique capable of generating unstructured triangular elements with optimal sizes, shapes, and polynomial approximation degrees. Several numerical experiments are presented to demonstrate the accuracy, efficiency, and robustness of the numerical method presented here.     

Interpolation error-based a posteriori error estimation for two-point boundary value problems and parabolic equations in one space dimension

Summary. I derive a posteriori error estimates for two-point boundary value problems and parabolic equations in one dimension based on interpolation error estimates. The interpolation error estimates are obtained from an extension of the error formula for the Lagrange interpolating polynomial in the case of symmetrically-spaced interpolation points. From this formula pointwise and seminorm a priori estimates of the interpolation error are derived. The interpolant in conjunction with the a priori estimates is used to obtain asymptotically exact a posteriori error estimates of the interpolation error. These a posteriori error estimates are extended to linear two-point boundary problems and parabolic equations. Computational results demonstrate the convergence of a posteriori error estimates and their effectiveness when combined with an hp-adaptive code for solving parabolic systems. Received April 17, 2000 / Revised version received September 25, 2000 / Published online May 30, 2001     

A posteriori error control for finite element approximations of elliptic eigenvalue problems

We develop a new approach to a posteriori error estimation for Galerkin finite element approximations of symmetric and nonsymmetric elliptic eigenvalue problems. The idea is to embed the eigenvalue approximation into the general framework of Galerkin methods for nonlinear variational equations. In this context residual-based a posteriori error representations are available with explicitly given remainder terms. The careful evaluation of these error representations for the concrete situation of an eigenvalue problem results in a posteriori error estimates for the approximations of eigenvalues as well as eigenfunctions. These suggest local error indicators that are used in the mesh refinement process.     

hp‐FEM for second moments of elliptic PDEs with stochastic data. II: Exponential convergence for stationary singular covariance functions

We prove exponential rates of convergence of a class of hp Galerkin Finite Element approximations of solutions to a model tensor nonhypoelliptic equation in the unit square □ = (0, 1)² which exhibit singularities on ?□ and on the diagonal Δ = {(x, y) ∈ □ : x = y}, but are otherwise analytic in □. As we explained in the first part (Pentenrieder and Schwab, Research Report, Seminar for Applied Mathematics, 2010) of this work, such problems arise as deterministic second moment equations of linear, second order elliptic operator equations Au = f with Gaussian random field data f. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2012     

A posteriori error estimates of hp spectral element methods for optimal control problems with L2-norm state constraint

In this paper, we investigate a distributed optimal control problem governed by elliptic partial differential equations with L²-norm constraint on the state variable. Firstly, the control problem is approximated by hp spectral element methods, which combines the advantages of the finite element methods with spectral methods; then, the optimality conditions of continuous system and discrete system are presented, respectively. Next, hp a posteriori error estimates are derived for the coupled state and control approximation. In the end, a projection gradient iterative algorithm is given, which solves the optimal control problems efficiently. Numerical experiments are carried out to confirm that the numerical results are in good agreement with the theoretical results.

     

Optimal error estimates for the hp-version interior penalty discontinuous Galerkin finite element method

We consider the hp-version interior penalty discontinuous Galerkinfinite-element method (hp-DGFEM) for second-order linear reaction–diffusionequations. To the best of our knowledge, the sharpest knownerror bounds for the hp-DGFEM are due to Rivière et al.(1999,Comput. Geosci., 3, 337–360) and Houston et al.(2002,SIAM J. Numer. Anal., 99, 2133–2163). These are optimalwith respect to the meshsize h but suboptimal with respect tothe polynomial degree p by half an order of p. We present improvederror bounds in the energy norm, by introducing a new functionspace framework. More specifically, assuming that the solutionsbelong element-wise to an augmented Sobolev space, we deducefully hp-optimal error bounds.     

Finite element theory on curved domains with applications to discontinuous Galerkin finite element methods

In this paper we provide key estimates used in the stability and error analysis of discontinuous Galerkin finite element methods (DGFEMs) on domains with curved boundaries. In particular, we review trace estimates, inverse estimates, discrete Poincaré–Friedrichs' inequalities, and optimal interpolation estimates in noninteger Hilbert–Sobolev norms, that are well known in the case of polytopal domains. We also prove curvature bounds for curved simplices, which does not seem to be present in the existing literature, even in the polytopal setting, since polytopal domains have piecewise zero curvature. We demonstrate the value of these estimates, by analyzing the IPDG method for the Poisson problem, introduced by Douglas and Dupont, and by analyzing a variant of the hp-DGFEM for the biharmonic problem introduced by Mozolevski and Süli. In both cases we prove stability estimates and optimal a priori error estimates. Numerical results are provided, validating the proven error estimates.