Analysis of Dual Consistency for Discontinuous Galerkin Discretizations of Source Terms

Abstract

The effects of dual consistency on discontinuous Galerkin (DG) discretizations of solution and solution gradient dependent source terms are examined. Two common discretizations are analyzed: the standard weighting technique for source terms and the mixed formulation. It is shown that if the source term depends on the first derivative of the solution, the standard weighting technique leads to a dual inconsistent scheme. A straightforward procedure for correcting this dual inconsistency and arriving at a dual consistent discretization is demonstrated. The mixed formulation, where the solution gradient in the source term is replaced by an additional variable that is solved for simultaneously with the state, leads to an asymptotically dual consistent discretization. A priori error estimates are derived to reveal the effect of dual inconsistent discretization on computed functional outputs. Combined with bounds on the dual consistency error, these estimates show that for a dual consistent discretization or the asymptotically dual consistent discretization resulting from the mixed formulation, O(h2p) convergence can be shown for linear problems and linear outputs. For similar but dual inconsistent schemes, only O(hp) can be shown. Numerical results for a one-dimensional test problem confirm that the dual consistent and asymptotically dual consistent schemes achieve higher asymptotic convergence rates with grid refinement than a similar dual inconsistent scheme for both the primal and adjoint solutions as well as a simple functional output.