# An efficient neural-network and finite-difference hybrid method for elliptic interface problems with applications

arXiv: https://arxiv.org/abs/2210.05523 · HTML: https://mlsanity.com/paper/2210.05523/

## Key facts

- Cited by 2 papers in the map (data through 14 Oct 2026).
- Builds on 4 of its 29 references that are in the map.
- First posted on arXiv 11 Oct 2022.

## Builds on (most cited first, 4 of 4)

- [Automatic differentiation in machine learning: a survey](https://mlsanity.com/paper/1502.05767/) (2015)
- [The Deep Ritz method: A deep learning-based numerical algorithm for solving variational problems](https://mlsanity.com/paper/1710.00211/) (2017)
- [A Discontinuity Capturing Shallow Neural Network for Elliptic Interface Problems](https://mlsanity.com/paper/2106.05587/) (2021)
- [A Shallow Ritz Method for Elliptic Problems with Singular Sources](https://mlsanity.com/paper/2107.12013/) (2021)

## Built on by (most cited first, 2 of 2)

- [The Hard-Constraint PINNs for Interface Optimal Control Problems](https://mlsanity.com/paper/2308.06709/) (2023)
- [Dirichlet-Neumann learning algorithm for solving elliptic interface problems](https://mlsanity.com/paper/2301.07361/) (2023)

## Abstract (arXiv)

A new and efficient neural-network and finite-difference hybrid method is developed for solving Poisson equation in a regular domain with jump discontinuities on embedded irregular interfaces. Since the solution has low regularity across the interface, when applying finite difference discretization to this problem, an additional treatment accounting for the jump discontinuities must be employed. Here, we aim to elevate such an extra effort to ease our implementation by machine learning methodology. The key idea is to decompose the solution into singular and regular parts. The neural network learning machinery incorporating the given jump conditions finds the singular solution, while the standard five-point Laplacian discretization is used to obtain the regular solution with associated boundary conditions. Regardless of the interface geometry, these two tasks only require supervised learning for function approximation and a fast direct solver for Poisson equation, making the hybrid method easy to implement and efficient. The two- and three-dimensional numerical results show that the present hybrid method preserves second-order accuracy for the solution and its derivatives, and it is comparable with the traditional immersed interface method in the literature. As an application, we solve the Stokes equations with singular forces to demonstrate the robustness of the present method.