Duality based error estimation in the presence of discontinuities

摘要

Goal-oriented mesh adaptation, in particular using the dual-weighted residual (DWR) method, is known in many cases to produce very efficient meshes. For obtaining such meshes the (numerical) solution of an adjoint problem is needed to weight the residuals appropriately with respect to their relevance for the overall error. For hyperbolic problems already the weak primal problem requires in general an additional entropy condition to assert uniqueness of solutions; this difficulty is also reflected when considering adjoints to hyperbolic problems involving discontinuities where again an additional requirement (reversibility) is needed to select appropriate solutions.Within this article, an approach to the DWR method for hyperbolic problems based on an artificial viscosity approximation is proposed. It is discussed why the proposed method provides a well-posed dual problem, while a direct, formal, application of the dual problem does not. Moreover, we will discuss a further, novel, approach in which the forward problem need not be modified, thus allowing for an unchanged forward solution. The latter procedure introduces an additional residual term in the error estimation, accounting for the inconsistency between primal and dual problem.Finally, the effectivity of the extended error estimator, assessing the global error by a suitable functional of interest, is tested numerically; and the advantage over a formal estimator approach is demonstrated. (C) 2019 IMACS. Published by Elsevier B.V. All rights reserved.