Numerical methods for systems of highly oscillatory ordinary differential equations

Current research made contribution to the numerical analysis of highly oscillatory ordinary differential equations. Highly oscillatory functions appear to be at the forefront of the research in numerical analysis. In this work we developed efficient numerical algorithms for solving highly oscillatory differential equations. The main important achievements are: to the contrary of classical methods, our numerical methods share the feature that asymptotically the approximation to the exact solution improves as the frequency of oscillation grows; also our methods are computationally feasible and as such do not require fine partition of the integration interval. In this work we show that our methods introduce better accuracy of approximation as compared with the state of the art solvers in Matlab and Maple. This thesis presents methods for efficient numerical approximation of linear and non-linear systems of highly oscillatory ordinary differential equations. Phenomena of high oscillation is considered a major computational problem occurring in Fourier analysis, computational harmonic analysis, quantum mechanics, electrodynamics and fluid dynamics. Classical methods based on Gaussian quadrature fail to approximate oscillatory integrals. In this work we introduce numerical methods which share the remarkable feature that the accuracy of approximation improves as the frequency of oscillation increases. Asymptotically, our methods depend on inverse powers of the frequency of oscillation, turning the major computational problem into an advantage. Evolving ideas from the stationary phase method, we first apply the asymptotic method to solve highly oscillatory linear systems of differential equations. The asymptotic method provides a background for our next, the Filon-type method, which is highly accurate and requires computation of moments. We also introduce two novel methods. The first method, we call it the FM method, is a combination of Magnus approach and the Filon-type method, to solve matrix exponential. The second method, we call it the WRF method, a combination of the Filon-type method and the waveform relaxation methods, for solving highly oscillatory non-linear systems. Finally, completing the theory, we show that the Filon-type method can be replaced by a less accurate but moment free Levin-type method.

Name dc:type.qualificationname

Doctor of Philosophy (PhD)

Level dc:type.qualificationlevel

Doctoral

Grantor dc:publisher.institution

University of Cambridge

Year dc:date.issued

2010

Author dc:creator

• Khanamiryan, Marianna

Advisor dc:contributor.advisor

• Iserles, Arieh

• Numerical analysis of differential equations
• Highly oscillatory ordinary differential equations
• Asymptotic methods
• Filon quadrature rules
• Levin method
• Lie groups methods

Licence

• https://www.repository.cam.ac.uk/bitstreams/f304c3ae-912c-46ce-80d3-34beaef9ae88/download
• https://www.rioxx.net/licenses/all-rights-reserved/

Language dc:language

eng

DOI dc:identifier.doi

https://doi.org/10.17863/CAM.16094

OAI identifier oai:identifier

oai:www.repository.cam.ac.uk:1810/226323