Numerical error analysis for nonsymmetric interior penalty discontinuous Galerkin method of Cahn–Hilliard equation

In this paper, we derive a theoretical analysis of nonsymmetric interior penalty discontinuous Galerkin methods for solving the Cahn–Hilliard equation. We prove unconditional unique solvability of the discrete system and derive stability bounds with a generalized chemical energy density. Convergence of the method is obtained by optimal a priori error estimates. Our analysis is valid for both symmetric and nonsymmetric versions of the discontinuous Galerkin formulation.     

Analysis of the parareal approach based on discontinuous Galerkin method for time-dependent Stokes equations

This paper analyzes a parareal approach based on discontinuous Galerkin (DG) method for the time-dependent Stokes equations. A class of primal discontinuous Galerkin methods, namely variations of interior penalty methods, are adopted for the spatial discretization in the parareal algorithm (we call it parareal DG algorithm). We study three discontinuous Galerkin methods for the time-dependent Stokes equations, and the optimal continuous in time error estimates for the velocities and pressure are derived. Based on these error estimates, the proposed parareal DG algorithm is proved to be unconditionally stable and bounded by the error of discontinuous Galerkin discretization after a finite number of iterations. Finally, some numerical experiments are conducted which confirm our theoretical results, meanwhile, the efficiency of the parareal DG algorithm can be seen through a parallel experiment.     

Local and pointwise error estimates of the local discontinuous Galerkin method applied to the Stokes problem

We prove local and pointwise error estimates for the local discontinuous Galerkin method applied to the Stokes problem in two and three dimensions. By using techniques originally developed by A. Schatz [Math. Comp., 67 (1998), 877-899] to prove pointwise estimates for the Laplace equation, we prove optimal weighted pointwise estimates for both the velocity and the pressure for domains with smooth boundaries.

     

Error estimates for a discontinuous Galerkin method with interior penalties 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 schemes are formulated. These schemes are symmetric. Hp‐version error estimates are analyzed for these schemes. For the semi‐discrete time scheme a priori L^∞(H¹) error estimate is derived and similarly, l^∞(H¹) and l²(H¹) for the fully‐discrete time schemes. These results indicate that spatial rates in H¹ and time truncation errors in L² are optimal. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Technical Note: A note on the selection of the penalty parameter for discontinuous Galerkin finite element schemes

We obtain a computable lower bound on the value of the interior penalty parameters sufficient for the existence of a unique discontinuous Galerkin finite element approximation of a second order elliptic problem. The bound obtained is valid for meshes containing an arbitrary number of hanging nodes and elements of arbitrary nonuniform polynomial order. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

A divergence-free weak Galerkin method for quasi-Newtonian Stokes flows

This paper proposes a weak Galerkin finite element method to solve incompressible quasi-Newtonian Stokes equations. We use piecewise polynomials of degrees k + 1(k 0) and k for the velocity and pressure in the interior of elements, respectively, and piecewise polynomials of degrees k and k + 1 for the boundary parts of the velocity and pressure, respectively. Wellposedness of the discrete scheme is established. The method yields a globally divergence-free velocity approximation. Optimal priori error estimates are derived for the velocity gradient and pressure approximations. Numerical results are provided to confirm the theoretical results.     

A selective immersed discontinuous Galerkin method for elliptic interface problems

This article proposes a selective immersed discontinuous Galerkin method based on bilinear immersed finite elements (IFE) for solving second‐order elliptic interface problems. This method applies the discontinuous Galerkin formulation wherever selected, such as those elements around an interface or a singular source, but the regular Galerkin formulation everywhere else. A selective bilinear IFE space is constructed and applied to the selective immersed discontinuous Galerkin method based on either the symmetric or nonsymmetric interior penalty discontinuous Galerkin formulation. The new method can solve an interface problem by a rectangular mesh with local mesh refinement independent of the interface even if its geometry is nontrivial. Meanwhile, if desired, its computational cost can be maintained very close to that of the standard Galerkin IFE method. It is shown that the selective bilinear IFE space has the optimal approximation capability expected from piecewise bilinear polynomials. Numerical examples are provided to demonstrate features of this method, including the effectiveness of local mesh refinement around the interface and the sensitivity to the penalty parameters. Copyright © 2013 John Wiley & Sons, Ltd.     

A discontinuous Galerkin method for the Rosenau equation

A priori error estimates for the Rosenau equation, which is a K-dV like Rosenau equation modelled to describe the dynamics of dense discrete systems, have been studied by one of the authors. But since a priori error bounds contain the unknown solution and its derivatives, it is not effective to control error bounds with only a given step size. Thus we need to estimate a posteriori errors in order to control accuracy of approximate solutions using variable step sizes. A posteriori error estimates of the Rosenau equation are obtained by a discontinuous Galerkin method and the stability analysis is discussed for the dual problem. Numerical results on a posteriori error and wave propagation are given, which are obtained by using various spatial and temporal meshes controlled automatically by a posteriori error.     

Connections between discontinuous Galerkin and nonconforming finite element methods for the Stokes equations

We study a discontinuous Galerkin finite element method (DGFEM) for the Stokes equations with a weak stabilization of the viscous term. We prove that, as the stabilization parameter γ tends to infinity, the solution converges at speed γ^?1 to the solution of some stable and well‐known nonconforming finite element methods (NCFEM) for the Stokes equations. In addition, we show that an a posteriori error estimator for the DGFEM‐solution based on the reconstruction of a locally conservative H(div, Ω)‐tensor tends at the same speed to a classical a posteriori error estimator for the NCFEM‐solution. These results can be used to affirm the robustness of the DGFEM‐method and also underline the close relationship between the two approaches. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Fully discrete interior penalty discontinuous Galerkin methods for nonlinear parabolic equations

In this article, we investigate interior penalty discontinuous Galerkin (IPDG) methods for solving a class of two‐dimensional nonlinear parabolic equations. For semi‐discrete IPDG schemes on a quasi‐uniform family of meshes, we obtain a priori bounds on solutions measured in the L² norm and in the broken Sobolev norm. The fully discrete IPDG schemes considered are based on the approximation by forward Euler difference in time and broken Sobolev space. Under a restriction related to the mesh size and time step, an hp ‐version of an a priori l^∞(L²) and l²(H¹) error estimate is derived and numerical experiments are presented.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 288–311, 2012     

An Eulerian‐Lagrangian discontinuous Galerkin method for transient advection‐diffusion equations

We develop an Eulerian‐Lagrangian discontinuous Galerkin method for time‐dependent advection‐diffusion equations. The derived scheme has combined advantages of Eulerian‐Lagrangian methods and discontinuous Galerkin methods. The scheme does not contain any undetermined problem‐dependent parameter. An optimal‐order error estimate and superconvergence estimate is derived. Numerical experiments are presented, which verify the theoretical estimates.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Local discontinuous Galerkin approximations to fractional Bagley-Torvik equation

In this research, numerical approximation to fractional Bagley-Torvik equation as an important model arising in fluid mechanics is investigated. Our discretization algorithm is based on the local discontinuous Galerkin (LDG) schemes along with using the natural upwind fluxes, which enables us to solve the model problem element by element. This means that we require to solve a low-order system of equations in each subinterval, hence avoiding the need for a full global solution. The proposed schemes are tested on a range of initial- and boundary-value problems including a variable coefficient example, a nonsmooth problem, and some oscillatory test cases with exact solutions. Also, the validation of the proposed methods was compared with those obtained available existing computational procedures. Overall, it was found that LDG methods indicated highly satisfactory performance with comparatively lower degree of polynomials and number of elements compared with other numerical models.     

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.     

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 C0 interior penalty discontinuous Galerkin method for fourth‐order total variation flow. I: Derivation of the method and numerical results

We consider the numerical solution of a fourth‐order total variation flow problem representing surface relaxation below the roughening temperature. Based on a regularization and scaling of the nonlinear fourth‐order parabolic equation, we perform an implicit discretization in time and a C⁰ Interior Penalty Discontinuous Galerkin (C⁰IPDG) discretization in space. The C⁰IPDG approximation can be derived from a mixed formulation involving numerical flux functions where an appropriate choice of the flux functions allows to eliminate the discrete dual variable. The fully discrete problem can be interpreted as a parameter dependent nonlinear system with the discrete time as a parameter. It is solved by a predictor corrector continuation strategy featuring an adaptive choice of the time step sizes. A documentation of numerical results is provided illustrating the performance of the C⁰IPDG method and the predictor corrector continuation strategy. The existence and uniqueness of a solution of the C⁰IPDG method will be shown in the second part of this paper.     

Solvability of the fractional Volterra‐Fredholm integro differential equation by hybridizable discontinuous Galerkin method

This article proves the existence and uniqueness of the solution obtained by the hybridizable discontinuous Galerkin (HDG) method of the fractional Volterra‐Fredholm integro differential equation. The method based on local solvers and transmission condition is applied to the equation using two auxiliary variables. The form of the equation is amenable for achieving the solvability criteria of the problem according to the HDG method. We also calculate a numerical solution of the problem whose exact solution is taken as a smooth or fractional function. This results in a tridiagonal, symmetric, and positive definite stiffness matrix.     

The discontinuous Galerkin method for discretely observed Asian options

Asian options represent an important subclass of the path-dependent contracts that are identified by payoff depending on the average of the underlying asset prices over the prespecified period of option lifetime. Commonly, this average is observed at discrete dates, and also, early exercise features can be admitted. As a result, analytical pricing formulae are not always available. Therefore, some form of a numerical approximation is essential for efficient option valuation. In this paper, we study a PDE model for pricing discretely observed arithmetic Asian options with fixed as well as floating strike for both European and American exercise features. The pricing equation for such options is similar to the Black-Scholes equation with 1 underlying asset, and the corresponding average appears only in the jump conditions across the sampling dates. The objective of the paper is to present the comprehensive methodological concept that forms and improves the valuation process. We employ a robust numerical procedure based on the discontinuous Galerkin approach arising from the piecewise polynomial generally discontinuous approximations. This technique enables a simple treatment of discrete sampling by incorporation of jump conditions at each monitoring date. Moreover, an American early exercise constraint is directly handled as an additional nonlinear source term in the pricing equation. The proposed solving procedure is accompanied by an empirical study with practical results compared to reference values.     

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

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

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.     

On discontinuous Galerkin method and semiregular family of triangulations

Discretization of second order elliptic partial differential equations by discontinuous Galerkin method often results in numerical schemes with penalties. In this paper we analyze these penalized schemes in the context of quite general triangular meshes satisfying only a semiregularity assumption. A new (modified) penalty term is presented and theoretical properties are proven together with illustrative numerical results. This work is a part of the research project MSM 0021620839 financed by MSMT and was partly supported by the project No. 201/04/1503 of the Grant Agency of the Czech Republic.