Skip to main content
SpringerLink
Log in

Adaptive residual refinement in an RBF finite difference scheme for 2D time-dependent problems

  • Published:
Computational and Applied Mathematics Aims and scope Submit manuscript

Abstract

The current study introduces a fast and accurate method based on a novel adaptive meshfree technique to solve time-dependent partial differential equations with solutions representing high gradients or quick changes in several local areas of the domain. Utilizing uniform grids for these problems is prohibitive computationally, since the solution reaches singularity. This study aims to suggest an adaptive strategy to produce a suitable and cost-effective irregular node refinement. For this purpose, a dynamic algorithm is proposed that finds areas with quick changes and applies a local node adaptive approach merely in those nearly singular areas. Additionally, within this algorithm, unlike the Kansa technique, the radial basis function collocation technique was mixed with a finite difference scheme. According to this approach, in place of using the adaptive algorithm on the complete domain of the problem, it can be used only on time steps. Therefore, we need only to solve small systems of linear equations on each time step instead of large systems on the entire domain. Besides performing a stability analysis of the numerical scheme, the new algorithm is tested on parabolic (heat equation) and hyperbolic (wave equation) PDEs over regular and irregular two-dimensional domains. The attained results prove the accuracy and effectiveness of the proposed technique. Especially, our computational method is able to reduce the nodes in the domain with no impairment in terms of accuracy, thus turning out to be effective in the localization of oscillations owing to sharp gradients in the solution.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Algorithm 1
Algorithm 2
Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12
Fig. 13
Fig. 14
Fig. 15
Fig. 16
Fig. 17
Fig. 18
Fig. 19

Similar content being viewed by others

Data Availability

Data are available upon reasonable request.

References

  • Ahmad M, Siraj-ul-Islam (2018) Meshless analysis of parabolic interface problems. Eng Anal Bound Elem 94:134–152

    MathSciNet  Google Scholar 

  • Ahmad M, Siraj-ul-Islam, Ullah B (2020) Local radial basis function collocation method for Stokes equations with interface conditions. Eng Anal Bound Elem 119:246–256

    MathSciNet  Google Scholar 

  • Ames WF (2014) Numerical methods for partial differential equations. Academic Press, London

    Google Scholar 

  • Arney DC, Flaherty JE (1990) An adaptive mesh-moving and local refinement method for time-dependent partial differential equations. ACM Trans Math Softw 16:48–71

    MathSciNet  Google Scholar 

  • Bell J, Berger M, Saltzman J, Welcome M (1994) Three-dimensional adaptive mesh refinement for hyperbolic conservation laws. SIAM J Sci Comput 15(1):127–138

    MathSciNet  Google Scholar 

  • Boyd JP (2005) Trouble with Gegenbauer reconstruction for defeating Gibbs phenomenon: runge phenomenon in the diagonal limit of Gegenbauer polynomial approximations. J Comput Phys 204:253–264

    MathSciNet  ADS  Google Scholar 

  • Boztosun I, Charafi A (2002) An analysis of the linear advection-diffusion equation using mesh-free and mesh-dependent methods. Eng Anal Bound Elem 26:889–895

    Google Scholar 

  • Buhmann MD (2003) Radial basis functions: theory and implementations, vol 12. Cambridge monographs on applied and computational mathematics. Cambridge University Press, Cambridge

    Google Scholar 

  • Buhmann M, Dyn N (1993) Spectral convergence of multiquadric interpolation. Proc Edinb Math Soc 36(2):319–333

    MathSciNet  Google Scholar 

  • Burgarelli D, Kischinhevsky M, Biezuner RJ (2006) A new adaptive mesh refinement strategy for numerically solving evolutionary PDE’s. J Comput Appl Math 196(1):115–131

    MathSciNet  ADS  Google Scholar 

  • Carlson RE, Foley T (1991) The parameter \(R^2\) in multiquatric interpolation. Comput Math Appl 21:29–42

    MathSciNet  Google Scholar 

  • Cavoretto R (2015) A numerical algorithm for multidimensional modeling of scattered data points. Comput Appl Math 34:65–80

    MathSciNet  Google Scholar 

  • Cavoretto R (2021) Adaptive radial basis function partition of unity interpolation: a bivariate algorithm for unstructured data. J Sci Comput 87:41

    MathSciNet  Google Scholar 

  • Cavoretto R (2022) Adaptive LOOCV-based kernel methods for solving time-dependent BVPs. Appl Math Comput 429:127228

    MathSciNet  Google Scholar 

  • Cavoretto R, De Rossi A (2019) Adaptive meshless refinement schemes for RBF-PUM collocation. Appl Math Lett 90:131–138

    MathSciNet  Google Scholar 

  • Cavoretto R, De Rossi A (2020a) Error indicators and refinement strategies for solving Poisson problems through a RBF partition of unity collocation scheme. Appl Math Comput 369:124824

    MathSciNet  Google Scholar 

  • Cavoretto R, De Rossi A (2020b) A two-stage adaptive scheme based on RBF collocation for solving elliptic PDEs. Comput Math Appl 79:3206–3222

    MathSciNet  Google Scholar 

  • Cavoretto R, De Rossi A (2023) An adaptive residual sub-sampling algorithm for kernel interpolation based on maximum likelihood estimations. J Comput Appl Math 418:114658

    MathSciNet  Google Scholar 

  • Cavoretto R, De Rossi A, Erb W (2021) Partition of unity methods for signal processing on graphs. J Fourier Anal Appl 27:66

    MathSciNet  Google Scholar 

  • Çayan S, Özhan BB, Sezer M (2022) An adaptive approach for solving fourth-order partial differential equations: algorithm and applications to engineering models. Comput Appl Math 41(8):1–17

    MathSciNet  Google Scholar 

  • Chengi AH-D, Golberg MA, Kansa EJ, Zammito G (2003) Exponential convergence and H-c multiquadric collocation method for partial differential equations. Numer Methods Partial Differ Equ 19(5):571–594

    MathSciNet  Google Scholar 

  • Chinchapatnam PP, Djidjeli K, Nair PB (2006) Unsymmetric and symmetric meshless schemes for the unsteady convection-diffusion equation. Comput Methods Appl Mech Eng 195:2432–2453

    MathSciNet  ADS  Google Scholar 

  • Davis SF, Flaherty JE (1982) An adaptive finite element method for initial-boundary value problems for partial differential equations. SIAM J Sci Stat Comput 3(1):6–22

    MathSciNet  Google Scholar 

  • Dehghan M, Shokri A (2007) A numerical method for two-dimensional Schrödinger equation using collocation radial basis functions. Comput Math Appl 54:136–146

    MathSciNet  Google Scholar 

  • Dolejsi V (2015) Anisotropic hp-adaptive discontinuous Galerkin method for the numerical solution of time dependent PDEs. Appl Math Comput 267:682–697

    MathSciNet  Google Scholar 

  • Driscoll TA, Heryudono ARH (2007) Adaptive residual subsampling for radial basis function interpolation collocations problems. Comput Math Appl 53:927–939

    MathSciNet  Google Scholar 

  • Esmaeilbeigi M, Chatrabgoun O (2019) An efficient method based on RBFs for multilayer data interpolation with application in air pollution data analysis. Comput Appl Math 38(4):1–20

    MathSciNet  Google Scholar 

  • Esmaeilbeigi M, Garmanjani G (2016) A shift-adaptive meshfree method for solving a class of initial-boundary value problems with moving boundaries in one-dimensional domain. Numer Methods Partial Differ Equ 32(6):1622–1646

    MathSciNet  Google Scholar 

  • Esmaeilbeigi M, Garmanjani G (2017) Gaussian radial basis function interpolant for the different data sites and basis centers. Calcolo 54(1):155–166

    MathSciNet  Google Scholar 

  • Esmaeilbeigi M, Hosseini MM (2012) Dynamic node adaptive strategy for nearly singular problems on large domains. Eng Anal Bound Elem 36:1311–1321

    MathSciNet  Google Scholar 

  • Esmaeilbeigi M, Hosseini MM, Mohyud-Din ST (2011) A new approach of the radial basis functions method for telegraph equations. Int J Phys Sci 6(6):1517–1527

    Google Scholar 

  • Esmaeilbeigi M, Paripour M, Garmanjani G (2018) Approximate solution of the fuzzy fractional Bagley–Torvik equation by the RBF collocation method. Comput Methods Differ Equ 6(2):186–214

    MathSciNet  Google Scholar 

  • Fasshauer GE (2007) Meshfree approximation methods with Matlab. World Scientific, Singapore

    Google Scholar 

  • Fasshauer GE (2011) Positive definite kernels: past, present and future. Dolomit Res Notes Approx 4:21–63

    Google Scholar 

  • Fornberg B, Zuev J (2007) The Runge phenomenon and spatially variable shape parameters in RBF interpolation. Comput Math Appl 54:379–398

    MathSciNet  Google Scholar 

  • Franke R (1982) Scattered data interpolation: test of some methods. Math Comput 38:181–200

    MathSciNet  Google Scholar 

  • Garmanjani G, Cavoretto R, Esmaeilbeigi M (2018) A RBF partition of unity collocation method based on finite difference for initial-boundary value problems. Comput Math Appl 75(11):4066–4090

    MathSciNet  Google Scholar 

  • Ge L, Sun T (2022) An adaptive hp-version stochastic Galerkin method for constrained optimal control problem governed by random reaction diffusion equations. Comput Appl Math 41(3):1–30

    MathSciNet  Google Scholar 

  • Golbabai A, Mohebianfar E, Rabiei H (2015) On the new variable shape parameter strategies for radial basis functions. Comput Appl Math 34:691–704

    MathSciNet  Google Scholar 

  • Hon YC, Schaback R, Zhou X (2003) An adaptive greedy algorithm for solving large RBF collocation problems. Numer Algorithms 32(1):13–25

    MathSciNet  ADS  Google Scholar 

  • Hussain M, Haq S (2020) A hybrid radial basis functions collocation technique to numerically solve fractional advection-diffusion models. Numer Methods Partial Differ Equ. https://doi.org/10.1002/num.22472

    Article  MathSciNet  Google Scholar 

  • Jannesari Z, Tatari M (2020) An adaptive strategy for solving convection dominated diffusion equation. Comput Appl Math 39(2):1–15

    MathSciNet  Google Scholar 

  • Lefloch PG (2000) Hyperbolic systems of conservation laws. Oxford University Press, Oxford

    Google Scholar 

  • Lu J, Nie Y (2021) A collocation method based on localized radial basis functions with reproducibility for nonlocal diffusion models. Comput Appl Math 40(8):1–23

    MathSciNet  Google Scholar 

  • Martinez FMB (2008) Meshless methods for elliptic and free-boundary problems. PhD thesis

  • Sarra SA (2005) Adaptive radial basis function methods for time dependent partial differential equations. Appl Numer Math 54(1):79–94

    MathSciNet  Google Scholar 

  • Singh V, Siraj-ul-Islam, Mohanty RK (2019) Local meshless method for convection dominated steady and unsteady partial differential equations. Eng Comput 35:803–812

    Google Scholar 

  • Siraj-ul-Islam, Khan MN, Hussain I, Ahmad I, Ahmad H (2021) A local meshless method for the numerical solution of space-dependent inverse heat problems. Math Methods Appl Sci 44(4):3066–3079

    MathSciNet  Google Scholar 

  • Ullah Z, Ullah B, Khan W, Siraj-ul-Islam (2022) Proportional topology optimisation with maximum entropy-based meshless method for minimum compliance and stress constrained problems. Eng Comput 38(6):5541–5561

    Google Scholar 

  • Wendland H (1995) Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. Adv Comput Math 4:389–396

    MathSciNet  Google Scholar 

  • Wendland H (2005) Scattered data approximation, vol 17. Cambridge monographs on applied and computational mathematics. Cambridge University Press, Cambridge

    Google Scholar 

  • Wu Z (2004) Dynamically knots setting in meshless method for solving time dependent propagation equation. Comput Methods Appl Mech Eng 193(12–14):1221–1229

    MathSciNet  ADS  Google Scholar 

Download references

Acknowledgements

The authors sincerely thank the anonymous referee(s) for the valuable comments and suggestions, which enabled to significantly improve the quality of this paper. The third author acknowledges financial support from the INdAM–GNCS Project “Computational methods for kernel-based approximation and its applications”, code CUP\(\_\)E55F22000270001. Moreover, this work was partially supported by the Department of Mathematics “G. Peano” of the University of Torino through 2020 project “Mathematical methods in computational sciences”. This research has been accomplished within the RITA (Research ITalian network on Approximation) and the UMI Group TAA (Approximation Theory and Applications).

Author information

Authors and Affiliations

  1. Department of Mathematics, Malayer University, Malayer, 65719-95863, Iran

    G. Garmanjani & M. Esmaeilbeigi

  2. Department of Mathematics “G. Peano”, University of Torino, via Carlo Alberto 10, 10123, Torino, Italy

    R. Cavoretto

  3. Member of the INdAM Research group GNCS, Rome, Italy

    R. Cavoretto

Authors
  1. G. Garmanjani
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. M. Esmaeilbeigi
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. R. Cavoretto
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to M. Esmaeilbeigi.

Ethics declarations

Conflict of interest

The authors declare that they have no any competing interests.

Additional information

Communicated by Fabricio Simeoni de Sousa.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Garmanjani, G., Esmaeilbeigi, M. & Cavoretto, R. Adaptive residual refinement in an RBF finite difference scheme for 2D time-dependent problems. Comp. Appl. Math. 43, 39 (2024). https://doi.org/10.1007/s40314-023-02541-1

Download citation

  • Received:

  • Revised:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s40314-023-02541-1

Keywords

Mathematics Subject Classification

Search

Navigation