Abstract:In this paper, we present continuous-stage partitioned Runge-Kutta (csPRK) methods for energy-preserving integration of Hamiltonian systems. A sufficient condition for the energy preservation of the csPRK methods is derived. It is shown that the presented condition contains the existing condition for energy-preserving continuous-stage Runge-Kutta methods as a special case. A noticeable and interesting result is that when we use the simplifying assumptions of order conditions and the normalized shifted Legendre polynomials for constructing high-order energy-preserving csPRK methods, both the Butcher "weight" coefficients $B_\tau$ and $\widehat{B}_\tau$ must be equal to $1$. As illustrative examples, new energy-preserving integrators are acquired by virtue of the presented condition, and for the sake of verifying our theoretical results, some numerical experiments are reported.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve the energy - preservation problem in the numerical integration of Hamiltonian systems. Specifically, the author proposes the continuous - stage partitioned Runge - Kutta (csPRK) method for energy - preserving numerical integration of Hamiltonian systems. The following are the core problems of this paper:
1. **Energy - preservation problem**:
- Hamiltonian systems are of great significance in the fields of physics and engineering, and one of their characteristics is the conservation of total energy. However, traditional numerical methods may lead to energy drift during long - time integration, affecting the accuracy of calculation results.
- The goal of the paper is to develop a new numerical method - the continuous - stage partitioned Runge - Kutta (csPRK) method - to ensure that the energy of the system remains unchanged during the numerical integration process.
2. **Limitations of existing methods**:
- Existing energy - preserving methods (such as the continuous - stage Runge - Kutta method) have certain limitations, especially when dealing with high - order energy - preserving problems.
- The paper extends the existing energy - preserving methods by introducing new conditions and construction methods, making them more effectively applicable to high - dimensional Hamiltonian systems.
3. **Combination of theory and application**:
- The paper not only provides sufficient theoretical conditions but also verifies the effectiveness of these conditions through specific numerical experiments.
- The author shows how to use simplified assumptions and normalized shifted Legendre polynomials to construct high - order energy - preserving csPRK methods and proves the superiority of these methods in practical applications.
### Markdown representation of formulas
Some of the key formulas involved in the paper are as follows:
- **Basic equations of Hamiltonian systems**:
\[
\dot{p}=-\nabla_q H(p, q), \quad \dot{q}=\nabla_p H(p, q), \quad p(t_0) = p_0\in\mathbb{R}^d, \quad q(t_0)=q_0\in\mathbb{R}^d
\]
- **Definition of the csPRK method**:
\[
Y_\tau = y_0+h\int_0^1 A_{\tau,\sigma} f(t_0 + C_\sigma h, Y_\sigma, Z_\sigma)\,d\sigma,
\]
\[
Z_\tau = z_0+h\int_0^1 \hat{A}_{\tau,\sigma} g(t_0+\hat{C}_\sigma h, Y_\sigma, Z_\sigma)\,d\sigma,
\]
\[
y_1 = y_0+h\int_0^1 B_\tau f(t_0 + C_\tau h, Y_\tau, Z_\tau)\,d\tau,
\]
\[
z_1 = z_0+h\int_0^1 \hat{B}_\tau g(t_0+\hat{C}_\tau h, Y_\tau, Z_\tau)\,d\tau.
\]
- **Sufficient conditions for energy preservation**:
\[
A_{0,\sigma}=0, \quad A_{1,\sigma}=B_\sigma, \quad \forall\sigma\in[0,1],
\]
\[
\hat{A}_{0,\sigma}=0, \quad \hat{A}_{1,\sigma}=\hat{B}_\sigma, \quad \forall\sigma\in[0,1],
\]
\[
\frac{\partial}{\partial\tau}A_{\tau,\sigma}=\frac{\partial}{\partial\