Efficient spectral-Galerkin method for the two-dimensional Gray-Scott model
DOI: https://doi.org/10.3846/mma.2026.24728Abstract
The Gray-Scott (GS) model is a nonlinear reaction diffusion system widely used in applied sciences. This paper delves into the investigation of a spectral method for the GS model with homogeneous Neumann boundary conditions. The proposed numerical scheme integrates a spectral approach based on generalized Jacobi polynomials for spatial discretization with the two-step backward differentiation formula (BDF2) for temporal discretization. We prove the boundedness, generalized stability, and convergence of this new method. Extensive numerical results demonstrate the efficiency of the new proposed scheme and provide numerical validation of the theoretical analysis. The key advantages of our new approach are twofold: (i) it utilizes generalized Jacobi polynomials with indices α = β = −2, which naturally satisfy the boundary conditions and thereby not only simplify the theoretical analysis but also yield a well-conditioned discrete system that enhances computational efficiency; and (ii) the numerical errors exhibit exponential decay in space.
Keywords:
spectral method, Gray-Scott model, generalized Jacobi polynomials, nonlinear system, reaction-diffusion systemHow to Cite
Share
License
Copyright (c) 2026 The Author(s). Published by Vilnius Gediminas Technical University.
This work is licensed under a Creative Commons Attribution 4.0 International License.
References
M. Abbaszadeh and M. Dehghan. A reduced order finite difference method for solving space-fractional reaction-diffusion systems: The Gray-Scott model. Eur. Phys. J. Plus, 134:620, 2019. https://doi.org/10.1140/epjp/i2019-12951-0
M. Abbaszadeh, M. Dehghan and I.M. Navon. A POD reduced-order model based on spectral Galerkin method for solving the space-fractional Gray-Scott model with error estimate. Eng. Comput., 38:2245–2268, 2022. https://doi.org/10.1007/s00366-020-01195-5
A.A.A. Amin and D.S. Mashat. Analysis of Gray-Scott’s model numerically. Amer. J. Comput. Math., 11:273–288, 2021. https://doi.org/10.4236/ajcm.2021.114018
R. Castelli. Rigorous computation of non-uniform patterns for the 2-dimensional Gray-Scott reaction-diffusion equation. Acta Appl. Math., 151:27–52, 2017. https://doi.org/10.1007/s10440-017-0101-x
M. Dehghan and M. Abbaszadeh. Numerical study of three-dimensional Turing patterns using a meshless method based on moving Kriging element free Galerkin (EFG) approach. Comput. Math. Appl., 72:427–454, 2016. https://doi.org/10.1016/j.camwa.2016.04.038
I.R. Epstein and J.A. Pojman. An introduction to nonlinear chemical dynamics: oscillations, waves, patterns, and chaos. Oxford University Press, New York, 1998. https://doi.org/10.1093/oso/9780195096705.001.0001
P. Gray and S.K. Scott. Autocatalytic reaction in the isothermal, continuous stirred tank reactor: Oscillations and instabilities in the system a + 2b → 3b; b → c. Chem. Eng. Sci., 39(6):1087–1097, 1984. https://doi.org/10.1016/0009-2509(84)87017-7
B. Guo, J. Shen and L. Wang. Generalized Jacobi polynomials/functions and their applications. Appl. Numer. Math., 59(5):1011–1028, 2009. https://doi.org/10.1007/s10915-005-9055-7
B. Guo, J. Shen and C. Xu. Spectral and pseudospectral approximations using Hermite functions: application to the Dirac equation. Adv. Comput. Math., 19:35–55, 2003. https://doi.org/10.1023/A:1022892132249
B. Guo, T. Sun and C. Zhang. Jacobi and Laguerre quasi-orthogonal approximation and related interpolations. Math. Comp., 82(281):413–441, 2013. https://doi.org/10.1090/S0025-5718-2012-02614-7
W. Hundsdorfer and J.G. Verwer. Numerical Solution of Time-Dependent Advection-Diffusion-Reaction Equations. Springer, Berlin, Heidelberg, 2010. https://doi.org/10.1007/978-3-662-09017-6
T. Kolokolnikov, M.J. Ward and J. Wei. Zigzag and breakup instabilities of stripes and rings in the two-dimensional Gray–Scott model. Stud. Appl. Math., 116:35–95, 2006. https://doi.org/10.1111/j.1365-2966.2005.0333.x
S. Kondo and T. Miura. Reaction-diffusion model as a framework for understanding biological pattern formation. Science, 329:1616–1620, 2010. https://doi.org/10.1126/science.1179047
J.S. McGough and K. Riley. Pattern formation in the Gray–Scott model. Nonlinear Anal. Real World Appl., 5:105–121, 2004. https://doi.org/10.1016/S1468-1218(03)00020-8
A. Munteanu and R.V. Solé. Pattern formation in noisy selfreplicating spots. Int. J. Bifurc. Chaos, 16(12):3679–3685, 2006. https://doi.org/10.1142/S0218127406017063
J.D. Murray. Mathematical Biology II: Spatial Models and Biomedical Applications. Springer, New York, NY, 2003. https://doi.org/10.1007/b98869
Y. Nishiura and D. Ueyama. A skeleton structure of self-replicating dynamics. Physica D, 130:73–104, 1999. https://doi.org/10.1016/S0167-2789(99)00010-X
S. Singh, R.C. Mittal, S.R. Thottoli and M. Tamsir. High-fidelity simulations for Turing pattern formation in multi-dimensional Gray–Scott reaction-diffusion system. Appl. Math. Comput., 452:128079, 2023. https://doi.org/10.1016/j.amc.2023.128079
T. Sun and B. Guo. Generalized Jacobi approximation in multiple dimensions and its applications. J. Sci. Comput., 55:327–350, 2013. https://doi.org/10.1007/s10915-012-9633-4
S.A.M. Tonekaboni and A. Shademani. Mathematical investigation of two dimensional pattern formation. Int. J. Biol. Res., 2(1):1–5, 2014. https://doi.org/10.14419/ijbr.v2i1.1502
L. Wang and J. Shen. Error analysis for mapped Jacobi spectral methods. J. Sci. Comput., 24(2):183–218, 2005. https://doi.org/10.1007/s10915-004-4613-y
W. Wang, Y. Lin, F. Yang, L. Zhang and Y. Tan. Numerical study of pattern formation in an extended Gray-Scott model. Commun. Nonlinear Sci. Numer. Simul., 16:2016–2026, 2011. https://doi.org/10.1016/j.cnsns.2010.09.002
J. Wei and M. Winter. Existence and stability of multiple-spot solutions for the Gray–Scott model in r2. Physica D, 176:147–180, 2003. https://doi.org/10.1016/S0167-2789(02)00743-1
O.P. Yadav and R. Jiwari. A finite element approach for analysis and computational modelling of coupled reaction diffusion models. Numer. Methods Partial Differential Equations, 35:830–850, 2019. https://doi.org/10.1002/num.22328
P.A. Zegeling and H.P. Kok. Adaptive moving mesh computations for reaction–diffusion systems. J. Comput. Appl. Math., 168:519–528, 2004. https://doi.org/10.1016/j.cam.2003.06.013
View article in other formats
CrossMark check
Published
Issue
Section
Copyright
Copyright (c) 2026 The Author(s). Published by Vilnius Gediminas Technical University.
License
This work is licensed under a Creative Commons Attribution 4.0 International License.