NUMERICAL SOLUTIONS OF REACTION-DIFFUSION EQUATION SYSTEMS WITH TRIGONOMETRIC QUINTIC B-SPLINE COLLOCATION ALGORITHM

Abstract

In this study, trigonometric quintic B-spline collocation method is constructed for computing numerical solutions of the reaction-diffusion system (RDS). Schnakenberg, Gray-Scott and Brusselator models are special cases of reaction-diffusion systems considered as examples in this paper. Crank-Nicolson formulae is used for the time discretization of the generalized RDS and the nonlinear terms in time-discretized form of RDS are linearized using the Taylor expansion. The fully integration of the generalized system is carried out using the collocation method based on the trigonometric quintic B-splines. The method is tested on different problems to illustrate the accuracy. The error norms are calculated for the linear problem whereas the relative error is given for nonlinear problems. Both simple and easy B-spline algorithms are illustrated to give the solutions of RDS and also the graphical representation of the efficient solutions are presented for the nonlinear RDSs. Combination of the quintic B-splines and the collocation method is shown to present numerical solutions of the RDS successfully. With the presented method, it is possible to get approximate solutions as well as their derivatives up to an order of four on the problem domain.

Keywords

• Reaction-diffusion; Collocation; B-spline; Finite element method
• Brusselator
• Schnakenberg
• Gray-Scott

References

1. [1] Ruuth SJ. Implicit-explicit methods for reaction-diffusion problems in pattern formation. J Math Bio 1995; 34: 148-176.
2. [2] Madzvamuse A, Wathen AJ, Maini PK. A moving grid finite element method applied to a biological pattern generator. J Comp Phys 2003; 190: 478-500.
3. [3] Zegeling PA, Kok HP. Adaptive moving mesh computations for reaction-diffusion systems. J Comp App Math 2004; 168: 519-528.
4. [4] Ropp DL, Shadid JN, Ober CC. Studies of the accuracy of time integration methods for reaction-diffusion equations. J Comp Phys 2004; 194(2): 544-574.
5. [5] Ropp DL, Shadid JN. Stability of operator splitting methods for systems with indefinite operators: reaction-diffusion systems. J Comp Phys 2005; 203(2): 449-466.
6. [6] Chou CS, Zhang Y, Zhao R, Nie Q. Numerical Methods for Stiff Reaction-Diffusion Systems. Discrete and Continuous Dynamical Systems-Series B 2007;7: 515-525.
7. [7] Yadav OP, Jiwari R. A finite element approach for analysis and computational modelling of coupled reaction diffusion models. Numerical Methods for Partial Differential Equations 2019; 35 (2): 830-850.
8. [8] Mittal RC, Rohila R. Numerical Simulation of Reaction-Diffusion Systems by Modified Cubic B-spline Differential Quadrature Method. Chaos, Solitons and Fractals 2016; 92 (1): 1339-1351.

9. [9] Jiwari R, Singh S, Kumar A. Numerical simulation to capture the pattern formation of coupled reaction-diffusion models. Chaos, Solitons and Fractals 2017; 103: 422-439.
10. [10] Jiwari R, Gerisch A. A local radial basis function differential quadrature semidiscretisation technique for the simulation of time-dependent reaction-diffusion problems. Eng Comp 2016; 38 (6): 2666-2691.
11. [11] Schoenberg IJ. On Trigonometric Spline Interpolation. J Math Mech 1964;13: 795-826.
12. [12] Hamid NNA, Majid AA, Ismail AIM. Cubic Trigonometric B-spline Applied to Linear Two-Point Boundary Value Problems of Order Two. World Academy of Science, Engineering and Technology 2010; 47: 478-803.
13. [13] Gupta Y, Kumar MA Computer Based Numerical Method for Singular Boundary Value Problems. Int J Comp App 2011;1: 0975-8887.
14. [14] Zin SM, Abbas M, Majid AA, Ismail AIM. A New Trigonometric Spline Approach to Numerical Solution of Generalized Nonlinear Klien-Gordon Equation. PLoS one, 2014;9(5): e95774.
15. [15] Abbas M, Majid AA, Ismail AIM, Rashid A. Numerical Method Using Cubic Trigonometric B-Spline Technique for Nonclassical Diffusion Problems. Abstract and Applied Analysis, 2014; Article ID 849682.
16. [16] Yağmurlu NM, Karakaş AS. Numerical solutions of the equal width equation by trigonometric cubic B‐spline collocation method based on Rubin–Graves type linearization. Num Meth Part Diff Eqs, 2020; 36(5): 1170-1183.
17. [17] Nikolis A. Numerical Solutions of Ordinary Differential Equations with Quadratic Trigonometric Splines. App Maths E-Notes 2004;4: 142-149.
18. [18] Nikolis A, Seimenis I. Solving Dynamical Systems with Cubic Trigonometric Splines. App Maths E-notes, 2005; 5: 116-123.
19. [19] Sahin A. Numerical solutions of the reaction-diffusion equations with B-spline finite element method. PhD, Eskişehir Osmangazi University, Eskişehir, Turkey, 2009.
20. [20] Ersoy O, Dag I. Numerical Solutions of the Reaction Diffusion System by Using Exponential Cubic B-spline Collocation Algorithms. Open Phys 2015; 13: 414-427.
21. [21] Onarcan AT, Adar N, Dag I. Trigonometric cubic B-spline collocation algorithm for numerical solutions of reaction–diffusion equation systems. Comp App Math 2018; 37: 6848–6869.
22. [22] Tok-Onarcan A, Adar N, Dag I. Wave simulations of Gray-Scott reaction-diffusion system. Math Meth in the App Sci 2019; 19: 5566–5581.
23. [23] Kaur N, Joshi V. Numerical solution to the Gray-Scott Reaction-Diffusion equation using Hyperbolic B-spline. J Phys: Conference Series 2022: 2267.
24. [24] Hepson ÖE, Yiğit G, Allahviranloo T. Numerical simulations of reaction–diffusion systems in biological and chemical mechanisms with quartic-trigonometric B-splines. Comp App Math 2021; 40(4): 1-23.
25. [25] Hepson ÖE, Numerical solutions of the Gardner equation via trigonometric quintic B-spline collocation method, Sakarya Üniversitesi Fen Bilimleri Enstitüsü Dergisi 2018; 22 (6): 1576–1584.
26. [26] Chandrasekharan Nair L, Awasthi A. Quintic trigonometric spline based numerical scheme for nonlinear modified Burgers’ equation, Num Meth for Part Diff Eq 2019; 35(3): 1269–1289.
27. [27] Khater MMA, Nisar KS, Mohamed MS. Numerical investigation for the fractional nonlinear space-time telegraph equation via the trigonometric Quintic B-spline scheme. Math Meth in App Sci 2021; 44(6); 4598–4606.
28. [28] Alam MP, Kumar D, Khan A. Trigonometric quintic B-spline collocation method for singularly perturbed turning point boundary value problems. Int J Comp Math 2021;98(5): 1029–1048.
29. [29] Prigogine I, Lefever R. Symmetry Breaking Instabilities in Dissipative Systems. J Chem Phys 1968; 48:1695-1700.
30. [30] Schnakenberg J. Simple Chemical Reaction Systems with Limit Cycle Behavior. J Theo Bio 1979; 81: 389-400.
31. [31] Gray P, Scott SK. Autocatalytic Reactions in The Isothermal, Continuous Stirred Tank Reactor: Oscillations and Instabilities in the system A+2B 3B, B C. Chem Eng Sci, 1984;39: 1087-1097.
32. [32] Craster RV, Sassi R. Spectral Algorithms for Reaction-Diffusion Equations. Technical Report. Note del Polo 2006: No:99.
33. [33] Rubin SG, Graves Jr RA. A cubic spline approximation for problems in fluid mechanics. Technical Report, Nasa, Washington, 1975.