A unified study of continuous and discontinuous Galerkin methods

A unified study is presented in this paper for the design and analysis of different finite element methods(FEMs), including conforming and nonconforming FEMs, mixed FEMs, hybrid FEMs, discontinuous Galerkin(DG) methods, hybrid discontinuous Galerkin(HDG) methods and weak Galerkin(WG) methods.Both HDG and WG are shown to admit inf-sup conditions that hold uniformly with respect to both mesh and penalization parameters. In addition, by taking the limit of the stabilization parameters, a WG method is shown to converge to a mixed method whereas an HDG method is shown to converge to a primal method. Furthermore,a special class of DG methods, known as the mixed DG methods, is presented to fill a gap revealed in the unified framework.     

UNIFORM STABILITY AND ERROR ANALYSIS FOR SOME DISCONTINUOUS GALERKIN METHODS

In this paper,we provide a number of new estimates on the stability and convergence of both hybrid discontinuous Galerkin(HDG)and weak Galerkin(WG)methods.By using the standard Brezzi theory on mixed methods,we carefully define appropriate norms for the various discretization variables and then establish that the stability and error estimates hold uniformly with respect to stabilization and discretization parameters.As a result,by taking appropriate limit of the stabilization parameters,we show that the HDG method converges to a primal conforming method and the WG method converges to a mixed conforming method.     

A nonstandard multigrid method with flexible multiple semicoarsening for the numerical solution of the pressure equation in a Navier-Stokes solver

Numerical methods for the incompressible Reynolds-averaged Navier-Stokes equations discretized by finite difference techniques on collocated cell-centered structured grids are considered in this paper. A widespread solution method to solve the pressure-velocity coupling problem is to use a segregated approach, in which the computational work is deeply controlled by the solution of the pressure problem. This pressure equation is an elliptic partial differential equation with possibly discontinuous or anisotropic coeffficients. The resulting singular linear system needs efficient solution strategies especially for 3-dimensional applications. A robust method (close to MG-S [22,34]) combining multiple cell-centered semicoarsening strategies, matrix-independent transfer operators, Galerkin coarse grid approximation is therefore designed. This strategy is both evaluated as a solver or as a preconditioner for Krylov subspace methods on various 2- or 3-dimensional fluid flow problems. The robustness of this method is shown. This revised version was published online in June 2006 with corrections to the Cover Date.     

An efficient unified approach for the numerical solution of delay differential equations

In this paper we propose a new framework for designing a delay differential equation (DDE) solver which works with any supplied initial value problem (IVP) solver that is based on a standard step-by-step approach, such as Runge-Kutta or linear multi-step methods, and can provide dense output. This is done by treating a general DDE as a special example of a discontinuous IVP. Using this interpretation we develop an efficient technique to solve the resulting discontinuous IVP. We also give a more clear process for the numerical techniques used when solving the implicit equations that arise on a time step, such as when the underlying IVP solver is implicit or the delay vanishes. The new modular design for the resulting simulator we introduce, helps to accelerate the utilization of advances in the different components of an effective numerical method. Such components include the underlying discrete formula, the interpolant for dense output, the strategy for handling discontinuities and the iteration scheme for solving any implicit equations that arise.     

Fully discrete dynamic mesh discontinuous Galerkin methods for the Cahn-Hilliard equation of phase transition

Fully discrete discontinuous Galerkin methods with variable mesh- es in time are developed for the fourth order Cahn-Hilliard equation arising from phase transition in materials science. The methods are formulated and analyzed in both two and three dimensions, and are proved to give optimal order error bounds. This coupled with the flexibility of the methods demonstrates that the proposed discontinuous Galerkin methods indeed provide an efficient and viable alternative to the mixed finite element methods and nonconforming (plate) finite element methods for solving fourth order partial differential equations.

     

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.     

The discontinuous Galerkin method for fractional degenerate convection-diffusion equations

We propose and study discontinuous Galerkin methods for strongly degenerate convection-diffusion equations perturbed by a fractional diffusion (Lévy) operator. We prove various stability estimates along with convergence results toward properly defined (entropy) solutions of linear and nonlinear equations. Finally, the qualitative behavior of solutions of such equations are illustrated through numerical experiments.     

SUPERCONVERGENCE OF DISCONTINUOUS GALERKIN METHOD FOR NONSTATIONARY HYPERBOLIC EQUATION

AbstractFor the first order nonstationary hyperbolic equation taking the piecewise linear discontinuous Galerkin solver, we prove that under the uniform rectangular partition, such a discontinuous solver, after postprossesing, can have two and half approximative order which is half order higher than the optimal estimate by Lesaint and Raviart under the rectangular partition.     

Domain Decomposition Preconditioners for Discontinuous Galerkin Discretizations of Compressible Fluid Flows

In this article we consider the application of Schwarz-type domain decomposition preconditioners to the discontinuous Galerkin finite element approximation of the compressible Navier-Stokes equations. To discretize this system of conservation laws, we exploit the (adjoint consistent) symmetric version of the interior penalty discontinuous Galerkin finite element method. To define the necessary coarse-level solver required for the definition of the proposed preconditioner, we exploit ideas from composite finite element methods, which allow for the definition of finite element schemes on general meshes consisting of polygonal (agglomerated) elements. The practical performance of the proposed preconditioner is demonstrated for a series of viscous test cases in both two- and three-dimensions.     

Implementing a discontinuous Galerkin method for the compressible,inviscid Euler equations in the DUNE framework

A discontinuous Galerkin scheme was implemented in the DUNE framework to solve the compressible, inviscid Euler equations in a multi-dimensional Cartesian grid. It uses a HLLC Riemann solver for the numerical fluxes in the interfaces, a total variation bounded limiter to handle discontinuities, a positivity preserving limiter for near vacuum conditions, and adaptive mesh refinement (AMR). (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Simulation of elastic wave propagation in geological media: Intercomparison of three numerical methods

For wave propagation in heterogeneous media, we compare numerical results produced by grid-characteristic methods on structured rectangular and unstructured triangular meshes and by a discontinuous Galerkin method on unstructured triangular meshes as applied to the linear system of elasticity equations in the context of direct seismic exploration with an anticlinal trap model. It is shown that the resulting synthetic seismograms are in reasonable quantitative agreement. The grid-characteristic method on structured meshes requires more nodes for approximating curved boundaries, but it has a higher computation speed, which makes it preferable for the given class of problems.     

Application of hp-discontinuous Galerkin finite element methods to the rotating disk electrode problems in electrochemistry

This paper presents the interior penalty discontinuous Galerkin finite element methods (DGFEM) for solving the rotating disk electrode problems in electrochemistry. We present results for the simple E reaction mechanism (convection-diffusion equations), the EC’ reaction mechanism (reaction-convection-diffusion equation) and the ECE and EC₂E reaction mechanisms (linear and nonlinear systems of reaction-convection-diffusion equations, respectively). All problems will be in one dimension.     

Development of discontinuous Galerkin methods for Maxwell’s equations in metamaterials and perfectly matched layers

The discontinuous Galerkin method has proved to be an accurate and efficient way to numerically solve many differential equations. In this paper, we extend this method to solve the time-dependent Maxwell’s equations when metamaterials and perfectly matched layers are involved. Numerical results are presented to demonstrate that our method is not only simple to implement, but also quite effective in solving Maxwell’s equations in complex media.     

Numerical solution of the Burgers’ equation by local discontinuous Galerkin method

In this paper, the Burgers’ equation is transformed into the linear diffusion equation by using the Hopf–Cole transformation. The obtained linear diffusion equation is discretized in space by the local discontinuous Galerkin method. The temporal discretization is accomplished by the total variation diminishing Runge–Kutta method. Numerical solutions are compared with the exact solution and the numerical solutions obtained by Adomian’s decomposition method, finite difference method, B-spline finite element method and boundary element method. The results show that the local discontinuous Galerkin method is one of the most efficient methods for solving the Burgers’ equation. Even with small viscosity coefficient, it can get the satisfied solution.     

Coupling of compressible and incompressible flow regions using the multiple pressure variables approach

In many cases, multiphase flows are simulated on the basis of the incompressible Navier–Stokes equations. This assumption is valid as long as the density changes in the gas phase can be neglected. Yet, for certain technical applications such as fuel injection, this is no longer the case, and at least the gaseous phase has to be treated as a compressible fluid. In this paper, we consider the coupling of a compressible flow region to an incompressible one based on a splitting of the pressure into a thermodynamic and a hydrodynamic part. The compressible Euler equations are then connected to the Mach number zero limit equations in the other region. These limit equations can be solved analytically in one space dimension that allows to couple them to the solution of a half‐Riemann problem on the compressible side with the help of velocity and pressure jump conditions across the interface. At the interface location, the flux terms for the compressible flow solver are provided by the coupling algorithms. The coupling is demonstrated in a one‐dimensional framework by use of a discontinuous Galerkin scheme for compressible two‐phase flow with a sharp interface tracking via a ghost‐fluid type method. The coupling schemes are applied to two generic test cases. The computational results are compared with those obtained with the fully compressible two‐phase flow solver, where the Mach number zero limit is approached by a weakly compressible fluid. For all cases, we obtain a very good agreement between the coupling approaches and the fully compressible solver. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

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.     

Simulation of the interaction of light with 3-D metallic nanostructures using a proper orthogonal decomposition-Galerkin reduced-order discontinuous Galerkin time-domain method

In this artice, we report on a reduced-order model (ROM) based on the proper orthogonal decomposition (POD) technique for the system of 3-D time-domain Maxwell's equations coupled to a Drude dispersion model, which is employed to describe the interaction of light with nanometer scale metallic structures. By using the singular value decomposition (SVD) method, the POD basis vectors are extracted offline from the snapshots produced by a high order discontinuous Galerkin time-domain (DGTD) solver. With a Galerkin projection and a second order leap-frog (LF₂) time discretization, a discrete ROM is constructed. The stability condition of the ROM is then analyzed. In particular, when the boundary is a perfect electric conductor condition, the global energy of the ROM is bounded, which is consistent with the characteristics of global energy in the DGTD method. It is shown that the ROM based on Galerkin projection can maintain the stability characteristics of the original high dimensional model. Numerical experiments are presented to verify the accuracy, demonstrate the capabilities of the POD-based ROM and assess its efficiency for 3-D nanophotonic problems.     

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 line search filter algorithm with inexact step computations for equality constrained optimization

In this paper, a new line search filter algorithm for equality constrained optimization is presented. The approach belongs to the class of inexact Newton-like methods. It can also be regarded as an inexact version of generic sequential quadratic programming (SQP) methods. The trial step is obtained by truncatedly solving the primal-dual system based on any robust and efficient linear system solver. Practical termination tests for the linear system solver are established to ensure global convergence. Preliminary numerical results demonstrate the approach is potentially useful.