Energy norm error estimates of a hybrid high-order method for the linear parabolic integro-differential equations on general meshes
Provides a stable and convergent numerical scheme for a class of integro-differential equations, but the approach is incremental as it extends existing HHO methods to a specific problem type.
The paper develops a hybrid high-order method for linear parabolic integro-differential equations, achieving error estimates of order O(τ² + h^{k+1}) in discrete norms, validated by numerical experiments on polygonal meshes.
We are concerned in designing a suitable numerical scheme based on the equal-order hybrid high-order (HHO) method for the linear parabolic integro-differential equations. The spatial discretization is made using the equal-order HHO method and subsequently we perform the stability analysis of the corresponding semi-discrete scheme. The convergence results are presented in suitably defined Bochner norms for the semi-discrete problem. Then a second-order temporal discretization is implemented on the time domain using a Crank-Nicolson scheme where the memory term is approximated using a composite trapezoidal quadrature rule. The stability of the resultant complete discrete schemes are analyzed followed by derivation of the error estimates of order $\mathcal{O}(τ^{2}+h^{k+1})$, $k\ge 0$ is the degree of local polynomial approximation, in discrete $l^{2}(0,T;H^{1}(Ω))$ and $l^{\infty}(0,T;H^{1}(Ω))$ like norms. Numerical illustrations are performed on some polygonal meshes validating the theoretical estimates.