论文信息 - The NUMOL solution of time-dependent PDEs using DESI Runge-Kutta formulae

The NUMOL solution of time-dependent PDEs using DESI Runge-Kutta formulae

In a recent publication Butcher and Cash introduced a new class of L-stable implicit Runge-Kutta formulae in an attempt to eliminate the disadvantages of SIRK and DIRK formulae and to combine their advantages. These formulae, namely Diagonally Extended Singly Implicit formulae (DESI), have been investigated by the present author and have been implemented in an experimental code in a variable-stepsize/variable-order mode. Numerical experimentation on a large number of stiff systems has shown that the DESI code is much more efficient than STRIDE (which is an implementation of SIRK formulae) and competitive with BDF implementations. In this paper we investigate the efficiency of the DESI code on large dimension stiff systems with banded Jacobians, which arise from the discretization of PDEs with the Numerical Method of Lines (NUMOL). The results obtained by the DESI code are compared with the results obtained by the widely used BDF code LSODE. These results indicate that the current implementation of DESI performs satisfactorily and, in particular, is competitive with LSODE.

Michail Diamantakis

[1] Kevin Burrage,et al. An implementation of singly-implicit Runge-Kutta methods , 1980 .

[2] John D. Pryce,et al. Two FORTRAN packages for assessing initial value methods , 1987, TOMS.

[3] A. Hindmarsh,et al. Stiff ode slovers: a review of current and coming attractions , 1987 .

[4] A. Hindmarsh. Avoiding BDF stability barriers in the MOL solution of advection-dominated problems , 1995 .

[5] S. McKee. Two species diffusion in a sphere with inter-species linkage at the boundary , 1982 .

[6] A. Prothero,et al. On the stability and accuracy of one-step methods for solving stiff systems of ordinary differential equations , 1974 .

[7] J. H. Wilkinson. The algebraic eigenvalue problem , 1966 .

[8] John C. Butcher,et al. Towards efficient Runge-Kutta methods for stiff systems , 1990 .

[9] L. Shampine. ODE solvers and the method of lines , 1994 .

[10] R. Alexander. Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[11] T. E. Hull,et al. Comparing numerical methods for stiff systems of O.D.E:s , 1975 .

[12] John C. Butcher,et al. A Transformed implicit Runge-Kutta Method , 1979, JACM.

[13] J. R. Cash,et al. Two New Finite Difference Schemes for Parabolic Equations , 1984 .

[14] Kevin Burrage,et al. A special family of Runge-Kutta methods for solving stiff differential equations , 1978 .