| Abstract: |
| Nonlinear elliptic problems arise in many fields, including plasma physics, astrophysics, and optimal transport. In this talk, we present our recent work on 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-Ampere type, including the Dirichlet Monge-Ampere 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. |
|