Document Type


Publication Date



The current paper establishes the computational efficiency and accuracy of the RBF-FD method for large-scale geoscience modeling with comparisons to state-of-the-art methods as high-order discontinuous Galerkin and spherical harmonics, the latter using expansions with close to 300,000 bases. The test cases are demanding fluid flow problems on the sphere that exhibit numerical challenges, such as Gibbs phenomena, sharp gradients, and complex vortical dynamics with rapid energy transfer from large to small scales over short time periods. The computations were possible as well as very competitive due to the implementation of hyperviscosity on large RBF stencil sizes (corresponding roughly to 6th to 9th order methods) with up to O(105) nodes on the sphere. The RBF-FD method scaled as O(N) per time step, where N is the total number of nodes on the sphere. In Appendix A, guidelines are given on how to chose parameters when using RBF-FD to solve hyperbolic PDEs.

Copyright Statement

NOTICE: This is the author’s version of a work accepted for publication by Elsevier. Changes resulting from the publishing process, including peer review, editing, corrections, structural formatting and other quality control mechanisms, may not be reflected in this document. Changes may have been made to this work since it was submitted for publication. The definitive version has been published in Journal of Computational Physics, 2012. DOI: 10.1016/

Included in

Mathematics Commons