Abstract:Using a recent characterization of energy-preserving B-series, we derive the explicit conditions on the coefficients of a Runge-Kutta method that ensure energy preservation (for Hamiltonian systems) up to a given order in the step size, which we refer to as the pseudo-energy-preserving (PEP) order. We study explicit Runge-Kutta methods with PEP order higher than their classical order. We provide examples of such methods up to PEP order six, and test them on Hamiltonian ODE and PDE systems. We find that these methods behave similarly to exactly energy-conservative methods over moderate time intervals and exhibit significantly smaller errors, relative to other Runge-Kutta methods of the same order, for moderately long-time simulations.

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is to improve the energy - preserving performance of explicit Runge - Kutta (RK) methods in the numerical integration of Hamiltonian systems. Specifically, the author aims to derive the conditions that ensure the RK method preserves energy as high - order as possible for a given step - size order. This property is known as the Pseudo - Energy - Preserving (PEP) order. ### Main problems 1. **Importance of energy preservation**: For the numerical integration of Hamiltonian systems, energy preservation is crucial. Traditional RK methods cannot accurately preserve the energy of all Hamiltonian systems, but the energy error can be reduced by designing high - order energy - preserving methods under specific conditions. 2. **Definition of Pseudo - Energy - Preserving (PEP) order**: An RK method with a classical order \(p\) is said to have a PEP order \(q\) if the Hamiltonian error after one step is \(O(h^{q + 1})\), where \(q\geq p\). 3. **Objective**: To develop explicit RK methods with the maximum PEP order and study the performance of these methods in different Hamiltonian systems, especially in long - time simulations. ### Solutions - **B - series theory**: Use B - series theory to analyze and derive energy - preserving conditions. The B - series is a tool for representing the errors of numerical methods, and through it, the energy - preserving conditions of the RK method can be directly written. - **Condition derivation**: Based on the B - series theory, derive the algebraic conditions that the RK method satisfying the PEP order needs to meet. These conditions involve the coefficients of the RK method (such as the coefficients in the Butcher tableau), and become more complex as the PEP order increases. - **Method construction**: Search for explicit RK methods satisfying a specific PEP order through numerical optimization techniques (such as the fmincon function in MATLAB). For a given number of stages \(s\) and classical order \(p\), find the method with the highest feasible PEP order \(q\). ### Experimental verification - **Convergence test**: Verify the theoretical convergence and energy - preserving performance of the newly constructed PEP methods by solving some standard Hamiltonian systems (such as exponential - conserving systems, nonlinear oscillators, etc.). - **Long - term behavior comparison**: Compare the PEP methods with classical RK methods and pseudo - symplectic methods, especially in long - time simulations, and observe their error growth. ### Conclusion The paper shows that explicit RK methods with high PEP orders exhibit smaller energy errors in long - time simulations, and the error growth is nearly linear within a certain time, which is better than traditional RK methods. This indicates that although these methods cannot completely preserve energy, they can effectively control the energy error for a long time, thus providing more accurate numerical solutions. Through the above methods and experiments, the paper successfully solves the problem of how to achieve high - order energy preservation in explicit RK methods and provides new ideas for future numerical algorithm design.

Pseudo-Energy-Preserving Explicit Runge-Kutta Methods

ENERGY-PRESERVING CONTINUOUS-STAGE EXPONENTIAL RUNGE--KUTTA INTEGRATORS FOR EFFICIENTLY SOLVING HAMILTONIAN SYSTEMS

Energy-superconvergent Runge-Kutta Time Discretizations

Energy-preserving continuous-stage partitioned Runge-Kutta methods

Arbitrary high-order linearly implicit energy-preserving algorithms for Hamiltonian PDEs

Strong stability preserving explicit Runge-Kutta methods of maximal effective order

Exponential Runge-Kutta schemes for inhomogeneous Boltzmann equations with high order of accuracy

Energy Conservative Relaxation-Free Runge-Kutta Schemes

High order strong stability preserving multi-derivative implicit and IMEX Runge--Kutta methods with asymptotic preserving properties

Simplified explicit exponential Runge-Kutta methods without order reduction

Strong stability preserving two-step Runge-Kutta methods

Strong Stability Preserving Integrating Factor Runge-Kutta Methods

Explicit Runge Kutta Methods that Alleviate Order Reduction

Symplectic Exponential Runge-Kutta Methods for Solving Nonlinear Hamiltonian Systems

Fourth-Order Paired-Explicit Runge-Kutta Methods

Cost-reduction implicit exponential Runge-Kutta methods for highly oscillatory systems

Efficient implementation of partitioned stiff exponential Runge-Kutta methods

Explicit and implicit error inhibiting schemes with post-processing

Energy-preserving continuous-stage Runge-Kutta-Nyström methods

Many-Stage Optimal Stabilized Runge–Kutta Methods for Hyperbolic Partial Differential Equations

Arbitrarily high-order energy-preserving schemes for the Camassa-Holm equation