Skip to main content
Log in

A method for exponential propagation of large systems of stiff nonlinear differential equations

  • Published:
Journal of Scientific Computing Aims and scope Submit manuscript

Abstract

A new time integrator for large, stiff systems of linear and nonlinear coupled differential equations is described. For linear systems, the method consists of forming a small (5–15-term) Krylov space using the Jacobian of the system and carrying out exact exponential propagation within this space. Nonlinear corrections are incorporated via a convolution integral formalism; the integral is evaluated via approximate Krylov methods as well. Gains in efficiency ranging from factors of 2 to 30 are demonstrated for several test problems as compared to a forward Euler scheme and to the integration package LSODE.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Similar content being viewed by others

References

  • Arneodo, A., and Elezgaray, J. (1990).Phys. Lett. A. 143, 25.

    Google Scholar 

  • Arnoldi, W. E. (1951).Q. Appl. Math. 9, 17.

    Google Scholar 

  • Brown, P., and Hindmarsh, A. C. (1986).SIAM J. Num. Anal. 23, 610.

    Google Scholar 

  • Brown, P., and Hindmarsh, A. C. (1989).J. Appl. Math. Comp. 31, 40.

    Google Scholar 

  • Byrne, G. D., and Hindmarsh, A. C. (1987).J. Comput. Phys. 70, 1.

    Google Scholar 

  • Christodoulou, K. N. and Scriven, L. E. (1988).J. Sci. Comput. 3, 355.

    Google Scholar 

  • Cullum, J., and Willoughby, R. (1985).Lanczos Algorithms for Large Symmetric Eigenvalue Computations, Birkhauser, Boston.

    Google Scholar 

  • Gear, C. W. (1971).Numerical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, Englewood Cliffs, New Jersey.

    Google Scholar 

  • Gear, C. W., and Saad, Y. (1983).SIAM J. Sci. Stat. Comput. 4, 583.

    Google Scholar 

  • Goldhirsch, I., Orszag, S. A., and Maulik, B. K. (1987).J. Sci. Comput. 2, 33.

    Google Scholar 

  • Hindmarsh, A. C. (1983). In R. S. Steplemanet al. (eds.),Scientific Computing, North-Holland, Amsterdam, pp. 55–64.

    Google Scholar 

  • Hindmarsh, A. C., and Brown, P. (1987). In R. Vichnevetsky and R. S. Stepleman (eds.),Advances in Computer Methods for Partial Differential Equations — VI, IMACS, New Brunswick, pp. 355–362.

    Google Scholar 

  • Jahnke, W., Skaggs, W. E., and Winfree, A. T. (1989).J. Phys. Chem. 93, 740.

    Google Scholar 

  • Leforestier, C, Bisseling, R., Cerjan, C, Feit, M., Friesner, R., Guldberg, A., Hammerich, A., Jolicard, G., Karrlein, W., Meyer, H. D., Lipkin, N. Roncero, O., and Kosloff, R. (1989).J. Comput. Phys., submitted.

  • Park, T., and Light, J. C. (1986).J. Chem. Phys. 85, 5870.

    Google Scholar 

  • Saad, Y. (1980).Lin. Alg. Appl. 34, 269.

    Google Scholar 

Download references

Author information

Authors and Affiliations

Authors

Rights and permissions

Reprints and permissions

About this article

Cite this article

Friesner, R.A., Tuckerman, L.S., Dornblaser, B.C. et al. A method for exponential propagation of large systems of stiff nonlinear differential equations. J Sci Comput 4, 327–354 (1989). https://doi.org/10.1007/BF01060992

Download citation

  • Received:

  • Issue Date:

  • DOI: https://doi.org/10.1007/BF01060992

Key words

Navigation