Goal Oriented Time Adaptivity Using Local Error Estimates

Abstract

We consider initial value problems (IVPs) where we are interested in a quantity of interest (QoI) that is the integral in time of a functional of the solution. For these, we analyze goal oriented time adaptive methods that use only local error estimates. A local error estimate and timestep controller for step-wise contributions to the QoI are derived. We prove convergence of the error in the QoI for tolerance to zero under a controllability assumption. By analyzing global error propagation with respect to the QoI, we can identify possible issues and make performance predictions. Numerical tests verify these results. We compare performance with classical local error based time-adaptivity and a posteriori based adaptivity using the dual-weighted residual (DWR) method. For dissipative problems, local error based methods show better performance than DWR and the goal oriented method shows good results in most examples, with significant speedups in some cases.

A time adaptive multirate Dirichlet–Neumann waveform relaxation method for heterogeneous coupled heat equations

“…As can be seen in the numerical results in Section 8.3, this does not lead to a loss of order. Other options for constructing partitioned time‐integration methods based on SDIRK2 are discussed in [33, Chap. 5].…”

Development of a Balanced Adaptive Time-Stepping Strategy Based on an Implicit JFNK-DG Compressible Flow Solver

“…The accumulation of these errors, which we refer to as global errors, is more difficult to estimate and will not be considered in this paper. Studies considering the global errors can be found in reference [19]. We recall that the numerical solution at the time step n + 1 is calculated by ( 19)…”

“…Noventa et al developed a similar time step adaptation strategy in unsteady incompressible turbulent flow simulations [21], which is further tested in compressible simulations focusing on the comparison of different step size controllers [22]. The adaptive method based on local temporal error was further developed for goal oriented time step adaptation and compared with a method based on global error estimate [19]. However, most of them such as, the studies in references [4,5,21,22], consider the temporal discretization methods separately as an ODE solver using the method of lines but do not take into account specific physical properties of the problem and specific spatial discretization adopted.…”

Goal-Oriented Adaptive Space-Time Finite Element Methods for Regularized Parabolic P-Laplace Problems

“…[6,14,20,15] and the references therein. One possible way to overcome this requires checkpoint techniques as presented in [19]. In this paper, we instead treat time as just another variable and perform an all-at-once discretization of the complete space-time cylinder 𝑄.…”