A posteriori error estimates for discontinuous Galerkin methods of obstacle problems

We present a posteriori error analysis of discontinuous Galerkin methods for solving the obstacle problem, which is a representative elliptic variational inequality of the first kind. We derive reliable error estimators of the residual type. Efficiency of the estimators is theoretically explored and numerically confirmed.     

A posteriori estimates for the Bubble Stabilized Discontinuous Galerkin Method

In this paper, two reliable and efficient a posteriori error estimators for the Bubble Stabilized Discontinuous Galerkin (BSDG) method for diffusion-reaction problems in two and three dimensions are derived. The theory is followed by some numerical illustrations.     

一种推广的对流扩散方程的局部化间断Galerkin方法

研究求解一种产生于径向渗流问题的推广的对流扩散方程的局部化间断Galerkin方法,对一般非线性情形证明了方法的L^2稳定性;对线性情形证明了,当方法取有限元空间为κ次多项式空间时,数值解逼近的L^∞（0,T;L^2)模的误差阶为κ。     

A posteriori error analysis of the discontinuous finite element methods for first order hyperbolic problems

In this paper, we consider the a posteriori error analysis of discontinuous Galerkin finite element methods for the steady and nonsteady first order hyperbolic problems with inflow boundary conditions. We establish several residual-based a posteriori error estimators which provide global upper bounds and a local lower bound on the error. Further, for nonsteady problem, we construct a fully discrete discontinuous finite element scheme and derive the a posteriori error estimators which yield global upper bound on the error in time and space. Our a posteriori error analysis is based on the mesh-dependent a priori estimates for the first order hyperbolic problems. These a posteriori error analysis results can be applied to develop the adaptive discontinuous finite element methods.     

A Posteriori Error Estimation for the Dual Mixed Finite Element Method of the Stokes Problem

The paper presents an a posteriori error estimator of residual type for the stationary Stokes problem in a possibly multiply-connected polygonal domain of the plane, using the dual mixed finite element method (FEM). We prove lower and upper error bounds with the explicit dependence of the viscosity parameter and without any regularity assumption on the solution.     

A discontinuous Galerkin method for the Cahn-Hilliard equation

The Cahn-Hilliard equation is modeled to describe the dynamics of phase separation in glass and polymer systems. A priori error estimates for the Cahn-Hilliard equation have been studied by the authors. In order to control accuracy of approximate solutions, a posteriori error estimation of the Cahn-Hilliard equation is obtained by discontinuous Galerkin method.     

Simple error estimators for the Galerkin BEM for some hypersingular integral equation in 2D

A posteriori error estimation is an important tool for reliable and efficient Galerkin boundary element computations. For hypersingular integral equations in 2D with a positive-order Sobolev space, we analyse the mathematical relation between the (h???h/2)-error estimator from [S. Ferraz-Leite and D. Praetorius, Simple a posteriori error estimators for the h-version of the boundary element method, Computing 83 (2008), pp. 135–162], the two-level error estimator from [M. Maischak, P. Mund, and E. Stephan, Adaptive multilevel BEM for acoustic scattering, 585 Comput. Methods Appl. Mech. Eng. 150 (1997), pp. 351–367], and the averaging error estimator from [C. Carstensen and D. Praetorius, Averaging techniques for the a posteriori bem error control for a hypersingular integral equation in two dimensions, SIAM J. Sci. Comput. 29 (2007), pp. 782–810]. All of these a posteriori error estimators are simple in the following sense: first, the numerical analysis can be done within the same mathematical framework, namely localization techniques for the energy norm. Second, there is almost no implementational overhead for the realization.     

Tikhonov regularization method for a backward problem for the time-fractional diffusion equation

This paper is devoted to solve a backward problem for a time-fractional diffusion equation with variable coefficients in a general bounded domain by the Tikhonov regularization method. Based on the eigenfunction expansion of the solution, the backward problem for searching the initial data is changed to solve a Fredholm integral equation of the first kind. The conditional stability for the backward problem is obtained. We use the Tikhonov regularization method to deal with the integral equation and obtain the series expression of solution. Furthermore, the convergence rates for the Tikhonov regularized solution can be proved by using an a priori regularization parameter choice rule and an a posteriori regularization parameter choice rule. Two numerical examples in one-dimensional and two-dimensional cases respectively are investigated. Numerical results show that the proposed method is effective and stable.     

A second-order splitting combined with orthogonal cubic spline collocation method for the Rosenau equation

A second-order splitting method is applied to a KdV-like Rosenau equation in one space variable. Then an orthogonal cubic spline collocation procedure is employed to approximate the resulting system. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index 1. Error estimates in L² and L^∞ norms have been obtained for the semidiscrete approximations. For the temporal discretization, the time integrator RADAU5 is used for the resulting system. Some numerical experiments have been conducted to validate the theoretical results and to confirm the qualitative behaviors of the Rosenau equation. Finally, orthogonal cubic spline collocation method is directly applied to BBM (Benjamin–Bona–Mahony) and BBMB (Benjamin–Bona–Mahony–Burgers) equations and the well-known decay estimates are demonstrated for the computed solution. © 1998 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 14: 695–716, 1998     

A combined mixed finite element method and local discontinuous Galerkin method for miscible displacement problem in porous media

A combined method consisting of the mixed finite element method for flow and the local discontinuous Galerkin method for transport is introduced for the one-dimensional coupled system of incompressible miscible displacement problem. Optimal error estimates in L∞(0,T;L2) for concentration c,in L2(0,T;L2)for cxand L∞(0,T;L2) for velocity u are derived. The main technical difficulties in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method,the nonlinearity,and the coupling of the models. Numerical experiments are performed to verify the theoretical results. Finally,we apply this method to the one-dimensional compressible miscible displacement problem and give the numerical experiments to confirm the efficiency of the scheme.     

A posteriori error estimates of discontinuous Galerkin methods for non-standard Volterra integro-differential equations

^** Email: jingtang{at}lsec.cc.ac.cn^*** Email: hermann{at}math.mun.ca In this paper we establish a posteriori error estimates forthe discontinuous Galerkin (DG) method applied to linear, semilinearand non-standard (non-linear) Volterra integro-differentialequations. We also present an analysis of the DG method withquadrature for the memory term. Numerical experiments basedon three integro-differential equations are used to illustratevarious aspects of the error analysis.     

Decay and scattering of small solutions for Rosenau equations

We study the long-time behavior of small solutions of the Cauchy problem for a Rosenau equation. For a class of nonlinearity of the perturbation, the global small solution was obtained, and the decay and scattering for small amplitude solution are established.     

Adaptive meshless Galerkin boundary node methods for hypersingular integral equations

Adaptive refinement techniques are developed in this paper for the meshless Galerkin boundary node method for hypersingular boundary integral equations. Two types of error estimators are derived. One is a perturbation error estimator that is formulated based on the difference between numerical solutions obtained using two consecutive nodal arrangements. The other is a projection error estimator that is formulated based on the difference between the numerical solution itself and its projection. These error estimators are proven to have an upper and a lower bound by the constant multiples of the exact error in the energy norm. A localization scheme is presented to accomodate the non-local property of hypersingular integral operators for the needed computable local error indicators. The convergence of the adaptive meshless techniques is verified theoretically. To confirm the theoretical results and to show the efficiency of the adaptive techniques, numerical examples in 2D and 3D with high singularities are provided.     

Error estimates for discontinuous Galerkin method for nonlinear parabolic equations

We consider the nonlinear parabolic partial differential equations. We construct a discontinuous Galerkin approximation using a penalty term and obtain an optimal L^∞(L²) error estimate.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

A POSTERIORI ERROR ESTIMATE OF OPTIMAL CONTROL PROBLEM OF PDE WITH INTEGRAL CONSTRAINT FOR STATE

In this paper, we study adaptive finite element discretization schemes for an optimal control problem governed by elliptic PDE with an integral constraint for the state. We derive the equivalent a posteriori error estimator for the finite element approximation, which particularly suits adaptive multi-meshes to capture different singularities of the control and the state. Numerical examples are presented to demonstrate the efficiency of a posteriori error estimator and to confirm the theoretical results.     

A discontinuous functional differential equation

In this paper we give a sufficient condition for the existence of maximal and minimal solutions to a discontinuous functional differential equation x′(t)=f(t,x(t),x(·)), x(0)=0. We apply the result to establish an existence theorem for the Darboux problem for a partial differential equation.     

Equilibrated error estimators for discontinuous Galerkin methods

We consider some diffusion problems in domains of ?^d, d = 2 or 3 approximated by a discontinuous Galerkin method with polynomials of any degree. We propose a new a posteriori error estimator based on H(div)‐conforming elements. It is shown that this estimator gives rise to an upper bound where the constant is one up to higher order terms. The lower bound is also established with a constant depending on the aspect ratio of the mesh, the dependence with respect to the coefficients being also traced. The reliability and efficiency of the proposed estimator is confirmed by some numerical tests. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Interior penalty discontinuous Galerkin method for Maxwell's equations: Energy norm error estimates

We develop the symmetric interior penalty discontinuous Galerkin (DG) method for the time-dependent Maxwell equations in second-order form. We derive optimal a priori error estimates in the energy norm for smooth solutions. We also consider the case of low-regularity solutions that have singularities in space.