An improved radial basis-pseudospectral method with hybrid Gaussian-cubic kernels
For researchers using pseudospectral methods, this work offers a practical improvement to overcome conditioning issues, though it is an incremental modification of existing RBF techniques.
The paper proposes a hybrid radial basis function combining cubic splines and Gaussian kernels to improve conditioning and scalability of radial basis-pseudospectral methods. The approach reduces ill-conditioning and maintains stability and accuracy for small shape parameters, outperforming pure Gaussian and cubic kernels in solving Poisson, Helmholtz, and Transport equations.
While pseudospectral (PS) methods can feature very high accuracy, they tend to be severely limited in terms of geometric flexibility. Application of global radial basis functions overcomes this, however at the expense of problematic conditioning (1) in their most accurate flat basis function regime, and (2) when problem sizes are scaled up to become of practical interest. The present study considers a strategy to improve on these two issues by means of using hybrid radial basis functions that combine cubic splines with Gaussian kernels. The parameters, controlling Gaussian and cubic kernels in the hybrid RBF, are selected using global particle swarm optimization. The proposed approach has been tested with radial basis-pseudospectral method for numerical approximation of Poisson, Helmholtz, and Transport equation. It was observed that the proposed approach significantly reduces the ill-conditioning problem in the RBF-PS method, at the same time, it preserves the stability and accuracy for very small shape parameters. The eigenvalue spectra of the coefficient matrices in the improved algorithm were found to be stable even at large degrees of freedom, which mimic those obtained in pseudospectral approach. Also, numerical experiments suggest that the hybrid kernel performs significantly better than both pure Gaussian and pure cubic kernels.