Abstract:This paper analyzes two eXtended finite element methods (XFEMs) for linear quadratic optimal control problems governed by Poisson equation in non-convex domains. We follow the variational discretization concept to discretize the continuous problems, and apply an XFEM with a cut-off function and a classic XFEM with a fixed enrichment area to discretize the state and co-state equations. Optimal error estimates are derived for the state, co-state and control. Numerical results confirm our theoretical results.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: in the linear - quadratic optimal control problem governed by the Poisson equation in non - convex domains, how to effectively improve the accuracy of the finite element method. Specifically, due to the singularities in non - convex domains, the weak solutions are usually not in the \(H^2(\Omega)\) space, which will reduce the accuracy of the standard finite element method. For this reason, the author introduced the extended finite element method (XFEM) to deal with the singularities in non - convex domains and improve the accuracy of the numerical solutions.
### Main research contents:
1. **Problem background**:
- Consider a linear - quadratic optimal control problem governed by the Poisson equation. The goal is to minimize the functional \(J(y, u)\), where \(y\) is the state variable and \(u\) is the control variable.
- This problem is constrained by the Poisson equation, and the control variable \(u\) satisfies certain upper and lower bound constraints.
2. **Challenges**:
- In non - convex domains, due to the existence of boundary angles, the regularity of the solution is low, and the convergence speed of the standard finite element method will be affected.
- Traditional improvement methods such as graded meshes and singular function methods have problems such as complex construction or high computational cost in practical applications.
3. **Solutions**:
- Introduce the extended finite element method (XFEM). By adding basis functions characterizing singularities to the standard finite element space, the complexity of body - fitted meshes is avoided.
- Use the concept of variational discretization to discretize the continuous problem. Apply XFEM with cutoff functions and XFEM with classical fixed enrichment regions to discretize the state equation and the co - state equation respectively.
4. **Theoretical results**:
- Derive the optimal error estimates of the state, co - state and control variables.
- Verify the theoretical results through numerical experiments, and show the effectiveness and superiority of the proposed method.
### Mathematical model:
The optimal control problem has the form:
\[
\min J(y, u)=\frac{1}{2}\int_{\Omega}(y - y_d)^2\,dx+\frac{\alpha}{2}\int_{\Omega}u^2\,dx
\]
Subject to the following Poisson equation constraints:
\[
\begin{cases}
-\Delta y = u + f&\text{in }\Omega,\\
y = 0&\text{on }\partial\Omega,
\end{cases}
\]
And the control variable \(u\) satisfies the constraint:
\[
u_0\leq u\leq u_1\quad\text{a.e. on }\Omega.
\]
### Conclusion:
By introducing the extended finite element method and combining the concept of variational discretization, this paper successfully solves the numerical difficulties caused by low regularity in the linear - quadratic optimal control problem governed by the Poisson equation in non - convex domains, and improves the accuracy and convergence speed of the numerical solutions.