Analysis of trigonometric implicit Runge-Kutta methods

Using generalized collocation techniques based on fitting functions that are trigonometric (rather than algebraic as in classical integrators), we develop a new class of multistage, one-step, variable stepsize, and variable coefficients implicit Runge-Kutta methods to solve oscillatory ODE problems. The coefficients of the methods are functions of the frequency and the stepsize. We refer to this class as trigonometric implicit Runge-Kutta (TIRK) methods. They integrate an equation exactly if its solution is a trigonometric polynomial with a known frequency. We characterize the order and A-stability of the methods and establish results similar to that of classical algebraic collocation RK methods.

[1]  J. M. Franco.  Exponentially fitted explicit Runge-Kutta-Nyström methods , 2004 .

[2]  Kazufumi Ozawa,et al.  A functional fitting Runge-Kutta method with variable coefficients , 2001 .

[3]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[4]  L. Kramarz.  Stability of collocation methods for the numerical solution ofy″=f (x,y) , 1980 .

[5]  Kazufumi Ozawa.  A functionally fitted three-stage explicit singly diagonally implicit Runge-Kutta method , 2005 .

[6]  J. M. Franco.  An embedded pair of exponentially fitted explicit Runge-Kutta methods , 2002 .

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

[8]  Kevin Burrage,et al.  Parallel and sequential methods for ordinary differential equations , 1995, Numerical analysis and scientific computation.

[9]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[10]  Simeon Ola Fatunla,et al.  Numerical integrators for stiff and highly oscillatory differential equations , 1980 .

[11]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[12]  Nguyen Huu Cong,et al.  A-stable diagonally implicit Runge-Kutta-Nyström methods for parallel computers , 1993, Numerical Algorithms.

[13]  John P. Coleman,et al.  Mixed collocation methods for y ′′ =f x,y , 2000 .

[14]  K. Ozawa.  A four-stage implicit Runge-Kutta-Nyström method with variable coefficients for solving periodic initial value problems , 1999 .

[15]  G. Vanden Berghe,et al.  Optimal implicit exponentially-fitted Runge–Kutta methods , 2001 .

[16]  Nguyen Huu Cong,et al.  Stability of collocation-based Runge-Kutta-Nyström methods , 1991 .

[17]  Nguyen Huu Cong,et al.  On functionally-fitted Runge–Kutta methods , 2006 .

[18]  T. E. Hull,et al.  Comparing Numerical Methods for Ordinary Differential Equations , 1972 .

[19]  Beatrice Paternoster,et al.  Runge-Kutta(-Nystro¨m) methods for ODEs with periodic solutions based on trigonometric polynomials , 1998 .

[20]  R. Kress.  Numerical Analysis , 1998 .

[21]  B. M. Fulk.  MATH , 1992 .

[22]  W. Gautschi.  Numerical integration of ordinary differential equations based on trigonometric polynomials , 1961 .

[23]  H. De Meyer,et al.  Frequency determination and step-length control for exponentially-fitted Runge---Kutta methods , 2001 .

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

[25]  Beatrice Paternoster,et al.  Stability regions of one step mixed collocation methods for y″=f(x,y) , 2005 .

[26]  K. Ozawa.  A functional fitting Runge-Kutta-Nyström method with variable coefficients , 2002 .

[27]  F. Krogh,et al.  Solving Ordinary Differential Equations , 2019, Programming for Computations - Python.

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