An energy‐preserving Crank–Nicolson Galerkin method for Hamiltonian partial differential equations

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

[1]  G. Quispel,et al.  A new class of energy-preserving numerical integration methods , 2008 .

[2]  Jie Shen,et al.  A Fourth-Order Time-Splitting Laguerre-Hermite Pseudospectral Method for Bose-Einstein Condensates , 2005, SIAM J. Sci. Comput..

[3]  Abdul-Qayyum M. Khaliq,et al.  Higher order exponential time differencing scheme for system of coupled nonlinear Schrödinger equations , 2014, Appl. Math. Comput..

[4]  Christophe Besse,et al.  Order Estimates in Time of Splitting Methods for the Nonlinear Schrödinger Equation , 2002, SIAM J. Numer. Anal..

[5]  G. Akrivis,et al.  On fully discrete Galerkin methods of second-order temporal accuracy for the nonlinear Schrödinger equation , 1991 .

[6]  Qin Sheng,et al.  Solving the Generalized Nonlinear Schrödinger Equation via Quartic Spline Approximation , 2001 .

[7]  Takayasu Matsuo,et al.  High-order schemes for conservative or dissipative systems , 2003 .

[8]  Georgios Akrivis,et al.  On optimal order error estimates for the nonlinear Schro¨dinger equation , 1993 .

[9]  M. S. Ismail.  Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method , 2008, Math. Comput. Simul..

[10]  R. McLachlan.  Symplectic integration of Hamiltonian wave equations , 1993 .

[11]  Yulong Xing,et al.  Fourth Order Exponential Time Differencing Method with Local Discontinuous Galerkin Approximation for Coupled Nonlinear Schrödinger Equations , 2015 .

[12]  Tingchun Wang,et al.  Fourth-order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions , 2013, J. Comput. Phys..

[13]  Abdul-Qayyum M. Khaliq,et al.  Exponential time differencing Crank-Nicolson method with a quartic spline approximation for nonlinear Schrödinger equations , 2014, Appl. Math. Comput..

[14]  L. Brugnano,et al.  A simple framework for the derivation and analysis of effective one-step methods for ODEs , 2010, Appl. Math. Comput..

[15]  Mechthild Thalhammer,et al.  High-Order Exponential Operator Splitting Methods for Time-Dependent Schrödinger Equations , 2008, SIAM J. Numer. Anal..

[16]  S. Reich,et al.  Multi-symplectic spectral discretizations for the Zakharov–Kuznetsov and shallow water equations , 2001 .

[17]  K. Feng.  Difference schemes for Hamiltonian formalism and symplectic geometry , 1986 .

[18]  O. Gonzalez.  Time integration and discrete Hamiltonian systems , 1996 .

[19]  Zhi-Zhong Sun,et al.  On the L∞ convergence of a difference scheme for coupled nonlinear Schrödinger equations , 2010, Comput. Math. Appl..

[20]  Jialin Hong,et al.  Compact and efficient conservative schemes for coupled nonlinear Schrödinger equations , 2015 .

[21]  Luigi Brugnano,et al.  Energy conservation issues in the numerical solution of the semilinear wave equation , 2014, Appl. Math. Comput..