Optimal a priori error estimates for the -version of the local discontinuous Galerkin method for convection-diffusion problems

We study the convergence properties of the -version of the local discontinuous Galerkin finite element method for convection-diffusion problems; we consider a model problem in a one-dimensional space domain. We allow arbitrary meshes and polynomial degree distributions and obtain upper bounds for the energy norm of the error which are explicit in the mesh-width , in the polynomial degree , and in the regularity of the exact solution. We identify a special numerical flux for which the estimates are optimal in both and . The theoretical results are confirmed in a series of numerical examples.

     

Spatial and spectral superconvergence of discontinuous Galerkin method for hyperbolic problems

In this paper, we analyze the spatial and spectral superconvergence properties of one-dimensional hyperbolic conservation law by a discontinuous Galerkin (DG) method. The analyses combine classical mathematical arguments with MATLAB experiments. Some properties of the DG schemes are discovered using discrete Fourier analyses: superconvergence of the numerical wave numbers, Radau structure of the X spatial error.     

On discontinuous Galerkin finite element method for singularly perturbed delay differential equations

A nonsymmetric discontinuous Galerkin FEM with interior penalties has been applied to one-dimensional singularly perturbed problem with a constant negative shift. Using higher order polynomials on Shishkin-type layer-adapted meshes, a robust convergence has been proved in the corresponding energy norm. Numerical experiments support theoretical findings.     

用Fourier Galerkin方法求解间断特征值问题

本文我们考虑一个带间断系数的特征值问题,使用Fourier G a lerk in方法求解.数值试验表明特征值收敛速度达到三阶,而特征函数收敛速度达到2.5阶.表明此方法对间断系数问题非常有效.     

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

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

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 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.     

On the coupling of finite volume and discontinuous Galerkin method for elliptic problems

The coupling of cell-centered finite volume method with primal discontinuous Galerkin method is introduced in this paper for elliptic problems. Convergence of the method with respect to the mesh size is proved. Numerical examples confirm the theoretical rates of convergence. Advantages of the coupled scheme are shown for problems with discontinuous coefficients or anisotropic diffusion matrix.     

A high-order non-conforming discontinuous Galerkin method for time-domain electromagnetics

In this paper, we discuss the formulation, stability and validation of a high-order non-dissipative discontinuous Galerkin (DG) method for solving Maxwell’s equations on non-conforming simplex meshes. The proposed method combines a centered approximation for the numerical fluxes at inter element boundaries, with either a second-order or a fourth-order leap-frog time integration scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary-level hanging nodes. The method is proved to be stable and conserves a discrete counterpart of the electromagnetic energy for metallic cavities. Numerical experiments with high-order elements show the potential of the method.     

Adaptive discontinuous Galerkin methods in multiwavelets bases

We use a multiwavelet basis with the Discontinuous Galerkin (DG) method to produce a multi-scale DG method. We apply this Multiwavelet DG method to convection and convection-diffusion problems in multiple dimensions. Merging the DG method with multiwavelets allows the adaptivity in the DG method to be resolved through manipulation of multiwavelet coefficients rather than grid manipulation. Additionally, the Multiwavelet DG method is tested on non-linear equations in one dimension and on the cubed sphere.     

A priori error estimates for interior penalty discontinuous Galerkin method applied to nonlinear Sobolev equations

A discontinuous Galerkin method with interior penalties is presented for nonlinear Sobolev equations. A semi-discrete and a family of Fully-discrete time approximate scheme are formulated. These schemes can be symmetric or nonsymmetric. Hp-version error estimates are analyzed for these schemes. Just because of a damping term ·(b(u)u_t) included in Sobolev equation, which is the distinct character different from parabolic equation, special test functions are chosen to deal with this term. Finally, a priori L^∞(H¹) error estimate is derived for the semi-discrete time scheme and similarly, l^∞(H¹) and l²(H¹) for the Fully-discrete time schemes. These results also indicate that spatial rates in H¹ and time truncation errors in L² are optimal.     

A uniformly convergent continuous-discontinuous Galerkin method for singularly perturbed problems of convection-diffusion type

In this paper, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving singularly perturbed convection-diffusion problems. On Shishkin mesh with linear elements, a rate O(N^-1lnN) in an associated norm is established, where N is the number of elements. Numerical experiments complement the theoretical results. Moreover, a rate O(N^-2ln²N) in a discrete L^∞ norm, and O(N^-2) in L² norm, are observed numerically on the Shishkin mesh.     

BDDC methods for discontinuous Galerkin discretization of elliptic problems

A discontinuous Galerkin (DG) discretization of Dirichlet problem for second-order elliptic equations with discontinuous coefficients in 2-D is considered. For this discretization, balancing domain decomposition with constraints (BDDC) algorithms are designed and analyzed as an additive Schwarz method (ASM). The coarse and local problems are defined using special partitions of unity and edge constraints. Under certain assumptions on the coefficients and the mesh sizes across ∂Ω_i, where the Ω_i are disjoint subregions of the original region Ω, a condition number estimate C(1+max_ilog(H_i/h_i))² is established with C independent of h_i, H_i and the jumps of the coefficients. The algorithms are well suited for parallel computations and can be straightforwardly extended to the 3-D problems. Results of numerical tests are included which confirm the theoretical results and the necessity of the imposed assumptions.     

A superconvergent LDG-hybridizable Galerkin method for second-order elliptic problems

We identify and study an LDG-hybridizable Galerkin method, which is not an LDG method, for second-order elliptic problems in several space dimensions with remarkable convergence properties. Unlike all other known discontinuous Galerkin methods using polynomials of degree for both the potential as well as the flux, the order of convergence in of both unknowns is . Moreover, both the approximate potential as well as its numerical trace superconverge in -like norms, to suitably chosen projections of the potential, with order . This allows the application of element-by-element postprocessing of the approximate solution which provides an approximation of the potential converging with order in . The method can be thought to be in between the hybridized version of the Raviart-Thomas and that of the Brezzi-Douglas-Marini mixed methods.

     

Locking-free adaptive discontinuous Galerkin FEM for linear elasticity problems

An adaptive discontinuous Galerkin finite element method for linear elasticity problems is presented. We develop an a posteriori error estimate and prove its robustness with respect to nearly incompressible materials (absence of volume locking). Furthermore, we present some numerical experiments which illustrate the performance of the scheme on adaptively refined meshes.

     

Analysis of a DG method for singularly perturbed convection-diffusion problems

In this article, we studied a discontinuous Galerkin finite element method for convection-diffusion-reaction problems with singular perturbation. Our approach is highly flexible by allowing the use of discontinuous approximating function on polytopal mesh without imposing extra conditions on the convection coefficient. A priori error estimate is devised in a suitable energy norm on general polytopal mesh. Numerical examples are provided.     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

Convergence analysis of the LDG method applied to singularly perturbed problems

Considering a two‐dimensional singularly perturbed convection–diffusion problem with exponential boundary layers, we analyze the local discontinuous Galerkin (DG) method that uses piecewise bilinear polynomials on Shishkin mesh. A convergence rate O(N^‐1 lnN) in a DG‐norm is established under the regularity assumptions, while the total number of mesh points is O(N²). The rate of convergence is uniformly valid with respect to the singular perturbation parameter ε. Numerical experiments indicate that the theoretical error estimate is sharp. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

A method of asymptotic constructions with improved accuracy for a quasilinear singularly perturbed parabolic convection-diffusion equation

The Dirichlet problem on an interval for quasilinear singularly perturbed parabolic convection-diffusion equation is considered. The higher order derivative of the equation is multiplied by a parameter ε that takes any values from the half-open interval (0, 1]. For this type of linear problems, the order of the ε-uniform convergence (with respect to x and t) for the well-known schemes is not higher than unity (in the maximum norm). For the boundary value problem under consideration, grid approximations are constructed that converge ε-uniformly at the rate of O(N ^?2ln² N + N ^?2 ₀), where N + 1 and N ₀ + 1 are the numbers of the mesh points with respect to x and t, respectively. On the x axis, piecewise uniform meshes that condense in the boundary layer are used. If the parameter value is small compared to the effective step of the spatial grid, the domain decomposition method is used, which is motivated by “asymptotic constructions.” Monotone approximations of “auxiliary” subproblems describing the main terms of the asymptotic expansion of the solution outside a neighborhood of the boundary layer neighborhood are used. In the neighborhood of the boundary layer (of the width O(ε ln N)) the first derivative with respect to x is approximated by the central difference derivative. These subproblems are successively solved in the subdomains on uniform grids. If the parameter values are not sufficiently small (compared to the effective step of the mesh with respect to x), the classical implicit difference schemes approximating the first derivative with respect to x by the central difference derivative are applied. To improve the accuracy in t, the defect correction technique is used. Notice that the calculation of the solution of the constructed difference scheme (the scheme based on the method of asymptotic constructions) can be considerably simplified for sufficiently small values of the parameter ε.