An efficient RBF-FD method using polyharmonic splines alongside polynomials for the numerical solution of two-dimensional PDEs held on irregular domains and subject to Dirichlet and Robin boundary conditions

Document Type : Research Paper

Authors

1 Department of Applied Mathematics, Imam Khomeini International University, Qazvin, 34149-16818, Iran

2 Department of Applied Mathematics, Imam Khomeini International University, Qazvin, 34148-96818, Iran

Abstract

In the present paper, the relatively new method of Radial Basis Function-Generated Finite Difference (RBF-FD) is used to solve a class of Partial Differential Equations (PDEs) with Dirichlet and Robin boundary conditions. For this approximation, Polyharmonic Splines (PHS) are used alongside Polynomials. This combination has many benefits. On the other hand, Polyharmonic Splines have no shape parameter and therefore relieve us of the hassle of calculating the optimal shape parameter. As the first problem, a two-dimensional Poisson equation with the Dirichlet boundary condition is investigated in various domains. Then, an elliptic PDE with the Robin boundary condition is solved by the proposed method. The results of numerical studies indicate the excellent efficiency, accuracy and high speed of the method, while for these studies, very fluctuating and special test functions have been used.

Keywords

[1] S. Abbasbandy, E. Shivanian, K. Hammood AL-Jizani, and S.N. Atluri, Pseudospectral meshless radial point interpolation for generalized biharmonic equation subject to simply supported and clamped boundary conditions, Engin. Anal. Boundary Elements 125 (2021), 23–32.
[2] M. Abbaszadeh, M. Dehghan, Simulation flows with multiple phases and components via the radial basis functions finite difference (RBF-FD) procedure: Shan-Chen model, Engineering Analysis with Boundary Elements, 119 (2020) 151-161.
[3] P. Assari, H. Adibi, and M. Dehghan, A meshless discrete Galerkin (MDG) method for the numerical solution of integral equations with logarithmic kernels, J. Comput. Appl. Math. 267 (2014), 160–181.
[4] P. Assari and M. Dehghan, A meshless discrete collocation method for the numerical solution of singular logarithmic boundary integral equations utilizing radial basis functions, Appl. Math. Comput. 315 (2017), 424–444.
[5] A.G. Atta, W.M. Abd-Elhameed, G.M. Moatimid, and Y.H. Youssri, Modal shifted fifth-kind Chebyshev tau integral approach for solving heat conduction equation, Fractal Fractional 6 (2022), no. 11, 619.
[6] A.G. Atta and Y.H. Youssri, Advanced shifted first-kind Chebyshev collocation approach for solving the nonlinear time-fractional partial integro-differential equation with a weakly singular kernel, Comp. Appl. Math. 41 (2022), no. 8, 381.
[7] A.G. Atta, W.M. Abd-Elhameed, G.M. Moatimid, and Y.H.Youssri, A fast Galerkin approach for solving the fractional Rayleigh-Stokes problem via sixth-kind Chebyshev polynomials, Mathematics 10 (2022), no. 11, 1843.
[8] A.G. Atta, W.M. Abd-Elhameed, and Y.H. Youssri, Shifted fifth-kind Chebyshev polynomials Galerkin-based procedure for treating fractional diffusion-wave equation, Int. J. Modern Phys. 33 (2022), no. 8, 2250102.
[9] G.A. Barnett, A robust RBF-FD formulation based on polyharmonic splines and polynomials, Ph.D. thesis, University of Colorado, 2015.
[10] G. Barton, Elements of Green’s Functions and Propagation: Potentials, Diffusion, and Waves, Oxford University Press, New York, 1989.
[11] V. Bayona, M. Moscoso, M. Carretero, and M. Kindelan, RBF-FD formulas and convergence properties, J. Comput. Phys. 229 (2010), 8281–8295.
[12] V. Bayona, M. Moscoso, and M. Kindelan, Optimal constant shape parameter for multiquadric based RBF-FD method, J. Comput. Phys. 230 (2011), 7384–7399.
[13] V. Bayona, M. Moscoso, and M. Kindelan, Optimal variable shape parameter for multiquadric based RBF-FD method, J. Comput. Phys. 231 (2012), 2466–2481.
[14] V. Bayona, M. Moscoso, and M. Kindelan, Gaussian RBF-FD weights and its corresponding local truncation errors, Engin. Anal. Boundary Elements 36 (2012), 1361–1369.
[15] V. Bayona, N. Flyer, B. Fornberg, and G.A. Barnett, On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs, J. Comput. Phys. 332 (2017), 257–273.
[16] V. Bayona, An insight into RBF-FD approximations augmented with polynomials, Comput. Math. Appl. 77 (2019), 2337–2353.
[17] V. Bayona, N. Flyer, and B. Fornberg, On the role of polynomials in RBF-FD approximations: III. Behavior near domain boundaries, J. Comput. Phys. 380 (2019), 378–399.
[18] M. Dehghan and V. Mohammadi, A numerical scheme based on radial basis function finite difference (RBF-FD) technique for solving the high-dimensional nonlinear Schrodinger equations using an explicit time discretization: Runge–Kutta method, Compu. Phys. Commun. 217 (2017), 23–34.
[19] H. Ding, C. Shu, D. B. Tang, Error estimates of local multiquadric-based differential quadrature (LMQDQ) method through numerical experiments, Int. J. Numer. Methods Eng. 63 (2005), 1513–1529.
[20] W. Fang, Y. Wang, and Y. Xu, An implementation of fast wavelet Galerkin methods for integral equations of the second kind, J. Sci. Comput. 20 (2004), 277–302.
[21] G.E. Fasshauer, Meshfree Approximation Methods with MATLAB, Interdisciplinary Mathematical Sciences 6, World Scientific Publishers, Singapore, 2007.
[22] G.E. Fasshauer and J.G. Zhang, On choosing “optimal” shape parameters for RBF approximation, Numer. Alg. 45 (2007), 345–368.
[23] N. Flyer, B. Fornberg, V. Bayona, and G.A. Barnett, On the role of polynomials in RBF-FD approximations: I. Interpolation and accuracy, J. Comput. Phys. 321 (2016), 21–38.
[24] N. Flyer, G.A. Barnett, and L.J. Wicker, Enhancing finite differences with radial basis functions: Experiments on the Navier-Stokes equations, J. Comput. Phys. 316 (2016), 39–62.
[25] B. Fornberg and G. Wright, Stable computation of multiquadric interpolants for all values of the shape parameter, Comput. Math. Appl. 48 (2004), 853–867.
[26] B. Fornberg and C. Piret, On choosing a radial basis function and a shape parameter when solving a convective PDE on a sphere, J. Comput. Phys. 227 (2008), 2758–2780.
[27] B. Fornberg and N. Flyer, A Primer on Radial Basis Functions with Applications to the Geosciences, SIAM. Philadelphia, 2015.
[28] B. Fornberg and N. Flyer, Solving PDEs with radial basis functions, Acta Numerica 24 (2015), 215–258.
[29] D. Gunderman, N. Flyer, and B. Fornberg, Transport schemes in spherical geometries using spline-based RBF-FD with polynomials, J. Comput. Phys. 408 (2020), 109256.
[30] R.L. Hardy, Multiquadric equations of topography and other irregular surfaces, J. Geophy. Res. 76 (1971), 1905-1915.
[31] J. Helsing and A. Karlsson, An explicit kernel-split panel-based Nystrom scheme for integral equations on axially symmetric surfaces, J. Comput. Phys. 272 (2014), 686–703.
[32] E.J. Kansa, Multiquadrics—a scattered data approximation scheme with applications to computational fluid dynamics I: Surface approximations and partial derivative estimates, Comput. Math. Appl. 19 (1990), 127–145.
[33] E.J. Kansa, Multiquadrics-a scattered data approximation scheme with applications to computational fluid dynamics II: Solutions to parabolic, hyperbolic and elliptic partial differential equations, Comput. Math. Appl. 19 (1990), 147–161.
[34] H.R. Khodabandehlo, E. Shivanian, and S. Abbasbandy, Numerical solution of nonlinear delay differential equations of fractional variable-order using a novel shifted Jacobi operational matrix, Engin. Comput. 38 (2022), 2593–2607.
[35] G. Liu and Y. Gu, An Introduction to Meshfree Methods and Their Programming, Springer, 2005.
[36] R. Nittka, Elliptic and parabolic problems with Robin boundary conditions on Lipschitz domains, Ph.D. Thesis, Universitat Ulm, 2010.
[37] D.T. Oanh, O. Davydov, and H.X. Phu, Adaptive RBF-FD method for elliptic problems with point singularities in 2D, Appl. Math. Comput. 313 (2017), 474–497.
[38] Y. Qiao, S. Zhai, and X. Feng, RBF-FD method for the high dimensional time fractional convection-diffusion equation, Int. Commun. Heat Mass Transfer 89 (2017), 230–240.
[39] A. Rahimi, C.A.E. Shivanian, and S. Abbasbandy, Analysis of new RBF-FD weights, calculated based on inverse quadratic functions, J. Math. 2022 (2022), Article ID 3718132, 7 pages.
[40] J. Rashidinia and M.N. Rasoulizadeh, Numerical methods based on radial basis function-generated finite difference (RBF-FD) for solution of GKdVB equation, Wave Motion. 90 (2019), 152–167.
[41] D. Rostamy, M. Emamjome, and S. Abbasbandy, A meshless technique based on the pseudospectral radial basis functions method for solving the two-dimensional hyperbolic telegraph equation, Eur. Phys. J. Plus 132 (2017), 1–11.
[42] R. Schaback, Error estimates and condition numbers for radial basis function interpolants, Adv. Comput. Math. 3 (1995), 251–264.
[43] V. Shankar, G.B. Wright, R.M. Kirby, and A.L. Fogelson, A radial basis function (RBF)- finite difference (FD) method for diffusion and reaction-diffusion equations on surfaces, J. Sci. Comput. 63 (2015), 745–768.
[44] E. Shivanian, Pseudospectral meshless radial point hermit interpolation versus pseudospectral meshless radial point interpolation, Int. J. Comput. Meth. 17 (2020), no. 7, 1950023.
[45] E. Shivanian, A. Rahimi, and M. Hosseini, Meshless local radial point interpolation to three-dimensional wave equation with Neumann’s boundary conditions, Int. J. Comput. Math. 93 (2016), no. 12, 2124–2140.
[46] C. Shu, H. Ding, K.S. Yeo, Local radial basis function-based differential quadrature method and its application to solve two-dimensional incompressible Navier-Stokes equations, Comput. Meth. Appl. Mech. Eng. 192 (2003), 941–954.
[47] F. Soleymani, M. Barfeie, and F.K. Haghani, Inverse multi-quadric RBF for computing the weights of FD method: Application to American options, Commun. Nonlinear Sci. Numer. Simul. 64 (2018) 74-88.
[48] M. Tillenius, E. Larsson, E. Lehto, and N. Flyer, A scalable RBF-FD method for atmospheric flow, J. Comput. Phys. 298 (2015), 406–422.
[49] A.I. Tolstykh, On using RBF-based differencing formulas for unstructured and mixed structured-unstructured grid calculations, Proc. 16th IMACSWorld Cong. 228, 2000, pp. 4606–4624.
[50] A.I. Tolstykh and D. A. Shirobokov, On using radial basis functions in a ”finite difference mode” with applications to elasticity problems, Comput. Mech. 33 (2003), 68–79.
[51] I. Tominec, P.F. Villard, E. Larsson, V. Bayona, and N. Cacciani, An unfitted radial basis function generated finite difference method applied to thoracic diaphragm simulations, J. Comput. Phys. 469 (2022), 111496.
[52] G.B. Wright, Radial basis function interpolation: Numerical and analytical developments, Ph.D. thesis, University of Colorado, Boulder, CO, 2003.
[53] G.B. Wright and B. Fornberg, Scattered node compact finite difference-type formulas generated from radial basis functions, J. Comput. Phys. 212 (2006), 99–123.
[54] Y.H. Youssri, Two Fibonacci operational matrix pseudospectral schemes for nonlinear fractional Klein–Gordon equation, Int. J. Modern Phys. 33 (2022), no. 4, 2250049.
Volume 15, Issue 4
April 2024
Pages 173-184
  • Receive Date: 23 August 2022
  • Revise Date: 03 March 2023
  • Accept Date: 23 March 2023