Analytical and numerical solutions of the convection-diffusion-reaction equations applying the differential transformation method and the Crank-Nicolson method, along with stability analysis and truncation error analysis

Habibah, Ummu; Rohman, Zalfa Camilla; Dewanti, Yashinta Novena; Hananti, Arfi Nadhifa

doi:10.22034/cmde.2025.64919.2958

Analytical and numerical solutions of the convection-diffusion-reaction equations applying the differential transformation method and the Crank-Nicolson method, along with stability analysis and truncation error analysis

Document Type : Research Paper

Authors

Department of Mathematics, Faculty of Science, Technology, and Mathematics, Brawijaya University, Indonesia.

10.22034/cmde.2025.64919.2958

Abstract

This study presents a unified approach for solving convection-diffusion-reaction equations by integrating the Differential Transformation Method (DTM) for analytical approximations with the Crank-Nicolson numerical scheme. The DTM is employed to derive an analytical solution, while the Crank-Nicolson method is used to compute the numerical solution. The results demonstrate that the analytical solution obtained via DTM is identical to the exact solution. Furthermore, the stability of the Crank-Nicolson numerical scheme is assessed using Von-Neumann stability analysis, confirming that the method is unconditionally stable. The local truncation error is determined via Taylor series expansion to establish its order of accuracy. This analysis reveals that the Crank-Nicolson scheme for the convection-diffusion-reaction equation exhibits a local truncation error of order $O(h^2+k^2)$, ensuring a second-order accurate scheme. Numerical simulations are conducted for various parameter values to examine their impact on the solution. The simulation results demonstrate the gradual transport of the substance from high to low concentration regions, observed through the diminishing displacement of material along the $x$-axis. Further numerical experiments investigate the effects of different values of $h$ and $k$. The results indicate a direct correlation between decreasing values of $h$ and $k$ and a reduction in the average error, underscoring the method’s accuracy and efficiency.

Keywords

Main Subjects

Numerical approximation in context of PDEs

References

[1] D. Adak and E. Natarajan, A unified analysis of nonconforming virtual element methods for convection diffusion reaction problem, (2016).
[2] F. Ayaz, Solutions of the system of differential equations by differential transform method, Appl. Math. Comput.,147 (2004), 547–567.
[3] L. A. Albesa, M. B. Garcia, M. L. P. Ferrando, and S. C. C. Navarrete, A study about the solution of convection- diffusion-reaction equation with Danckwerts boundary conditions by analytical, method of lines, and Crank- Nicolson techniques, Math. Methods Appl. Sci., 46 (2021), 2133–2164.
[4] M. Al-Sawalha and M. Noorani, Application of the differential transformation method for the solution of the hyperchaotic Rössler system, Commun. Nonlinear Sci. Numer. Simul., 14 (2009), 1509–1514.
[5] S. Aydın and C. Çiftçi, The finite element method solution of variable diffusion coefficient convection-diffusion equations, AIP Conf. Proc., 1479 (2012), 1835–1838.
[6] J. Bonilla and S. Badia, Maximum-principle preserving space–time isogeometric analysis, Comput. Methods Appl. Mech. Engrg., 354 (2019), 422–440.
[7] M. Balàzsovà, M. Feistauer, and M. Vlasàk, Stability of the ALE space-time discontinuous Galerkin method for nonlinear convection-diffusion problems in time-dependent domains, ESAIM Math. Model. Numer. Anal., 52 (2018), 2327–2356.
[8] M. Benedetto, S. Berrone, A. Borio, S. Pieraccini, and S. Scialò, Order preserving SUPG stabilization for the virtual element formulation of advection–diffusion problems, Comput. Methods Appl. Mech. Engrg., 311 (2016), 18–40.
[9] S. Berrone, A. Borio, and G. Manzini, SUPG stabilization for the nonconforming virtual element method for advection–diffusion–reaction equations, Comput. Methods Appl. Mech. Engrg., 340 (2018), 500–529.
[10] R. David, P. Saucez, J. L. Vasel, and A. V. Wouwer, Modeling and numerical simulation of secondary settlers: A method of lines strategy, Water Res., 43 (2009), 319–330.
[11] D. A. Garzòn-Alvarado, C. H. Galeano, and J. M. Mantilla, Computational examples of reaction–convection–diffusion equations solution under the influence of fluid flow: First example, Appl. Math. Model., 36 (2012), 5029–5045.
[12] W. Hundsdorfer and J. Verwer, Numerical solution of time-dependent advection-diffusion-reaction equations, Springer Series in Computational Mathematics, Springer, Berlin, 2003.
[13] A. S. V. R. Kanth and K. Aruna, Differential transform method for solving the linear and nonlinear Klein-Gordon equation, Comput. Phys. Commun., 180 (2009), 708–711.
[14] Y. Khan, Solving certain classes of Lane-Emden type equation using the differential transformation method, Adv. Differ. Equ., 2012 (2012), 174.
[15] E. Khatizah, P. T. Karima, and D. I. Astuti, Aplikasi Metode Transformasi Diferensial pada Sistem Persamaan Diferensial, Jurnal Matematika dan Aplikasinya, 14(2) (2015), 1-8.
[16] S. A. Lima, M. Kamrujjaman, and M. S. Islam, Numerical solution of convection-diffusion-reaction equations by a finite element method with error correlation, AIP Adv., 11 (2021), 085225.
[17] J. Makungu, H. Haario, and W. C. Mahera, A generalized 1-dimensional particle transport method for convection- diffusion-reaction model, Afrika Mat., 23 (2012), 21–39.
[18] N. Mphephu, Numerical solution of 1-D convection-diffusion-reaction equation, Master’s thesis, University of Venda, African Institute for Mathematical Sciences, 2013.
[19] A. Niemi, N. Collier, and V. Calo, Automatically stable discontinuous Petrov–Galerkin methods for stationary transport problems: quasi-optimal test space norm, Comput. Math. Appl., 66 (2013), 2096–2113.
[20] H. Nayar and P. A. Phiri, The solution of linearised Korteweg-De Vries equation using the differential transform method, J. Adv. Math. Comput. Sci., 29 (2018), 1–10.
[21] S. Wu and J. Xu, Simplex-averaged finite element methods for H(grad), H(curl), and H(div) convection-diffusion problems, SIAM J. Numer. Anal., 58 (2020), 884–906.
[22] Q. Wang and Z. Zhou, SUPG-stabilized virtual element method for optimal control problem governed by a convection-dominated diffusion equation, Entropy, 23 (2021), 723.