An official website of the United States government
Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher.
Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?
Some links on this page may take you to non-federal websites. Their policies may differ from this site.
-
We consider a partial differential equation (PDE) approach to numerically solve the reflector antenna problem by solving an optimal transport problem on the unit sphere with cost function . At each point on the sphere, we replace the surface PDE with a generalized Monge–Ampère type equation posed on the local tangent plane. We then use a provably convergent finite difference scheme to approximate the solution and construct the reflector. The method is easily adapted to take into account highly nonsmooth data and solutions, which makes it particularly well adapted to real-world optics problems. Computational examples demonstrate the success of this method in computing reflectors for a range of challenging problems including discontinuous intensities and intensities supported on complicated geometries.more » « less
-
We consider the numerical construction of minimal Lagrangian graphs, which is related to recent applications in materials science, molecular engineering, and theoretical physics. It is known that this problem can be formulated as an additive eigenvalue problem for a fully nonlinear elliptic partial differential equation. We introduce and implement a two-step generalized finite difference method, which we prove converges to the solution of the eigenvalue problem. Numerical experiments validate this approach in a range of challenging settings. We further discuss the generalization of this new framework to Monge-Ampère type equations arising in optimal transport. This approach holds great promise for applications where the data does not naturally satisfy the mass balance condition, and for the design of numerical methods with improved stability properties.more » « less
