Explicit multi-symplectic extended leap-frog methods for Hamiltonian wave equations

被引:27
作者
Shi, Wei [2 ]
Wu, Xinyuan [2 ]
Xia, Jianlin [1 ]
机构
[1] Purdue Univ, Dept Math, W Lafayette, IN 47907 USA
[2] Nanjing Univ, State Key Lab Novel Software Technol, Dept Math, Nanjing 210093, Jiangsu, Peoples R China
基金
美国国家科学基金会;
关键词
Multi-symplectic integrators; Hamiltonian wave equations; Extended leap-frog methods; Conservation laws; Dispersive properties; TRIGONOMETRICALLY-FITTED FORMULAS; RUNGE-KUTTA METHODS; NUMEROV-TYPE METHOD; NUMERICAL-INTEGRATION; NYSTROM METHODS; HIGH-ORDER; CONSERVATION; PDES;
D O I
10.1016/j.jcp.2012.07.004
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
In this paper, we study the integration of Hamiltonian wave equations whose solutions have oscillatory behaviors in time and/or space. We are mainly concerned with the research for multi-symplectic extended Runge-Kutta-Nystrom (ERKN) discretizations and the corresponding discrete conservation laws. We first show that the discretizations to the Hamiltonian wave equations using two symplectic ERKN methods in space and time respectively lead to an explicit multi-symplectic integrator (Eleap-frogI). Then we derive another multi-symplectic discretization using a symplectic ERKN method in time and a symplectic partitioned Runge-Kutta method, which is equivalent to the well-known Stormer-Verlet method in space (Eleap-frogII). These two new multi-symplectic schemes are extensions of the leap-frog method. The numerical stability and dispersive properties of the new schemes are analyzed. Numerical experiments with comparisons are presented, where the two new explicit multi-symplectic methods and the leap-frog method are applied to the linear wave equation and the Sine-Gordon equation. The numerical results confirm the superior performance and some significant advantages of our new integrators in the sense of structure preservation. (C) 2012 Elsevier Inc. All rights reserved.
引用
收藏
页码:7671 / 7694
页数:24
相关论文
共 45 条
[1]  
[Anonymous], FRONT PROSPECTS CONT
[2]   Multisymplectic box schemes and the Korteweg-de Vries equation [J].
Ascher, UM ;
McLachlan, RI .
APPLIED NUMERICAL MATHEMATICS, 2004, 48 (3-4) :255-269
[4]   Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity [J].
Bridges, TJ ;
Reich, S .
PHYSICS LETTERS A, 2001, 284 (4-5) :184-193
[5]   A geometric formulation of the conservation of wave action and its implications for signature and the classification of instabilities [J].
Bridges, TJ .
PROCEEDINGS OF THE ROYAL SOCIETY A-MATHEMATICAL PHYSICAL AND ENGINEERING SCIENCES, 1997, 453 (1962) :1365-1395
[6]   Trigonometrically fitted explicit Numerov-type method for periodic IVPs with two frequencies [J].
Fang, Yonglei ;
Song, Yongzhong ;
Wu, Xinyuan .
COMPUTER PHYSICS COMMUNICATIONS, 2008, 179 (11) :801-811
[7]   A trigonometrically fitted explicit Numerov-type method for second-order initial value problems with oscillating solutions [J].
Fang, Yonglei ;
Wu, Xinyuan .
APPLIED NUMERICAL MATHEMATICS, 2008, 58 (03) :341-351
[8]   Runge-Kutta-Nystrom methods adapted to the numerical integration of perturbed oscillators [J].
Franco, JM .
COMPUTER PHYSICS COMMUNICATIONS, 2002, 147 (03) :770-787
[9]   Linear PDEs and numerical methods that preserve a multisymplectic conservation law [J].
Frank, J ;
Moore, BE ;
Reich, S .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2006, 28 (01) :260-277
[10]   A new family of Runge-Kutta type methods for the numerical integration of perturbed oscillators [J].
González, AB ;
Martín, P ;
Farto, JM .
NUMERISCHE MATHEMATIK, 1999, 82 (04) :635-646