A Characteristic Mapping Method for Tracer Transport on the Sphere

Abstract

A semi-Lagrangian Characteristic Mapping method for the solution of the tracer transport equations on the sphere is presented. The method solves for the solution operator of the equations by approximating the inverse of the diffeomorphism generated by a given velocity field. The evolution of any tracer and mass density can then be computed via pullback with this map. We present a spatial discretization of the manifold-valued map using a projection-based approach with spherical spline interpolation. The numerical scheme yields $C^1$ continuity for the map and global second-order accuracy for the solution of the tracer transport equations. Error estimates are provided and supported by convergence tests involving solid body rotation, moving vortices, deformational, and compressible flows. Additionally, we illustrate some features of computing the solution operator using a numerical mixing test and the transport of a fractal set in a complex flow environment.

arXiv commentUpdated accepted manuscript

References in corpus (7)

• The "Yin-Yang Grid": An Overset Grid in Spherical Geometry
• A gradient-augmented level set method with an optimally local, coherent advection scheme
• Geometric integration on spheres and some interesting applications
• Mesh-free Semi-Lagrangian Methods for Transport on a Sphere Using Radial Basis Functions
• Jet schemes for advection problems
• The remapped particle-mesh advection scheme
• A Characteristic Mapping Method for the two-dimensional incompressible Euler equations

Cited by in corpus (2)

• A Characteristic Mapping Method for Vlasov-Poisson with Extreme Resolution Properties
• A Characteristic Mapping Method with Source Terms: Applications to Ideal Magnetohydrodynamics