Comparison of some recent numerical methods for initial-value problems for stiff ordinary differential equations

  • Authors:
  • Shirley Abelman;Kailash C. Patidar

  • Affiliations:
  • School of Computational and Applied Mathematics, University of the Witwatersrand, Johannesburg, Private Bag-3, Wits 2050, South Africa;Department of Mathematics and Applied Mathematics, University of the Western Cape, Private Bag X17, Bellville 7535, South Africa

  • Venue:
  • Computers & Mathematics with Applications
  • Year:
  • 2008

Quantified Score

Hi-index 0.09

Visualization

Abstract

We consider the combustion equation as one of the candidates from the class of stiff ordinary differential equations. A solution over a length of time that is inversely proportional to @d0 (where @d0 is a small disturbance of the pre-ignition state) is sought. This problem has a transient at the midpoint of the integration interval. The solution changes from being non-stiff to stiff, and afterwards becomes non-stiff again. We provide its asymptotic and numerical solution obtained via a variety of methods. Comparisons are made for the numerical results which we obtain with the MATLAB ode solvers (ode45, ode15s and ode23s) and some nonstandard finite difference methods. Results corresponding to standard finite difference method are also presented. Furthermore, the discussion on these approaches along with the others, provides several open problems for new and young researchers.