A posteriori error analysis for a continuous space-time finite element method for a hyperbolic integro-differential equation

Fardin Saedpanah

DOI: https://doi.org/10.48550/arXiv.1205.0159

2012-11-15

Abstract:An integro-differential equation of hyperbolic type, with mixed boundary conditions, is considered. A continuous space-time finite element method of degree one is formulated. A posteriori error representations based on space-time cells is presented such that it can be used for adaptive strategies based on dual weighted residual methods. A posteriori error estimates based on weighted global projections and local projections are also proved.

Numerical Analysis

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the problems related to the numerical solution of the **hyperbolic integro - differential equation** and its error estimation. Specifically, the paper mainly focuses on the following points:

1. **Problem background**:
- The paper studies an initial - boundary value problem in the form of:
\[
\ddot{u}(x,t)-\nabla\cdot\sigma(u;x,t)=f(x,t)\quad\text{in }\Omega\times(0,T),
\]
where \(u(x,t)\) is the displacement vector, \(f(x,t)\) and \(g(x,t)\) represent the volume load and surface load respectively, and the stress \(\sigma(t)\) is defined as:
\[
\sigma(t)=\sigma_0(t)-\int_0^tK(t - s)\sigma_0(s)\,ds,
\]
where \(\sigma_0(t)=2\mu_0\epsilon(t)+\lambda_0\text{Tr}(\epsilon(t))I\), and \(\epsilon\) is the strain tensor.

2. **Continuous space - time finite element method**:
- The paper proposes a numerical solution method based on the continuous space - time finite element method for solving the above - mentioned hyperbolic integro - differential equation. This method uses first - order polynomials for spatial and temporal discretization.

3. **A posteriori error analysis**:
- The main contribution of the paper is to propose **a posteriori error representations**, which can be used for adaptive strategies based on the dual weighted residual (DWR) method. The paper presents three different a posteriori error representations and proves two a posteriori error estimates based on global and local projections.

4. **Goal - oriented adaptive discretization**:
- By introducing the dual problem, the paper provides a framework that enables goal - oriented adaptive discretization through the dual weighted residual method. This helps to allocate computational resources more effectively during the calculation process and improve the accuracy of the numerical solution.

5. **Influence of the convolution kernel**:
- The paper also discusses the influence of the convolution kernel \(K(t)\), especially how to handle the problem of the amount of data that needs to be stored as the time step increases when the convolution kernel has weakly singular properties. The paper mentions using methods such as Lubich convolution quadrature to overcome this challenge.

### Summary

In general, this paper is dedicated to developing an efficient numerical method for solving hyperbolic integro - differential equations and provides a method for evaluating and improving the accuracy of numerical solutions through a posteriori error analysis. In addition, the paper also explores how to effectively manage computational resources in the presence of convolution terms.