Fast Operator-Splitting Methods for Nonlinear Elliptic Equations

Nonlinear elliptic problems arise in many fields, including plasma physics, astrophysics, and optimal transport. In this article, we propose a novel operator-splitting/finite element method for solving such problems. We begin by introducing an auxiliary function in a new way for a semilinear elliptic partial differential equation, leading to the development of a convergent operator-splitting/finite element scheme for this equation. The algorithm is then extended to fully nonlinear elliptic equations of the Monge-Ampère type, including the Dirichlet Monge-Ampère equation and Pucci’s equation. This is achieved by reformulating the fully nonlinear equations into forms analogous to the semilinear case, enabling the application of the proposed splitting algorithm. In our implementation, a mixed finite element method is used to approximate both the solution and its Hessian matrix. Numerical experiments show that the proposed method outperforms existing approaches in efficiency and accuracy, and can be readily applied to problems defined on domains with curved boundaries.