Error analysis for parabolic optimal control problems with measure data in a nonconvex polygonal domain

Pratibha Shakya

2024-03-11

Abstract:This paper considers the finite element approximation to parabolic optimal control problems with measure data in a nonconvex polygonal domain. Such problems usually possess low regularity in the state variable due to the presence of measure data and the nonconvex nature of the domain. The low regularity of the solution allows the finite element approximations to converge at lower orders. We prove the existence, uniqueness and regularity results for the solution to the control problem satisfying the first order optimality condition. For our error analysis we have used piecewise linear elements for the approximation of the state and co-state variables, whereas piecewise constant functions are employed to approximate the control variable. The temporal discretization is based on the implicit Euler scheme. We derive both a priori and a posteriori error bounds for the state, control and co-state variables. Numerical experiments are performed to validate the theoretical rates of convergence.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the error analysis in the finite - element approximation of parabolic optimal control problems with measure data in non - convex polygonal domains. Specifically, the paper focuses on the following aspects:

1. **Low - regularity problem**: Due to the existence of measure data and the non - convex nature of the domain in the problem, the solution of the state variable usually has low regularity. This causes the finite - element approximation to converge at a lower order.

2. **Existence, uniqueness and regularity results**: The paper proves the existence, uniqueness and regularity results of the solution of the control problem satisfying the first - order optimality conditions.

3. **Error estimates**: The paper derives a priori and a posteriori error bounds for the state variable, the control variable and the co - state variable. Specifically:
- Piecewise linear functions are used to approximate the state variable and the co - state variable.
- Piecewise constant functions are used to approximate the control variable.
- Time discretization is based on the implicit Euler method.

4. **Numerical experiments**: Numerical experiments are carried out to verify the theoretical convergence rate.

### Main contributions
- **Theoretical analysis**: It provides detailed a priori and a posteriori error estimates, which are of great significance for understanding and improving numerical methods.
- **Numerical verification**: The theoretical results are verified by numerical experiments, demonstrating the effectiveness of the method.

### Formula summary
- **Objective functional**:
\[
J(y, u)=\frac{1}{2}\int_{0}^{T}\|y - y_d\|_{L^{2}(\Omega)}^{2}dt+\frac{\Lambda}{2}\int_{0}^{T}\|u\|_{L^{2}(\Omega)}^{2}dt
\]

- **State equation**:
\[
\begin{cases}
\frac{\partial y}{\partial t}-\Delta y = \sigma\tau+u&\text{in }\Omega_T,\\
y = 0&\text{on }\Gamma_T,\\
y(\cdot,0)=y_0&\text{in }\Omega.
\end{cases}
\]

- **A priori error estimate**:
\[
\|y - Y_h(u)\|_{L^{2}(0,T;L^{2}(\Omega))}\leq C\left(h^{2s}k^{- 1/2}+k^{1/2}+h^{s}\right)\left(\|y_0\|_{L^{2}(\Omega)}+\|\sigma\|_{L^{\infty}(0,T;L^{2}(\Omega))}\|\tau\|_{B[0,T]}+\|u\|_{L^{2}(0,T;L^{2}(\Omega))}\right)
\]

- **A posteriori error estimate**:
\[
\|y - Y_h\|_{L^{2}(0,T;L^{2}(\Omega))}+\|u - u_\rho\|_{L^{2}(0,T;L^{2}(\Omega))}+\|\phi - \Phi_h\|_{L^{2}(0,T;L^{2}(\Omega))}\leq C h^{s}
\]

### Conclusion
This paper provides theoretical a priori and a posteriori error estimates through a detailed analysis of the finite - element approximation of parabolic optimal control problems with measure data in non - convex polygonal domains, and verifies the effectiveness of these estimates through numerical experiments. These results are of great significance for understanding and improving the numerical methods for such problems.