RBF-FD Formulas and Convergence Properties

The local RBF is becoming increasingly popular as an alternative to the global version that suffers from ill-conditioning. In this paper, we study analytically the convergence behavior of the local RBF method as a function of the number of nodes employed in the scheme, the nodal distance, and the shape parameter. We derive exact formulas for the first and second derivatives in one dimension, and for the Laplacian in two dimensions. Using these formulas we compute Taylor expansions for the error. From this analysis, we find that there is an optimal value of the shape parameter for which the error is minimum. This optimal parameter is independent of the nodal distance. Our theoretical results are corroborated by numerical experiments.This work has been supported by Spanish MECD Grants FIS2007-62673, FIS2008-04921 and by Madrid Autonomous Region Grant S2009-1597

A least-squares implicit RBF-FD closest point method and applications to PDEs on moving surfaces

The closest point method (Ruuth and Merriman, J. Comput. Phys. 227(3):1943-1961, [2008]) is an embedding method developed to solve a variety of partial differential equations (PDEs) on smooth surfaces, using a closest point representation of the surface and standard Cartesian grid methods in the embedding space. Recently, a closest point method with explicit time-stepping was proposed that uses finite differences derived from radial basis functions (RBF-FD). Here, we propose a least-squares implicit formulation of the closest point method to impose the constant-along-normal extension of the solution on the surface into the embedding space. Our proposed method is particularly flexible with respect to the choice of the computational grid in the embedding space. In particular, we may compute over a computational tube that contains problematic nodes. This fact enables us to combine the proposed method with the grid based particle method (Leung and Zhao, J. Comput. Phys. 228(8):2993-3024, [2009]) to obtain a numerical method for approximating PDEs on moving surfaces. We present a number of examples to illustrate the numerical convergence properties of our proposed method. Experiments for advection-diffusion equations and Cahn-Hilliard equations that are strongly coupled to the velocity of the surface are also presented

Pricing Financial Derivatives using Radial Basis Function generated Finite Differences with Polyharmonic Splines on Smoothly Varying Node Layouts

In this paper, we study the benefits of using polyharmonic splines and node layouts with smoothly varying density for developing robust and efficient radial basis function generated finite difference (RBF-FD) methods for pricing of financial derivatives. We present a significantly improved RBF-FD scheme and successfully apply it to two types of multidimensional partial differential equations in finance: a two-asset European call basket option under the Black--Scholes--Merton model, and a European call option under the Heston model. We also show that the performance of the improved method is equally high when it comes to pricing American options. By studying convergence, computational performance, and conditioning of the discrete systems, we show the superiority of the introduced approaches over previously used versions of the RBF-FD method in financial applications

A well-balanced meshless tsunami propagation and inundation model

We present a novel meshless tsunami propagation and inundation model. We discretize the nonlinear shallow-water equations using a well-balanced scheme relying on radial basis function based finite differences. The inundation model relies on radial basis function generated extrapolation from the wet points closest to the wet-dry interface into the dry region. Numerical results against standard one- and two-dimensional benchmarks are presented.Comment: 20 pages, 13 figure