Convergent approaches for the Dirichlet Monge amp\`ere problem

In this article, we introduce and study three numerical methods for the Dirichlet Monge Amp\`ere equation in two dimensions. The approaches consist in considering new equivalent problems. The latter are discretized by a wide stencil finite difference discretization and monotone schemes are obtained. Hence, we apply the Barles-Souganidis theory to prove the convergence of the schemes and the Damped Newtons method is used to compute the solutions of the schemes. Finally, some numerical results are illustrated.Comment: 14pages, 7figure

Discretization of the 3D Monge-Ampere operator, between Wide Stencils and Power Diagrams

We introduce a monotone (degenerate elliptic) discretization of the Monge-Ampere operator, on domains discretized on cartesian grids. The scheme is consistent provided the solution hessian condition number is uniformly bounded. Our approach enjoys the simplicity of the Wide Stencil method, but significantly improves its accuracy using ideas from discretizations of optimal transport based on power diagrams. We establish the global convergence of a damped Newton solver for the discrete system of equations. Numerical experiments, in three dimensions, illustrate the scheme efficiency

Eigenvalue problems for fully nonlinear elliptic partial differential equations with transport boundary conditions

Fully nonlinear elliptic partial differential equations (PDEs) arise in a number of applications. From mathematical finance to astrophysics, there is a great deal of interest in solving them. Eigenvalue problems for fully nonlinear PDEs with transport boundary conditions are of particular interest as alternative formulations of PDEs that require data to satisfy a solvability condition, which may not be known explicitly or may be polluted by noisy data. Nevertheless, these have not yet been well-explored in the literature. In this dissertation, a convergence framework for numerically solving eigenvalue problems for fully nonlinear PDEs is introduced. In addition, existing two-dimensional methods for nonlinear equations are extended to handle transport boundary conditions and eigenvalue problems. Finally, new techniques are designed to enable appropriate discretization of a large range of fully nonlinear three-dimensional equations

Local Discontinuous Galerkin Methods for One-Dimensional Second Order Fully Nonlinear Elliptic and Parabolic Equations

This paper is concerned with developing accurate and efficient nonstandard discontinuous Galerkin methods for fully nonlinear second order elliptic and parabolic partial differential equations (PDEs) in the case of one spatial dimension. The primary goal of the paper to develop a general framework for constructing high order local discontinuous Galerkin (LDG) methods for approximating viscosity solutions of these fully nonlinear PDEs which are merely continuous functions by definition. In order to capture discontinuities of the first order derivative ux of the solution u, two independent functions q- and q+ are introduced to approximate one-sided derivatives of u. Similarly, to capture the discontinuities of the second order derivative uxx, four independent functions p--,p-+,p+-, and p++ are used to approximate one-sided derivatives of q- and q+. The proposed LDG framework, which is based on a nonstandard mixed formulation of the underlying PDE, embeds a given fully nonlinear problem into a mostly linear system of equations where the given nonlinear differential operator must be replaced by a numerical operator which allows multiple value inputs of the first and second order derivatives ux and uxx. An easy to verify set of criteria for constructing “good” numerical operators is also proposed. It consists of consistency and generalized monotonicity. To ensure such a generalized monotonicity property, the crux of the construction is to introduce the numerical moment in the numerical operator, which plays a critical role in the proposed LDG framework. The generalized monotonicity gives the LDG methods the ability to select the viscosity solution among all possible solutions. The proposed framework extends a companion finite difference framework developed by Feng and Lewis (J Comp Appl Math 254:81–98, 2013) and allows for the approximation of fully nonlinear PDEs using high order polynomials and non-uniform meshes. Numerical experiments are also presented to demonstrate the accuracy, efficiency and utility of the proposed LDG methods

Mixed Interior Penalty Discontinuous Galerkin Methods for Fully Nonlinear Second Order Elliptic and Parabolic Equations in High Dimensions

This article is concerned with developing efficient discontinuous Galerkin methods for approximating viscosity (and classical) solutions of fully nonlinear second-order elliptic and parabolic partial differential equations (PDEs) including the Monge–Ampère equation and the Hamilton–Jacobi–Bellman equation. A general framework for constructing interior penalty discontinuous Galerkin (IP-DG) methods for these PDEs is presented. The key idea is to introduce multiple discrete Hessians for the viscosity solution as a means to characterize the behavior of the function. The PDE is rewritten in a mixed form composed of a single nonlinear equation paired with a system of linear equations that defines multiple Hessian approximations. To form the single nonlinear equation, the nonlinear PDE operator is replaced by the projection of a numerical operator into the discontinuous Galerkin test space. The numerical operator uses the multiple Hessian approximations to form a numerical moment which fulfills consistency and g-monotonicity requirements of the framework. The numerical moment will be used to design solvers that will be shown to help the IP-DG methods select the “correct” solution that corresponds to the unique viscosity solution. Numerical experiments are also presented to gauge the effectiveness and accuracy of the proposed mixed IP-DG methods

Monotone and Consistent discretization of the Monge-Ampere operator

We introduce a novel discretization of the Monge-Ampere operator, simultaneously consistent and degenerate elliptic, hence accurate and robust in applications. These properties are achieved by exploiting the arithmetic structure of the discrete domain, assumed to be a two dimensional cartesian grid. The construction of our scheme is simple, but its analysis relies on original tools seldom encountered in numerical analysis, such as the geometry of two dimensional lattices, and an arithmetic structure called the Stern-Brocot tree. Numerical experiments illustrate the method's efficiency

Monotone discretizations of levelset convex geometric PDEs

We introduce a novel algorithm that converges to level-set convex viscosity solutions of high-dimensional Hamilton-Jacobi equations. The algorithm is applicable to a broad class of curvature motion PDEs, as well as a recently developed Hamilton-Jacobi equation for the Tukey depth, which is a statistical depth measure of data points. A main contribution of our work is a new monotone scheme for approximating the direction of the gradient, which allows for monotone discretizations of pure partial derivatives in the direction of, and orthogonal to, the gradient. We provide a convergence analysis of the algorithm on both regular Cartesian grids and unstructured point clouds in any dimension and present numerical experiments that demonstrate the effectiveness of the algorithm in approximating solutions of the affine flow in two dimensions and the Tukey depth measure of high-dimensional datasets such as MNIST and FashionMNIST.Comment: 42 pages including reference

A Monge-Ampère-solver for free-form reflector design

In this article we present a method for the design of fully free-form reflectors for illumination systems. We derive an elliptic partial differential equation of the Monge-Ampère type for the surface of a reflector that converts an arbitrary parallel beam of light into a desired intensity output pattern. The differential equation has an unusual boundary condition known as the transport boundary condition. We find a convex or concave solution to the equation using a state of the art numerical method. The method uses a nonstandard discretization based on the diagonalization of the Hessian. The discretized system is solved using standard Newton iteration. The method was tested for a circular beam with uniform intensity, a street light, and a uniform beam that is transformed into a famous Dutch painting. The reflectors were verified using commercial ray tracing software

Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

An equation of Monge-Ampère type has, for the first time, been solved numerically on the surface of the sphere in order to generate optimally transported (OT) meshes, equidistributed with respect to a monitor function. Optimal transport generates meshes that keep the same connectivity as the original mesh, making them suitable for r-adaptive simulations, in which the equations of motion can be solved in a moving frame of reference in order to avoid mapping the solution between old and new meshes and to avoid load balancing problems on parallel computers. The semi-implicit solution of the Monge-Ampère type equation involves a new linearisation of the Hessian term, and exponential maps are used to map from old to new meshes on the sphere. The determinant of the Hessian is evaluated as the change in volume between old and new mesh cells, rather than using numerical approximations to the gradients. OT meshes are generated to compare with centroidal Voronoi tesselations on the sphere and are found to have advantages and disadvantages; OT equidistribution is more accurate, the number of iterations to convergence is independent of the mesh size, face skewness is reduced and the connectivity does not change. However anisotropy is higher and the OT meshes are non-orthogonal. It is shown that optimal transport on the sphere leads to meshes that do not tangle. However, tangling can be introduced by numerical errors in calculating the gradient of the mesh potential. Methods for alleviating this problem are explored. Finally, OT meshes are generated using observed precipitation as a monitor function, in order to demonstrate the potential power of the technique