The New Jersey Institute of Technology's
Electronic Theses & Dissertations Project

Title:	Numerical methods for optimal transport and optimal information transport on the sphere
Author:	Turnquist, Axel G. R.
View Online:	njit-etd2022-027 (xiv, 196 pages ~ 9.1 MB pdf)
Department:	Department of Mathematical Sciences
Degree:	Doctor of Philosophy
Program:	Mathematical Sciences
Document Type:	Dissertation
Advisory Committee:	Hamfeldt, Brittany Froese (Committee chair) Boubendir, Yassine (Committee member) Muratov, Cyrill B. (Committee member) Shirokoff, David (Committee member) Lai, Rongjie (Committee member)
Date:	2022-05
Keywords:	Diffeomorphic mapping Differential geometry Numerical analysis Optics Optimal transport Partial differential equations
Availability:	Unrestricted
Abstract:	The primary contribution of this dissertation is in developing and analyzing efficient, provably convergent numerical schemes for solving fully nonlinear elliptic partial differential equation arising from Optimal Transport on the sphere, and then applying and adapting the methods to two specific engineering applications: the reflector antenna problem and the moving mesh methods problem. For these types of nonlinear partial differential equations, many numerical studies have been done in recent years, the vast majority in subsets of Euclidean space. In this dissertation, the first major goal is to develop convergent schemes for the sphere. However, another goal of this dissertation is application-centered, that is evaluating whether the partial differential equation techniques using Optimal Transport are actually the best methods for solving such problems. The reflector antenna is an optics inverse problem where one finds the shape of a reflector surface in order to refocus light into a prescribed far-field output intensity. This problem can be solved using Optimal Transport. The moving mesh methods problem is an adaptive mesh technique where one redistributes the density of the vertices of a mesh without tangling the edges connecting the vertices. Both Optimal Transport and Optimal Information Transport approaches can be used in solving this problem. The Monge Problem of Optimal Transport is concerned with computing the “optimal” mapping between two probability distributions. This actually can define a Riemannian distance between probability measures in a probability space. An-other choice of Riemannian metric on this space, the infinite-dimensional Fisher-Rao metric, gives an “information geometric” structure to the space of probability measures. It turns out that a simple partial differential equation can be solved for a mapping that relates to the underlying information geometry given by the Fisher-Rao metric. Solving for such an “information geometric” mapping is known as Optimal Information Transport. In this dissertation, a convergence framework is first established for com-puting the solution to the partial differential equation formulation of Optimal Transport on the sphere. This convergence framework uses geodesic normal coor-dinates to perform computations in local tangent planes. The numerical scheme also has a control on the Lipschitz constant of the discrete solution, which allows a convergence theorem for consistent and monotone discretizations to be proved in the absence of a comparison principle for the partial differential equation. Then, a finite-difference scheme for the partial differential equation formulation of Opti-mal Transport on the sphere is constructed which satisfies the hypotheses of the convergence theorem. An explicit formula for the mixed Hessian term is derived for two different cost functions. In order to construct a monotone discretization, discrete Laplacian terms are carefully added into the scheme. Current work has established convergence rates for solutions of monotone discretizations of linear elliptic partial differential equations on compact 2D manifolds without boundary. The goal is to then generalize these linearized arguments for the Optimal Transport case on the sphere. Computations are performed for the reflector antenna problem. Other ad hoc schemes exist for computing the reflector antenna problem, but the proposed scheme is the most efficient provably convergent scheme. Further adaptations are made that allow for the scheme to deal with non-smooth cases more explicitly. For the moving mesh methods problem, a comparison of computations via Optimal Transport and Optimal Information Transport is performed for the sphere using provably convergent monotone schemes for both computations. These comparisons show the merits of using Optimal Information Transport for some challenging computations. Optimal Information Transport also seems like a natural generalization to other compact 2D surfaces beyond the sphere.