The numerical solution of the nonlinear system of stiff differential equations by the modified matrix-exponential method

Fattahi, Mohammad; Matinfar, Mashallah

doi:10.22075/ijnaa.2023.30796.4496

The numerical solution of the nonlinear system of stiff differential equations by the modified matrix-exponential method

Document Type : Research Paper

Authors

Mohammad Fattahi

Mashallah Matinfar

Department of Mathematics, Science of Mathematics Faculty, University of Mazandaran, Babolsar, Iran

10.22075/ijnaa.2023.30796.4496

Abstract

In this paper, the modified matrix exponential (MME) method under zero-order hold (ZOH) assumption, is applied to solve systems of stiff ordinary differential equations. Some examples are given to illustrate the accuracy and effectiveness of the method. We compare our results with results obtained by matrix exponential (ME) method and by the Matlab ode23 solver.

Keywords

Modified matrix exponential

Matrix exponential

Abel equation of the second kind

Nonlinear differential equations

Jacobian matrix

[1] A. Abdi, S.A. Hosseini and H. Podhaisky, Adaptive linear barycentric rational finite differences method for stiff ODEs, J. Comput. Appl. Math. 357 (2019), 204–214.

[2] F. Allffi-Pentini, V.D. Fonzo and V. Paris, A novel algorithm for the numerical integration of systems of ordinary differential equations arising in chemical problem, J. Math. Chem. 33 (2003), 1–15.

[3] M. Calvo, S. Gonzalez-Pinto and J.I. Montijano, Runge–Kutta methods for the numerical solution of stiff semilinear system, BIT Numer. Math. 40 (2000), 611–639.

[4] M.S.H. Chowdhury, I. Hashim and Md. Alal Hosen, Solving linear and non-linear stiff system of ordinary differential equations by multistage Adomian decomposition method, Proc. The Third Intl. Conf. Adv. Appl. Sci. Environ. Technol. 125 (2015), 79–82.

[5] M. Fattahi and M. Matinfar, The numerical solution of the second kind of Abel equations by the modified matrix-exponential method, Int. J. Nonlinear Anal. Appl. In Press, (2023), 1–7.

[6] G.Y. Kulikov, R. Weiner and G.M. Amiraliyev, Doubly quasi-consistent fixed-step size numerical integration of stiff ordinary differential equations with implicit two-step peer methods, J. Comput. Appl. Math. 340 (2018), 256–275.

[7] U. Lepik, Haar wavelet method for solving stiff differential equations, Math. Model. Anal. 14 (2009), 79–82.

[8] S. Li, Canonical Euler splitting method for nonlinear composite stiff evolution equations, Appl. Math. Comput. 289 (2016), 220–236.

[9] Y. Park and K.T. Chong, The numerical solution of the point kinetics equation using the matrix exponential method, Ann. Nuclear Energy 55 (2013), 42–48.