• Journal of information and communication convergence engineering
• /
• v.12 no.1
• /
• pp.39-45
• /
• 2014

A numerical method for solving Hamiltonian equations is said to be symplectic if it preserves the symplectic structure associated with the equations. Various symplectic methods are widely used in many fields of science and technology. A symplectic method preserves an approximate Hamiltonian perturbed from the original Hamiltonian. It theoretically supports the effectiveness of symplectic methods for long-term integration. Although it is also related to long-term integration, numerical stability of symplectic methods have received little attention. In this paper, we consider explicit symplectic methods defined for Hamiltonian equations with Hamiltonians of the special form, and study their numerical stability using the harmonic oscillator as a test equation. We propose a new stability criterion and clarify the stability of some existing methods that are visually based on the criterion. We also derive a new method that is better than the existing methods with respect to a Courant-Friedrichs-Lewy condition for hyperbolic equations; this new method is tested through a numerical experiment with a nonlinear wave equation.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.173-182
• /
• 2014

The mixed convolved action principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation, convolution, and fractional calculus. In this paper, its potential in the development of numerical methods for transient problems in various dynamical systems when adopting temporally second order approximation is investigated. For this, the classical single-degree-of-freedom linear elastic dynamical systems are primarily considered to investigate computational characteristics of the developed algorithms. For the undamped system, all the developed algorithms are symplectic with respect to the time step. For the damped system, they are shown to be accurate with good convergence characteristics.