High-order accurate thin layer approximations for time-domain electromagnetics, Part II: Transmission layers

We continue the development of high-order accurate thin layer approximations for time-domain electromagnetics and focus in this paper on a new family of models for thin transmission layers. The thin transmission layer approximations are valid for general isotropic materials and certain types of anisotropic materials. The models also allow the inclusion of smoothly curved layers. Both dielectric and magnetic materials can be considered. These models are non-trivial and we discuss their formulation, properties, and implementation in the context of discontinuous Galerkin methods which emerge as being particularly well suited for this family of models. The range of validity, accuracy, and stability of the models and numerical approximations is demonstrated through one- and two-dimensional examples.     

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.     

A robust discontinuous Galerkin method for solving convection-diffusion problems

In this paper, a new DG method was designed to solve the model problem of the one-dimensional singularly-perturbed convection-diffusion equation. With some special chosen numerical traces, the existence and uniqueness of the DG solution is provided. The superconvergent points inside each element are observed. Particularly, the 2p ＋ 1-order superconvergence and even uniform superconvergence under layer-adapted mesh are observed numerically.     

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.     

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.     

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.     

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.     

A discontinuous Galerkin finite element method for time dependent partial differential equations with higher order derivatives

In this paper, we develop a new discontinuous Galerkin (DG) finite element method for solving time dependent partial differential equations (PDEs) with higher order spatial derivatives. Unlike the traditional local discontinuous Galerkin (LDG) method, the method in this paper can be applied without introducing any auxiliary variables or rewriting the original equation into a larger system. Stability is ensured by a careful choice of interface numerical fluxes. The method can be designed for quite general nonlinear PDEs and we prove stability and give error estimates for a few representative classes of PDEs up to fifth order. Numerical examples show that our scheme attains the optimal -th order of accuracy when using piecewise -th degree polynomials, under the condition that is greater than or equal to the order of the equation.

     

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.     

A Fourier spectral-discontinuous Galerkin method for time-dependent 3-D Schrödinger–Poisson equations with discontinuous potentials

In this paper, we propose a high order Fourier spectral-discontinuous Galerkin method for time-dependent Schrödinger–Poisson equations in 3-D spaces. The Fourier spectral Galerkin method is used for the two periodic transverse directions and a high order discontinuous Galerkin method for the longitudinal propagation direction. Such a combination results in a diagonal form for the differential operators along the transverse directions and a flexible method to handle the discontinuous potentials present in quantum heterojunction and supperlattice structures. As the derivative matrices are required for various time integration schemes such as the exponential time differencing and Crank Nicholson methods, explicit derivative matrices of the discontinuous Galerkin method of various orders are derived. Numerical results, using the proposed method with various time integration schemes, are provided to validate the method.     

Wavelets and adaptive grids for the discontinuous Galerkin method

In this paper, space adaptivity is introduced to control the error in the numerical solution of hyperbolic systems of conservation laws. The reference numerical scheme is a new version of the discontinuous Galerkin method, which uses an implicit diffusive term in the direction of the streamlines, for stability purposes. The decision whether to refine or to unrefine the grid in a certain location is taken according to the magnitude of wavelet coefficients, which are indicators of local smoothness of the numerical solution. Numerical solutions of the nonlinear Euler equations illustrate the efficiency of the method.     

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.     

The local discontinuous Galerkin method for the Oseen equations

We introduce and analyze the local discontinuous Galerkin method for the Oseen equations of incompressible fluid flow. For a class of shape-regular meshes with hanging nodes, we derive optimal a priori estimates for the errors in the velocity and the pressure in - and negative-order norms. Numerical experiments are presented which verify these theoretical results and show that the method performs well for a wide range of Reynolds numbers.

     

High-order discontinuous element-based schemes for the inviscid shallow water equations: Spectral multidomain penalty and discontinuous Galerkin methods

Two commonly used types of high-order-accuracy element-based schemes, collocation-based spectral multidomain penalty methods (SMPM) and nodal discontinuous Galerkin methods (DGM), are compared in the framework of the inviscid shallow water equations. Differences and similarities in formulation are identified, with the primary difference being the dissipative term in the Rusanov form of the numerical flux for the DGM that provides additional numerical stability; however, it should be emphasized that to arrive at this equivalence between SMPM and DGM requires making specific choices in the construction of both methods; these choices are addressed. In general, both methods offer a multitude of choices in the penalty terms used to introduce boundary conditions and stabilize the numerical solution. The resulting specialized class of SMPM and DGM are then applied to a suite of six commonly considered geophysical flow test cases, three linear and three non-linear; we also include results for a classical continuous Galerkin (i.e., spectral element) method for comparison. Both the analysis and numerical experiments show that the SMPM and DGM are essentially identical; both methods can be shown to be equivalent for very special choices of quadrature rules and Riemann solvers in the DGM along with special choices in the type of penalty term in the SMPM. Although we only focus our studies on the inviscid shallow water equations the results presented should be applicable to other systems of nonlinear hyperbolic equations (such as the compressible Euler equations) and extendable to the compressible and incompressible Navier-Stokes equations, where viscous terms are included.     

A unified discontinuous Galerkin framework for time integration

We introduce a new discontinuous Galerkin approach for time integration. On the basis of the method of weighted residual, numerical quadratures are employed in the finite element time discretization to account for general nonlinear ordinary differential equations. Many different conditions, including explicit, implicit, and symplectic conditions, are enforced for the test functions in the variational analysis to obtain desirable features of the resulting time‐stepping scheme. The proposed discontinuous Galerkin approach provides a unified framework to derive various time‐stepping schemes, such as low‐order one‐step methods, Runge–Kutta methods, and multistep methods. On the basis of the proposed framework, several explicit Runge–Kutta methods of different orders are constructed. The derivation of symplectic Runge–Kutta methods has also been realized. The proposed framework allows the optimization of new schemes in terms of several characteristics, such as accuracy, sparseness, and stability. The accuracy optimization is performed on the basis of an analytical form of the error estimation function for a linear test initial value problem. Schemes with higher formal order of accuracy are found to provide more accurate solutions. We have also explored the optimization potential of sparseness, which is related to the general compressive sensing in signal/imaging processing. Two critical dimensions of the stability region, that is, maximal intervals along the imaginary and negative real axes, are employed as the criteria for stability optimization. This gives the largest Courant–Friedrichs–Lewy time steps in solving hyperbolic and parabolic partial differential equations, respectively. Numerical experiments are conducted to validate the optimized time‐stepping schemes. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

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.

     

Uniqueness of numerical solutions to nonlinear parabolic equations by a fully implicit discontinuous Galerkin method

We prove uniqueness of numerical solutions to nonlinear parabolic equations approximated by a fully implicit interior penalty discontinuous Galerkin (IPDG) method, with a mesh-independent constraint on time step.     

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.