On choosing a radial basis function and a shape parameter when solving a convective PDE on a sphere

详细信息    查看全文

• 作者：Bengt Fornberg ; Cé ; cile Piret
• 关键词：Radial basis functions ; RBF ; Shape parameter ; Wendland functions ; Sphere ; Spherical harmonics ; RBF-QR method
• 刊名：Journal of Computational Physics
• 出版年：2008
• 出版时间：20 February 2008
• 年：2008
• 卷：227
• 期：5
• 页码：2758-2780
• 全文大小：1399 K

文摘

Radial basis function (RBF) approximations have been used for some time to interpolate data on a sphere (as well as on many other types of domains). Their ability to solve, to spectral accuracy, convection-type PDEs over a sphere has been demonstrated only very recently. In such applications, there are two main choices that have to be made: (i) which type of radial function to use, and (ii) what value to choose for their shape parameter (denoted by method=retrieve&_udi=B6WHY-4R68NMY-2&_mathId=mml62&_user=10&_cdi=6863&_rdoc=3&_acct=C000050221&_version=1&_userid=10&md5=dc1790d1b4a50e4df3935d71ccec79d5"" title=""Click to view the MathML source"">ε, and with flat basis functions – stretched out in the radial direction – corresponding to method=retrieve&_udi=B6WHY-4R68NMY-2&_mathId=mml63&_user=10&_cdi=6863&_rdoc=3&_acct=C000050221&_version=1&_userid=10&md5=8b8a6680f76ffcabbd3bb52c5db94238"" title=""Click to view the MathML source"">ε=0). The recent RBF-QR algorithm has made it practical to compute stably also for small values of method=retrieve&_udi=B6WHY-4R68NMY-2&_mathId=mml64&_user=10&_cdi=6863&_rdoc=3&_acct=C000050221&_version=1&_userid=10&md5=193264dc12225d598326792620cc7e29"" title=""Click to view the MathML source"">ε. Results from solving a convective-type PDE on a sphere are compared here for many choices of radial functions over the complete range of method=retrieve&_udi=B6WHY-4R68NMY-2&_mathId=mml65&_user=10&_cdi=6863&_rdoc=3&_acct=C000050221&_version=1&_userid=10&md5=533e2385de600b85b24aef915742e056"" title=""Click to view the MathML source"">ε-values (from very large down to the limit of method=retrieve&_udi=B6WHY-4R68NMY-2&_mathId=mml66&_user=10&_cdi=6863&_rdoc=3&_acct=C000050221&_version=1&_userid=10&md5=3131eca712d9d91db916cece51eb9409"" title=""Click to view the MathML source"">ε→0). The results are analyzed with a methodology that has similarities to the customary Fourier analysis in equispaced 1-D periodic settings. In particular, we find that high accuracy can be maintained also over very long time integrations. We furthermore gain insights into why RBFs sometimes offer higher accuracy than spherical harmonics (since the latter arise as an often non-optimal special case of the former). Anticipated future application areas for RBF-based methods in spherical geometries include weather and climate modeling.