Abstract:Neural operators, which learn mappings between the function spaces, have been applied to solve boundary value problems in various ways, including learning mappings from the space of the forcing terms to the space of the solutions with the substantial requirements of data pairs. In this work, we present a data-free neural operator integrated with physics, which learns the Green kernel directly. Our method proceeds in three steps: 1. The governing equations for the Green's function are reformulated into an interface problem, where the delta Dirac function is removed; 2. The interface problem is embedded in a lifted space of higher-dimension to handle the jump in the derivative, but still solved on a two-dimensional surface without additional sampling cost; 3. Deep neural networks are employed to address the curse of dimensionality caused by this lifting operation. The approximate Green's function obtained through our approach is then used to construct preconditioners for the linear systems allowed by its mathematical properties. Furthermore, the spectral bias of it revealed through both theoretical analysis and numerical validation contrasts with the smoothing effects of traditional iterative solvers, which motivates us to propose a hybrid iterative method that combines these two solvers. Numerical experiments demonstrate the effectiveness of our approximate Green's function in accelerating iterative methods, proving fast convergence for solving indefinite problems even involving discontinuous coefficients.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: how to use neural networks to learn Green's function to accelerate the solution of indefinite boundary value problems. Specifically, the author proposes a data - free neural operator method, which directly learns the Green kernel by combining physical information, thereby constructing a pre - conditioner for linear systems, and proposes a hybrid iterative method to accelerate the solution of indefinite problems with discontinuous coefficients.
### Main problems and solutions
1. **Problem background**:
- Indefinite boundary value problems are very important in numerical simulation, but traditional iterative methods such as Jacobi, Gauss - Seidel, etc. face challenges when dealing with indefinite systems, especially when the coefficient function has a jump at the interface.
- Neural operators have been applied to solve boundary value problems, but most methods rely on a large number of data pairs and are difficult to handle high - dimensional problems.
2. **Proposed method**:
- **Data - free neural operator**: The author proposes a neural operator without data pairs that directly learns the Green kernel. This method is divided into three steps:
1. Reformulate the governing equation of the Green function as an interface problem and remove the Dirac delta function.
2. Embed the interface problem into a higher - dimensional space. By introducing auxiliary variables to handle the jump of derivatives, but still keep solving on a two - dimensional surface without increasing the sampling cost.
3. Use a deep neural network to solve the curse - of - dimensionality problem caused by the dimension - raising operation.
3. **Application and effect**:
- Use the obtained approximate Green function to construct a pre - conditioner to accelerate the GMRES method.
- Propose a hybrid iterative method that combines the neural Green function and the damped Jacobi method to achieve fast convergence of all frequency components.
- Numerical experiments show that this method performs well in solving indefinite problems (including problems with discontinuous coefficients) and can converge quickly.
### Formula summary
- **Integral form of Green's function**:
\[
u(x)=\int_{\Omega} G(x, y) f(y) \, dy
\]
where \( G(x, y) \) is the Green function and \( f(y) \) is the source term.
- **Equation after dimension - raising**:
\[
\mathcal{L}_y \bar{G}(x, y, r(x, y)) = 0, \quad y \in \Omega\setminus\Gamma
\]
\[
\bar{G}(x, y, r(x, y)) = 0, \quad y \in \partial\Omega
\]
\[
2\frac{\partial\bar{G}(x, y, r(x, y))}{\partial r}=-c^{-1}(x), \quad y \in \Gamma
\]
- **Loss function**:
\[
L(\theta)=L_{\Omega\setminus\Gamma}(\bar{G})+\beta_{\partial\Omega} L_{\partial\Omega}(\bar{G})+\beta_{\Gamma} L_{\Gamma}(\bar{G})+\beta_{\text{sym}} L_{\text{sym}}^{\Omega}(\bar{G})
\]
Through this method, the author successfully overcomes the limitations of traditional methods in dealing with indefinite boundary value problems and demonstrates its effectiveness in complex geometries and high - dimensional problems.