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] Ruuth SJ. Implicit-explicit methods for reaction-diffusion problems in pattern formation. J Math Bio 1995; 34: 148-176.
- [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] Zegeling PA, Kok HP. Adaptive moving mesh computations for reaction-diffusion systems. J Comp App Math 2004; 168: 519-528.
- [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] 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] 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] 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] 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] 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] 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.