Back to Search
Start Over
Efficient method for solving highly oscillatory ordinary differential equations with applications to physical systems
- Source :
- Physical Review Research. 2
- Publication Year :
- 2020
- Publisher :
- American Physical Society (APS), 2020.
-
Abstract
- We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schr\"odinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications.<br />Comment: Physical Review Research, accepted. 23 pages, 15 figures. The associated code is available online at https://github.com/fruzsinaagocs/oscode
- Subjects :
- Cosmology and Nongalactic Astrophysics (astro-ph.CO)
Oscillation
Numerical analysis
Physical system
FOS: Physical sciences
Numerical Analysis (math.NA)
Computational Physics (physics.comp-ph)
Ordinary differential equation
FOS: Mathematics
Applied mathematics
Mathematics - Numerical Analysis
Astrophysics - Instrumentation and Methods for Astrophysics
Instrumentation and Methods for Astrophysics (astro-ph.IM)
Physics - Computational Physics
Astrophysics - Cosmology and Nongalactic Astrophysics
Mathematics
Subjects
Details
- ISSN :
- 26431564
- Volume :
- 2
- Database :
- OpenAIRE
- Journal :
- Physical Review Research
- Accession number :
- edsair.doi.dedup.....e61a37f5664ed85715e178b72376a89f