Amplitude-shape method for solving partial differential equations of chemical kinetics Original Research Article

The original amplitude-shape method [N. Parumasur, Amplitude-Shape Method for the Numerical Solution of Ordinary Differential Equations, Ph.D. Thesis, Department of Mathematics, University of Natal, Durban, 1997; Quaestiones Math. 24 (2001) 69] for solving ordinary differential equations consisted in transforming certain systems of stiff ordinary differential equations into two subsystems consisting of a rapidly changing amplitude function and slowly varying shape vector, respectively. In this paper we consider a new formulation of the amplitude and shape equations, which is particularly well suited for certain classes of problems occurring in chemical kinetics. The resulting equations are treated with a partitioned method consisting of the Rosenbrock method and the trapezoidal rule. It is shown that the use of the proposed new version of the amplitude-shape method leads to a substantial gain in computational efficiency in comparison with the conventional implicit or semi-implicit schemes employed in stiff solvers.