2017/07/15 by Edward Boey, Boey, Edward, Yves Bourgault +3
Computer Science · Engineering · Physics and Astronomy · #65M50 #65M60 #Advanced Mathematical Modeling in Engineering #Advanced Numerical Methods in Computational Mathematics #FOS: Mathematics #Model Reduction and Neural Networks #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.1707.04787
openalex publication_date 2017/07/15 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/01
A residual error estimator is proposed for the energy norm of the error for a scalar reaction-diffusion problem and for the monodomain model used in cardiac electrophysiology. The problem is discretized using P1 finite elements in space, and the backward difference formula of second order (BDF2) in time. The estimator for space makes use of anisotropic interpolation estimates, assuming only minimal regularity. Reliability of the estimator is proven under certain mild assumptions on the convergence of the approximate solution. The monodomain model couples a nonlinear parabolic partial differential equation (PDE) with an ordinary differential equation (ODE) and this setting presents challenges theoretically as well as numerically. A space-time adaptation algorithm is proposed to control the global error, using a non-Euclidean metric for mesh adaptation and a simple method to adjust the time step. Numerical examples are used to verify the reliability and efficiency of the estimator, and to test the adaptive algorithm. The potential gains in efficiency of the proposed algorithm compared to methods using uniform meshes is discussed.