Abstract:In this paper we analyze a homogeneous parabolic problem with initial data in the space of regular Borel measures. The problem is discretized in time with a discontinuous Galerkin scheme of arbitrary degree and in space with continuous finite elements of orders one or two. We show parabolic smoothing results for the continuous, semidiscrete and fully discrete problems. Our main results are interior $L^\infty$ error estimates for the evaluation at the endtime, in cases where the initial data is supported in a subdomain. In order to obtain these, we additionally show interior $L^\infty$ error estimates for $L^2$ initial data and quadratic finite elements, which extends the corresponding result previously established by the authors for linear finite elements.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the fully discrete pointwise smoothing error estimate for homogeneous parabolic problems with initial values of measures. Specifically, the author focuses on how to evaluate the error between the numerical solution and the exact solution when the initial value condition is a regular Borel measure supported in certain sub - regions.
### Specific description of the problem
The paper studies the following homogeneous parabolic equation:
\[
\begin{cases}
\partial_t v - \Delta v = 0 & \text{in } I\times\Omega, \\
v = 0 & \text{on } I\times\partial\Omega, \\
v(0)=v_0 & \text{in }\Omega,
\end{cases}
\]
where $\Omega\subset\mathbb{R}^N$ is a bounded, convex polygonal or polyhedral region, $I=(0, T]$ is a bounded time interval, the initial condition $v_0$ is a regular Borel measure $v_0\in M(\Omega)$, and is supported in a sub - domain $\Omega_0\subset\Omega$, for example, the Dirac delta function in linear combination $v_0 = \sum_j\beta_j\delta_{x_j}$.
### Research objectives
The main objective of the author is to establish a fully discrete error estimate in the form of:
\[
\|(v - v_{kh})(T)\|_{L^\infty(\Omega_0)}\leq C(\Omega_0, T)\left(k^{2r + 1}+\ell_kh^{s + 1}\right)\|v_0\|_{M(\Omega)},
\]
where $r\geq0$ is the order of time discretization, $s = 1,2$ is the order of space discretization, and $\ell_k$ is a logarithmic term depending on the mesh size $h$ and the maximum time step $k$.
### Main contributions
1. **Fully discrete error estimate**: For discontinuous Galerkin time discretization of any degree and first - order or second - order continuous finite - element space discretization, the author has established a fully discrete error estimate.
2. **Internal $L^\infty$ error estimate**: For $L^2$ initial - value data and quadratic finite elements, the author has also extended previous results to obtain an internal $L^\infty$ error estimate.
3. **No additional smoothness assumptions required**: For the piecewise linear case ($s = 1$), no additional smoothness assumptions are required; while for the quadratic element case ($s = 2$), some additional smoothness assumptions are required, and these assumptions are available on rectangular domains.
### Application background
These results are particularly useful for estimating the initial values of parabolic problems from final - time observations, especially when the final - time observations are made at only a finite number of points, which is more relevant in practical applications.