Title:Error Analysis on a Novel Class of Exponential Integrators with Local Linear Extension Techniques for Highly Oscillatory ODEs

Abstract:This paper studies a class of non-autonomous highly oscillatory ordinary differential equations (ODEs) featuring a linear component inversely proportional to a small parameter $\varepsilon$ with purely imaginary eigenvalues, alongside an $\varepsilon$-independent nonlinear component. When $0<\varepsilon\ll 1$, the rapidly oscillatory solution constrains the step size selection and numerical accuracy, resulting in significant computational challenges. Motivated by linearization through introducing auxiliary polynomial variables, a new class of explicit exponential integrators (EIs) has recently been developed. The methods do not require the linear part to be diagonal or with all eigenvalues to be integer multiples of a fixed value - a general assumption in multiscale methods - and attain arbitrarily high convergence order without any order conditions. The main contribution of this work is to establish a rigorous error analysis for the new class of methods. To do this, we first demonstrate the equivalence between the high-dimensional system and the original problem by employing algebraic techniques. Building upon these fundamental results, we prove that the numerical schemes have a uniform convergence order of $O(h^{k+1})$ for the solution when using at most $k$-degree auxiliary polynomial variables with time step sizes smaller than $\varepsilon$. For larger step sizes under the bounded oscillatory energy condition, the methods achieve a convergence order of $O(\varepsilon h^k)$ for the solution. These theoretical results are further applied to second-order oscillatory equations, yielding improved uniform accuracy with respect to $\varepsilon$. Finally, numerical experiments confirm the optimality of the derived error estimates.