Abstract
A fourth-order three-stage symplectic integrator similar to the second-order Störmer—Verlet method has been proposed and used before [Chin. Phys. Lett. 28 (2011) 070201; Eur. Phys. J. Plus 126 (2011) 73]. Continuing the work initiated in the publications, we investigate the numerical performance of the integrator applied to a one-dimensional wave equation, which is expressed as a discrete Hamiltonian system with a fourth-order central difference approximation to a second-order partial derivative with respect to the space variable. It is shown that the Störmer—Verlet-like scheme has a larger numerical stable zone than either the Störmer—Verlet method or the fourth-order Forest—Ruth symplectic algorithm, and its numerical errors in the discrete Hamiltonian and numerical solution are also smaller.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.