Uniformly accurate multiscale time integrators for highly oscillatory second order differential equations

Abstract

In this paper, two multiscale time integrators (MTIs), motivated from two types of multiscale decomposition by either frequency or frequency and amplitude, are proposed and analyzed for solving highly oscillatory second order differential equations with a dimensionless parameter $0<\epsilon \le 1$. In fact, the solution to this equation propagates waves with wavelength at $O({\epsilon }^2)$ when $0<\epsilon \ll 1$, which brings significantly numerical burdens in practical computation. We rigorously establish two independent error bounds for the two MTIs at $O({\tau }^2/{\epsilon }^2)$ and $O({\epsilon }^2)$ for $\epsilon \in (0,1]$ with $\tau >0$ as step size, which imply that the two MTIs converge uniformly with linear convergence rate at $O(\tau )$ for $\epsilon \in (0,1]$ and optimally with quadratic convergence rate at $O({\tau }^2)$ in the regimes when either $\epsilon =O(1)$ or $0<\epsilon \le \tau$. Thus the meshing strategy requirement (or $\epsilon$-scalability) of the two MTIs is $\tau =O(1)$ for $0<\epsilon \ll 1$, which is significantly improved from $\tau =O({\epsilon }^3)$ and $\tau =O({\epsilon }^2)$ requested by finite difference methods and exponential wave integrators to the equation, respectively. Extensive numerical tests and comparisons with those classical numerical integrators are reported, which gear towards better understanding on the convergence and resolution properties of the two MTIs. In addition, numerical results support the two error bounds very well.

39 pages and 1 figure

References in corpus (1)

• Uniformly accurate numerical schemes for highly oscillatory Klein-Gordon and nonlinear Schrödinger equations

Cited by in corpus (11)

• Comparison of numerical methods for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime
• A uniformly accurate (UA) multiscale time integrator Fourier pseoduspectral method for the Klein-Gordon-Schrodinger equations in the nonrelativistic limit regime
• Numerical methods and comparison for the Dirac equation in the nonrelativistic limit regime
• A uniformly accurate (UA) multiscale time integrator pseudospectral method for the Dirac equation in the nonrelativistic limit regime
• Hamiltonian Particle-in-Cell methods for Vlasov-Poisson equations
• Simply improved averaging for coupled oscillators and weakly nonlinear waves
• Uniformly accurate exponential-type integrators for klein-gordon equations with asymptotic convergence to classical splitting schemes in the nonlinear schroedinger limit
• Long time error analysis of the fourth-order compact finite difference methods for the nonlinear Klein-Gordon equation with weak nonlinearity
• Nonlinear Geometric Optics method based multi-scale numerical schemes for highly-oscillatory transport equations
• Uniform error bounds of an exponential wave integrator for the long-time dynamics of the nonlinear Klein-Gordon equation
• Uniformly accurate numerical schemes for the nonlinear Dirac equation in the nonrelativistic limit regime