A posteriori error estimates robust with respect to nonlinearities and final time.

Martin Vohralík: Thursday, 25th May at 11:00


A posteriori estimates enable us to certify the error committed in a numerical simulation. In particular, the equilibrated flux reconstruction technique yields a guaranteed error upper bound, where the flux obtained by local postprocessing is of independent interest since it is always locally conservative. In this talk, we tailor this methodology to model nonlinear and time-dependent problems to obtain estimates that are robust, i.e., of quality independent of the strength of the nonlinearities and the final time. These estimates include and build on common iterative linearization schemes such as Zarantonello, Picard, Newton, or M- and L-ones. We first consider steady problems and conceive two settings: we either augment the energy difference by the discretization error of the current linearization step, or we design iteration-dependent norms that feature weights given by the current iterate. We then turn to unsteady problems. Here we first consider the linear heat equation and finally move to the Richards one, which is doubly nonlinear and exhibits both parabolic–hyperbolic and parabolic–elliptic degeneracies. Robustness with respect to the final time and local efficiency in both time and space are addressed here. Numerical experiments illustrate the theoretical findings all along the presentation.

Comments are closed.