Abstract:In the optimization of convex domains under a PDE constraint numerical difficulties arise in the approximation of convex domains in $\mathbb{R}^3$. Previous research used a restriction to rotationally symmetric domains to reduce shape optimization problems to a two-dimensional setting. In the current research, two approaches for the approximation in $\mathbb{R}^3$ are considered. First, a notion of discrete convexity allows for a nearly convex approximation with polyhedral domains. An alternative approach is based on the recent observation that higher order finite elements can approximate convex functions conformally. As a second approach these results are used to approximate optimal convex domains with isoparametric convex domains. The proposed algorithms were tested on shape optimization problems constrained by a Poisson equation and both algorithms achieved similar results.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the difficulties encountered in the numerical approximation of convex domains in three - dimensional space $\mathbb{R}^3$, especially in shape optimization problems under partial differential equation (PDE) constraints. Specifically, the paper focuses on how to effectively approximate and optimize three - dimensional convex domains while maintaining convexity constraints.
#### Main problems
1. **Numerical approximation of convex domains**: In high - dimensional spaces, especially in $\mathbb{R}^3$, the numerical approximation of convex domains is very difficult. The traditional linear finite - element method cannot guarantee a consistent approximation of convex functions, which leads to the inaccuracy of numerical solutions.
2. **Limitations of rotation - symmetric domains**: Previous studies simplified the problem to a two - dimensional problem by restricting it to rotation - symmetric domains. However, this restriction is generally unreasonable, and its impact is difficult to evaluate.
3. **Handling of convexity constraints**: In order to handle convexity constraints in three - dimensional space, the paper proposes two methods:
- **Discrete convexity**: A method based on discrete geometry that allows approximation using polyhedral domains.
- **Isoparametric convexity**: Using high - order finite - element methods to approximate convex functions and further applying them to the approximation of convex domains.
#### Specific form of the optimization problem
The paper considers the following shape optimization problem with PDE constraints:
\[
\begin{cases}
\text{Minimize} & \int_{\Omega} j(x, u(x), \nabla u(x)) \, dx \\
\text{w.r.t.} & \Omega \subset \mathbb{R}^3, \quad u \in H_0^1(\Omega) \\
\text{s.t.} & -\Delta u = f \quad \text{in } \Omega \\
& \Omega \subset Q \quad \text{(convex and open)}
\end{cases}
\]
where $Q \subset \mathbb{R}^3$ is a bounded, open and convex region, $j: Q \times \mathbb{R} \times \mathbb{R}^3 \to \mathbb{R}$ is a suitable Carathéodory function, and $f \in L^2(Q)$.
### Solutions
1. **Discrete convexity method**: A discrete convexity condition is defined, which is suitable for the P1 finite - element method. Through this method, almost convex polyhedral domains can be constructed for approximation.
2. **Isoparametric convexity method**: High - order finite - element methods (such as quadratic isoparametric elements) are used to approximate convex functions and further applied to the approximation of convex domains. This method can maintain convexity constraints more accurately.
### Experimental results
The paper verifies the effectiveness of these two methods through numerical experiments. The experimental results show that both methods can achieve similar results when dealing with shape optimization problems under Poisson equation constraints, and can achieve effective optimization while maintaining convexity constraints.
### Summary
The main contribution of this paper is to propose and verify two effective methods for handling convexity constraints in three - dimensional space, thereby solving the problems encountered by traditional methods in approximating convex domains in high - dimensional spaces. These methods not only improve the accuracy of numerical solutions but also provide new ideas and tools for subsequent research.