Discontinuous Galerkin method

From Wikipedia, the free encyclopedia
Jump to: navigation, search

Discontinuous Galerkin methods (DG methods) in mathematics form a class of numerical methods for solving differential equations. They combine features of the finite element and the finite volume framework and have been successfully applied to hyperbolic, elliptic, parabolic and mixed form problems arising from a wide range of applications. DG methods have in particular received considerable interest for problems with a dominant first-order part, e.g. in electrodynamics, fluid mechanics and plasma physics.

Discontinuous Galerkin methods were first proposed and analyzed in the early 1970s as a technique to numerically solve partial differential equations. In 1973 Reed and Hill introduced a DG method to solve the hyperbolic neutron transport equation.

The origin of the DG method for elliptic problems cannot be traced back to a single publication as features such as jump penalization in the modern sense were developed gradually. However, among the early influential contributors were Babuška, J.-L. Lions, Nitsche and Zlamal. DG methods for elliptic problems were already developed in a paper by Baker in the setting of 4th order equations in 1977. A more complete account of the historical development and an introduction to DG methods for elliptic problems is given in a publication by Arnold, Brezzi, Cockburn and Marini. A number of research directions and challenges on DG methods are collected in the proceedings volume edited by Cockburn, Karniadakis and Shu.


Much like the continuous Galerkin (CG) method, the discontinuous Galerkin (DG) method is a finite element method formulated relative to a weak formulation of a particular model system. Unlike traditional CG methods that are conforming, the DG method works over a trial space of functions that are only piecewise continuous, and thus often comprise more inclusive function spaces than the finite-dimensional inner product subspaces utilized in conforming methods.

As an example, consider the continuity equation arising in electromagnetic theory:

 \frac{\partial \rho} {\partial t} + \nabla \cdot \mathbf{J} = 0.

Now consider the finite-dimensional space of discontinuous piecewise polynomial functions over the spatial domain \Omega restricted to a discrete triangulation \Omega_{h}, written as

 S_{h}^{p}(\Omega_{h})=\{v_{|\Omega_{e_{i}}}\in P^{p}(\Omega_{e_{i}}), \ \ \forall\Omega_{e_{i}}\in \Omega_{h}\}

for P^{p}(\Omega_{e_{i}}) the space of degree \leq p polynomials over element \Omega_{e_{i}} indexed by i. Then for finite element shape functions N_{j}\in P^{p} the solution is represented by

\rho_{h} = \sum_{j=1}^{dofs}\rho_{j}^{i}(t)N_{j}^{i}(\boldsymbol{x}), \quad \forall\boldsymbol{x}\in\Omega_{e_{i}}.

Then similarly choosing a test function

\phi_{h}(\boldsymbol{x})=\sum_{j=1}^{dofs}\phi_{j}^{i}N_{j}^{i}(\boldsymbol{x}), \quad \forall \boldsymbol{x}\in\Omega_{e_{i}},

multiplying the continuity equation by \phi_h and integrating by parts in space, the semidiscrete DG formulation becomes:

 \frac{d}{dt}\int_{\Omega_{e_{i}}}\rho_{h}\phi_{h}d\boldsymbol{x} + \int_{\partial\Omega_{e_{i}}} \phi_{h}\mathbf{J}_{h}\cdot\boldsymbol{n} d\boldsymbol{x} = \int_{\Omega_{e_{i}}}\mathbf{J}_{h}\cdot\nabla\phi_{h}d\boldsymbol{x}.

See also[edit]


  • D.N. Arnold, F. Brezzi, B. Cockburn and L.D. Marini, Unified analysis of discontinuous Galerkin methods for elliptic problems, SIAM J. Numer. Anal. 39(5):1749-1779, 2002.
  • G. Baker, Finite element methods for elliptic equations using nonconforming elements, Math. Comp. 31 (1977), no. 137, 45–59.
  • B. Cockburn, G. E. Karniadakis and C.-W. Shu (eds.), Discontinuous Galerkin methods. Theory, computation and applications, Lecture Notes in Computational Science and Engineering, 11. Springer-Verlag, Berlin, 2000.
  • D.A. Di Pietro and A. Ern, Mathematical Aspects of Discontinuous Galerkin Methods. Mathématiques et Applications, Vol. 69, Springer-Verlag, Berlin, 2011.
  • J.S. Hesthaven and T. Warburton, Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications. Springer Texts in Applied Mathematics 54. Springer Verlag, New York, 2008.
  • CFD Wiki http://www.cfd-online.com/Wiki/Discontinuous_Galerkin
  • W.H. Reed and T.R. Hill, Triangular mesh methods for the neutron transport equation, Tech. Report LA-UR-73-479, Los Alamos Scientific Laboratory, 1973.