A finite element method for the Monge-Ampère equation with transport boundary conditions

Abstract

We address the numerical solution via Galerkin type methods of the Monge-Ampère equation with transport boundary conditions arising in optimal mass transport, geometric optics and computational mesh or grid movement techniques. This fully nonlinear elliptic problem admits a linearisation via a Newton-Raphson iteration, which leads to an oblique derivative boundary value problem for elliptic equations in nondivergence form. We discretise these by employing the nonvariational finite element method, which lead to empirically observed optimal convergence rates, provided recovery techinques are used to approximate the gradient and the Hessian of the unknown functions. We provide extensive numerical testing to illustrate the strengths of our approach and the potential applications in optics and mesh movement.

2 (x6) figures

References in corpus (2)

• The Geometry of r-adaptive meshes generated using Optimal Transport Methods
• A DGFEM for Uniformly Elliptic Two Dimensional Oblique Boundary Value Problems

Cited by in corpus (3)

• Discontinuous Galerkin and $C^0$-IP finite element approximation of periodic Hamilton--Jacobi--Bellman--Isaacs problems with application to numerical homogenization
• Finite element theory on curved domains with applications to DGFEMs
• Adaptive $C^0$ interior penalty methods for Hamilton-Jacobi-Bellman equations with Cordes coefficients