Space-time discontinuous Galerkin discretizations for linear first-order hyperbolic evolution systems

We introduce a space-time discretization for linear first-order hyperbolic evolution systems using a discontinuous Galerkin approximation in space and a Petrov-Galerkin scheme in time. We show well-posedness and convergence of the discrete system. Then we introduce an adaptive strategy based on goal-oriented dual-weighted error estimation. The full space-time linear system is solved with a parallel multilevel preconditioner. Numerical experiments for the linear transport equation and the Maxwell equation in 2D underline the effciency of the overall adaptive solution process

Space-time discontinuous Galerkin discretizations for linear first-order hyperbolic evolution systems. Revised March 2016

We introduce a space-time discretization for linear first-order hyperbolic evolution systems using a discontinuous Galerkin approximation in space and a Petrov-Galerkin scheme in time. We show well-posedness and convergence of the discrete system. Then we introduce an adaptive strategy based on goal-oriented dual-weighted error estimation. The full space-time linear system is solved with a parallel multilevel preconditioner. Numerical experiments for the linear transport equation and the Maxwell equation in 2D underline the effciency of the overall adaptive solution process

A unified error analysis for spatial discretizations of wave-type equations with applications to dynamic boundary conditions

This thesis provides a unified framework for the error analysis of non-conforming space discretizations of linear wave equations in time-domain, which can be cast as symmetric hyperbolic systems or second-order wave equations. Such problems can be written as first-order evolution equations in Hilbert spaces with linear monotone operators. We employ semigroup theory for the well-posedness analysis and to obtain stability estimates for the space discretizations. To compare the finite dimensional approximations with the original solution, we use the concept of a lift from the discrete to the continuous space. Time integration with the Crank–Nicolson method is analyzed. In this framework, we derive a priori error bounds for the abstract space semi-discretization in terms of interpolation and discretization errors. These error bounds yield previously unkown convergence rates for isoparametric finite element discretizations of wave equations with dynamic boundary conditions in smooth domains. Moreover, our results allow to consider already investigated space discretizations in a unified way. Here it successfully reproduces known error bounds. Among the examples which we dicuss in this thesis are discontinuous Galerkin discretizations of Maxwell’s equations and finite elements with mass lumping for the scalar wave equation

A posteriori analysis of fully discrete method of lines DG schemes for systems of conservation laws

We present reliable a posteriori estimators for some fully discrete schemes applied to nonlinear systems of hyperbolic conservation laws in one space dimension with strictly convex entropy. The schemes are based on a method of lines approach combining discontinuous Galerkin spatial discretization with single- or multi-step methods in time. The construction of the estimators requires a reconstruction in time for which we present a very general framework first for odes and then apply the approach to conservation laws. The reconstruction does not depend on the actual method used for evolving the solution in time. Most importantly it covers in addition to implicit methods also the wide range of explicit methods typically used to solve conservation laws. For the spatial discretization, we allow for standard choices of numerical fluxes. We use reconstructions of the discrete solution together with the relative entropy stability framework, which leads to error control in the case of smooth solutions. We study under which conditions on the numerical flux the estimate is of optimal order pre-shock. While the estimator we derive is computable and valid post-shock for fixed meshsize, it will blow up as the meshsize tends to zero. This is due to a breakdown of the relative entropy framework when discontinuities develop. We conclude with some numerical benchmarking to test the robustness of the derived estimator