Stephen Metcalfe

Andrea Cangiani, Emmanuil H. Georgoulis, Irene Kyza et al.

This work is concerned with the development of a space-time adaptive numerical method, based on a rigorous a posteriori error bound, for a semilinear convection-diffusion problem which may exhibit blow-up in finite time. More specifically, a posteriori error bounds are derived in the $L^{\infty}(L^2)+L^2(H^1)$-type norm for a first order in time implicit-explicit (IMEX) interior penalty discontinuous Galerkin (dG) in space discretization of the problem, although the theory presented is directly applicable to the case of conforming finite element approximations in space. The choice of the discretization in time is made based on a careful analysis of adaptive time stepping methods for ODEs that exhibit finite time blow-up. The new adaptive algorithm is shown to accurately estimate the blow-up time of a number of problems, including one which exhibits regional blow-up.

Irene Kyza, Stephen Metcalfe, Thomas Wihler

This work is concerned with the derivation of an a posteriori error estimator for Galerkin approximations to nonlinear initial value problems with an emphasis on finite-time existence in the context of blow-up. The stucture of the derived estimator leads naturally to the development of both h and hp versions of an adaptive algorithm designed to approximate the blow-up time. The adaptive algorithms are then applied in a series of numerical experiments, and the rate of convergence to the blow-up time is investigated.

Andrea Cangiani, Emmanuil H. Georgoulis, Stephen Metcalfe

This work is concerned with the derivation of a robust a posteriori error estimator for a discontinuous Galerkin method discretisation of linear non-stationary convection-diffusion initial/boundary value problems and with the implementation of a corresponding adaptive algorithm. More specifically, we derive a posteriori bounds for the error in the $L^2(H^1)$-type norm for an interior penalty discontinuous Galerkin (dG) discretisation in space and a backward Euler discretisation in time. An important feature of the estimator is robustness with respect to the Péclet number of the problem which is verified in practice by a series of numerical experiments. Finally, an adaptive algorithm is proposed utilising the error estimator. Optimal rate of convergence of the adaptive algorithm is observed in a number of test problems.

Irene Kyza, Stephen Metcalfe

This work is concerned with the development of a space-time adaptive numerical method, based on a rigorous a posteriori error bound, for the semilinear heat equation with a general local Lipschitz reaction term whose solution may blow-up in finite time. More specifically, conditional a posteriori error bounds are derived in the $L^{\infty}L^{\infty}$ norm for a first order in time, implicit-explicit (IMEX), conforming finite element method in space discretization of the problem. Numerical experiments applied to both blow-up and non blow-up cases highlight the generality of our approach and complement the theoretical results.