Numerical methods for systems of highly oscillatory ordinary differential equations
University of Cambridge
Department of Applied Mathematics and Theoretical Physics
Doctor of Philosophy (PhD)
MetadataShow full item record
Khanamiryan, M. (2010). Numerical methods for systems of highly oscillatory ordinary differential equations (Doctoral thesis). https://doi.org/10.17863/CAM.16094
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.
Numerical analysis of differential equations, Highly oscillatory ordinary differential equations, Asymptotic methods, Filon quadrature rules, Levin method, Lie groups methods
The work was supported by Trinity College, University of Cambridge.
This record's DOI: https://doi.org/10.17863/CAM.16094