This study presents a comprehensive analysis and comparison of two numerical approaches for handling the nonlinear Lotka-Volterra interaction and diffusion system: the novel Taylor series method (NTSM) and the Galerkin finite element method (GFEM). We develop a semi-analytical NTSM method by iteratively obtaining coefficients from the governing partial differential equations and initial and boundary conditions to construct high-order Taylor series in space and time. The method is validated by comparison with an exact analytic solution, demonstrating high accuracy at the real-solution level (achieving errors below machine precision (approximately $ 10^{-16} $) in the tested scenarios) and rapid convergence, attributed to the analytic nature of the underlying solution. Furthermore, the use of linear (hierarchical) basis functions and Galerkin weighting with GFEM provides robust and physically acceptable approximations that rapidly achieve the stable equilibrium state of the system. Using Maple 18 and MATLAB, numerical tests indicate that the GFEM is robust and suitable for complex geometries and extended dynamics; however, the NTSM method significantly outperforms it in terms of accuracy and convergence speed for smooth and well-defined boundary problems in straightforward analytical problems. The results obtained show the accuracy and efficiency of NTSM in solving analytical and numerical problems, while GFEM is characterized by its strength in complex domains. Theoretical convergence analysis confirms an $ O(h^2) $ order convergence for GFEM and an exponential convergence for NTSM in the smooth test cases examined in this study, thus supporting the numerical results.
Citation: Badran Jasim Salim. Numerical solution of nonlinear reaction-diffusion systems by a novel Taylor series method versus the Galerkin finite element method[J]. AIMS Mathematics, 2026, 11(7): 21781-21798. doi: 10.3934/math.2026881
This study presents a comprehensive analysis and comparison of two numerical approaches for handling the nonlinear Lotka-Volterra interaction and diffusion system: the novel Taylor series method (NTSM) and the Galerkin finite element method (GFEM). We develop a semi-analytical NTSM method by iteratively obtaining coefficients from the governing partial differential equations and initial and boundary conditions to construct high-order Taylor series in space and time. The method is validated by comparison with an exact analytic solution, demonstrating high accuracy at the real-solution level (achieving errors below machine precision (approximately $ 10^{-16} $) in the tested scenarios) and rapid convergence, attributed to the analytic nature of the underlying solution. Furthermore, the use of linear (hierarchical) basis functions and Galerkin weighting with GFEM provides robust and physically acceptable approximations that rapidly achieve the stable equilibrium state of the system. Using Maple 18 and MATLAB, numerical tests indicate that the GFEM is robust and suitable for complex geometries and extended dynamics; however, the NTSM method significantly outperforms it in terms of accuracy and convergence speed for smooth and well-defined boundary problems in straightforward analytical problems. The results obtained show the accuracy and efficiency of NTSM in solving analytical and numerical problems, while GFEM is characterized by its strength in complex domains. Theoretical convergence analysis confirms an $ O(h^2) $ order convergence for GFEM and an exponential convergence for NTSM in the smooth test cases examined in this study, thus supporting the numerical results.
| [1] | P. E. Allaire, Basics of the finite element method, Brown Publishers, 1985. |
| [2] |
J. H. He, L. Verma, B. Pandit, A. K. Verma, R. P. Agarwal, A new Taylor series based numerical method: simple, and promising, J. Appl. Comput. Mech., 9 (2023), 1122–1134. https://doi.org/10.22055/jacm.2023.43228.4040 doi: 10.22055/jacm.2023.43228.4040
|
| [3] | V. A. Kostin, A. V. Kostin, B. J. Salim, On the well-posedness of the Cauchy problem for the generalized telegraph equations, Bull. South Ural State Univ., 7 (2014), 50–59. |
| [4] |
E. G. Tsega, A numerical solution of three-dimensional unsteady state heat conduction equation, Int. J. Math. Model. Comput., 11 (2021), 49–60. https://doi.org/10.30495/ijm2c.2021.684809 doi: 10.30495/ijm2c.2021.684809
|
| [5] |
B. J. Salim, D. M. Merie, Numerical and analytical solutions for a nonlinear system, Int. J. Differ. Equations, 2026 (2026), 9569337. https://doi.org/10.1155/ijde/9569337 doi: 10.1155/ijde/9569337
|
| [6] |
B. J. Salim, A numerical solution by finite volume method to solving three-dimensional system, Int. J. Anal. Appl., 24 (2026), 81. https://doi.org/10.28924/2291-8639-24-2026-81 doi: 10.28924/2291-8639-24-2026-81
|
| [7] |
Y. Lou, S. Martínez, P. Poláčik, Loops and branches of coexistence states in Lotka-Volterra competition model, J. Differ. Equations, 230 (2006), 720–742. https://doi.org/10.1016/j.jde.2006.04.005 doi: 10.1016/j.jde.2006.04.005
|
| [8] | R. L. Burden, J. D. Faires, Numerical analysis, 2 Eds., PWS-Kent, 1985. |
| [9] |
B. J. Salim, O. A. Jasim, CAS wavelets method to solve reaction-diffusion system and compare it with (G.F.E) method, Period. Tche Quim., 17 (2020), 1110–1123. https://doi.org/10.52571/PTQ.V17.N35.2020.91_SALIM_PGS_1110_1123.PDF doi: 10.52571/PTQ.V17.N35.2020.91_SALIM_PGS_1110_1123.PDF
|
| [10] |
M. D. Jasim, Using the centrodal mean to building parallel predictor corrector Runge Kutta method, AIP Conf. Proc., 2398 (2022), 0060004. https://doi.org/10.1063/5.0093381 doi: 10.1063/5.0093381
|
| [11] |
D. Mathale, P. G. Dlamini, M. Khumalo, Compact finite difference relaxation method for chaotic and hyperchaotic initial value systems, Comput. Appl. Math., 37 (2018), 5187–5202. https://doi.org/10.1007/s40314-018-0624-4 doi: 10.1007/s40314-018-0624-4
|
| [12] |
İ. Avcı, N. I. Mahmudov, Numerical solutions for multi-term fractional order differential equations with fractional Taylor operational matrix of fractional integration, Mathematics, 8 (2020), 96. https://doi.org/10.3390/math8010096 doi: 10.3390/math8010096
|
| [13] | S. S. Behzadi, T. Allahviranloo, An efficient numerical technique for solving the fractional Lotka-Volterra model, Comput. Methods Differ. Equations, 8 (2020), 249–264. |
| [14] |
J. Lu, Y. Liu, C. W. Shu, An oscillation-free discontinuous Galerkin method for scalar hyperbolic conservation laws, SIAM J. Numer. Anal., 59 (2021), 1299–1324. https://doi.org/10.1137/20M1354192 doi: 10.1137/20M1354192
|
| [15] | Z. Li, J. Li, A finite element analysis of traveling wave solutions for diffusive Lotka-Volterra systems, Commun. Nonlinear Sci. Numer. Simul., 75 (2019), 1–14. |
| [16] | V. Thomée, Finite element methods, In: C. Cancès, P. Omens, Finite element methods, Springer, 2018, 77–103. |
| [17] |
S. Paul, S. P. Mondal, P. Bhattacharya, Numerical solution of Lotka Volterra prey predator model by using Runge–Kutta–Fehlberg method and Laplace adomian decomposition method, Alex. Eng. J., 55 (2016), 613–617. https://doi.org/10.1016/j.aej.2015.12.026 doi: 10.1016/j.aej.2015.12.026
|
| [18] | G. Psihoyios, Recent advances in the numerical solution of nonlinear dynamical systems, ICNAAM, 2019. |
| [19] | C. Wei, Y. Li, Event-triggered control for stability of nonlinear stochastic delay systems driven by fractional Brownian motion, IAENG Int. J. Appl. Math., 55 (2025), 4393–4397. |
| [20] | G. S. Britto, A. X. G. Mohamed, A. A. Saraswathi, D. Iruthayaraj, Multi-series expansions via generalized alpha difference operators, IAENG Int. J. Appl. Math., 55 (2025), 3611–3618. |
| [21] |
T. Liu, H. Zhang, S. Wang, A new high-order compact CN-ADI scheme on graded meshes for three-dimensional nonlinear PIDEs with multiple weakly singular kernels, Appl. Math. Lett., 171 (2025), 109697. https://doi.org/10.1016/j.aml.2025.109697 doi: 10.1016/j.aml.2025.109697
|
| [22] |
Z. Chen, H. Zhang, H. Chen, ADI compact difference scheme for the two-dimensional integro-differential equation with two fractional Riemann–Liouville integral kernels, Fractal Fract., 8 (2024), 707. https://doi.org/10.3390/fractalfract8120707 doi: 10.3390/fractalfract8120707
|
| [23] |
X. Yang, Y. Shi, H. Zhang, Z. Zhang, Analysis of a new time-space two-grid method for the two-dimensional nonlinear nonlocal mobile/immobile transport model, J. Sci. Comput., 107 (2026), 96. https://doi.org/10.1007/s10915-026-03314-8 doi: 10.1007/s10915-026-03314-8
|
| [24] |
S. H. Rasouli, Extension of the Lotka-Volterra competition model, Hacettepe J. Math. Stat., 53 (2024), 1401–1407. https://doi.org/10.15672/hujms.1315963 doi: 10.15672/hujms.1315963
|
| [25] |
F. Tian, X. Tang, T. Xu, J. Yang, L. Li, A hybrid adaptive finite element phase-field method for quasi-static and dynamic brittle fracture, Int. J. Numer. Methods Eng., 120 (2019), 1108–1125. https://doi.org/10.1002/nme.6172 doi: 10.1002/nme.6172
|
| [26] | G. Strang, G. J. Fix, An analysis of the finite element method, Wellesley-Cambridge Press, 2008. https://doi.org/10.1137/1.9780980232707 |
| [27] | C. T. Kelley, Iterative methods for linear and nonlinear equations, Society for Industrial and Applied Mathematics, 1995. https://doi.org/10.1137/1.9781611970944 |
| [28] | J. P. Boyd, Chebyshev and fourier spectral methods, 2 Eds., Dover Publications, Inc., 2001. |
| [29] | K. W. Morton, D. F. Mayers, Numerical solution of partial differential equations, 2 Eds., Cambridge University Press, 2005. https://doi.org/10.1017/CBO9780511812248 |
| [30] | R. C. Smith, Uncertainty quantification: theory, implementation, and applications, Philadelphia, 2014. |
Figures(5) / Tables(7)
Badran Jasim Salim. Numerical solution of nonlinear reaction-diffusion systems by a novel Taylor series method versus the Galerkin finite element method[J]. AIMS Mathematics, 2026, 11(7): 21781-21798. doi: 10.3934/math.2026881
DownLoad: