Abstract
A semidiscretization based method for solving Hamiltonian partial differential equations is proposed in this article. Our key idea consists of two approaches. First, the underlying equation is discretized in space via a selected finite element method and the Hamiltonian PDE can thus be casted to Hamiltonian ODEs based on the weak formulation of the system. Second, the resulting ordinary differential system is solved by an energy-preserving integrator. The relay leads to a fully discretized and energy-preserved scheme. This strategy is fully realized for solving a nonlinear Schrodinger equation through a combination of the Galerkin discretization in space and a Crank–Nicolson scheme in time. The order of convergence of our new method is O ( τ 2 + h 2 ) if the discrete L2-norm is employed. An error estimate is acquired and analyzed without grid ratio restrictions. Numerical examples are given to further illustrate the conservation and convergence of the energy-preserving scheme constructed.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1485–1504, 2016
Sign-in/Register to access full text options 
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: Numerical Methods for Partial Differential Equations
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.
