The extended discrete gradient method is an extension of traditional discrete gradient method, which is specially designed to solve oscillatory Hamiltonian systems efficiently while preserving their energy exactly. In...The extended discrete gradient method is an extension of traditional discrete gradient method, which is specially designed to solve oscillatory Hamiltonian systems efficiently while preserving their energy exactly. In this paper, based on the extended discrete gradient method, we present an efficient approach to devising novel schemes for numerically solving conservative (dissipative) nonlinear wave partial differential equations. The new scheme can preserve the energy exactly for conservative wave equations. With a minor remedy to the extended discrete gradient method, the new scheme is applicable to dissipative wave equations. Moreover, it can preserve the dissipation structure for the dissipative wave equation as well. Another important property of the new scheme is that it is linearly-fitted, which guarantees much fast convergence for the fixed-point iteration which is required by an energy-preserving integrator. The efficiency of the new scheme is demonstrated by some numerical examples.展开更多
文摘The extended discrete gradient method is an extension of traditional discrete gradient method, which is specially designed to solve oscillatory Hamiltonian systems efficiently while preserving their energy exactly. In this paper, based on the extended discrete gradient method, we present an efficient approach to devising novel schemes for numerically solving conservative (dissipative) nonlinear wave partial differential equations. The new scheme can preserve the energy exactly for conservative wave equations. With a minor remedy to the extended discrete gradient method, the new scheme is applicable to dissipative wave equations. Moreover, it can preserve the dissipation structure for the dissipative wave equation as well. Another important property of the new scheme is that it is linearly-fitted, which guarantees much fast convergence for the fixed-point iteration which is required by an energy-preserving integrator. The efficiency of the new scheme is demonstrated by some numerical examples.