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

DOI10.1016/J.JCP.2015.12.018zbMATH Open1351.86008arXiv1512.02935OpenAlexW2189647597MaRDI QIDQ729449FDOQ729449

Authors: Hilary Weller, Mike Cullen, P. A. Browne, C. J. Budd 

Publication date: 20 December 2016

Published in: Journal of Computational Physics (Search for Journal in Brave)

Abstract: An equation of Monge-Amp`ere 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`ere 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.

Full work available at URL: https://arxiv.org/abs/1512.02935

Recommendations

• Optimal-transport -- based mesh adaptivity on the plane and sphere using finite elements
• Moving Mesh Generation Using the Parabolic Monge–Ampère Equation
• The scaling and skewness of optimally transported meshes on the sphere
• The geometry of r-adaptive meshes generated using optimal transport methods
• An optimal robust equidistribution method for two-dimensional grid adaptation based on Monge-Kantorovich optimization

zbMATH Keywords

modellingoptimal transportrefinementadaptivemesh generationmeshatmosphereMonge-Ampére

Mathematics Subject Classification ID

Mesh generation, refinement, and adaptive methods for boundary value problems involving PDEs (65N50) Computational methods for problems pertaining to geophysics (86-08) Meteorology and atmospheric physics (86A10)

Cites Work

• Polar factorization and monotone rearrangement of vector‐valued functions
• Title not available (Why is that?)
• Polar factorization of maps on Riemannian manifolds
• Variational mesh adaptation: Isotropy and equidistribution
• Mixed finite elements for numerical weather prediction
• Adaptive mesh refinement for hyperbolic partial differential equations
• Numerical representation of geostrophic modes on arbitrarily structured C-grids
• Centroidal Voronoi Tessellations: Applications and Algorithms
• Convergent finite difference solvers for viscosity solutions of the elliptic Monge-Ampère equation in dimensions two and higher
• A numerical method for the elliptic Monge-Ampère equation with transport boundary conditions
• Fast three dimensional r-adaptive mesh redistribution
• Adaptivity with moving grids
• Modelling atmospheric flows with adaptive moving meshes
• An efficient numerical algorithm for the \(L^{2}\) optimal transport problem with periodic densities
• Moving mesh methods in multiple dimensions based on harmonic maps
• The geometry of r-adaptive meshes generated using optimal transport methods
• An augmented Lagrangian approach to the numerical solution of the Dirichlet problem for the elliptic Monge-Ampère equation in two dimensions
• Parabolic Monge–Ampère methods for blow-up problems in several spatial dimensions
• Two numerical methods for the elliptic Monge-Ampère equation
• Wide stencil finite difference schemes for the elliptic Monge-Ampère equation and functions of the eigenvalues of the Hessian
• Numerical methods for fully nonlinear elliptic equations of the Monge-Ampère type
• Mixed finite element methods for the fully nonlinear Monge-Ampère equation based on the vanishing moment method
• The Monge-Ampère trajectory correction for semi-Lagrangian schemes

Cited In (22)

• A continuation multiple shooting method for Wasserstein geodesic equation
• A surface moving mesh method based on equidistribution and alignment
• Mesh optimization for singular axisymmetric harmonic maps from the disc into the sphere
• The scaling and skewness of optimally transported meshes on the sphere
• Hybridizable discontinuous Galerkin methods for the two-dimensional Monge-Ampère equation
• An adaptive moving mesh method for forced curve shortening flow
• A convergent finite difference method for optimal transport on the sphere
• Conservation with moving meshes over orography
• Multi-scale hydro-morphodynamic modelling using mesh movement methods
• A convergence framework for optimal transport on the sphere
• Numerical homogenization by an adaptive Fourier spectral method on non-uniform grids using optimal transport
• A Volumetric approach to Monge's optimal transport on surfaces
• Bandwidth-based mesh adaptation in multiple dimensions
• Sliced optimal transport on the sphere
• Optimal-transport -- based mesh adaptivity on the plane and sphere using finite elements
• Well balanced residual distribution for the ALE spherical shallow water equations on moving adaptive meshes
• Moving mesh method for simulating high-dimensional time dependent PDEs with fast propagating shock waves
• The Sinkhorn algorithm, parabolic optimal transport and geometric Monge-Ampère equations
• The geometry of r-adaptive meshes generated using optimal transport methods
• A fast dynamic smooth adaptive meshing scheme with applications to compressible flow
• Adaptive mesh methods on compact manifolds via optimal transport and optimal information transport
• An adaptive conservative moving mesh method

This page was built for publication: Mesh adaptation on the sphere using optimal transport and the numerical solution of a Monge-Ampère type equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q729449)