Optimal-transport-based mesh adaptivity on the plane and sphere using finite elements

In moving mesh methods, the underlying mesh is dynamically adapted without changing the connectivity of the mesh. We specifically consider the generation of meshes which are adapted to a scalar monitor function through equidistribution. Together with an optimal transport condition, this leads to a Monge–Ampere equation ` for a scalar mesh potential. We adapt an existing finite element scheme for the standard Monge–Ampere ` equation to this mesh generation problem; this is a mixed finite element scheme, in which an extra discrete variable is introduced to represent the Hessian matrix of second derivatives. The problem we consider has additional nonlinearities over the basic Monge–Ampere equation due to the implicit dependence of the monitor func- ` tion on the resulting mesh. We also derive an equivalent Monge–Ampere-like equa- ` tion for generating meshes on the sphere. The finite element scheme is extended to the sphere, and we provide numerical examples. All numerical experiments are performed using the open-source finite element framework Firedrake