A neural network approach to learning solutions of a class of elliptic variational inequalities

Amal Alphonse,Michael Hintermüller,Alexander Kister,Chin Hang Lun,Clemens Sirotenko
2024-11-28
Abstract:We develop a weak adversarial approach to solving obstacle problems using neural networks. By employing (generalised) regularised gap functions and their properties we rewrite the obstacle problem (which is an elliptic variational inequality) as a minmax problem, providing a natural formulation amenable to learning. Our approach, in contrast to much of the literature, does not require the elliptic operator to be symmetric. We provide an error analysis for suitable discretisations of the continuous problem, estimating in particular the approximation and statistical errors. Parametrising the solution and test function as neural networks, we apply a modified gradient descent ascent algorithm to treat the problem and conclude the paper with various examples and experiments. Our solution algorithm is in particular able to easily handle obstacle problems that feature biactivity (or lack of strict complementarity), a situation that poses difficulty for traditional numerical methods.
Optimization and Control,Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: **Using the neural network method to solve the obstacle problem of a class of elliptic variational inequalities**. Specifically, the authors developed a weak adversarial approach to solve the obstacle problem through neural networks. The obstacle problem is a class of elliptic variational inequality problems, usually expressed as: \[ \text{find } u \in K: \langle Au - f, u - v \rangle_{V^*, V} \leq 0 \quad \forall v \in K, \] where: - \( V = H^1(\Omega) \) is the standard Sobolev space defined on the bounded Lipschitz domain \( \Omega \subset \mathbb{R}^n \). - \( \langle \cdot, \cdot \rangle_{V^*, V} \) represents the dual pairing between \( V \) and its topological dual \( V^* \). - The constraint set \( K \) is defined as: \[ K = \left\{ u \in H^1(\Omega) \mid u \geq \psi \text{ in } \Omega, u = h \text{ on } \partial \Omega \right\}, \] where \( h \in H^{1/2}(\partial \Omega) \) is the given boundary data, \( \psi \in H^1(\Omega) \) is the obstacle function satisfying \( \psi \leq h \), and \( f \in L^2(\Omega) \) is the given source term. - The operator \( A: K \subset V \to V^* \) is assumed to be Lipschitz continuous and coercive, that is, there exist constants \( C_a, C_b > 0 \) such that: \[ \|Au - Av\|_{V^*} \leq C_b \|u - v\|_V \quad \forall u, v \in K, \] \[ \langle Au - Av, u - v \rangle_{V^*, V} \geq C_a \|u - v\|^2_V \quad \forall u, v \in K. \] To solve this problem, the authors adopted the following steps: 1. **Rewrite the obstacle problem as a minimax problem**: By introducing the regularized gap function, the obstacle problem is transformed into a minimax optimization problem. This method is applicable even when the elliptic operator is not symmetric. 2. **Use neural network to parameterize the solution and test functions**: Parameterize the solution \( u \) and the test function \( v \) as neural networks, and apply the modified gradient descent - ascent algorithm to numerically solve the problem. 3. **Provide error analysis**: Estimate the error of the appropriate discretization of the continuous problem, with particular attention to the approximation error and the statistical error. 4. **Handle the biactivity problem**: This method can easily handle cases with biactivity or lack of strict complementarity, which is a challenge for traditional numerical methods. In summary, the main goal of this paper is to use the neural network method to solve the obstacle problem of elliptic variational inequalities, and it provides detailed theoretical analysis, discretization methods and numerical experiment results.