Multistep numerical methods for the integration of oscillatory problems in several frequencies

Authors:
F. García-Alonso;J. A. Reyes;J. M. Ferrándiz;J. Vigo-Aguiar
Affiliations:
Departamento de Matemática Aplicada, Universidad de Alicante, Apartado 99, E-03080 Alicante, Spain;Departamento de Matemática Aplicada, Universidad de Alicante, Apartado 99, E-03080 Alicante, Spain;Departamento de Matemática Aplicada, Universidad de Alicante, Apartado 99, E-03080 Alicante, Spain;Dept. of Applied Mathematics, University of Salamanca, Plaza de los Caídos s/n (37008) Salamanca, Spain
Venue:
Advances in Engineering Software
Year:
2009

Citing 6
Cited 2

Numerical methods for ordinary differential systems: the initial value problem

Numerical methods for ordinary differential systems: the initial value problem
A new numerical method for the integration of highly oscillatory second-order ordinary differential equations

Selected papers of the sixth conference on Numerical Treatment of Differential Equations
Multistep Numerical Methods Based on the Scheifele G-Functions with Application to Satellite Dynamics

SIAM Journal on Numerical Analysis
A General Procedure For the Adaptation of Multistep Algorithms to the Integration of Oscillatory Problems

SIAM Journal on Numerical Analysis
Higher-order variable-step algorithms adapted to the accurate numerical integration of perturbed oscillators

Computers in Physics
The various order explicit multistep exponential fitting for systems of ordinary differential equations

Journal of Computational and Applied Mathematics

An algorithm for exact integration of some forced and damped oscillatory problems, based in the τ-functions

Advances in Engineering Software
A new approach for multistep numerical methods in several frequencies for perturbed oscillators

Advances in Engineering Software

Quantified Score

Hi-index	0.00

Visualization

Abstract

The perturbed harmonic oscillator appear frequently in the mathematical modelling of many problems in physics and engineering. The harmonic oscillator has a special purpose in Astrodynamics, because the Kunstaanheimo-Stiefel (KS) and Burdet-Ferrandiz (BF) transformations reduce the Kepler problem to harmonic oscillators. A new multi-step methods of numerical integration are introduced that generalize SMF ones. They are defined for arbitrary order and have similar properties to the former methods. Modified methods allowing step variations, whose coefficients are computed from relations of recurrence, are derived, what considerably improve the implementation of the algorithms. These methods are based in a sequence of analytical @f-functions dependent on two parameters @a and @b that generalizes the Scheifele's G-functions and that, under wide hypothesis, allow us obtain the solution of harmonic oscillator. In this paper a new methodology is generated to solve the problem which the @f-functions series create regarding the calculus of recurrence relations transforming the method into a multistep scheme. The methodology is implemented in a computational algorithm, which lets us solve in a general way the problems of physics and engineering which are modalized by means of the study of the harmonic oscillators. Numerical examples already used by other authors are presented. They show how the new developed methods in this paper may compete in accuracy or efficiency with other well-reputed algorithms.