On the finite volume element method

Summary The finite volume element method (FVE) is a discretization technique for partial differential equations. It uses a volume integral formulation of the problem with a finite partitioning set of volumes to discretize the equations, then restricts the admissible functions to a finite element space to discretize the solution. this paper develops discretization error estimates for general selfadjoint elliptic boundary value problems with FVE based on triangulations with linear finite element spaces and a general type of control volume. We establishO(h) estimates of the error in a discreteH ¹ semi-norm. Under an additional assumption of local uniformity of the triangulation the estimate is improved toO(h ²). Results on the effects of numerical integration are also included.This research was sponsored in part by the Air Force Office of Scientific Research under grant number AFOSR-86-0126 and the National Science Foundation under grant number DMS-8704169. This work was performed while the author was at the University of Colorado at Denver     

Higher order fully discrete scheme combined withH 1-Galerkin mixed finite element method for semilinear reaction-diffusion equations

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by anH ¹-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) ofindex one. Apriori error estimates for semidiscrete scheme are derived for both differential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.     

Two-grid methods for –P1 mixed finite element approximation of general elliptic optimal control problems with low regularity

In this paper, we present a two-grid mixed finite element scheme for distributed optimal control governed by general elliptic equations. –P₁ mixed finite elements are used for the discretization of the state and co-state variables, whereas piecewise constant function is used to approximate the control variable. We first use a new approach to obtain the superclose property between the centroid interpolation and the numerical solution of the optimal control u with order h² under the low regularity. Based on the superclose property, we derive the optimal a priori error estimates. Then, using a postprocessing projection operator, we get a second-order superconvergent result for the control u. Next, we construct a two-grid mixed finite element scheme and analyze a priori error estimates. In the two-grid scheme, the solution of the elliptic optimal control problem on a fine grid is reduced to the solution of the elliptic optimal control problem on a much coarser grid and the solution of a linear algebraic system on the fine grid and the resulting solution still maintains an asymptotically optimal accuracy. Finally, a numerical example is presented to verify the theoretical results.     

A posteriori error estimates for mixed finite element approximations of parabolic problems

We derive residual based a posteriori error estimates for parabolic problems on mixed form solved using Raviart–Thomas–Nedelec finite elements in space and backward Euler in time. The error norm considered is the flux part of the energy, i.e. weighted L ²(Ω) norm integrated over time. In order to get an optimal order bound, an elementwise computable post-processed approximation of the scalar variable needs to be used. This is a common technique used for elliptic problems. The final bound consists of terms, capturing the spatial discretization error and the time discretization error and can be used to drive an adaptive algorithm.     

Asymptotic expansions andL ∞-error estimates for mixed finite element methods for second order elliptic problems

Summary Asymptotic expansions for mixed finite element approximations of the second order elliptic problem are derived and Richardson extrapolation can be applied to increase the accuracy of the approximations. A new procedure, which is called the error corrected method, is presented as a further application of the asymptotic error expansion for the first order BDM approximation of the scalar field. The key point in deriving the asymptotic expansions for the error is an establishment ofL ¹-error estimates for mixed finite element approximations for the regularized Green's functions. As another application of theL ¹-error estimates for the regularized Green's functions, we shall present maximum norm error estimates for mixed finite element methods for second order elliptic problems.     

On mixed finite element methods for linear elastodynamics

Summary We construct and analyze finite element methods for approximating the equations of linear elastodynamics, using mixed elements for the discretization of the spatial variables. We consider two different mixed formulations for the problem and analyze semidiscrete and up to fourth-order in time fully discrete approximations.L ² optimal-order error estimates are proved for the approximations of displacement and stress.Work supported in part by the Hellenic State Scholarship Foundation     

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)     

A posteriori error analysis of nonconforming finite volume elements for general second‐order elliptic PDEs

In this article, we study the a posteriori H¹ and L² error estimates for Crouzeix‐Raviart nonconforming finite volume element discretization of general second‐order elliptic problems in ?². The error estimators yield global upper and local lower bounds. Finally, numerical experiments are performed to illustrate the theoretical findings. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Mixed Discontinuous Galerkin Time-Stepping Method for Linear Parabolic Optimal Control Problems

In this paper, we discuss the mixed discontinuous Galerkin (DG) finite element approximation to linear parabolic optimal control problems. For the state variables and the co-state variables, the discontinuous finite element method is used for the time discretization and the Raviart-Thomas mixed finite element method is used for the space discretization. We do not discretize the space of admissible control but implicitly utilize the relation between co-state and control for the discretization of the control. We derive a priori error estimates for the lowest order mixed DG finite element approximation. Moveover, for the element of arbitrary order in space and time, we derive a posteriori $L^2(0, T ;L^2(Ω))$ error estimates for the scalar functions, assuming that only the underlying mesh is static. Finally, we present an example to confirm the theoretical result on a priori error estimates.     

Analysis of the Fictitious Domain Method with an L 2-Penalty for Elliptic Problems

The L ²-penalty fictitious domain method is based on a reformulation of the original problem in a larger simple-shaped domain by introducing a discontinuous reaction term with a penalty parameter ε > 0. We first derive regularity results and some a priori estimates and then prove several error estimates. We also give several error estimates for discretization problems by the finite element and finite volume methods.     

Convergence of finite element method for linear second‐order wave equations with discontinuous coefficients

A finite element method is proposed and analyzed for hyperbolic problems with discontinuous coefficients. The main emphasize is given on the convergence of such method. Due to low global regularity of the solutions, the error analysis of the standard finite element method is difficult to adopt for such problems. For a practical finite element discretization, optimal error estimates in L^∞(L²) and L^∞(H¹) norms are established for continuous time discretization. Further, a fully discrete scheme based on a symmetric difference approximation is considered, and optimal order convergence in L^∞(H¹) norm is established. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Error estimates of the finite element method for the exterior Helmholtz problem with a modified DtN boundary condition

A priori error estimates in the H¹- and L²-norms are established for the finite element method applied to the exterior Helmholtz problem, with modified Dirichlet-to-Neumann (MDtN) boundary condition. The error estimates include the effect of truncation of the MDtN boundary condition as well as that of discretization of the finite element method. The error estimate in the L²-norm is sharper than that obtained by the author [D. Koyama, Error estimates of the DtN finite element method for the exterior Helmholtz problem, J. Comput. Appl. Math. 200 (1) (2007) 21-31] for the truncated DtN boundary condition.     

On the Convergence of Finite Element Method for Second Order Elliptic Interface Problems

The purpose of this paper is to study the convergence of finite element approximation to the exact solution of general self-adjoint elliptic equations with discontinuous coefficients. Due to low global regularity of the solution, it is difficult to achieve optimal order of convergence with classical finite element methods [Numer. Math. 1998; 79:175–202]. In this paper, an isoparametric type of discretization is used to prove optimal order error estimates in L ² and H ¹ norms when the global regularity of the solution is low. The interface is assumed to be of arbitrary shape and is smooth for our purpose. Further, for the purpose of numerical computations, we discuss the effect of numerical quadrature on finite element solution, and the related optimal order estimates are also established.     

Mortar element methods for parabolic problems

In this article a standard mortar finite element method and a mortar element method with Lagrange multiplier are used for spatial discretization of a class of parabolic initial‐boundary value problems. Optimal error estimates in L^∞(L²) and L^∞(H¹)‐norms for semidiscrete methods for both the cases are established. The key feature that we have adopted here is to introduce a modified elliptic projection. In the standard mortar element method, a completely discrete scheme using backward Euler scheme is discussed and optimal error estimates are derived. The results of numerical experiments support the theoretical results obtained in this article. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Error Estimation for Discontinuous Galerkin Solutions of Two-Dimensional Hyperbolic Problems

We analyze the discretization errors of discontinuous Galerkin solutions of steady two-dimensional hyperbolic conservation laws on unstructured meshes. We show that the leading term of the error on each element is a linear combination of orthogonal polynomials of degrees p and p+1. We further show that there is a strong superconvergence property at the outflow edge(s) of each element where the average discretization error converges as O(h ^2p+1) compared to a global rate of O(h ^p+1). Our analyses apply to both linear and nonlinear conservation laws with smooth solutions. We show how to use our theory to construct efficient and asymptotically exact a posteriori discretization error estimates and we apply these to some examples.     

A posteriori error estimation with finite element methods of lines for one-dimensional parabolic systems

Summary Consider the solution of one-dimensional linear initial-boundary value problems by a finite element method of lines using a piecewiseP ^th -degree polynomial basis. A posteriori estimates of the discretization error are obtained as the solutions of either local parabolic or local elliptic finite element problems using piecewise polynomial corrections of degreep+1 that vanish at element ends. Error estimates computed in this manner are shown to converge in energy under mesh refinement to the exact finite element discretization error. Computational results indicate that the error estimates are robust over a wide range of mesh spacings and polynomial degrees and are, furthermore, applicable in situations that are not supported by the analysis.This research was partially supported by the U.S. Air Force Office of Scientific Research, Air Force Systems Command, USAF, under Grant Number AFOSR 90-0194; by the U.S. Army Research Office under Contract Number DAAL03-91-G-0215; and by the National Science Foundation under Institutional Infrastructure Grant Number CDA-8805910     

Guaranteed velocity error control for the pseudostress approximation of the Stokes equations

The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L². Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016     

Improved L2 and H1 error estimates for the Hessian discretization method

The Hessian discretization method (HDM) for fourth-order linear elliptic equations provides a unified convergence analysis framework based on three properties namely coercivity, consistency, and limit-conformity. Some examples that fit in this approach include conforming and nonconforming finite element methods (ncFEMs), finite volume methods (FVMs) and methods based on gradient recovery operators. A generic error estimate has been established in L², H¹, and H²-like norms in literature. In this paper, we establish improved L² and H¹ error estimates in the framework of HDM and illustrate it on various schemes. Since an improved L² estimate is not expected in general for FVM, a modified FVM is designed by changing the quadrature of the source term and a superconvergence result is proved for this modified FVM. In addition to the Adini ncFEM, in this paper, we show that the Morley ncFEM is an example of HDM. Numerical results that justify the theoretical results are also presented.     

Local mesh refinement for the discretization of Neumann boundary control problems on polyhedra

This paper deals with optimal control problems constrained by linear elliptic partial differential equations. The case where the right‐hand side of the Neumann boundary is controlled, is studied. The variational discretization concept for these problems is applied, and discretization error estimates are derived. On polyhedral domains, one has to deal with edge and corner singularities, which reduce the convergence rate of the discrete solutions, that is, one cannot expect convergence order two for linear finite elements on quasi‐uniform meshes in general. As a remedy, a local mesh refinement strategy is presented, and a priori bounds for the refinement parameters are derived such that convergence with optimal rate is guaranteed. As a by‐product, finite element error estimates in the H¹(Ω)‐norm, L²(Ω)‐norm and L²(Γ)‐norm for the boundary value problem are obtained, where the latter one turned out to be the main challenge. Copyright © 2015 John Wiley & Sons, Ltd.     

Low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions

In this paper, we consider low‐order stabilized finite element methods for the unsteady Stokes/Navier‐Stokes equations with friction boundary conditions. The time discretization is based on the Euler implicit scheme, and the spatial discretization is based on the low‐order element (P₁−P₁ or P₁−P₀) for the approximation of the velocity and pressure. Moreover, some error estimates for the numerical solution of fully discrete stabilized finite element scheme are obtained. Finally, numerical experiments are performed to confirm our theoretical results.