Premium
Semi‐Lagrangian advection on a spherical geodesic grid
Author(s) -
Carfora Maria Francesca
Publication year - 2007
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.1445
Subject(s) - discretization , geodesic , mathematics , grid , advection , interpolation (computer graphics) , triangulation , geometry , rotation (mathematics) , mathematical analysis , solving the geodesic equations , flow (mathematics) , grid method multiplication , surface (topology) , classical mechanics , physics , motion (physics) , thermodynamics
A simple and efficient numerical method for solving the advection equation on the spherical surface is presented. To overcome the well‐known ‘pole problem’ related to the polar singularity of spherical coordinates, the space discretization is performed on a geodesic grid derived by a uniform triangulation of the sphere; the time discretization uses a semi‐Lagrangian approach. These two choices, efficiently combined in a substepping procedure, allow us to easily determine the departure points of the characteristic lines, avoiding any computationally expensive tree‐search. Moreover, suitable interpolation procedures on such geodesic grid are presented and compared. The performance of the method in terms of accuracy and efficiency is assessed on two standard test cases: solid‐body rotation and a deformation flow. Copyright © 2007 John Wiley & Sons, Ltd.