Numerical Treatment of Non-local Integral Operators in the Framework of Evolutionary Equations
Sebastian Franz, Sascha Trostorff
TL;DR
This work extends the abstract evolutionary equation framework to include a non-local integral operator $T_K$ with kernel $K$, proving well-posedness and causality under a kernel-bound condition in a weighted Hilbert space $L^{2,\rho}(\mathbb{R};H)$. It then develops a space-time Galerkin discretisation (conforming finite elements in space and discontinuous Galerkin in time) and provides a rigorous convergence analysis, showing $\|U-U^{\tau,h}\|_{\rho} \le C(h^k + T^{1/2}\tau^{q+1})$ under suitable regularity and kernel assumptions. Numerical experiments with smooth and weakly singular kernels corroborate the theoretical rates and illustrate practical considerations for history-term computations and adaptive quadrature. The results lay a solid foundation for numerically solving Volterra/Fredholm-type integro-differential problems within the Picard framework, with potential extensions to time-derivative kernels and more general non-local operators.
Abstract
Using the theory of evolutionary equations, we consider abstract differential equations including non-local integral operators. After providing a condition for the well-posedness of the addressed equation we consider a numerical method of approximating its solution. We provide convergence proofs under conditions on the kernel of the integral operator and the solution and finish the paper with some simulation results.
