Search

In this paper, two novel classes of implicit exponential Runge-Kutta (ERK) methods are studied for solving highly oscillatory systems. First of all, we analyze the symplectic conditions of two kinds of exponential integrators, and present a first-order symplectic method. In order to solve highly oscillatory problems, the highly accurate implicit ERK integrators (up to order four) are formulated by comparing the Taylor expansions of numerical and exact solutions, it is shown that the order conditions of two new kinds of exponential methods are identical to the order conditions of classical Runge-Kutta (RK) methods. Moreover, we investigate the linear stability properties of these exponential methods. Finally, numerical results not only present the long time energy preservation of the first-order symplectic method, but also illustrate the accuracy and efficiency of these formulated methods in comparison with standard ERK methods.

In this paper, two new families of fourth-order explicit exponential Runge--Kutta (ERK) methods with four stages are studied for solving first-order differential systems $y'(t)+My(t)=f(y(t))$. By comparing the Taylor series of the exact solution, the order conditions of these ERK methods are derived, which are exactly identical to the order conditions of explicit Runge--Kutta methods, and these ERK methods reduce to classical Runge--Kutta methods once $M\rightarrow \mathbf{0}$. Moreover, we analyze the stability properties and the convergence of the new methods. Several numerical examples are implemented to illustrate the accuracy and efficiency of these ERK methods by comparison with standard exponential integrators.

In this paper, one-stage explicit trigonometric integrators for solving quasilinear wave equations are formulated and studied. For solving wave equations, we first introduce trigonometric integrators as the semidiscretization in time and then consider a spectral Galerkin method for the discretization in space. We show that one-stage explicit trigonometric integrators in time have second-order convergence and the result is also true for the fullydiscrete scheme without requiring any CFL-type coupling of the discretization parameters. The results are proved by using energy techniques, which are widely applied in the numerical analysis of methods for partial differential equations.

The study of the long time conservation for numerical methods poses interesting and challenging questions from the point of view of geometric integration. In this paper, we analyze the long time energy and magnetic moment conservations of two-step symmetric methods for charged-particle dynamics. A two-step symmetric method is proposed and its long time behaviour is shown not only in a normal magnetic feld but also in a strong magnetic feld. The approaches to dealing with these two cases are based on the backward error analysis and modulated Fourier expansion, respectively. It is obtained from the analysis that the method has better long time conservations than the variational method which was researched recently in the literature.

In this paper, we investigate the long-time near-conservations of energy and kinetic energy by the widely used exponential integrators to highly oscillatory conservative systems. The modulated Fourier expansions of two kinds of exponential integrators have been constructed and the long-time numerical conservations of energy and kinetic energy are obtained by deriving two almost-invariants of the expansions. Practical examples of the methods are given and the theoretical results are confirmed and demonstrated by a numerical experiment.

In this paper, a novel class of exponential Fourier collocation methods (EFCMs) is presented for solving systems of first-order ordinary differential equations. These so-called exponential Fourier collocation methods are based on the variation-of-constants formula, incorporating a local Fourier expansion of the underlying problem with collocation methods. We discuss in detail the connections of EFCMs with trigonometric Fourier collocation methods (TFCMs), the well-known Hamiltonian Boundary Value Methods (HBVMs), Gauss methods and Radau IIA methods. It turns out that the novel EFCMs are an essential extension of these existing methods. We also analyse the accuracy in preserving the quadratic invariants and the Hamiltonian energy when the underlying system is a Hamiltonian system. Other properties of EFCMs including the order of approximations and the convergence of fixed-point iterations are investigated as well. The analysis given in this paper proves further that EFCMs can achieve arbitrarily high order in a routine manner which allows us to construct higher-order methods for solving systems of first-order ordinary differential equations conveniently. We also derive a practical fourth-order EFCM denoted by EFCM(2,2) as an illustrative example. The numerical experiments using EFCM(2,2) are implemented in comparison with an existing fourth-order HBVM, an energy-preserving collocation method and a fourth-order exponential integrator in the literature. The numerical results demonstrate the remarkable efficiency and robustness of the novel EFCM(2,2).

In this paper we discuss the efficient implementation of RKN-type Fourier collocation methods, which are used when solving second-order differential equations. The proposed implementation relies on an alternative formulation of the methods and the blended formulation. The features and effectiveness of the implementation are confirmed by the performance of the methods on two numerical tests.