Scheifele two-step methods for perturbed oscillators

Two-step methods specially adapted to the numerical integration of perturbed oscillators are obtained. The formulation of the methods is based on a refinement of classical Taylor expansions due to Scheifele [G. Scheifele, On the numerical integration of perturbed linear oscillating systems, Z. Angew. Math. Phys. 22 (1971) 186-210]. The key property is that those algorithms are able to integrate exactly harmonic oscillators with frequency @w. The methods depend on a parameter @[email protected], where h is the stepsize. Based on the B2-series theory of Coleman [J.P. Coleman, Order conditions for a class of two-step methods for y^''=f(x,y), IMA J. Numer. Anal. 23 (2003) 197-220] we derive the order conditions of this new type of method. The linear stability and phase properties are examined. The theory is illustrated with some fourth- and fifth-order explicit schemes. Numerical results carried out on an assortment of test problems (such as the integration of the orbital motion of earth satellites) show the relevance of the theory.

[1]  Generalization of the Störmer method for perturbed oscillators without explicit first derivatives , 1999 .

[2]  J. M. Franco.  A class of explicit two-step hybrid methods for second-order IVPs , 2006 .

[3]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[4]  J. M. Franco.  New methods for oscillatory systems based on ARKN methods , 2006 .

[5]  Liviu Gr. Ixaru,et al.  P-stability and exponential-fitting methods for y″″ = f(x, y) , 1996 .

[6]  D. G. Bettis,et al.  Stabilization of Cowell's method , 1969 .

[7]  Hans Van de Vyver.  A PHASE-FITTED AND AMPLIFICATION-FITTED EXPLICIT TWO-STEP HYBRID METHOD FOR SECOND-ORDER PERIODIC INITIAL VALUE PROBLEMS , 2006 .

[8]  M. M. Chawla,et al.  An explicit sixth-order method with phase-lag of order eight for y ″= f ( t , y ) , 1987 .

[9]  Hans Van de Vyver.  Phase-fitted and amplification-fitted two-step hybrid methods for y˝= f ( x,y ) , 2007 .

[10]  J. M. Franco.  A 5(3) pair of explicit ARKN methods for the numerical integration of perturbed oscillators , 2003 .

[11]  J. M. Franco.  Runge–Kutta–Nyström methods adapted to the numerical integration of perturbed oscillators , 2002 .

[12]  Ch. Tsitouras,et al.  Explicit Numerov Type Methods with Reduced Number of Stages , 2003 .

[13]  J. Ferrándiz,et al.  A general canonical transformation increasing the number of variables with application to the two-body problem , 1987 .

[14]  J. Ferrándiz,et al.  Increased accuracy of computations in the main satellite problem through linearization methods , 1992 .

[15]  P. Henrici.  Discrete Variable Methods in Ordinary Differential Equations , 1962 .

[16]  Ben P. Sommeijer,et al.  Explicit Runge-Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions , 1987 .

[17]  T. E. Simos,et al.  Symmetric Eighth Algebraic Order Methods with Minimal Phase-Lag for the Numerical Solution of the Schrödinger Equation , 2002 .

[18]  Pablo Martín,et al.  A new family of Runge–Kutta type methods for the numerical integration of perturbed oscillators , 1999, Numerische Mathematik.

[19]  John P. Coleman,et al.  Order conditions for a class of two‐step methods for y″ = f (x, y) , 2003 .

[20]  John P. Coleman.  Numerical Methods for y″ =f(x, y) via Rational Approximations for the Cosine , 1989 .

[21]  Gerhard Scheifele.  On numerical integration of perturbed linear oscillating systems , 1971 .

[22]  J. M. Franco.  Stability of explicit ARKN methods for perturbed oscillators , 2005 .

[23]  T. E. Simos,et al.  Controlling the error growth in long–term numerical integration of perturbed oscillations in one or several frequencies , 2004, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[24]  Linda R. Petzold,et al.  Numerical solution of highly oscillatory ordinary differential equations , 1997, Acta Numerica.

[25]  L.Gr. Ixaru,et al.  Numerov method maximally adapted to the Schro¨dinger equation , 1987 .

[26]  T. E. Simos,et al.  ON THE CONSTRUCTION OF EFFICIENT METHODS FOR SECOND ORDER IVPS WITH OSCILLATING SOLUTION , 2001 .

[27]  José-Miguel Farto,et al.  An algorithm for the systematic construction of solutions to perturbed problems , 1998 .

[28]  Hans Van de Vyver.  An adapted explicit hybrid method of Numerov-type for the numerical integration of perturbed oscillators , 2007, Appl. Math. Comput..

[29]  J. Lambert,et al.  Symmetric Multistip Methods for Periodic Initial Value Problems , 1976 .

[30]  José M. Ferrándiz,et al.  Multistep Numerical Methods Based on the Scheifele G -Functions with Application to Satellite Dynamics , 1997 .

[31]  Tom E. Simos.  Explicit eighth order methods for the numerical integration of initial-value problems with periodic or oscillating solutions , 1999 .

[32]  M. M. Chawla,et al.  Numerov made explicit has better stability , 1984 .