Numerical methods for fully nonlinear degenerate diffusions

Abstract

We propose finite difference methods for degenerate fully nonlinear elliptic equations and prove the convergence of the schemes. Our focus is on the pure equation and a related free boundary problem of transmission type. The cornerstone of our argument is a regularisation procedure. It decouples the degeneracy term from the elliptic operator driving the diffusion process. In the free boundary setting, a two-parameters regularisation entails new, genuine difficulties. To bypass them, we resort to the intrinsic properties of the regularised problem and a fixed-point type of argument. We present numerical experiments supporting our theoretical results. Our methods are flexible, and our approach can be extended to a broader class of non-variational problems.