Numerical integrators for stiff and highly oscillatory differential equations
HTML articles powered by AMS MathViewer
 by Simeon Ola Fatunla PDF
 Math. Comp. 34 (1980), 373390 Request permission
Abstract:
Some Lstable fourthorder explicit onestep numerical integration formulas which require no matrix inversion are proposed to cope effectively with systems of ordinary differential equations with large Lipschitz constants (including those having highly oscillatory solutions). The implicit integration procedure proposed in Fatunla [11] is further developed to handle a larger class of stiff systems as well as those with highly oscillatory solutions. The same pair of nonlinear equations as in [11] is solved for the stiffness/oscillatory parameters. However, the nonlinear systems are transformed into linear forms and an efficient computational procedure is developed to obtain these parameters. The new schemes compare favorably with the backward differentiation formula (DIFSUB) of Gear [13], [14] and the blended linear multistep methods of Skeel and Kong [24], and the symmetric multistep methods of Lambert and Watson [17].References

V. AMDURSKY & A. ZIV, On Numerical Treatment of Stiff, Highly Oscillatory Systems, IBM Technical Report 015, IBM Israel Scientific Center, 1974.
V. AMDURSKY & A. ZIV, On the Numerical Solution of Stiff Linear Systems of the Oscillatory Type, IBM Technical Report 032, IBM Israel Scientific Center, 1975.
V. AMDURSKY & A. ZIV, The Numerical Treatment of Linear Highly Oscillatory ODE Systems by Reduction to Nonoscillatory Types, IBM Technical Report 039, IBM Israel Scientific Center, 1976.
G. BJUREL, G. DAHLQUIST, B. LINDBERG, S. LINDEN & L. ODEN, Survey of Stiff Ordinary Differential Equations, Computer Science Report NA 70.11, Royal Institute of Technology, Stockholm, Sweden, 1970.
 J. C. Butcher, Implicit RungeKutta processes, Math. Comp. 18 (1964), 50–64. MR 159424, DOI 10.1090/S00255718196401594249
 Germund G. Dahlquist, A special stability problem for linear multistep methods, Nordisk Tidskr. Informationsbehandling (BIT) 3 (1963), 27–43. MR 170477, DOI 10.1007/bf01963532
 W. H. Enright, Second derivative multistep methods for stiff ordinary differential equations, SIAM J. Numer. Anal. 11 (1974), 321–331. MR 351083, DOI 10.1137/0711029 W. H. ENRIGHT, T. E. HULL & B. LINDBERG, “Comparing numerical methods for stiff systems of ODEs,” BIT, v. 15, 1975, pp. 1048. S. O. FATUNLA, “A new algorithm for numerical solutions of ODEs,” Comput. Math. Appl., v. 2, 1976, pp. 247253.
 Simeon O. Fatunla, A variable order onestep scheme for numerical solution of ordinary differential equations, Comput. Math. Appl. 4 (1978), no. 1, 33–41. MR 501922, DOI 10.1016/08981221(78)900056
 Simeon Ola. Fatunla, An implicit twopoint numerical integration formula for linear and nonlinear stiff systems of ordinary differential equations, Math. Comp. 32 (1978), no. 141, 1–11. MR 474830, DOI 10.1090/S00255718197804748300
 Walter Gautschi, Numerical integration of ordinary differential equations based on trigonometric polynomials, Numer. Math. 3 (1961), 381–397. MR 138200, DOI 10.1007/BF01386037
 C. W. Gear, The automatic integration of stiff ordinary differential equations. , Information Processing 68 (Proc. IFIP Congress, Edinburgh, 1968) NorthHolland, Amsterdam, 1969, pp. 187–193. MR 0260180 C. W. GEAR, “Algorithm 407: DIFSUB for solution of ordinary differential equations,” Comm. ACM, v. 14, 1971, pp. 185190.
 L. W. Jackson and Surender Kumar Kenue, A fourth order exponentially fitted method, SIAM J. Numer. Anal. 11 (1974), 965–978. MR 362926, DOI 10.1137/0711075
 J. D. Lambert, Nonlinear methods for stiff systems of ordinary differential equations, Conference on the Numerical Solution of Differential Equations (Univ. of Dundee, Dundee, 1973) Lecture Notes in Math., Vol. 363, Springer, Berlin, 1974, pp. 75–88. MR 0426436
 J. D. Lambert and I. A. Watson, Symmetric multistep methods for periodic initial value problems, J. Inst. Math. Appl. 18 (1976), no. 2, 189–202. MR 431691
 J. Douglas Lawson, Generalized RungeKutta processes for stable systems with large Lipschitz constants, SIAM J. Numer. Anal. 4 (1967), 372–380. MR 221759, DOI 10.1137/0704033
 Bengt Lindberg, On smoothing and extrapolation for the trapezoidal rule, Nordisk Tidskr. Informationsbehandling (BIT) 11 (1971), 29–52. MR 281356, DOI 10.1007/bf01935326 W. LINIGER & R. A. WILLOUGHBY, Efficient Numerical Integration of Stiff Systems of Ordinary Differential Equations, IBM Research Report RC 1970, IBM, Yorktown Heights, New York, 1969.
 W. L. Miranker and G. Wahba, An averaging method for the stiff highly oscillatory problem, Math. Comp. 30 (1976), no. 135, 383–399. MR 423817, DOI 10.1090/S00255718197604238170
 W. L. Miranker, M. van Veldhuizen, and G. Wahba, Two methods for the stiff highly oscillatory problem, Topics in numerical analysis, III (Proc. Roy. Irish Acad. Conf., Trinity Coll., Dublin, 1976) Academic Press, London, 1977, pp. 257–273. MR 0657229
 W. L. Miranker and M. van Veldhuizen, The method of envelopes, Math. Comp. 32 (1978), no. 142, 453–496. MR 494952, DOI 10.1090/S00255718197804949528
 Robert D. Skeel and Antony K. Kong, Blended linear multistep methods, ACM Trans. Math. Software 3 (1977), no. 4, 326–345. MR 461922, DOI 10.1145/355759.355762
 Arthur David Snider and Gary Charles Fleming, Approximation by aliasing with application to “Certaine” stiff differential equations, Math. Comp. 28 (1974), 465–473. MR 343637, DOI 10.1090/S00255718197403436373
 E. Stiefel and D. G. Bettis, Stabilization of Cowell’s method, Numer. Math. 13 (1969), 154–175. MR 263250, DOI 10.1007/BF02163234
Additional Information
 © Copyright 1980 American Mathematical Society
 Journal: Math. Comp. 34 (1980), 373390
 MSC: Primary 65L05
 DOI: https://doi.org/10.1090/S0025571819800559191X
 MathSciNet review: 559191