Abstract:The construction of modified equations is an important step in the backward error analysis of symplectic integrator for Hamiltonian systems. In the context of partial differential equations, the standard construction leads to modified equations with increasingly high frequencies which increase the regularity requirements on the analysis. In this paper, we consider the next order modified equations for the implicit midpoint rule applied to the semilinear wave equation to give a proof-of-concept of a new construction which works directly with the variational principle. We show that a carefully chosen change of coordinates yields a modified system which inherits its analytical properties from the original wave equation. Our method systematically exploits additional degrees of freedom by modifying the symplectic structure and the Hamiltonian together.
What problem does this paper attempt to address?
This paper attempts to solve the problem of high - frequency components brought by the standard modified equation construction method when applying variational integrators to Hamiltonian partial differential equations, especially semi - linear wave equations. Specifically:
1. **Limitations of the standard construction method**: For Hamiltonian partial differential equations, the standard modified equation construction method will lead to an asymptotic series expansion containing arbitrarily high - order frequency terms, which increases the regularity requirements in the analysis, and when the time step is close to the CFL limit with the space step, these high - frequency components will make the numerical solution unstable.
2. **Objectives of the new construction method**: This paper proposes a new construction method that avoids introducing high - frequency components by directly using the variational principle. This method makes the modified system inherit the analytical properties of the original wave equation through an appropriate transformation of coordinates, and will not introduce frequency components beyond the scale of the original problem.
3. **Specific application scenarios**: Taking the application of the implicit midpoint rule to the semi - linear wave equation as an example, the article shows how to construct the modified equation by choosing an appropriate coordinate transformation, thus proving the effectiveness of the new method. The results show that the new modified equation not only does not introduce additional high - frequency components, but also is locally well - posed in the energy space and can maintain good numerical stability without depending on the time - step parameter.
In summary, this paper aims to solve the problem of high - frequency components introduced by existing methods when dealing with Hamiltonian partial differential equations, and proposes a new construction method based on the variational principle to ensure that the modified equation retains the physical characteristics of the original system and improves the stability and accuracy of the numerical solution.
### Formula summary
- **Original semi - linear wave equation**:
\[
\ddot{u}=\partial_{xx}u + f(u)
\]
- **Discrete Lagrangian**:
\[
L(u_k, u_{k + 1};h)=h\mathcal{L}\left(\frac{u_k + u_{k + 1}}{2},\frac{u_{k + 1}-u_k}{h}\right)
\]
- **Modified Lagrangian**:
\[
\mathcal{L}_{\text{mod}}(u,\dot{u},\ddot{u};h)=\mathcal{L}(u,\dot{u})-h^2\int_{S^1}\left(\frac{1}{24}\ddot{u}^2-\frac{1}{8}(\partial_x\dot{u})^2+\frac{1}{8}f'(u)\dot{u}^2\right)dx
\]
- **Modified Hamiltonian**:
\[
H_{\text{mod}}=H + h^2\int_{S^1}\left(\frac{1}{12}(\partial_x\dot{u})^2-\frac{1}{12}f'(u)\dot{u}^2\right)dx
\]
- **Operator \(K(u)\)**:
\[
K(u)v=\left(1-\frac{h^2}{6}f'(u)-\frac{h^2}{6}\partial_{xx}\right)v
\]
These formulas show how the new method effectively eliminates high - frequency components while maintaining the physical characteristics of the original system, thereby improving the stability and accuracy of the numerical solution.