A kernel-free boundary integral method for reaction-diffusion equations

i

This article has been updated. A Correction to this article was published on .

i

The Original Article was published on

i

The Original Article was published on

i

This article has been updated

i

This article was retracted on

i

This article has been updated

Abstract

This paper was based on a kernel-free boundary integral (KFBI) method for solving the reaction-diffusion equation. The KFBI method serves as a general elliptic solvers for boundary value problems in an irregular problem domain. Unlike traditional boundary integral methods, the KFBI method avoids complicated direct integral calculations. Instead, a Cartesian grid-based five-point compact difference scheme was used to discretize the equivalent simple interface problem, whose solution is the integral involved in the corresponding boundary integral equations (BIEs). The resulting linear system was treated with a fast Fourier transform (FFT)-based elliptic solver, and the BIEs were iteratively solved by the generalized minimal residual (GMRES) method. The first step in solving the reaction-diffusion equation was to discretize the time variable with a two-stage second-order semi-implicit Runge-Kutta (SIRK) method, which transforms the problem into a spatial modified Helmholtz equation in each time step and can be solved by the KFBI method later. The proposed algorithm had second-order accuracy in both time and space even for small diffusion problems, and the computational work was roughly proportional to the number of grid nodes in the Cartesian grid due to the fast elliptic solver used. Numerical results verified the stability, efficiency, and accuracy of the method.

Keywords

kernel-free boundary integral method / Cartesian grid-based method / reaction-diffusion equation / two-step semi-implicit Runge-Kutta method

CHANGE HISTORY

A Correction to this paper has been published:
https://doi.org

The Retraction Note was updated to reflect reception of statements of disagreement from all authors.

This article has been retracted.Please see the Retraction Notice for more detail:
https://doi.org

References

[1] R. Bastiaansen, A. Doelman, The dynamics of disappearing pulses in a singularly perturbed reaction-diffusion system with parameters that vary in time and space, Physica D, 388 (2019), 45-72. https://doi.org/10.1016/j.physd.2018.09.003

[2] M. Chirilus-Bruckner, A. Doelman, P. van Heijster, J. D. M. Rademacher, Butterfly catastrophe for fronts in a three-component reaction-diffusion system, J. Nonlinear Sci., 25 (2015), 87-129. https://doi.org/10.1007/s00332-014-9222-9

[3] D. Gomez, J. Wei, Z. Yang, Multi-spike solutions to the one-dimensional subcritical fractional Schnakenberg system, Physica D, 448 (2023), 133720. https://doi.org/10.1016/j.physd.2023.133720

[4] O. Jaibi, A. Doelman, M. Chirilus-Bruckner, E. Meron, The existence of localized vegetation patterns in a systematically reduced model for dryland vegetation, Physica D, 412 (2020), 132637. https://doi.org/10.1016/j.physd.2020.132637

[5] A. Iuorio, F. Veerman, The influence of autotoxicity on the dynamics of vegetation spots, Physica D, 427 (2021), 133015. https://doi.org/10.1016/j.physd.2021.133015

[6] P. Carter, A. Doelman, K. Lilly, E. Obermayer, S. Rao, Criteria for the (in)stability of planar interfaces in singularly perturbed 2-component reaction-diffusion equations, Physica D, 444 (2023), 133596. https://doi.org/10.1016/j.physd.2022.133596

[7] L. K. Bieniasz, M. Vynnycky, S. McKee, Integral equation-based simulation of transient experiments for an EC2 mechanism: Beyond the steady state simplification, Electrochim. Acta, 428 (2022), 140896. https://doi.org/10.1016/j.electacta.2022.140896

[8] V. Lucarini, L. Serdukova, G. Margazoglou, Lévy noise versus Gaussian-noise-induced transitions in the Ghil-Sellers energy balance model, Nonlinear Processes Geophys., 29 (2022), 183-205. https://doi.org/10.5194/npg-29-183-2022

[9] A. Gupta, A. Kaushik, M. Sharma, A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction-diffusion problems with Robin-boundary conditions, Numer. Methods Partial Differ. Equations, 39 (2023), 1220-1250. https://doi.org/10.1002/num.22931

[10] J. Zhang, X. Liu, Convergence and supercloseness in a balanced norm of finite element methods on Bakhvalov-type meshes for reaction-diffusion problems, J. Sci. Comput., 88 (2021), 27. https://doi.org/10.1007/s10915-021-01542-8

[11] A. Kaushik, V. Kumar, M. Sharma, N. Sharma, A modified graded mesh and higher order finite element method for singularly perturbed reaction-diffusion problems, Math. Comput. Simul., 185 (2021), 486-496. https://doi.org/10.1016/j.matcom.2021.01.006

[12] S. I. Ei, M. Kuwamura, Y. Morita, A variational approach to singular perturbation problems in reaction-diffusion systems, Physica D, 207 (2005), 171-219. https://doi.org/10.1016/j.physd.2005.05.020

[13] K. Aarthika, R. Shiromani, V. Shanthi, A higher-order finite difference method for two-dimensional singularly perturbed reaction-diffusion with source-term-discontinuous problem, Comput. Math. Appl., 118 (2022), 56-73. https://doi.org/10.1016/j.camwa.2022.04.016

[14] X. Cai, D. L. Cai, R. Q. Wu, K. H. Xie, High accuracy non-equidistant method for singular perturbation reaction-diffusion problem, Appl. Math. Mech., 30 (2009), 175-182. https://doi.org/10.1007/s10483-009-0205-8

[15] S. C. S. Rao, A. K. Chaturvedi, Analysis and implementation of a computational technique for a coupled system of two singularly perturbed parabolic semilinear reaction-diffusion equations having discontinuous source terms, Commun. Nonlinear Sci. Numer. Simul., 108 (2022), 106232. https://doi.org/10.1016/j.cnsns.2021.106232

[16] X. Meng, M. Stynes, Balanced-norm and energy-norm error analyses for a backward Euler/FEM solving a singularly perturbed parabolic reaction-diffusion problem, J. Sci. Comput., 92 (2022), 67. https://doi.org/10.1007/s10915-022-01931-7

[17] M. Faustmann, J. M. Melenk, Robust exponential convergence of hp-FEM in balanced norms for singularly perturbed reaction-diffusion problems: Corner domains, Comput. Math. Appl., 74 (2017), 1576-1589. https://doi.org/10.1016/j.camwa.2017.03.015

[18] L. P. Franca, A. L. Madureira, F. Valentin, Towards multiscale functions: Enriching finite element spaces with local but not bubble-like functions, Comput. Methods Appl. Mech. Eng., 194 (2005), 3006-3021. https://doi.org/10.1016/j.cma.2004.07.029

[19] A. Secer, I. Onder, M. Ozisik, Sinc-Galerkin method for solving system of singular perturbed reaction-diffusion problems, Sigma J. Eng. Nat. Sci., 39 (2021), 203-212. https://doi.org/10.14744/sigma.2021.00010

[20] R. Lin, A discontinuous Galerkin least-squares finite element method for solving coupled singularly perturbed reaction-diffusion equations, J. Comput. Appl. Math., 307 (2016), 134-142. https://doi.org/10.1016/j.cam.2016.02.052

[21] X. Ma, J. Zhang, Supercloseness in a balanced norm of the NIPG method on Shishkin mesh for a reaction diffusion problem, Appl. Math. Comput., 444 (2023), 127828. https://doi.org/10.1016/j.amc.2022.127828

[22] T. Valanarasu, N. Ramanujan, Asymptotic initial-value method for singularly-perturbed boundary-value problems for second-order ordinary differential equations, J. Optim. Theory Appl., 116 (2003), 167-182. https://doi.org/10.1023/A:1022118420907

[23] J. Zhang, L. Ge, J. Kouatchou, A two colorable fourth-order compact difference scheme and parallel iterative solution of the 3D convection diffusion equation, Math. Comput. Simul., 54 (2000), 65-80. https://doi.org/10.1016/S0378-4754(00)00205-6

[24] K. Pan, D. He, H. Hu, An extrapolation cascadic multigrid method combined with a fourth-order compact scheme for 3D Poisson equation, J. Sci. Comput., 70 (2017), 1180-1203. https://doi.org/10.1007/s10915-016-0275-9

[25] J. Zhang, An explicit fourth-order compact finite difference scheme for three-dimensional convection-diffusion equation, Commun. Numer. Methods Eng., 14 (1998), 209-218. https://doi.org/10.1002/(SICI)1099-0887(199803)14:3%3C209::AID-CNM139%3E3.0.CO;2-P

[26] J. Kouatchou, J. Zhang, Optimal injection operator and high order schemes for multigrid solution of 3D Poisson equation, Int. J. Comput. Math., 76 (2000), 173-190. https://doi.org/10.1080/00207160008805018

[27] J. Zhang, Fast and high accuracy multigrid solution of the three dimensional Poisson equation, J. Comput. Phys., 143 (1998), 449-461. https://doi.org/10.1006/jcph.1998.5982

[28] M. M. Gupta, J. Zhang, High accuracy multigrid solution of the 3D convection-diffusion equation, Appl. Math. Comput., 113 (2000), 249-274. https://doi.org/10.1016/S0096-3003(99)00085-5

[29] M. M. Gupta, J. Kouatchou, Symbolic derivation of finite difference approximations for the three-dimensional Poisson equation, Numer. Methods Partial Differ. Equations, 14 (1998), 593-606. https://doi.org/10.1002/(SICI)1098-2426(199809)14:5%3C593::AID-NUM4%3E3.0.CO;2-D

[30] W. F. Spotz, G. F. Carey, A high-order compact formulation for the 3D Poisson equation, Numer. Methods Partial Differ. Equations, 12 (1996), 235-243. https://doi.org/10.1002/(SICI)1098-2426(199603)12:2%3C235::AID-NUM6%3E3.0.CO;2-R

[31] Y. Kyei, J. P. Roop, G. Tang, A family of sixth-order compact finite-difference schemes for the three-dimensional Poisson equation, Adv. Numer. Anal., 2010 (2010), 352174. https://doi.org/10.1155/2010/352174

[32] Y. Ge, Multigrid method and fourth-order compact difference discretization scheme with unequal meshsizes for 3D Poisson equation, J. Comput. Phys., 229 (2010), 6381-6391. https://doi.org/10.1016/j.jcp.2010.04.048

[33] Y. Ge, F. Cao, J. Zhang, A transformation-free HOC scheme and multigrid method for solving the 3D Poisson equation on nonuniform grids, J. Comput. Phys., 234 (2013), 199-216. https://doi.org/10.1016/j.jcp.2012.09.034

[34] A. Mayo, A. Greenbaum, Fourth order accurate evaluation of integrals in potential theory on exterior 3D regions, J. Comput. Phys., 220 (2007), 900-914. https://doi.org/10.1016/j.jcp.2006.05.042

[35] G. Saldanha, Single cell high order difference schemes for Poisson’s equation in three variables, Appl. Math. Comput., 186 (2007), 548-557. https://doi.org/10.1016/j.amc.2006.07.126

[36] R. F. Boisvert, A fourth-order-accurate Fourier method for the Helmholtz equation in three dimensions, ACM Trans. Math. Software, 13 (1987), 221-234. https://doi.org/10.1145/29380.29863

[37] W. Hackbusch, Z. P. Nowak, On the fast matrix multiplication in the boundary element method by panel clustering, Numer. Math., 54 (1989), 463-491. https://doi.org/10.1007/BF01396324

[38] Y. Wang, J. Zhang, Fast and robust sixth-order multigrid computation for the three-dimensional convection-diffusion equation, J. Comput. Appl. Math., 234 (2010), 3496-3506. https://doi.org/10.1016/j.cam.2010.05.022

[39] C. Bacuta, D. Hayes, J. Jacavage, Efficient discretization and preconditioning of the singularly perturbed reaction-diffusion problem, Comput. Math. Appl., 109 (2022), 270-279. https://doi.org/10.1016/j.camwa.2022.01.031

[40] H. Chen, Y. Su, B. D. Shizgal, A direct spectral collocation Poisson solver in polar and cylindrical coordinates, J. Comput. Phys., 160 (2000), 453-469. https://doi.org/10.1006/jcph.2000.6461

[41] G. Biros, L. Ying, D. Zorin, A fast solver for the Stokes equations with distributed forces in complex geometries, J. Comput. Phys., 193 (2004), 317-348. https://doi.org/10.1016/j.jcp.2003.08.011

[42] A. Greenbaum, A. Mayo, Rapid parallel evaluation of integrals in potential theory on general three-dimensional regions, J. Comput. Phys., 145 (1998), 731-742. https://doi.org/10.1006/jcph.1998.6048

[43] L. Ying, G. Biros, D. Zorin, A high-order 3D boundary integral equation solver for elliptic PDEs in smooth domains, J. Comput. Phys., 219 (2006), 247-275. https://doi.org/10.1016/j.jcp.2006.03.021

[44] A. N. Marques, J. C. Nave, R. R. Rosales, A correction function method for Poisson problems with interface jump conditions, J. Comput. Phys., 230 (2011), 7567-7597. https://doi.org/10.1016/j.jcp.2011.06.014

[45] J. R. Phillips, J. K. White, A precorrected-FFT method for electrostatic analysis of complicated 3D structures, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst., 16 (1997), 1059-1072. https://doi.org/10.1109/43.662670

[46] W. Ying, C. S. Henriquez, A Kernel-free boundary integral method for elliptic boundary value problems, J. Comput. Phys., 227 (2007), 1046-1074. https://doi.org/10.1016/j.jcp.2007.08.021

[47] E. Braverman, M. Israeli, A. Averbuch, L. Vozovoi, A fast 3D Poisson solver of arbitrary order accuracy, J. Comput. Phys., 144 (1998), 109-136. https://doi.org/10.1006/jcph.1998.6001

[48] W. Ying, W. C. Wang, A Kernel-free boundary integral method for implicitly defined surfaces, J. Comput. Phys., 252 (2013), 606-624. https://doi.org/10.1016/j.jcp.2013.06.019

[49] Y. Xie, W. Ying, A fourth-order Kernel-free boundary integral method for the modified Helmholtz equation, J. Sci. Comput., 78 (2019), 1632-1658. https://doi.org/10.1007/s10915-018-0821-8

[50] A. N. Marques, J. C. Nave, R. R. Rosales, High order solution of Poisson problems with piecewise constant coefficients and interface jumps, J. Comput. Phys., 335 (2017), 497-515. https://doi.org/10.1016/j.jcp.2017.01.029

[51] O. Steinbach, Numerical Approximation Methods for Elliptic Boundary Value Problems: Finite and Boundary Elements , Springer, 2007. https://doi.org/10.1007/978-0-387-68805-3

[52] A. Frangi, A. di Gioia, Multipole BEM for the evaluation of damping forces on MEMS, Comput. Mech., 37 (2005), 24-31. https://doi.org/10.1007/s00466-005-0694-1

[53] A. A. Samarskii, The Theory of Difference Schemes , CRC Press, 2001. https://doi.org/10.1201/9780203908518

[54] X. Feng, H. Song, T. Tang, J. Yang, Nonlinear stability of the implicit-explicit methods for the Allen-Cahn equation, Inverse Probl. Imaging, 7 (2013), 679-695. https://doi.org/10.3934/ipi.2013.7.679