S. Saberi,L. Zhang,A. Vogel
Abstract:The weak imposition of essential boundary conditions is an integral aspect of unfitted finite element methods, where the physical boundary does not in general coincide with the computational domain. In this regard, the symmetric Nitsche's method is a powerful technique that preserves the symmetry and variational consistency of the unmodified weak formulation. The stabilization parameter in Nitsche's method plays a crucial role in the stability of the resultant formulation, whose estimation is computationally intensive and dependent on the particular cut configuration using the conventional eigenvalue-based approach. In this work, we employ as model problem the finite cell method in which the need for the generation of a boundary-conforming mesh is circumvented by embedding the physical domain in a, typically regular, background mesh. We propose a data-driven estimate based on machine learning methods for the estimation of the stabilization parameter in Nitsche's method that offers an efficient constant-complexity alternative to the eigenvalue-based approach independent of the cut configuration. It is shown, using numerical benchmarks, that the proposed method can estimate the stabilization parameter accurately and is by far more computationally efficient. The data-driven estimate can be integrated into existing numerical codes with minimal modifications and thanks to the wide adoption of accelerators such as GPUs by machine learning frameworks, can be used with virtually no extra implementation cost on GPU devices, further increasing the potential for computational gains over the conventional eigenvalue-based estimate.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the computational complexity problem in the estimation of the stabilization parameter in the Nitsche method. Specifically, the Nitsche method is a powerful technique for weakly imposing essential boundary conditions, especially suitable for unfitted finite element methods, where the physical boundary usually does not coincide with the computational domain. The stabilization parameter plays a crucial role in the Nitsche method, ensuring the stability of the resulting formula. However, the traditional eigenvalue - based method for estimating the stabilization parameter is computationally intensive and depends on the specific cutting configuration.
To solve this problem, the authors propose a data - driven method based on machine learning to estimate the stabilization parameter in the Nitsche method. This method has the following advantages:
1. **Efficiency**: Compared with the traditional eigenvalue method, the data - driven method has a constant complexity (O(1)), independent of the cutting configuration.
2. **Accuracy**: Numerical benchmark tests show that the data - driven method can accurately estimate the stabilization parameter.
3. **Easy integration**: It can be easily integrated into existing numerical codes with almost no additional implementation cost, especially when using GPU accelerators.
### Specific content
- **Background and motivation**:
- Unfitted finite element methods (such as the finite cell method) avoid the need to generate boundary - conforming meshes, but require weakly imposing boundary conditions.
- The Nitsche method maintains symmetry and variational consistency when weakly imposing boundary conditions, but the choice of its stabilization parameter is crucial.
- **Problem description**:
- The estimation of the stabilization parameter is usually done by solving the generalized eigenvalue problem, which is very expensive computationally.
- The computational complexity of the traditional method increases as the proportion of the cut part decreases, resulting in inefficiency.
- **Solution**:
- A data - driven method based on machine learning is proposed to estimate the stabilization parameter.
- The neural network is used to approximate the stabilization parameter, reducing the computational complexity and improving the efficiency.
- **Experimental verification**:
- Numerical experiments are carried out using the Poisson equation as a model problem.
- The results show that the data - driven method is superior to the traditional method in both accuracy and computational efficiency.
### Formula representation
The formulas involved in the paper include:
- The strong form of the Poisson equation:
\[
-\nabla^2 u = f \quad \text{in} \; \Omega,
\]
\[
u = g \quad \text{on} \; \Gamma_D,
\]
\[
\nabla u \cdot n = h \quad \text{on} \; \Gamma_N.
\]
- The weak form:
\[
\int_{\Omega} \nabla v \cdot \nabla u \, dx - \int_{\Gamma_D} v (\nabla u \cdot n) \, ds = \int_{\Omega} v f \, dx + \int_{\Gamma_N} v h \, ds.
\]
- The weak form of the Nitsche method:
\[
\int_{\Omega_e} \alpha \nabla v \cdot \nabla u \, dx - \int_{\Gamma_D} v (\nabla u \cdot n) \, ds - \int_{\Gamma_D} (u - g) (\nabla v \cdot n) \, ds + \int_{\Gamma_D} \lambda v (u - g) \, ds = \int_{\Omega_e} \alpha v f \, dx + \int_{\Gamma_N} v h \, ds.
\]
- The generalized eigenvalue problem: