In applied mathematics, discontinuous Galerkin methods (DG methods) 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 firstorder 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, Joachim Nitsche and Miloš Zlámal. DG methods for elliptic problems were already developed in a paper by Garth 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.
YouTube Encyclopedic

1/5Views:2 9811 0453111 240601

✪ "Discontinuous Galerkin Methods for Hyerbolic PDEs: 1"  Olindo Zanotti

✪ "Discontinuous Galerkin Methods for Hyperbolic PDEs: 2"  Olindo Zanotti

✪ Weight Adjusted Discontinuous Galerkin methods for Acoustic and Elastic Wave Propagation

✪ John Sipe 3

✪ DOE CSGF 2012: Toward Adaptive SpaceTime Discontinuous Galerkin Discretizations of Unsteady Aero...
Transcription
Contents
Overview
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 finitedimensional inner product subspaces utilized in conforming methods.
As an example, consider the continuity equation for a scalar unknown in a spatial domain without "sources" or "sinks" :
where is the flux of .
Now consider the finitedimensional space of discontinuous piecewise polynomial functions over the spatial domain restricted to a discrete triangulation , written as
for the space of polynomials with degrees less than or equal to over element indexed by . Then for finite element shape functions the solution is represented by
Then similarly choosing a test function
multiplying the continuity equation by and integrating by parts in space, the semidiscrete DG formulation becomes:
Scalar hyperbolic conservation law
A scalar hyperbolic conservation law is of the form
where one tries to solve for the unknown scalar function , and the functions are typically given.
Space discretization
The space will be discretized as
Furthermore, we need the following definitions
Basis for function space
We derive the basis representation for the function space of our solution . The function space is defined as
where denotes the restriction of onto the interval , and denotes the space of polynomials of maximal degree . The index should show the relation to an underlying discretization given by . Note here that is not uniquely defined at the intersection points .
At first we make use of a specific polynomial basis on the interval , the Legendre polynomials , i.e.,
Note especially the orthogonality relations
Transformation onto the interval , and normalization is achieved by functions
which fulfill the orthonormality relation
Transformation onto an interval is given by
which fulfill
For normalization we define , and for normalization we define , s.t.
Finally, we can define the basis representation of our solutions
Note here, that is not defined at the interface positions.
Besides, prism bases are employed for planarlike structures, and are capable for 2D/3D hybridation.
DGscheme
The conservation law is transformed into its weak form by multiplying with test functions, and integration over test intervals
By using partial integration one is left with
The fluxes at the interfaces are approximated by numerical fluxes with
where denotes the left and righthand sided limits. Finally, the DGScheme can be written as
Scalar elliptic equation
A scalar elliptic equation is of the form
This equation is the steadystate heat equation, where is the temperature. Space discretization is the same as above. We recall that the interval is partitioned into intervals of length .
We introduce jump and average of functions at the node :
The interior penalty discontinuous Galerkin (IPDG) method is: find satisfying
where the bilinear forms and are
and
The linear forms and are
and
The penalty parameter is a positive constant. Increasing its value will reduce the jumps in the discontinuous solution. The term is chosen to be equal to for the symmetric interior penalty Galerkin method; it is equal to for the nonsymmetric interior penalty Galerkin method.
See also
References
 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.
 A. Cangiani, Z. Dong, E.H. Georgoulis, and P. Houston, hpVersion Discontinuous Galerkin Methods on Polygonal and Polyhedral Meshes, SpringerBriefs in Mathematics, (December 2017).
 W. Mai, J. Hu, P. Li, and H. Zhao, “An efficient and stable 2D/3D hybrid discontinuous Galerkin timedomain analysis with adaptive criterion for arbitrarily shaped antipads in dispersive parallelplate pair,” IEEE Trans. Microw. Theory Techn., vol. 65, no. 10, pp. 3671–3681, Oct. 2017.
 W. Mai et al., “A straightforward updating criterion for 2D/3D hybrid discontinuous Galerkin timedomain method controlling comparative error,” IEEE Trans. Microw. Theory Techn., vol. 66, no. 4, pp. 1713–1722, Apr. 2018.
 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. SpringerVerlag, Berlin, 2000.
 P. Lesaint, and P. A. Raviart. "On a finite element method for solving the neutron transport equation." Mathematical aspects of finite elements in partial differential equations 33 (1974): 89–123.
 D.A. Di Pietro and A. Ern, Mathematical Aspects of Discontinuous Galerkin Methods. Mathématiques et Applications, Vol. 69, SpringerVerlag, 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.
 B. Rivière, Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations: Theory and Implementation. SIAM Frontiers in Applied Mathematics, 2008.
 CFD Wiki http://www.cfdonline.com/Wiki/Discontinuous_Galerkin
 W.H. Reed and T.R. Hill, Triangular mesh methods for the neutron transport equation, Tech. Report LAUR73–479, Los Alamos Scientific Laboratory, 1973.