Abstract:Oscillatory second order linear ordinary differential equations arise in many scientific calculations. Because the running times of standard solvers increase linearly with frequency when they are applied to such problems, a variety of specialized methods, most of them quite complicated, have been proposed. Here, we point out that one of the simplest approaches not only works, but yields a scheme for solving oscillatory second order linear ordinary differential equations which is significantly faster than current state-of-the-art techniques. Our method, which operates by constructing a slowly varying phase function representing a basis of solutions of the differential equation, runs in time independent of the frequency and can be applied to second order equations whose solutions are oscillatory in some regions and slowly varying in others. In the high-frequency regime, our algorithm discretizes the nonlinear Riccati equation satisfied by the derivative of the phase function via a Chebyshev spectral collocation method and applies the Newton-Kantorovich method to the resulting system of nonlinear algebraic equations. We prove that the iterates converge quadratically to a nonoscillatory solution of the Riccati equation. The quadratic convergence of the Newton-Kantorovich method and the simple form of the linearized equations ensure that this procedure is extremely efficient. Our algorithm then extends the slowly varying phase function calculated in the high-frequency regime throughout the solution domain by solving a certain third order linear ordinary differential equation related to the Riccati equation. We describe the results of numerical experiments showing that our algorithm is orders of magnitude faster than existing schemes, including the modified Magnus method [18], the current state-of-the-art approach [7] and the recently introduced ARDC method [1].

Solving time-varying periodic coefficient Riccati differential equations via Fourier series and precise integration method

NUMERICAL APPROACH FOR A SYSTEM OF SECOND KIND VOLTERRA INTEGRAL EQUATIONS IN MAGNETO-ELECTRO-ELASTIC DYNAMIC PROBLEMS

A numerical method for computing time-periodic solutions in dissipative wave systems

Fourier-Gegenbauer Pseudospectral Method for Solving Time-Dependent One-Dimensional Fractional Partial Differential Equations with Variable Coefficients and Periodic Solutions

A High Precision Direct Integration Method for Periodically Time-Varying Linear System

Taylor-Fourier integration

Fourier series (based) multiscale method for computational analysis in science and engineering: III. Fourier series multiscale method for linear differential equation with constant coefficients

On periodic approximate solutions of ordinary differential equations

Series solutions of non-linear Riccati differential equations with fractional order

A high-order integral equation-based solver for the time-dependent Schrodinger equation

An efficient numerical algorithm for solving nonlinear fractional Volterra integro-differential equation

Exact solutions and dynamic properties of a nonlinear fourth-order time-fractional partial differential equation

Numerical solution of nonlinear Volterra integro-differential equations of fractional order by the reproducing kernel method

High-accuracy power series solutions with arbitrarily large radius of convergence for fractional nonlinear differential equations

An Efficient Finite Difference/hermite-Galerkin Spectral Method for Time-Fractional Coupled Sine-Gordon Equations on Multidimensional Unbounded Domains and Its Application in Numerical Simulations of Vector Solitons

A Multi-Domain Spectral Method for Time-Fractional Differential Equations

Existence of Positive Definite Solution to Periodic Riccati Differential Equation

Backward differentiation formula method and random forest method to solve continuous‐time differential Riccati equations

An accelerated frequency-independent solver for oscillatory differential equations

Power Series Method for Time-Integration in Computational Acoustics

An algebraic method for finding a series of exact solutions to integrable and nonintegrable nonlinear evolutions equations