Long-term analysis of exponential integrators for highly oscillatory conservative systems

Bin Wang,Jiyong Li,Yonglei Fang

DOI: https://doi.org/10.48550/arXiv.1809.07268

2018-09-20

Abstract: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.

Numerical Analysis

What problem does this paper attempt to address?

### The problems the paper attempts to solve

This paper aims to study the approximate conservation properties of energy and kinetic energy in the long - term when exponential integrators are used to solve highly oscillatory conservative systems. Specifically:

1. **Problem background**:
- The paper focuses on the long - term analysis of highly oscillatory conservative systems in the form of \(y'(t) = Q\nabla H(y(t))\).
- The energy \(H(y)\) and kinetic energy \(K(y)\) of the system are defined as follows respectively:
\[
H(y)=\frac{1}{2}y^{\top}\left(\frac{1}{\epsilon}M\right)y + V(y)
\]
\[
K(y)=\frac{1}{2}y^{\top}\left(\frac{1}{\epsilon}M\right)y
\]
where \(Q\) is a \(d\times d\) skew - symmetric matrix, \(M\) is a \(d\times d\) symmetric real matrix, \(V:\mathbb{R}^d\rightarrow\mathbb{R}\) is a differentiable function, and \(\epsilon\) is a small parameter.

2. **Research objectives**:
- By constructing the modulated Fourier expansions of two types of exponential integrators, derive the approximate conservation properties of energy and kinetic energy of these methods in the long - term.
- Provide specific numerical examples to verify the theoretical results.

3. **Innovative points**:
- This is the first time to use the modulated Fourier expansion technique to conduct a long - term analysis of exponential integrators for Hamiltonian ordinary differential equations.
- It is proved that the symmetric EI - T and symplectic EI - O integrators can well preserve energy and kinetic energy in the long - term.

### Main contributions

- **Theoretical analysis**: Through the modulated Fourier expansion technique, the conservation properties of energy and kinetic energy of exponential integrators in the long - term are proved.
- **Numerical experiments**: Specific numerical experiments are provided to verify the correctness of the theoretical analysis.
- **Method generalization**: It provides a basis and reference for the long - term analysis of other types of exponential integrators (such as explicit exponential integrators).

### Conclusion

Through the study of the long - term energy and kinetic energy conservation properties of exponential integrators, this paper provides an important theoretical basis for understanding and applying these numerical methods. In particular, the superior performance of the symmetric EI - T and symplectic EI - O integrators in long - term simulations has been confirmed.